ALBERT

Description: 〈span〉〈div〉ABSTRACT〈/div〉Late Cretaceous–to–present-day mixed carbonate–clastic deposition along the Nicaraguan platform, western Caribbean Sea, has evolved from a tectonically controlled, rifted upper Eocene shallow–to–deep-marine carbonate–siliciclastic shelf to an upper Miocene–to–present-day tectonically stable shallow-marine carbonate platform and passive margin. By integrating subsurface data of 287 two-dimensional seismic lines and 27 wells, we interpret the Cenozoic stratigraphic sequence as 3 cycles of transgression and regression beginning with an upper Eocene rhodolitic–algal carbonate shelf that interfingered with marginal siliciclastic sediments derived from exposed areas of Central America bordering the margin to the west. During the middle Eocene, a carbonate platform was established with both rimmed reefs and isolated patch reefs. A late Eocene forced regression produced widespread erosion and subaerial exposure across much of the platform and was recorded by a regional unconformity. The Oligocene–upper Miocene sedimentary record includes a southeastward prograding delta of the proto-Coco river, which drained the emergent area of what is now northern Nicaragua. The late Miocene–to–present-day period marks a period of strong subsidence with the development of small pinnacle reefs. We describe favorable petroleum system elements of the Nicaraguan platform that include (1) Eocene fossiliferous limestone source rocks documented as thermally mature in vintage exploration wells and seen as active gas chimneys emanating from inferred carbonate reservoirs; (2) upper–to–middle Eocene reservoirs in patch and pinnacle reefs, middle Eocene calcareous slumps, and Oligocene fluvial-deltaic facies documented in wells; and (3) regional seal intervals that consist of both regional unconformities and Eocene–Oligocene intraformational shale.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Fault damage zones may significantly affect subsurface fluid migration and the development of unconventional resources. Most analyses of fault damage zones are based on direct field observations, and we expand these analyses to the subsurface by investigating the damage zone structure of an approximately 32-km (∼10〈sup〉5〈/sup〉-ft)-long right-lateral strike-slip fault in Oklahoma. We used the three-dimensional (3-D) seismic attribute of coherence to first define its regional and background levels, and then we evaluated the damage zone dimensions at multiple sites. We found damage zone thickness of approximately 1600 m (∼5300 ft) at a segment that is dominated by subsidiary faults, and it is slightly thicker at a segment with a pull-apart basin. The damage zone intensity decays exponentially with distance from the fault core, in agreement with field observations and distribution of seismic events. The coherence map displays a strong asymmetry of the damage zone between the two sides of the 3-D fault, which is related to the subsidiary structures of the fault zone. We discuss the effects of heterogeneous stress field on damage zone evolution through the detected subsidiary structures. It appears that seismic coherence is an effective tool for subsurface characterization of fault damage zones.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Some fault zones leak vertically to the ground surface or seafloor, whereas most others remain naturally sealed. Understanding the factors that cause this leakage is essential for predicting and preventing such leakage for both conventional reservoir development and subsurface CO〈sub〉2〈/sub〉 storage. This study, a comparison of leaking and nonleaking natural CO〈sub〉2〈/sub〉 gas accumulations, provides such constraints. We compare and contrast trap configurations, fluid pressures, and stress states for several natural CO〈sub〉2〈/sub〉 accumulations from the Colorado Plateau. Extensive surface geologic data are integrated with subsurface data from a large suite of groundwater and hydrocarbon wells. Leakage of CO〈sub〉2〈/sub〉 is documented by geochemical surveys and the occurrence of extensive travertine deposits. The leakage occurs exclusively in fault fracture damage zones where the total fluid pressure reduces the minimum horizontal effective stress to approximately zero. These results are consistent with natural and accidentally induced fault seeps from some deep-water hydrocarbon reservoirs. These criteria can be used to evaluate the potential for fault zones to provide vertical leakage pathways and loss of fluid containment.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The three-dimensionally complex, highly progradational mixed siliciclastic–carbonate strata of the San Andres and Grayburg Formations have long been the backbone of conventional hydrocarbon reservoir production from the Permian Basin, and significant recovery continues via waterflooding and CO〈sub〉2〈/sub〉 injection. Besides, nonreservoir equivalents of these formations have recently taken increasing significance as produced water disposal targets. However, seismic-stratigraphic interpretations are challenged by complex internal shelfal-stratal geometries and numerous laterally continuous but vertically thin fluid barriers in overlying platforms. We built a three-dimensional (3-D) geocellular model of Guadalupian 8–13 high-frequency sequences (G8–G13 HFSs) and then conducted forward seismic modeling (35-Hz 0° phase). This allows investigations on the validity of applying conventional reflection-geometry–based interpretation to delineate the G9 HFS top and base, which can potentially serve as bounding/constraining surfaces for upper San Andres shelf–Grayburg platform reservoirs. This study contributes to 3-D modeling methodologies by introducing a query tree to select geostatistical methods for modeling dual-scale heterogeneities and by integrating data from diverse sources for seamless and realistic 3-D models. Our seismic-stratigraphic evaluation demonstrates that conventional reflection–geometry-based interpretation does not adequately resolve the G9 top and base; deviations from the geocellular model reach up to 80 m (260 ft) and are thus well beyond the maximum acceptable error limits of ±0.5 wavelength. We suggest improving conventional interpretations of the G9 base by selective interpolation or mixed-polarity event picking near the error-prone shelf margin and upper slope. Besides, instead of picking the highly discontinuous seismic peak as G9 top, bulk-shifting of a shallower trough horizon near actual G10 top should deliver a more accurate surface representing G9 top.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Paleogene shale of the Dongying depression, a continental basin in eastern China, is taken as the study subject to examine the microscopic features of lacustrine shale reservoirs in the oil window. This study shows that shale pores in this evolutionary stage are present at the micrometer to nanometer scale, but fractures commonly have extension distances at the millimeter scale. Pores and fractures can be divided into three types, namely, primary pores, secondary pores, and cracks. Primary pores commonly have good connectivity at shallow burial depth. With the increase of burial depth, primary porosity is reduced because of compaction and cementation. Secondary pores are important in shale, including dissolved pores inside grains and at grain edge, and dissolution pores inside the hybrid of organic matter (OM) and clay minerals, and evaporite minerals, including carbonates or sulfates. Types of cracks were observed: bedding fissures, dissolution fractures, and structural fractures. The development of bedding fissures is related to the deposition of shale laminae. The formation of dissolution fractures is related to acidic fluids, such as organic acids and hydrogen sulfide, whereas the formation of structural fractures is jointly controlled by fault development, fluid overpressure, and lithofacies. The pores and fractures in the oil window of lacustrine shale can store and channel oil and gas. The hybrid OM–clay–carbonate (sulfate) and the pores inside are important through the oil window. Moreover, the development of the pores depends not only on hydrocarbon generation but also on the interaction of hydrocarbons and organic acid dissolution. This finding has important significance in the accumulation of oil and gas in continental shales.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In the past, determination of rock properties using image analysis relied upon petrographic transmitted-light images, but with limited success because of a lack of resolution and restricted computer processing power. A new technique that employs confocal laser scanning microscopy (CLSM) can be considered complementary to laboratory measurements and applicable to several samples, saving time and money and requiring only a limited amount of rock sample for analysis. We have studied several types of rocks with CLSM and fluorescent dye–impregnated thin sections. The two-dimensional scans of each thin section images is an area of 12 mm〈sup〉2〈/sup〉, with a pixel size of 0.198 µm and were used to simulate capillary pressure curves for pore bodies and pore throats. The CLSM technique also enables three-dimensional (3-D) visualization of the rock porosity. The studied rock samples were taken from diverse oil and gas field reservoirs: case A, a conventional sandstone (15.1% porosity, 29.8 md permeability); case B, a tight sandstone (3.7%, 0.02 md); case C, an oolitic carbonate (9.6%, 0.1 md); case D, a rhodolithic algal carbonate (19.8%, 43.7 md); case E, dolomitized carbonate (17%, 21.7 md); and case F, a naturally fractured carbonate (2.4%, 0.6 md). Our results confirm that the CLSM technique can be applied to rocks of contrasting porosity and permeability to obtain computed synthetic capillary pressure curves faster than with conventional measurement methods. The technique quantifies different pore-body and pore-throat sizes and distributions, with the added ability to visualize 3-D porosity and to extract from thin section analysis petrologic properties.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Oil API gravity predictions using published basin modeling source rock (SR) reaction kinetics have displayed poor matches between modeled output and field observations because these kinetic models do not predict increasing API gravities with increasing maturity. Ideally, an SR kinetic model should use at least two liquid components of different densities, which are generated and expelled from the SR such that the API gravities are a consequence of relative mixing. Very few available kinetic models predict APIs with reasonable trends, but those are either not adjustable to calibrate to field observations or do not consider sorption, which is a necessary process when evaluating unconventional resources. Five new kinetics data sets are presented in this paper, each representing a standard SR type, which provide geologically reasonable API gravity trends and ranges. Each kinetic model uses two liquid pseudocomponents and two vapor pseudocomponents. The relative ratios between the pseudocomponents at full kerogen transformation are average ratios available from public and proprietary kinetic data sets. The primary generation follows published activation energies, including minor shifts, which allow peak generation to occur at lower activation energies for the heavier liquid pseudocomponent and at higher energies for the lighter one. This systematic shift of activation energies thus results in a constant change in API gravity as primary generation progresses. Additional in-SR sorption and secondary cracking schemes support the primary generated API gravity trends. The default ranges of API gravity for the new five kinetic models represent observed averages but can be adjusted easily.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Shale gas in the Sichuan Basin and its periphery potentially plays an important role in the world shale gas industry. An understanding of remigration and leakage from continuous shale reservoirs is very important for shale gas exploration, especially in the Sichuan Basin and its periphery. The shale gas accumulation models that relate to remigration and leakage were developed within the Wufeng and Longmaxi black shales in the Jiaoshiba and the Youyang blocks. First, a tectono-sedimentary history of the Wufeng and Longmaxi black shales in the Sichuan Basin and its periphery was developed based on the published literature. The history exhibits a continuous distribution of high-quality Wufeng and Longmaxi black shale, which is the foundation of the shale gas formation. Second, the shale gas remigration–accumulation model in the anticlines was clarified by using data collected from the shale gas fields in Jiaoshiba block. The shale gas model for the Jiaoshiba block was developed on the basis of a continuous shale reservoir distribution, differentiated structural deformation, and a gas self-sealed system. Third, the shale gas fault failure leakage model in the fault blocks and the erosion model in the residual areas were revealed based on the shale reservoir and shale gas content heterogeneity in the Youyang block. These two models were validated by available data including 13 two-dimensional seismic lines and 2 shale gas exploration vertical wells in the Youyang block. Shale gas areas with high gas resource and gas production rates in the anticlines were defined by the remigration–accumulation model. The fault failure leakage model was used to find shale gas with limited commercial potential, whereas commercial shale gas was largely lacking according to the erosion residual model. The study on remigration and leakage from continuous shale reservoirs in the Sichuan Basin and its periphery can be used to better understand and improve the exploration efforts based on resource preservation.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Hyperspectral systems that image drill core can capture detail mineralogical information at the millimeter scale and thus have the potential to enable investigators to characterize shale composition and heterogeneity, complementing the direct chemical and x-ray diffraction analysis of core samples and guiding detailed sampling. This method provides insight into petrophysical and geomechanical properties because these properties are significantly correlated to rock composition. We tested this approach on a continuous long core from the shale sequence of the Horn River Group in the Horn River Basin, British Columbia, sampled at a spacing of 1 m (40 in.) and analyzed for geochemical composition. These data enable the calibration of spectral imagery to rock composition and specifically predict total organic carbon (TOC) and the abundance of SiO〈sub〉2〈/sub〉, Al〈sub〉2〈/sub〉O〈sub〉3〈/sub〉, K〈sub〉2〈/sub〉O, and CaO. We then imaged nine samples from the Woodford Shale from the Permian Basin, Texas, for a blind test to assess the predictive models. The models were then used to predict TOC and geochemical data over detailed imagery of 300 m (984 ft) of Horn River Group shale core and portray their spatial variability downhole as images and profiles. In its simplest form, hyperspectral imagery can be enhanced to highlight fabric in shale core that otherwise is difficult to visualize because of low brightness. In addition, we show that spectral imagery of shale can also be processed to either convey mineralogical (quartz, clay, and carbonate) or geochemical information. The resulting views can readily be used to guide the selection of samples and may provide tools for scaling reservoir properties from individual plugs to reservoir volumes.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Considerable attention has been directed to the Devonian Horn River Formation in western Canada with respect to geochemical evaluation of gas-generation and storage potential. Because organic geochemical analyses are not always useful for characterizing the type and amount of original organic matter, we surmise the original kerogen type and original hydrogen index (HIo) and subsequently estimate a reliable original total organic carbon (TOCo) based on a combination of inorganic and organic geochemical data. Productivity (SiO〈sub〉2〈/sub〉 and Ba) and terrestrial input (Al〈sub〉2〈/sub〉O〈sub〉3〈/sub〉, Hf, Nb, and Zr) proxies are used to estimate original kerogen types, which suggest that the Evie and Muskwa Members formed under conditions of high productivity and minor terrestrial input. These members also formed under reducing conditions, as indicated by the redox proxies (Mo, U, and Th/U). Under such conditions, primarily type II kerogen was preserved.By considering the fraction of biogenic silica, the estimated HIo values (400–500 mg hydrocarbon/g total organic carbon [TOC]) for the middle Otter Park Member are lower than that for Evie and Muskwa Members and higher than the upper and lower Otter Park Member. The stronger correlation between TOCo and trace elements suggests that HIo is useful for reconstructing the coherent variation in TOCo. Based on the original kerogen type and TOCo, the gas-generation and storage potentials of the Evie, middle Otter Park, and Muskwa Members are higher than those of other members. The source-rock potential is excellent for the Evie Member with an approximately 75% difference between TOCo and measured present-day TOC.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The global Precambrian–Cambrian system includes an important series of hydrocarbon-bearing strata. However, because rocks of this age are typically deeply buried, few petroleum exploration breakthroughs have been made, and the presence of source rocks remains somewhat controversial. Recently, commercial condensate and gas were discovered from the deep (∼6900 m [∼22,600 ft]) Zhongshen 1C (ZS1C) exploratory well drilled in the Tazhong uplift of the Tarim Basin, China, leading to renewed interest in the development of Cambrian source rocks in the basin. On the basis of outcrop reconnaissance and sample testing from around the Tarim Basin, we show that a set of high-quality source rocks were developed within the lower Cambrian Yuertusi Formation (Є〈sub〉1〈/sub〉y), at the base of the lower Cambrian. These rocks are black shales and typically have a total organic carbon content between 2% and 6% but extending as high as 17% in selected regions. This marine sequence is 10–15 m (33–49 ft) thick in some outcrops along the margins of the basin. Seismic data indicate that these high-quality source rocks may cover an area as large as 260,000 km〈sup〉2〈/sup〉 (100,000 mi〈sup〉2〈/sup〉). Their main organic parent material was benthic multicellular algae. On the basis of high-temperature thermal simulations conducted on these source rocks, we show that the gas composition and carbon isotopes from the ZS1C well are similar to the products generated at a thermal evolution stage corresponding to a vitrinite reflectance of between 2.2% and 2.5%. Late-stage natural gas accumulated within these rocks over time. The δ〈sup〉34〈/sup〉S correlation of organic sulfur compounds in the condensate with Cambrian sulfates provides further evidence for a Є〈sub〉1〈/sub〉y source rock origin of the ZS1C condensate and gas. The Cambrian dolomites in association with a salt seal exhibit favorable geological conditions for large-scale hydrocarbon accumulation. A new set of deep exploration strata can, therefore, be developed, guiding future deep Cambrian hydrocarbon exploration in the Tarim Basin.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The effects of reservoir heterogeneity on the development of submarine channel fields are still poorly understood because of lack of direct evidence for fluid flow. This study uses integrated well logs and three-dimensional seismic data from the Niger Delta Basin to characterize the previously undocumented spatial distribution of shale units and permeability contrasts within a submarine channel system. Combining these data with four-dimensional (4-D) seismic data facilitates the exploration of the controls of reservoir heterogeneity on fluid flow during development. The results show that the studied submarine channel system consists of multiple vertically stacked channel complex sets (CCSs) from CCS1 (oldest) to CCS5 (youngest), which are separated from each other by continuous shale barriers. The CCS2–CCS4, which are located in the stratigraphic middle of the channel system, are the main development layers because of their higher permeabilities and lower permeability contrasts. The 4-D seismic responses validate that the presence of shale barriers between vertically adjacent CCSs can hinder the flow of fluids between CCSs. Fluid flow between vertically adjacent CCSs barely occurs except in localized erosional locations where the sand fills of different CCSs are vertically connected. Each CCS consists of multiple individual channels, which can be separated by inclined shale baffles if they laterally migrate in one direction. As the 4-D seismic responses demonstrate, such inclined shale baffles can hinder fluid flow between adjacent individual channels and help to form multiple narrow flow paths in map view. The absence of inclined shale baffles also produces prominent permeability contrasts within each CCS, which are characterized by relatively high–permeability zones that are parallel to the channel axis. Comparison of this permeability distribution and the 4-D seismic responses shows that injected water preferentially sweeps along relatively high–permeability zones, which can help to form single wide flow paths with higher sweep efficiency or single narrow flow paths with lower sweep efficiency.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Source-to-sink analyses show that northern Gulf of Mexico (GOM) Wilcox Group siliciclastic deep-water systems are linked to transport of sediments from the Laramide tectonic belt into the deep basin. Less is known, however, about southern GOM sedimentation. New drilling and discoveries in the Mexican deep water have generated considerable interest since the opening of Mexico to international exploration. To investigate Paleogene deposition in Mexico’s offshore areas, a three-phased approach was employed: (1) seismic mapping of deep-water depocenters, (2) regional stratigraphic analysis of potential basin entry points, and (3) prediction of submarine-fan dimensions using empirical scaling relationships. Isochore and structural mapping of the Wilcox depocenters used available well and seismic data. Potential basin entry points were identified by evaluation of Wilcox fluvial–deltaic systems and tectonic elements. Empirical scaling relationships previously established between fluvial and deep-water segments provide first-order predictions of submarine-fan dimensions.Paleogene Wilcox source-to-sink systems of the greater GOM basin change north to south as a function of varied tectonics and sedimentary accommodation. The United States sector was a passive margin: continental-scale drainage systems fed a broad, gently dipping shelf. By contrast, the southern GOM basin was a tectonically active margin: smaller-scale fluvial systems sourced from the Hidalgoan uplands flowed directly into foreland basins located on the slope. Results presented here indicate that several systems rimming the southern GOM were able to effectively transfer sediment from the mountain belt into the basin. Regional observations and semiquantitative predictions of fan dimensions provide a context for future detailed work based on new well and seismic information.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Sequence stratigraphy based on wire-line logs, cores, and outcrops is entering its fourth decade of mainstream usage in industry and academia. The technique has proved to be an invaluable tool for improving stratigraphic analyses in both clastic and carbonate settings. Here we present a simple quantitative technique to support sequence stratigraphic interpretations in clastic shallow marine systems. The technique uses two pieces of data that are readily available from every subsurface field or outcrop study: (1) parasequence thickness (T) and (2) parasequence sandstone fraction (SF). The key assumptions are that parasequence thickness can be used as a proxy for accommodation at the time of deposition and parasequence sandstone fraction can be used as a proxy for sediment supply. This means that quantitative proxies for rates of accommodation development and sediment supply can be acquired from wire-line logs, cores, and outcrop data. Vertical trends in parasequence thickness divided by sandstone fraction (T/SF) approximate trends expected in systems tracts for changes in ratios of rate of accommodation development to rate of sediment supply. The technique, termed “TSF analysis,” can also be applied at lower-order sequence and composite sequence scales. It provides a quantitative and objective methodology for determining rank and order of sequence stratigraphic surfaces and units. Absolute T/SF values can be used to determine shoreline, stacked shoreline, and shelf-margin trajectories. Four case studies are presented, which demonstrate the robustness of the technique across a range of different data sets. Implications and potential future applications of TSF analyses are discussed.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Eagle Ford Formation has attracted considerable industry attention as a self-sourced unconventional shale reservoir. The productive interval in the Eagle Ford Formation is the transgressive systems tract, which contains parasequences whose lithologic content varies upward with increasing proportions of limestones. Optimum success in both exploration and production depends on the adequate characterization of fracture systems as a function of lithology. The outcrops present along US Highway 90 in Val Verde and Terrell Counties, Texas, provide considerable insight into the regional natural fracture system of the Eagle Ford Formation. Fracture-orientation analysis reveals two sets of conjugate hybrid shear fractures and two sets of regional fractures. Abutting relationships suggest that hybrid shear fractures formed first, followed by the thoroughgoing northeast-striking fracture set, and finally by a northwest-striking set, which tends to be confined to individual mechanical units. The orientation of these fractures suggests that they formed during post-Laramide stress relaxation and progressive exhumation. Spacing-frequency distribution analysis of the fracture population reveals a mature hypersaturated fracture system that likely formed at depth by overburden load and/or fluid pressure near maximum burial. Our results indicate that the Eagle Ford Formation displays a well-developed fracture network regionally distributed in the Val Verde Basin, and likely present in the productive Eagle Ford play. These observations provide evidence for pathways and vertical connectivity for potential fluid pathways throughout the Eagle Ford Formation.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉We document a novel approach to balanced three-dimensional structural restoration based on an adaptation of the GeoChron model. Conventionally, the GeoChron model defines a transformation of a geological model to a flattened space (U-V-T), with paleogeographic coordinates defined by the horizontal axes (U-V) and geologic time on the vertical axis (T). In our new balanced structural restoration scheme, the complete stratigraphy is restored using a transformation constrained only by the datum horizon. Scaling the vertical “T” axis to depth in a manner that preserves volume or layer thickness results in a geometric restoration that approximately minimizes strain globally. This restoration provides a geometrically plausible representation of the geologic structure at the time when the datum horizon was deposited. Restoration is independent of mechanical rock properties and is thus most applicable to regions in which mechanical rock properties are approximately homogeneous. Restoration kinematics may be constrained by growth strata if present.We validate the method with kinematic forward models and a laboratory sandbox model and apply it to two natural examples to demonstrate its capabilities for model validation and palinspastic restoration.We identify four criteria for assessing the validity of a structural model using the results of restoration: (1) anomalous fault throw, (2) timing of fault activity, (3) fault compliance, and (4) restoration strain. Analysis of the sandbox results and limitations of volume conservation derived from uncertainties in compaction states suggest accuracy of the method to be in the 5%–20% range.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The pore structure of shale has a significant effect on hydrocarbon migration and the long-term gas supply of shale gas wells. The present study investigates the spontaneous imbibition characteristics to evaluate the pore connectivity and wettability of marine Longmaxi shale samples from the southeastern Chongqing area and continental Yanchang shale samples from the Ordos Basin. The pore-size distribution obtained from N〈sub〉2〈/sub〉 adsorption and mercury intrusion porosimetry, field emission–scanning electron microscopy, and focused ion beam–scanning electron microscopy photos are used to interpret the imbibition behaviors. Our results show that the difference in dominant pore type between marine and continental samples, which is dominated by thermal maturity, controls on their imbibition behaviors as well as their wettability. Organic matter (OM) pores within Yanchang samples are poorly developed because of their low thermal maturity, and a large amount of water-wet inorganic pores are preserved in these samples because of relatively weak compaction. Oil-wet OM pores are well developed in Longmaxi samples with higher thermal maturity, and inorganic pores have been largely eliminated because of strong compaction. The low pore connectivity to water for both the Longmaxi and Yanchang samples is indicated by the low water imbibition slopes. Furthermore, the more oil-wet property of the Longmaxi samples and more water-wet characteristics of the Yanchang samples are obtained by comparing the directional water/oil imbibition slopes. In addition, the positive meaning of quartz in the protection of pore spaces is found in both the Longmaxi and the Yanchang samples used in this study.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Multiple natural gas fields have been discovered in the Baiyun depression and the adjoining Panyu lower uplift in the Pearl River Mouth Basin, northern South China Sea. The natural gases are associated with condensate and are characterized by relatively heavy carbon isotopes, with methane and ethane δ〈sup〉13〈/sup〉C values ranging from –44.2‰ to –33.6‰ and –30.0‰ to –25.4‰, respectively. Nearly all methane and ethane are derived from oil-prone type II kerogen in the Wenchang Formation source rock, whereas the heavy hydrocarbon gases (propane, butanes, and pentanes) are derived from both the Wenchang and Enping (type III kerogen) Formations, based on an integrated comparison of carbon isotopic compositions of the natural gases, typical type I/II and type III kerogen-derived gases, and the Enping and Wenchang kerogens. The gases from the eastern parts of the Baiyun depression and the Panyu lower uplift mainly originate from secondary oil cracking and primary kerogen cracking, respectively. The gases from the northern slope of the Baiyun depression are a mixture of oil-cracking and kerogen-cracking gases. Both oil-cracking and kerogen-cracking gases were mainly generated from the Wenchang Formation source rock in the maturity range of 1.5%–2.5% vitrinite reflectance, with a corresponding present-day depth range of 5400–6500 m (17,700–21,300 ft). The apparent contribution of the Wenchang Formation to the discovered gas accumulations demonstrates that it is the most important source rock in the area, instead of the Enping Formation. The search for more gas derived from oil cracking will be the next natural gas exploration direction in the Baiyun depression.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The northern Appalachian Basin depocenter of Pennsylvania represents one of the most economically important hydrocarbon-producing areas in the United States, yet the thermal conditions that promoted hydrocarbon formation within the basin are only marginally constrained. The prolific coal, oil, and natural gas fields of Pennsylvania are the direct result of thermal maturation of once deeply buried organic-rich sediment. Understanding how, why, and where thermal maturation occurred in the Appalachian Basin requires high-quality heat flow and thermal conductivity measurements, as well as paleotemperature estimates and basin modeling. To improve the understanding of heat flow, we present, to our knowledge, the first direct measurements of (1) thermal conductivity on Devonian core samples and (2) equilibrium temperature versus depth logs for the northern Appalachian Basin depocenter. Results from three well sites demonstrate that heat flow is conductive and nearly uniform, averaging 34 ± 2.5 mW/m〈sup〉2〈/sup〉, with an average thermal gradient of 29 ± 4°C/km. The new heat-flow measurements are significantly lower (30%–50% less) than previously published estimates that used nonequilibrium bottomhole temperature values and empirically derived thermal conductivity estimates. Our analysis indicates that previous studies correctly estimated the regional thermal gradient using bottomhole temperatures but overestimated heat flow in this region by as much as 50% because of inaccurate extrapolation of thermal conductivity. The results highlight the importance of directly measuring thermal conductivity to accurately quantify heat flow in deep sedimentary basins. Ultimately, additional paleotemperature data are necessary to improve our understanding of Appalachian Basin thermal evolution.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Shale samples of the Marcellus Shale from a well drilled in northeastern Pennsylvania were used to study diagenetic effects on the mineral and organic matter and their impact on petrophysical response. We analyzed an interval of high gamma ray and anomalously low electrical resistivity from a high thermal maturity (mean maximum vitrinite reflectance 〉 4%) part of the shale‐gas play. A suite of microanalytical techniques was used to study features of the shale down to the nanoscale and assess the level of thermal alteration of the mineral and organic phases.The samples are organic rich, with total organic carbon contents of 3–7 wt. %; the vast majority of the organic matter was identified as highly porous pyrobitumen. Matrix porosity is also present, especially within the clay aggregates and at the interface between rigid clasts and clay minerals.Mineral- and organic-based thermal maturity indices suggest that during burial the sediment had been exposed to temperatures as high as 285°C (545°F). Under these conditions, the residual, migrated organic matter assumed a partially crystalline habit as confirmed by the identification of turbostratic structures via electron microscopy imaging. Experimental dielectric measurements on organic matter–rich samples confirm that the anomalous electrical properties observed in the wire-line logs can be ascribed to the presence of an electrically conductive interconnected network of partially graphitized organic matter. The preservation of porosity suggests that this organic network can contribute not only to the electrical properties but also to the gas flow properties within the Marcellus Shale.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Understanding natural fracture networks in the subsurface is highly challenging, as direct one-dimensional borehole data are unable to reflect their spatial complexity, and three-dimensional seismic data are limited in spatial resolution to resolve individual meter-scale fractures.Here, we present a prototype workflow for automated fracture detection along horizontal scan lines using terrestrial light detection and ranging (t-LIDAR). Data are derived from a kilometer-scale Pennsylvanian (locally upper Carboniferous) reservoir outcrop analog in the Lower Saxony Basin, northwestern Germany. The workflow allows the t-LIDAR data to be integrated into conventional reservoir-modeling software for characterizing natural fracture networks with regard to orientation and spatial distribution. The analysis outlines the lateral reorientation of fractures from a west–southwest/east–northeast strike, near a normal fault with approximately 600 m (∼1970 ft) displacement, toward an east–west strike away from the fault. Fracture corridors, 10–20 m (33–66 ft) wide, are present in unfaulted rocks with an average fracture density of 3.4–3.9 m〈sup〉−1〈/sup〉 (11.2–12.8 ft〈sup〉−1〈/sup〉). A reservoir-scale digital outcrop model was constructed as a basis for data integration. The fracture detection and analysis serve as input for a stochastically modeled discrete fracture network, demonstrating the transferability of the derived data into standard hydrocarbon exploration-and-production-industry approaches.The presented t-LIDAR workflow provides a powerful tool for quantitative spatial analysis of outcrop analogs, in terms of natural fracture network characterization, and enriches classical outcrop investigation techniques. This study may contribute to a better application of outcrop analog data to naturally fractured reservoirs in the subsurface, reducing uncertainties in the characterization of this reservoir type at depth.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Using cores, well logs, and other borehole data, the results of this study show that the shallow-water lacustrine delta has its own unique depositional characteristics of the third member of Oligocene Dongying Formation (Ed〈sub〉3〈/sub〉) in the Baxian sag, Bohai Bay Basin, eastern China. During the Ed〈sub〉3〈/sub〉 stage, the rift–thermal basin subsidence transition stage, the paleoslope was divided into multilevel slopes by faults along the Wen’an slope with slope angles from approximately 0.19° to 2.02°. The paleogeographic conditions, low-discharge channel, and low accommodation controlled the sedimentary characteristics. The distributions of the shallow-water delta system were controlled by multilevel flexure slopes. The delta plain was distributed on the first- and second-level slope belts, and the delta front was distributed on the third-level slope belt. The high-sinuosity fluvial channel of the delta plain was the dominant facies in the whole shallow-water delta. Most sand was deposited in these channels along the second-level slope belt. Therefore, not enough sand was present to be transported into the lake (shallow water) to form mouth bars in the delta front. Therefore, mouth bars of the shallow-water delta front were few, and the sand beds were thin. Additionally, no more sand was available to be supplied right along to deep lake, the lacustrine basin was small, and there was insufficient accommodation and sand to develop a subaqueous fan in the delta front.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Three-dimensional reservoir modeling is an important aspect to determine the heterogeneity of organic-rich shale reservoirs, an area of study that continues to be explored and refined. A large proportion of data acquired from horizontal wells causes issues in the structural and property modeling for shale reservoirs. Since horizontal wells are designed to drill into a specific, narrow zone, their horizontal section tends to parallel or nearly parallel formation surfaces. As a result, formation surfaces have a much more complex spatial location relationship with horizontal wellbores than with vertical wellbores. The existing algorithms are not good at addressing this issue during structural modeling. The major problem of using horizontal well data in property modeling is the biased data set because their horizontal section tends to stay within a narrow zone. The property distribution feature estimated from this biased data set, as a significant, default input of geostatistical simulation algorithms, causes the constructed property models to deviate away from the real case in the subsurface. A method to infer more formation tops in pseudovertical wells according to a series of assumptions was developed to provide more constraint points for structural modeling within the areas of the horizontal well section. To use the biased database from horizontal wells, distribution function and trend model methods were developed for continuous property modeling, and percentage and probability trend models were developed for discrete property modeling. The Longmaxi–Wufeng shale in the Fuling gas field of Sichuan Basin was used as an example to express and verify these methods.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉An integrated approach to detect new areas of potential interest associated with stratigraphic traps in mature basins is presented. The study was carried out in the Middle Magdalena Valley basin, Colombia. The workflow integrates outcrop and subsurface interpretations of facies, activity of faults, and distribution of depocenters and paleocurrents and makes use of them to construct a three-dimensional exploration-scale geocellular facies model of the basin. The outcrop and well log sedimentological analysis distinguished facies associations of alluvial fan, overbank, floodplain, and channel fill, the last one constituting the reservoir rock. The seismic analysis showed that tectonic activity was coeval with the deposition of the productive units in the basin and that the activity ended earlier (before the middle Miocene) along the western margin than along the eastern margin. Paleogeographic reconstructions depict transverse and longitudinal fluvial systems, alluvial fans adjacent to the active basin margins, and floodplain facies dominating the structural highs and the southwestern depositional limit. These reconstructions provided statistical data (lateral variograms) to construct the model. The exploration-scale facies model depicts the complete structure of the basin in three dimensions and the gross distribution of the reservoir and seal rocks. The predictive capability of the model was evaluated positively, and the model was employed to detect zones of high channel fill facies probability that form bodies that are isolated or that terminate upward in pinchouts or are truncated by a fault. Our approach can prove helpful in improving general exploration workflows in similar settings.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Knowledge of in situ stress distribution is fundamental for coalbed methane production; however, it is poorly understood in the eastern Yunnan region, South China. In this study, the horizontal maximum (〈span〉S〈/span〉〈sub〉〈span〉Hmax〈/span〉〈/sub〉) and minimum (〈span〉S〈/span〉〈sub〉〈span〉hmin〈/span〉〈/sub〉) principal stress and vertical stress (〈span〉S〈/span〉〈sub〉〈span〉v〈/span〉〈/sub〉) were systematically analyzed for the first time. The results indicated that the magnitudes of 〈span〉S〈/span〉〈sub〉〈span〉Hmax〈/span〉〈/sub〉, 〈span〉S〈/span〉〈sub〉〈span〉hmin〈/span〉〈/sub〉, and 〈span〉S〈/span〉〈sub〉〈span〉v〈/span〉〈/sub〉 showed positive correlations with burial depth. In general, three types of in situ stress fields were determined: (1) 〈span〉S〈/span〉〈sub〉〈span〉Hmax〈/span〉〈/sub〉 〉 〈span〉S〈/span〉〈sub〉〈span〉v〈/span〉〈/sub〉 〉 〈span〉S〈/span〉〈sub〉〈span〉hmin〈/span〉〈/sub〉 in shallow layers with burial depths less than approximately 600 m (∼1970 ft) below ground level (bgl), indicating a dominant strike-slip faulting stress regime; (2) in medium layers approximately 600–800 m (∼1970–2625 ft) bgl, the in situ stress state followed multiple relationships, suggesting that the in situ stress regime was transformed; and (3) 〈span〉S〈/span〉〈sub〉〈span〉v〈/span〉〈/sub〉 〉 〈span〉S〈/span〉〈sub〉〈span〉Hmax〈/span〉〈/sub〉 〉〈span〉S〈/span〉〈sub〉〈span〉hmin〈/span〉〈/sub〉 in deep layers with burial depths greater than approximately 800 m (∼2625 ft) bgl, indicating a dominant normal faulting stress regime. Coal permeabilities obtained from injection–falloff well tests showed that they were widely distributed, and no obvious relationships were found between coal permeability and effective in situ stress magnitude. In the study area, the development and orientation of previously generated natural fractures combined with the present-day in situ stress distribution controlled the permeability in coal reservoirs. Differential stress and presence of natural fractures significantly affected the geometry and pattern of hydraulic fractures. In addition, in the eastern Yunnan region, locations with relatively deep depths in vertical wells and approximately west–northwest/east–southeast-trending horizontal wells suffered high potential of borehole instability because of the high differential stress.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In carbonate rock reservoirs, spatial distribution models and elastic properties are complex because of diagenetic processes and mineralogical composition, which together directly interfere with variations in pore shape and interconnectivity. The main objective of this paper is to propose a workflow to aid in three-dimensional quantitative carbonate reservoir characterization of the Quissamã Formation (Macaé Group) in the Pampo field of the Campos Basin, offshore Brazil. Model-based seismic inversion, sequential Gaussian simulation with cokriging for porosity modeling, and truncated Gaussian simulation with trend for facies modeling were used to characterize the carbonate reservoirs. Our results show that the carbonate platform is located between the upper Aptian and lower Albian seismic surfaces. Interpretation of a new surface, called the intra-Albian, was possible via acoustic-impedance (AI) analysis. Our workflow facilitated identification of low AI, high porosity, and best facies areas in structural highs where the most productive wells have been drilled. Facies modeling suggests that intercalation of facies with high and low porosities is connected to shallowing-upward cycles. Finally, several debris facies with low AI and high porosities were identified in an area that could be targeted for new exploration.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Thick marine shales occur in the Wufeng Formation and Longmaxi Formation of Sanquan Town of the Nanchuan District, Chongqing Municipality, which is located on the southeast margin of the Sichuan Basin. However, few details of the characteristics of the Wufeng–Longmaxi shales in this area have been reported. In this study, a well approximately 100 m (∼328 ft) deep was drilled. A high-quality shale (total organic carbon [TOC] 〉2.0 wt. %, clay 〈40%) interval that was approximately 24 m (∼79 ft) thick with an average TOC value of 3.0 wt. % mainly occurs in the Ordovician Wufeng Formation (Katian and Hirnantian) and base of the Silurian Longmaxi Formation (Rhuddanian). Shales with higher TOC values commonly have a higher porosity and specific surface area. Tectonic movements may also have been very important factors that influenced the petrophysical properties of the shales. For example, a detachment layer that resulted from complex tectonic movements is extensive in the Wufeng Formation. The cracks and microcracks in the detachment layer can result in good pore connectivity. Consequently, the detachment layer can be an effective migration pathway. The Longmaxi–Wufeng shales of Sanquan Town are also compared with those of the famous Jiaoye 1 well in the Jiaoshiba shale gas field in the eastern Sichuan Basin. Although the shales in Sanquan Town have considerable shale gas generation potential, the shale gas resource potential in Sanquan Town is probably poor because the escape of shale gas may be accelerated by the detachment layer in destroyed anticlines and synclines.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The ability to accurately predict the probability of fluid migration from depth through existing wells based on known well properties, such as age and depth, would be enormously helpful in understanding how migration pathways develop and the identification of potential migration without extensive field tests. The presence of fluid pathways is an important environmental issue because such pathways allow gas, either naturally occurring methane or sequestered CO〈sub〉2〈/sub〉, to move into the atmosphere. In this paper, we explore the ability of various predictive models to forecast gas migration at existing wells in Alberta, Canada, based upon the characteristics of existing deep wells. Alberta was selected as a case study because of the availability of data in an area that has required wells to be tested for pathway development after rig release since 1995. Wells that do not demonstrate pathway development require no further testing until the well is abandoned. We show that accurately predicting fluid migration requires detailed information on well construction, production, and fluid properties, and even then, the models considered in this study misclassify a large number of wells. This suggests other factors may contribute to pathway formation. Of the models investigated, random forests provide the best results on this data set, correctly identifying 78% of the wells used.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉This study aims to decipher the groundwater status of the parts of Tigray area, Ethiopia using an integrated methodology of remote sensing and geographic information systems (GIS). Digitized vector maps of the study area, that is, geology, land use and/or cover, and drainage, were generated and converted to raster data. The theme weight and class weights were assigned to the raster maps of the respective parameters. Weight age to the layers was assigned using an analytical hierarchy process and further overlay analysis was carried out in the ArcGIS environment to decipher the groundwater resources of the study area.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The present contribution aims for a characterization of microstructure and pore-space distribution of upper Visean Rudov beds, considered the main source rock for conventional oil deposits in the Ukrainian Dneiper–Donets Basin and a prospect for unconventional hydrocarbon production in recent years. Broad ion beam–scanning electron microscopy (SEM) mapping revealed a remarkably heterogeneous microstructure controlled by diagenetic precipitates (Fe/Mg carbonates, albite). Formation of these precipitates is likely triggered by organic matter decomposition and represents an important influencing factor for overall porosity and permeability. Furthermore, shale diagenesis also influences mechanical properties, as suggested by nanoindentation tests. The SEM-visible organic matter porosity is restricted to solid bitumen; although pores less than 2–3 nm in vitrinites of overmature samples are indicated by focused ion beam–SEM results, they cannot be resolved clearly by this method. Pore generation in solid bitumen that likely formed in situ in primary amorphous organic matter already starts at the early oil window in samples from the basinal oil-prone organofacies, whereas most porous solid bitumen at peak oil maturity was interpreted as relicts of primary oil migration, representing an earlier oil phase that predominantly accumulated in quartz-rich layers and became nanoporous during secondary cracking. In the terrestrially dominated transitional to marginal organofacies, pore generation in pyrobitumen resulting from gas generation occurs significantly later and is less intense. Formation of authigenic clay and carbonate minerals within pyrobitumen is likely related to organic acids formed during bitumen decomposition and implies the presence of an aqueous phase even in pores that are apparently filled exclusively with solid bitumen.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Ahdeb oil field is located in the Mesopotamian Basin of central Iraq within a northwest–southeast-trending anticline. Seven oil-bearing layers exist in the eastern area in the field, but there is only one oil-bearing layer in the western area. This study reveals that the reservoir filling process resulted from the difference in the elements in the petroleum system, the oil generation and migration process, and the formation of the structural trap. Most oils in the field, with pristane/phytane 〈 1 and a high relative abundance of hopanes exceeding C〈sub〉30〈/sub〉, were generated from the Upper Jurassic–Lower Cretaceous Chia Gara Formation, whereas some oils were generated from the Lower Cretaceous Ratawi and Zubair Formations. The mid-Upper Cretaceous reservoirs in the field are composed of lime grainstones, packstones, and wackestones.The main oil accumulation occurred during the Maastrichtian, coinciding with peak oil generation from the Chia Gara Formation with a 50% transformation ratio from organic matter to oil. The reservoirs of the eastern structural trap in the field were filled with large amounts of medium to heavy oils. After the formation of two structural traps in the western area in the mid-Miocene, oils pre-existing in the second layer of the Khasib Formation in the east began migrating toward the structural traps in the west during the late Miocene, as verified by relatively higher 1-/4-methylcarbazole and 1,8-/2,7-dimethycarbazole ratios of oils in the west than that in the east and residual solid bitumen in the east. The strike-slip fault might also have restricted oil or gas migration during the Miocene, limiting oil accumulation in the west.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The success of hydraulic fracturing and increasing use of basin-modeling packages drive the need to understand the effects of hydrocarbon (HC) generation on the mechanical properties of source rocks. A better understanding of relationships among geological, geochemical, and geomechanical parameters can potentially reduce the uncertainties associated with conventional and unconventional prospect evaluation.We present a simulation of microcrack growth based on a three-dimensional source-rock system. Upon thermal maturation, the kerogen transforms into lighter products, most of which are HCs. The generated products exert excessive pore pressure to the system resulting from the effect of volume expansion; this pressure is released through the expansion of pore space and formation of microcracks. Using linear elasticity and linear elastic fracture mechanics, our model calculates microcrack sizes (surface areas, lengths, apertures, and volumes) and the amount of overpressure throughout the maturation process. We validated this model with experimental data from 〈a href="https://pubs.geoscienceworld.org/aapgbull#b20"〉Kobchenko et al. (2011)〈/a〉, and performed sensitivity analysis for both laboratory and geological settings. Much larger microcracks are generated in laboratory settings compared to the subsurface because of the lack of overburden, resulting in secondary porosity over 100 times larger than the original organic porosity and crack lengths obtaining millimeter scale. In contrast, microcracks are much smaller in geological settings because of the presence of significant overburden and stiffer rock frames: the crack apertures are in the submicron regime with a crack length ranging from 100 to 300 μm. The formation of microcracks connects isolated microscale HC pockets, providing pathways for primary migration.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Studies of lacustrine carbonate rocks in continental rifts have received huge interest in recent years because of their great economic value in the South Atlantic. However, most existing facies and tectonosedimentary models for carbonate platforms are based on marine carbonate systems, whereas models for nonmarine systems are scarce. The main aim of this paper is to establish such models and to further our understanding of the hydrocarbon-bearing late synrift Lower Cretaceous carbonate successions of the Campos Basin, Brazil. This paper is based on a proximal to distal industrial data set of three-dimensional (3-D) seismic, cores, and well logs from the Coqueiros Formation (Coquina), southern Campos Basin. The dominant carbonate facies in the Coqueiros Formation are mollusk-rich grainstones, rudstones, and floatstones, which form the main reservoir facies. The 3-D seismic interpretations show an oblique extensional rift system, characterized by a series of grabens, half grabens, accommodation zones, and horsts oriented northeast–southwest to north–northeast-south–southwest. Three tectonic domains are recognized based on structural style, stretching factors, and subsidence rates as well as facies and different types of lacustrine carbonate platforms. Proximal rift margin areas are characterized by a series of half grabens with footwall and hanging-wall dip slopes of shallow lacustrine carbonates and fluviodeltaic mixed carbonate and siliciclastic deposits in marginal, hanging-wall basins. Central areas are carbonate rich with platforms established over horst blocks surrounded by deeper-water carbonate facies. Distal areas have the highest amount of stretching and subsidence and accumulate the thickest carbonate successions over a template of buried horsts and grabens. The entire carbonate succession underlies a thick layer of Aptian salt, which forms the seal to this prolific hydrocarbon system.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The uppermost Middle Triassic Leikoupo Formation in the western Sichuan Basin of China has recently been shown to host as much as 5.3 tcf (1.5 × 10〈sup〉12〈/sup〉 m〈sup〉3〈/sup〉) of natural gas resources. The reservoir rocks, composed mainly of microbially derived dolomudstone (e.g., thrombolites and stromatolites), are characterized by low porosity (〈8%) and permeability (〈0.001 to 10 md). The limestone is commonly tight and not of reservoir quality because of abundant meteoric calcite cementation, whereas the dolostone has various types of pores dominated by solution-enlarged pores and vugs, microbial framework pores, and micropores. Breccias are well developed in places, probably because of dissolution of underlying evaporites (e.g., anhydrite) by an influx of low-salinity fluids (e.g., freshwater and seawater) during an early burial stage. Early dolomitization created micropores in the dolomudstone, and subsequent diagenetic events were dominated by calcite, dolomite, quartz cementation, pyrite replacement, compaction, fracturing, and development of stylolites. Localized hydrothermal activity has been evidenced by high homogenization temperatures (∼160°C–200°C) obtained from fluid inclusions in fracture-filling cements. Bacterial sulfate reduction probably resulted in H〈sub〉2〈/sub〉S generation, pyrite precipitation, and solution-enlarged pore and vug formation, whereas part of the current H〈sub〉2〈/sub〉S in these reservoirs may have been sourced from thermochemical sulfate reduction or an underlying formation (e.g., the Feixiangguan Formation). Development of microfractures and associated micropores was probably the final diagenetic event, which improved pore interconnectivity. This study confirms the effect of diagenesis on the development of a microbial dolomudstone reservoir, which may be applicable to other similar microbial carbonate reservoirs elsewhere, for example, Middle Triassic sections of the Tethys region and offshore Brazil.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Clinoforms, the basic large-scale architectural form within which sediments are stored and eventually fed down depositional dip in clastic wedges, exist in many shapes and sizes. Understanding how they form, evolve, and degrade is critical to understanding how transport mechanisms affect the shelf margin and sediment partitioning and distribution, and their implications on the presence of a working petroleum system. The Neogene stratigraphic succession of the Taranaki Basin in New Zealand contains clinoform packages that display a variety of architectures well imaged on seismic data. Quantitative characterization of this interval was performed to unravel the processes by which clinoforms evolve under the influence of tectonic- and isostatic-driven subsidence, sea-level change, and sediment supply fluctuations. Nine different clinoform packages were identified on the basis of changes in their seismic stratigraphic characteristics. Two-dimensional stratigraphic forward modeling was used to conduct a sensitivity analysis and determine the relative importance of different geologic controls on their genesis. Our results show that during the early to late Pliocene, clinoform architectures were influenced by the opening of a back-arc rifting structure in the Taranaki Basin (northern graben), which controlled sediment redistribution and partitioning. At the same time, a drop in global sea level allowed sediment bypass to distal parts of the basin. During the late Pliocene, changes in the Australian–Pacific subduction zone forced rapid uplifting of the Southern Alps, generating a significant increase in sediment supply. Model simulations suggest that clinoform architectures during the late Pliocene were controlled by this increase in sediment supply and associated loading.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉A detailed, rock-based investigation of three Upper Cretaceous Eagle Ford Group cores situated behind, at, and downdip of the Lower Cretaceous Stuart City paleoreef-shelf margin in south Texas was conducted to understand stratigraphic, sedimentological, and geochemical relationships across this buried shelf margin. An understanding of how the Eagle Ford Group lithofacies vary across the paleoreef-shelf margin is currently lacking. We therefore examined a dip section of three cores across the antecedent shelf margin and delineated seven Eagle Ford lithofacies: (1) massive argillaceous mudstone, (2) massive to laminated foraminiferal lime wackestone, (3) radiolarian and foraminiferal dolomitic to lime packstone, (4) massive to bioturbated skeletal lime wackestone, (5) laminated foraminiferal lime packstone, (6) laminated inoceramid and foraminiferal lime grainstone, and (7) massive to bioturbated claystone. A basinward decrease in calcite from 60% to 48% is accompanied by an increase in clay minerals from 12% to 20%. The low-relief raised rim of the older, buried Stuart City paleoshelf margin may have acted as a barrier, dividing the Eagle Ford Group into two sedimentological systems: (1) a restricted drowned shelf to the north, and (2) an open-marine basinal setting to the south. The lower to upper Cenomanian Eagle Ford strata on the drowned shelf are cyclic and enriched in molybdenum, suggesting anoxic to euxinic water masses. The anoxic, open-marine, basinward strata are less cyclical and have a lower molybdenum (compared with the drowned shelf) content. Ash beds and gravity-flow deposits are rare south of the margin. A depositional model was constructed of the lower and upper Eagle Ford formations.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Predicting the lateral distribution of petroleum play elements (reservoirs, source rocks, and seals) requires basic understanding of regional basin evolution and depositional history. In remote areas where little data are available or where the basins have undergone episodes of tectonic deformation, this understanding relies on integrated analysis of the plate tectonic framework and the resulting paleogeography. The Arctic has experienced several episodes of tectonic deformation, which fundamentally changed the basin configuration and patterns of sediment routing. Here, we present a set of paleogeographic maps highlighting these changes during the Triassic–Paleogene. In the Triassic, the Arctic was characterized by a central restricted basin, which predominantly received clastic input from the Polar Urals and Arctic Canada. The Alaskan and Siberian passive margins received clastics from continent-scale drainage systems extending into the North American craton and the central Asian fold belt, respectively. In the Jurassic, the region was dominated by rifting as the central Arctic landmass rifted away from Laurentia. In the Early Cretaceous, the northern margin of the Barents Sea underwent regional uplift resulting in new provenance areas shedding sediments southward. Compression along the Pacific margin formed continuous topography and high sediment input to the Canada Basin during the Late Cretaceous. Regression in the Canada Basin continued in the Paleogene when major rift–tip deltas formed. This overview of Arctic paleogeography demonstrates the complexity of this overall data-poor area and shows the need for integrated, regional models to understand sediment routing and stratigraphic development in such areas.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉For both modeling and management of a reservoir, pathways to and through the seal into the overburden are of vital importance. Therefore, we suggest applying the presented structural modeling workflow that analyzes internal strain, elongation, and paleogeomorphology of the given volume. It is assumed that the magnitude of strain is a proxy for the intensity of subseismic scale fracturing. Zones of high strain may correlate with potential migration pathways. Because of the enhanced need for securing near-surface layer integrity when CO〈sub〉2〈/sub〉 storage is needed, an interpretation of three-dimensional (3-D) seismic data from the Cooperative Research Centre for Greenhouse Gas Technologies Otway site, Australia, was undertaken. The complete 3-D model was retrodeformed. Compaction- plus deformation-related strain was calculated for the whole volume. The strain distribution after 3-D restoration showed a tripartition of the study area, with the most deformation (30%–50%) in the southwest. Of 24 faults, 4 compartmentalize different zones of deformation. The paleomorphology of the seal formation is determined to tilt northward, presumably because of a much larger normal fault to the north. From horizontal extension analysis, it is evident that most deformation occurred before 66 Ma and stopped abruptly because of the production of oceanic crust in the Southern Ocean. Within the seal horizon, various high-strain zones and therefore subseismic pathways were determined. These zones range in width from 50 m (164 ft) up to 400 m (1312 ft) wide and do not simply follow fault traces, and—most importantly—none of them continue into the overburden. Such information is relevant for reservoir management and public communication and to safeguard near-surface ecologic assets.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In the last 30 yr, basin and petroleum system modeling (BPSM) has evolved into a large and diverse field encompassing a broad range of scientific disciplines. As BPSM is applied to an increasingly wide range of problems, what are, or should be, the future directions in the evolution of BPSM comes into question.To address this question, a survey was conducted at the AAPG Hedberg Research Conference on “The Future of Basin and Petroleum Systems Modeling,” held in Santa Barbara, California, April 3–8, 2016. To capture the full range of thoughts, participants were asked to list in priority order what they think are the three most important future directions in BPSM. The responses were collated into six general categories for analysis. The categorization process involved some qualitative judgements because some areas spanned several of the general areas.The results show that the most frequently cited directions are related to BPSM workflows, organizations, and processes. This category includes how modelers are used in an organization, how projects are executed, and how the results are interpreted and integrated.Migration modeling (primary and secondary) is the most frequently cited technical need. The results indicate that migration processes are not well understood and there are still substantial differences of thought about the processes involved and the best ways to model them.Some subjects, such as uncertainty and unconventionals, were mentioned in several of the general categories, whereas other subjects, such as increased functionality in the models, were only seldom mentioned.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Natural fractures are important storage spaces and fluid-flow channels in tight-oil sandstones. Intraformational open fractures are the major channels for fluid flow in tight-oil sandstones. Small faults may provide fluid-flow channels across different layers. According to analogous outcrops, cores, and borehole image logs, small faults and intraformational open fractures are developed in the tight-oil sandstones of the Upper Triassic Yanchang Formation in the southwestern Ordos Basin, China. Among them, high dip-angle intraformational open fractures are the most abundant. Northeast-southwest–trending fractures are the principal fractures for fluid flow because that is the present-day maximum horizontal compressive stress direction. Combined with production data, horizontal wells, striking normal to or at a large angle relative to the major flow pathways, are beneficial for tight-oil production improvement. Fractures with high dip angles are the main factor that influences initial oil production. Linkage and tip damage zones are more favorable for oil production improvement than wall damage zones. This study provides an example of natural fracture characterization and unravels fracture contributions to reservoir physical properties and oil production of tight-oil sandstones, which could provide a geological basis for oil exploration and development in tight sandstones.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Using recently acquired three-dimensional seismic data, we summarize typical patterns for seismic-based identification and stage analysis of sedimentary units in the Eocene succession of the southern slope-break belts of the Bozhong sag, Bohai Bay Basin, China. The sedimentary units in the study area are characterized by progradational reflectors and mound-shaped, bidirectional downlapping reflectors in dip and strike directions, respectively. Differential characteristics of a distinct sedimentary unit within one lobe are documented. The major provenance direction is defined and characterized by the largest dip angles of reflectors, the longest transport distance of sediments, and the thickest deposits in comparison to other dip directions—all recognized in this study and serving as typical characteristics for sedimentary unit identification and separation from the overlapped sedimentary complex. This study also summarizes diverse patterns—including collateral and prograding types—of sedimentary unit contact relationships and stage analysis along dip and strike directions. Collateral patterns are composed of three subtypes: superimposed, antithetic, and isolated. Three sedimentary units—S1, S2, and S3—are recognized in the study area. Summarized patterns of sedimentary unit contact relationships indicate that S1 was deposited earliest and S3 latest. The proposed patterns supplement seismic-based sedimentologic studies. This work may serve as a useful reference for sand-body characterization and stage analysis in other basins and similar areas.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Instead of using discrete values for properties that influence the volumetric calculation for recoverable reserves from the Middle Bakken, Pronghorn, and Three Forks reservoir rocks in the Williston Basin in North Dakota, an uncertainty-based assessment method was used. Various estimates have been published in the past that attempt to quantify recoverable reserves from the Bakken petroleum system. The Bakken–Three Forks trend is regarded as an unconventional tight oil play typical of a continuous-type basin-centered accumulation. However, production data reveal that areas are unequal and that certain regions stand out as sweet spots whereas others exhibit fairly high water cuts. This paper is based on 28 well models, which have been porosity-calibrated and adjusted for the prevalent thermal regime. The area of interest was delineated by geological parameters such as shale maturity and reservoir rock presence as well as existing production data. The purpose of this study is to use an uncertainty assessment method based on hundreds of basin model simulations that sample ranges of probable input parameters to quantify the recoverable reserves from the Bakken petroleum system in North Dakota. The results are displayed in reverse cumulative probability plots, tornado sensitivity charts, as well as in maps of the 10% chance, 50% chance (P50), 90% chance values. This means that there is an X% chance of success or an X probablity of realizing a certain amount of hydrocarbon. The P50 results of the uncertainty assessment indicate that approximately 4 billion bbl of oil and 3.6 tcf (102 billion m〈sup〉3〈/sup〉) of gas are recoverable from the Middle Bakken, Pronghorn, and Three Forks reservoir rocks in North Dakota. The Bakken–Three Forks trend appears to be an overcharged petroleum system, where the available pore space in reservoir rocks is the limiting factor for each accumulation.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Węglówka oil field is located in the outer Carpathians. The outer Carpathians are a region where hydrocarbons were discovered and exploited at the end of the nineteenth century in several dozen oil fields, which are relatively small. The Węglówka oil field is one of the largest in this region. In the 150 yr or so of hydrocarbon exploration in the area, more than 1 million t (〉1,237,000 tons [〉8,841,000 bbl]) of oil have been produced. Hydrocarbons are concentrated in Lower Cretaceous sandstones (Grodziszcze and Lgota sandstones) that form an anticline sealed by Upper Cretaceous marls called the Węglówka marls. These cap rocks are up to 600 m (2000 ft) thick. Because of the thrust-related exhumation, they were exposed at the surface and represent the youngest deposits in the region. The present work is focused on a detailed petrographic characterization of the Węglówka marls. This study allows petroleum geologists to better understand the evolution of porosity in these cap rocks and can serve as a foundation for the prediction of their sealing properties. The marls appear as a succession of interbedded red and green varieties, which occur in up to 2-m (6-ft)-thick beds. These beds are nonarenaceous, soft, and bioturbated. Grain size corresponds to approximately 80% clay and less than 20% silt fractions. X-ray diffraction (XRD) reveals that the marls contain, on average, 54% clay, 28% calcite, 16% quartz, up to 3% feldspars and, in red marls, 3% hematite. The XRD patterns of clay are typical of mixed-layer illite–smectite ([I–S]; 40% illite in I–S). The clay structures are dioctahedral with similar octahedral Mg and relatively high Fe〈sup〉3〈/sup〉〈sup〉+〈/sup〉 contents both in the red and green intervals. As revealed by standard petrography combined with high-resolution petrography performed through the use of a field emission scanning electron microscope, the marls have mudstone textures according to Dunham’s (1962) classification and are mostly composed of coccoliths and clay with rare nanoquartz. This rock may be considered an impure chalk. Sealing properties of the Węglówka marls are indicated by the specific surface area, porosity, pore size, and permeability, calculated using N〈sub〉2〈/sub〉 gas adsorption, helium, and mercury porosimetry. The sealing potential is postulated to result from a combination of the following: (1) origin of components (i.e., deposition of minute calcareous bioclasts and volcanic material as a source for clay); (2) oxygenated sedimentary environment (as a result of the presence of oxygen in the sediments, burrowing caused the rocks to be homogenized); and (3) tectonic-induced clogging of pore space because of reorganization of clay flakes (the rocks were strongly tectonically deformed, which resulted in reduction of porosity in clay aggregates).〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Calcite cementation has been identified as an active process in the Upper Triassic Yanchang Formation throughout its burial history and as a major diagenetic factor causing strong reservoir heterogeneities. The origins of calcite cements and their relevance to reservoir heterogeneities were investigated using a suite of petrographic and geochemical methods, including optical microscopy with fluorescence and cathodoluminescence, scanning and backscattered electron microscopy with energy-dispersive spectrometry, x-ray diffraction, x-ray fluorescence, electron probe microanalysis, quantitative evaluation of minerals by scanning electron microscopy, fluid inclusion analysis, and carbon and oxygen stable isotope analyses. The sandstones are compositionally immature with relatively high amounts of volcanic rock fragments. The two generations of calcite cements are Ca-I and Ca-II. The Ca-I calcites are distributed along the interface of sandstone and mudstone units and were formed during the Late Triassic to Early Jurassic at formation temperatures of approximately 90°C. The Ca-II calcite mainly developed in the lower part of the fining-upward sandstone units and was formed in the Late Jurassic at higher temperatures of approximately 110°C. The origins of calcite cements were constrained by geochemical and isotope measurements, fluid inclusion homogenization temperature, and in situ element analysis. The Ca-I calcite cement originated from dissolution of the lacustrine depositional carbonates in the interbedded mudstones and reprecipitation in the adjacent sandstones. The Ca-II calcite was mainly related to organic matter decarboxylation, with Ca〈sup〉2+〈/sup〉 having been provided internally by volcanic fragment alteration and plagioclase dissolution. Calcite cementation had caused strong reservoir heterogeneities in the Yanchang Formation tight sandstones. The Ca-I calcite cementation destroyed reservoir properties along the interface of sandstones and mudstones. The lower parts of the fining-upward sandstone units were tightly cemented by Ca-II calcite, although they originally had high porosity and permeability. The middle–upper parts of the fining-upward sandstone units contain less calcite cements and thus have better preserved reservoir pores because of oil emplacement inhibiting the calcite cementation processes.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Fuling shale gas field is located in a mountainous area, with well-developed underground rivers and karst caves. It also has a highly concentrated population, so the shale gas development in this field is faced with environmental protection problems. Combined with the characteristics of surface natural environment in the Fuling shale gas field and the features of shale gas development engineering, the main environmental issues encountered in the development of the Fuling shale gas field were analyzed. Studies on intensive land use, water conservation and protection, harmless use and disposal of oil-based drill cuttings, recycling of wastewater from drilling and fracturing, and green environment management mode for shale gas development were conducted, and the green development technology system suitable for the Fuling shale gas field was established. Field applications showed that, after applying the green development technology, the land occupation was reduced by 62.l%, the recycling rate of drilling and fracturing wastewater was up to 100%, the oil content of treated oil-based drill cuttings was less than 0.3%, and carbon dioxide emission was reduced by 64.47 × 10〈sup〉4〈/sup〉 t (1.41 × 10〈sup〉9〈/sup〉 lb). Thus, the goal of zero contamination was realized during shale gas field development. Research showed that the green and environmental protection development technology for the Fuling shale gas field has served as a valuable demonstration in the environmental protection in large-scale development of shale gas fields in China.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Increased oil and gas production in many areas has led to concerns over the effects these activities may be having on nearby groundwater quality. In this study, we determine the lateral and vertical extent of groundwater with less than 10,000 mg/L total dissolved solids near the Lost Hills–Belridge oil fields in northwestern Kern County, California, and document evidence of impacts by produced water disposal within the Tulare aquifer and overlying alluvium, the primary protected aquifers in the area.The depth at which groundwater salinity surpasses 10,000 mg/L ranges from 150 m (500 ft) in the northwestern part of the study area to 490–550 m (1600–1800 ft) in the south and east, respectively, as determined by geophysical log analysis and lab analysis of produced water samples. Comparison of logs from replacement wells with logs from their older counterparts shows relatively higher-resistivity intervals representing the vadose zone or fresher groundwater being replaced by intervals with much lower resistivity because of infiltration of brines from surface disposal ponds and injection of brines into disposal wells. The effect of the surface ponds is confined to the alluvial aquifer—the underlying Tulare aquifer is largely protected by a regional clay layer at the base of the alluvium. Sand layers affected by injection of produced waters in nearby disposal wells commonly exhibit log resistivity profiles that change from high resistivity in their upper parts to low resistivity near the base because of stratification by gravity segregation of the denser brines within each affected sand. The effects of produced water injection are mainly evident within the Tulare Formation and can be noted as far as 550 m (1800 ft) from the main group of disposal wells located along the east flank of South Belridge.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Recent oil discoveries in an Aptian–Cenomanian clinothem in Arctic Alaska demonstrate the potential for hundred-million- to billion-barrel oil accumulations in Nanushuk Formation topsets and Torok Formation foresets–bottomsets. Oil-prone source rocks and the clinothem are draped across the Barrow arch, a structural hinge between the Colville foreland basin and Beaufort Sea rifted margin. Stratigraphic traps lie in a favorable thermal maturity domain along multiple migration pathways across more than 30,000 km〈sup〉2〈/sup〉 (10,000 mi〈sup〉2〈/sup〉). Sediment from the Chukotkan orogen (Russia) filled the western Colville basin and spilled over the Beaufort rift shoulder, forming east- and north-facing shelf margins. Progradational shelf margin trajectories change abruptly to “sawtooth” trajectories at midclinothem, the result of reduction in sediment influx. Two stratigraphic trap types are inferred in Nanushuk basal topsets in the eastern part of the clinothem: (1) lowstand systems tracts, inferred to reflect forced regression, include a narrow, thick progradational stacking pattern perched on a sequence boundary on the upper slope; and (2) highstand-progradational systems tracts include a broad, thin wedge of shingled parasequences above a toplap surface. Both include stratigraphically isolated sandstone sealed by mudstone. Trap geometries in Torok foreset and bottomset facies in the same area include basin-floor fan, slope-apron, and slope-channel deposits that pinch out upslope and are sealed by mudstone. Significant potential exists for the discovery of additional oil accumulations in these stratigraphic trap types in the eastern part of the clinothem. Less potential may exist in the western part because reservoir–seal pairs may not be well developed.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Porosity is one of the most important rock properties in describing hydrocarbon reservoirs. Tests on core samples provide direct and representative porosity data, and the measurement of porosity at high confining pressures is recognized to correlate well with subsurface reservoir porosity. Whereas theoretical deductions of the changes and relationships of pressures, volumes, and compressibility suggest that porosity is reduced during the coring and lifting processes, the porosity measurement at elevated confining pressure does not evaluate original reservoir porosity. This theory is quantitatively validated by repeated laboratory experiments of loading and unloading on sandstone core samples. When the in situ confining pressure is approximately 30–35 MPa (∼4350–5076 psi), coring and lifting would cause a porosity reduction of approximately 1.2%–1.6%, and the porosity test under high confining stress results in further porosity loss. A revised approach in calculating reservoir porosity from cored samples is proposed and can have significant implications for reserve calculations, recovery factors, and geostatistical reservoir models. The study is important for both conventional and unconventional reservoirs because it discusses a fundamental mechanism of porosity change.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In this paper, high-resolution three-dimensional seismic data are used to interpret a transpressional salt tectonic structure in the Yingxiongling area, Qaidam Basin, China. The geometries of the salt structure and the Shizigou fault system that intersects it are precisely depicted. The Shizigou fault system is composed of suprasalt and subsalt components. The suprasalt component is a Y-shaped reverse fault, and the subsalt component is a complex flower structure. In previous studies, suprasalt and subsalt components were interpreted as two independent fault systems. This paper proposes instead that the suprasalt and subsalt faults are kinematically related and decoupled across the salt layer.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Miocene carbonate reservoirs in Central Luconia, offshore Sarawak, Malaysia, have been delivering gas for over 30 yr. In this paper, learnings from that period of production are used to understand the key drivers affecting flow during production and recovery optimization in existing fields as well as development decisions for new discoveries. The large data set, generated over more than 40 yr, was analyzed in a consistent manner through a holistic database, constrained by a stratigraphic framework, to allow reservoir units to be compared like-for-like (“integrated knowledge base” [IKB] concept). Carbonate reservoir heterogeneities impacting flow are grouped into “horizontal–heterogeneities”—argillaceous flooding layers and exposure-related karst—and “vertical–heterogeneities”—large-scale architectural elements, found especially along platform margins. Both types of heterogeneities control water ingress during production and influence the recovery mechanism. Argillaceous flooding layers can act as baffles, holding back water rise during production, or can form pressure compartments. Long-lived, fault-bounded reef margins, carbonate shoals, islands, and karsts can be vertical conduits for aquifer inflow. Platform shape and architecture impact column height and hence recovery efficiency. Additional drivers impacting recovery were found to be gas-column height, aquifer size and permeability, pressure connection to neighboring fields, and field development concepts. All drivers identified impact decisions throughout the field life, e.g., well count and design, intervention capabilities, evaluation and mitigation of early-water breakthrough, reservoir management, selecting enhanced recovery methods, and abandonment pressure. The IKB allowed to derive “big rules” on what matters for flow, which were used to decide on development strategies for greenfields in Central Luconia. The presented outcomes can be extrapolated to comparable carbonate systems, whereas the IKB approach can be adapted and applied to other mature basins and reservoir types where equally vast and historic data sets are awaiting to be used in the current era of digitalization.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Petroleum mobility in shale is closely correlated with the attributes of shale petroleum and pores; however, the relationship between these attributes is poorly understood. To characterize petroleum mobility in self-sourcing reservoirs, a suite of mature Eocene shales was selected and subjected to organic solvent extraction, and both the raw and solvent-treated samples were analyzed using pyrolysis, nitrogen adsorption, and x-ray diffraction. The results show that the pore surface area and pore volume of these shales are mainly controlled by their clay and quartz content rather than their organic matter (OM) content and are limited by the presence of carbonates. Correlations of soluble OM with pore surface area and volume after solvent extraction indicate that petroleum mobility of studied shales is initiated when the petroleum content reaches 0.70 wt. % of the rock and the pore diameter is over 12.1 nm. These thresholds are established in the studied area and should be similar for the self-sourcing reservoirs from similar sedimentary environments. This work proposes a method to reveal the thresholds of petroleum content and pore diameter for petroleum mobility in self-sourcing reservoirs, which is useful in the assessment of petroleum producibility and is of significance for unconventional petroleum exploration and exploitation.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉For oil-rich shales, current solvent extraction– and thermal extraction–based methods inaccurately measure hydrocarbon-filled porosity (〈span〉φ〈/span〉〈sub〉〈span〉HC〈/span〉〈/sub〉). Moreover, the hydrocarbon composition is not characterized by either method. Here, we show how open-system programmed thermal extraction and pyrolysis, LECO total organic carbon, Archimedes bulk density, and helium pycnometry measurements are integrated to calculate oil and gas pore volumes, characterize their composition, and estimate mobility. Use of a modified multiramp, slow-heating thermal extract, and pyrolysis temperature program further subdivides the 〈span〉φ〈/span〉〈sub〉〈span〉HC〈/span〉〈/sub〉. Saturate–aromatic–resin–asphaltene (SARA) separation and gas chromatography of solvent-extracted organic matter and thermally extracted oils are used to compositionally classify the 〈span〉φ〈/span〉〈sub〉〈span〉HC〈/span〉〈/sub〉. The segregated bulk compositions of gas- and oil-filled porosity measured via this method are shown to overlap and are broken into the following categories: gas-filled porosity (∼C〈sub〉1〈/sub〉–C〈sub〉14〈/sub〉), light oil–filled porosity (∼C〈sub〉6〈/sub〉–C〈sub〉36〈/sub〉), and heavy oil–filled porosity (∼C〈sub〉32〈/sub〉–C〈sub〉36〈/sub〉+). Furthermore, slow-heating multiramp thermal extraction can subdivide the light oil–filled porosity into four components capturing the C〈sub〉11〈/sub〉–C〈sub〉13〈/sub〉, C〈sub〉12〈/sub〉–C〈sub〉16〈/sub〉, C〈sub〉14〈/sub〉–C〈sub〉20〈/sub〉, and C〈sub〉17〈/sub〉–C〈sub〉36〈/sub〉 ranges of the extractable organic matter. Analysis of solvent-extracted oils by SARA identifies abundant saturates and aromatics in the light oil–filled porosity and abundant resins and asphaltenes in the heavy oil–filled porosity. Low-maturity shales can be dominated by heavy (C〈sub〉32〈/sub〉+) oils rich in asphaltene and resin fractions not observed in the produced fluid. The ratios of SARA components in the C〈sub〉15〈/sub〉+ fraction of produced fluid and core extract can be used to better estimate the potentially mobile 〈span〉φ〈/span〉〈sub〉〈span〉HC〈/span〉〈/sub〉.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Fuling gas field in Sichuan Basin, China, has produced greater than 1.5 × 10〈sup〉10〈/sup〉 m〈sup〉3〈/sup〉 (0.53 tcf) of natural gas from overmature Upper Ordovician Wufeng and lower Silurian Longmaxi shales. To systemically investigate the characteristics of wettability and connectivity and to understand the underlying causes of production behavior, we study five samples of Wufeng and Longmaxi shales with different total organic carbon contents and mineral compositions. Complementary approaches include mercury intrusion capillary pressure (MICP), contact angle measurement, spontaneous imbibition and saturated diffusion, and tracer (both nonsorbing and sorbing) migration mapped via laser ablation inductively coupled plasma mass spectrometry. According to measured contact angles and imbibition tests conducted on aqueous (deionized water and brine) and oleic (n-decane) phases, Wufeng and Longmaxi shales are strongly oil wet and moderately strong water wet. The lower boundary of estimated permeability obtained from n-decane imbibition can reach 137 nd, which is higher than the geometric mean permeability derived from the MICP method (5.5–68.8 nd). Effective diffusion coefficients of the Wufeng and Longmaxi shales are in the range of 10〈sup〉−13〈/sup〉 m〈sup〉2〈/sup〉/s (1.1 × 10〈sup〉−12〈/sup〉 ft〈sup〉2〈/sup〉/s). Tests of imbibition and saturated diffusion using tracer-containing brine show that concentrations of nanometer-sized tracers decrease rapidly (a factor of 〉10) over a migration distance of a few millimeters from the sample edge, suggesting the presence of poorly edge-connected water-wet pores. Sparsely connected hydrophilic pores, mixed wettability, and highly restricted pathways collectively contribute to the limited migration of nano-sized tracers, which probably results in the production behavior of initial steep decline and low overall recovery in the Fuling gas field.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The upper zone of the Lower Cretaceous Kharaib Formation (151–177 ft [46–54 m] thick in the studied wells) is a major oil reservoir in several giant oil fields. Wide variations in porosity and permeability of this zone have been shown to result from both the inhibition of burial cementation by oil in the crest of each field and localized cementation adjacent to stylolites, combined with the more subtle influence of widely varying depositional mud content and grain size. The present study examines these relationships in closer detail, using core and petrographic observations from two wells on the oil-filled crest and two wells on the water-filled flanks of a giant domal oil field.Although porosities are higher overall in the crestal cores, each well shows wide variations within each of seven main groupings of the samples by depositional texture. This heterogeneity results mainly from the distribution of clay, which is concentrated along depositional laminations and causes widely varying porosity losses in all textures by promoting stylolite development and associated calcite cementation. Higher clay abundance (and lower porosity) within the upper and lower 12–17 ft (4–5 m) of the reservoir reflects increased influx of siliciclastic fines across the epeiric Barremian carbonate platform immediately following and preceding, respectively, third-order falls in global sea level. Most (95%) of porosity-permeability data from the studied wells lie within Lucia rock-fabric class 3, showing distinct but relatively subtle differences between texture groups, whereas a subordinate part of the data from the upper, relatively mud-poor third of the reservoir plot at higher permeabilities. Development of a predictive model for the petrophysical heterogeneity of this example requires a combination of the following: (1) a diagenetic model for porosity controls; (2) the use of a modestly higher porosity-permeability transform (upper class 3) in the upper part of the reservoir than in the lower reservoir (lower class 3); and (3) a recognition of the scattered and widely varying occurrences of exceptionally high permeabilities in the upper reservoir.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In the Paleocene to Eocene Wilcox Group in the northern Gulf of Mexico, exploration targets are reaching into deep to ultradeep burial depths. At these great depths, reservoir quality (porosity and permeability) becomes an important risk factor in determining the chance of encountering an economic reservoir. Major controls on reservoir quality are pore types and abundances, pore-throat sizes, and pore network composition. These factors can be analyzed by integrating petrographic, core plug porosity and permeability, and mercury injection capillary pressure (MICP) analyses. The Wilcox sandstones are mostly lithic arkoses and feldspathic litharenites that contain primary interparticle pores, secondary dissolution pores, and micropores. However, these pore types evolve with depth and temperature. As temperature increases, the relative abundance of primary interparticle pores decreases, whereas the relative abundance of secondary dissolution pores and nano- to micropores increases. Associated with this evolution of pore networks with increasing temperature, there is a decrease in reservoir quality. This decrease in reservoir quality is caused by a transition to finer pore-throat sizes that correspond to changes in pore types. Petrographic analysis provides information on pore types, core plug porosity and permeability analysis provides information on volume of pores and effectiveness of flow, and MICP analysis provides information on pore-throat radius distribution. Through forecasting the pore network in the target temperature zone, a realistic porosity versus permeability transform can be selected to estimate permeability from wire-line log porosity.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Canning Basin is a largely unexposed and underexplored frontier basin, formed mostly in the Paleozoic. Geological knowledge of this basin is based predominantly on sparse regional “vintage” two-dimensional seismic and small three-dimensional (3-D) seismic surveys and less than 230 exploration wells. Following seismic interpretation, an integrated interpretation was completed on airborne gravity gradiometer (AGG), magnetic, seismic, well, and complementary data along the southwestern margin of the Fitzroy trough and Gregory subbasin. Seismic data were reinterpreted using AGG data to produce a better constrained geological model. A basement structure map, two intrasedimentary structure maps, and a formation distribution map were produced. The interpretation of seismic profiles, validated through 2.5-dimensional gravity gradiometer modeling, is essential to this workflow.Repeatedly reactivated west–northwest and northwest structural trends, inherited from Proterozoic orogenies, respectively delineate the Fitzroy trough and the Gregory subbasin with its northwestern structural extension into the Fitzroy trough, the Gregory subbasin trend. Subsidence occurred during two periods of extension. An asymmetric extensional system of the Fitzroy trough controlled Ordovician–Silurian deposition of the Carribuddy Group. Devonian–Carboniferous subsidence defines the Gregory subbasin trend. This Pillara extension reactivated structures in the east of the Fitzroy trough. Simultaneous activity of both extensional fault systems and growth faulting controlled the facies and thickness distribution of carbonates and clastics of the early Carboniferous Fairfield Group. The Meda and Fitzroy transpressional phases inverted faults of the Gregory subbasin trend and Fitzroy trough, producing prospects by structural interference.The improved understanding of tectono-stratigraphic relationships, including the 3-D distribution of carbonate reservoirs, benefited the planning of seismic surveys, prospect evaluation, drilling, and acreage relinquishment.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Thermal conductivity is a major influencing factor on subsurface conductive heat transport and resulting temperature distribution, which in turn is a key parameter in basin modeling. Basin modeling studies commonly use representative literature values of thermal conductivity despite their impact on modeling results. We introduce a workflow for quantifying the effect of uncertain thermal conductivity on subsurface temperature distribution and thus on basin modeling results and test this workflow on a two-dimensional generic model from the Nordkapp Basin; a prior ensemble of possible models is conditioned according to Bayes’ theorem to incorporate prior knowledge of temperature data. This conditional probability yields a posterior ensemble of temperature fields with a significantly reduced standard deviation. To verify our approach, we use five characteristic scenarios from the posterior ensemble for transient petroleum systems modeling. How considering uncertain thermal conductivity affects variance in hydrocarbon generation is assessed by modeling corresponding vitrinite reflectances (〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉).Temperature uncertainty increases with depth. It also increases with increasing offset from the salt diapirs, which can be associated with a large lateral heat-flow component in the complex tectonic environment of the Nordkapp Basin. The introduced workflow can reduce temperature uncertainty significantly, especially in regions with high prior uncertainty. The 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉 is very sensitive to changes in thermal conductivity because the onset depth of the gas window in the Nordkapp Basin may vary by up to 800 m (2600 ft) within the 95% confidence interval. This demonstrates the importance of quantification of the uncertainty in thermal conductivity on thermal basin modeling.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Yinggehai–Song Hong Basin has received a large amount of terrigenous sediment from different continental blocks since the Paleogene. The Yingdong slope, which is located on the eastern side of this basin, is an important potential gas province, but the provenance of the marine sediments in this area are poorly understood. The detrital zircon U-Pb geochronology of sedimentary rocks from the lower Miocene to Quaternary is examined in this study to investigate the temporal and spatial variations in provenance since the early Miocene. The U-Pb ages of detrital zircon range from 3078 to 30 Ma, suggesting that sediment input is derived from multiple sources. Detailed analyses of these components indicate that both the Red River and Hainan are likely the major sources of the sediments on the Yingdong slope, with additional minor contributions from central Vietnam (eastern Indochina block) and possibly the Songpan–Garze block. Variations in the dominant detrital zircon populations within stratigraphic successions display an increasing contribution from the Red River since the middle Miocene. This resulted from the progradation of the Red River Delta in the northern basin and may have also been influenced by regional surface uplift and associated climate changes in East Asia. This study shows that the Red River has had a relatively stable provenance since at least the early Miocene, indicating that any large-scale drainage capture of the Red River should have occurred before circa 23 Ma.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Jurassic black mudstone and coal beds in the central Junggar Basin, northwestern China, are the major source rocks for the basin with type II〈sub〉2〈/sub〉 and type III (gas-prone) kerogens. Widespread overpressures are developed in the Jurassic stratigraphic interval. Sonic and resistivity logs display strong characteristic responses of overpressure in the mudstones, with anomalously high acoustic traveltimes and low resistivity compared with the normally pressured mudstones. The overpressured Jurassic sediment sequences appear to have undergone normal compaction because the mudstones exhibit no anomalously low bulk density. The overpressured mudstones deviate from the normally pressured mudstones in density–effective vertical stress space. The overpressure in the Jurassic source rocks is, therefore, not caused by disequilibrium compaction. The overpressured Jurassic sandstone reservoirs are predominantly oil and gas saturated or oil bearing. The well-log responses of the overpressured mudstones and seismic velocity characteristics indicate that the top depth of the overpressure zone ranges from 3800 to 4600 m (12,500 to 15,100 ft), corresponding to formation temperatures of approximately 94°C to 111°C (∼201°F to 232°F), with estimated vitrinite reflectance values of 0.6% to 0.75%. The Jurassic source rocks with overpressure are capable of generating hydrocarban at present and are currently overpressured. All the evidence suggests that the overpressure in the Jurassic source rocks in the central Junggar Basin is caused by hydrocarbon (HC) generation. The overpressure evolution was modeled quantitatively in response to pressure changes caused by HC generation during basin evolution. The results indicate that multiple episodes of overpressure development and release occurred within the Jurassic source rocks, suggesting multiple episodes of HC expulsion. The timing and numbers of these episodes of HC expulsion were thus determined from the modeled overpressure evolution.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Thermal properties of rocks are essential parameters for investigating the geothermal regime of sedimentary basins, and they are also important factors in assessments of hydrocarbon and geothermal energy resources. The Tarim Basin, the largest basin located in the north of the Tibetan Plateau, northwestern China, has great hydrocarbon resource potential and is an ongoing target for industry exploration. However, the thermal properties of sedimentary rocks within the basin are yet to be systematically investigated at a basin scale, thereby limiting our understanding of the thermal regime in the basin. Here, we collected 101 samples of sedimentary rocks and measured their thermal properties. Our results show that the ranges (and means) of thermal conductivity, radiogenic heat production, and specific heat capacity are 1.08–5.35 W/mK (2.52 ± 0.99 W/mK), 0.03–3.24 μW/m〈sup〉3〈/sup〉 (1.24 ± 0.87 μW/m〈sup〉3〈/sup〉), and 0.75–1.10 kJ/(kg·°C) (0.87 ± 0.07 kJ/(kg·°C)), respectively. Volumetric heat capacity and thermal diffusivity at the temperature of 40°C range from 1.61 to 2.79 MJ/(m〈sup〉3〈/sup〉·K) (2.26 ± 0.25 MJ/[m〈sup〉3〈/sup〉·K]) and 0.44–2.95 × 10〈sup〉−6〈/sup〉 m〈sup〉2〈/sup〉/s ((1.12 ± 0.53) × 10〈sup〉−6〈/sup〉 m〈sup〉2〈/sup〉/s), respectively. The thermal properties vary considerably for different lithologies, even within the same lithotype, indicating that thermal properties alone cannot be used to distinguish lithology. Thermal conductivity increases with increased burial depth, density, and stratigraphic age, suggesting the dominant influence is porosity variation on thermal conductivity. Furthermore, a strong contrast in the thermal properties of rock salt and other sedimentary rocks perturbs the geothermal pattern, which should be taken into consideration when performing basin modeling.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The relationship between base metal deposits, especially Mississippi Valley–type (MVT) Pb–Zn deposits, and hydrocarbons is not well constrained. This is despite the fact that hydrocarbons generally occur in MVT deposits; the ores are emplaced in the same temperature range as hydrocarbon maturation and migration, and the deposits commonly occur in proximity to metal-rich black shales. Better understanding should lead to better exploration models for both hydrocarbons and MVT deposits. This connection is better understood with the help of Pb isotope patterns. Sphalerite Pb isotope compositions from the northern Arkansas and Tri-State mining districts and Woodford–Chattanooga and Fayetteville Shales were determined to assess the potential of shales as source rocks for the ore metals. The ores in both districts have a broad range of Pb isotope ratios and define linear trends, suggesting mixing of Pb from two distinct end members. Current results and previous depositional environment studies indicate the following: (1) shales deposited mainly under nonsulfidic anoxic conditions represent the less radiogenic end member, or (2) shales are the only source of ore metals. Given the array of organic molecules, each with their own thermochemical range, and the ways metals can be associated with them, the release of metals may cover varying ranges. Thus, the compositions of the released fluids would change through time and not have a single static composition, closely approximating the isotopic composition of the released metals at various times. Mineralization derived from a dynamically evolving fluid may show apparent end members, without the need to call on mixing of fluids from separate sources.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉A subaqueous clinoform system has been identified from high-quality three-dimensional seismic data from the northeast Exmouth Plateau, North West Shelf, Australia, and was interpreted as a shelf–slope–basin clinoformal component of a Jurassic fluviodeltaic system (the Legendre delta). Several geomorphological features associated with shelf-slope development and subsequent rift tectonics were identified, including (1) submarine channels at slope to basin floor; (2) gullies on the slope; (3) slumps on the shelf; and (4) canyons, canyon-derived gravity flow deposits, and a fan lobe developed in subsequent rift processes.The results of this study provide insights into the controlling factors on the sinuosity, degree of erosion, and sediment gravity flows of channels developed at slope to basin-floor settings, which shed light on the way fluvial sands were transported across the shelf and slope to the basin floor. The geometries and distributions of gravity flow deposits, if confirmed by drilling, may serve as an analog for reservoir prediction in the deep-water fluviodeltaic settings. The gullies on the slope were interpreted as a result of dilute, sheetlike flows. The slumps on the shelf were interpreted as a result of nonslope-related causes.The syntectonic canyons, the canyon-derived gravity flow deposits, and the fan lobe present vivid examples of the erosion and sedimentation processes during active rift tectonics and have significant implications for understanding the rift processes of the North West Shelf, Australia, as well as other rift-related basins.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Giant petroleum accumulations worldwide with burial depths more than 7000 m (〉23,000 ft) occur mostly in Mesozoic and Cenozoic reservoirs and yield predominantly natural gas. Recently, however, a giant oil accumulation with reservoir depths between 7000 m (23,000 ft) and 8000 m (26,000 ft) was discovered in the lower Paleozoic section in the southern part of the Halahatang region in the Tarim Basin, China. Petroleum sourced from lower Paleozoic rocks is contained in Ordovician karst fracture-cave reservoirs and sealed by Middle–Upper Ordovician limestones and mudstones. The newly discovered superdeep accumulation is among the deepest black single-phase oil accumulations worldwide and opens up new avenues for petroleum exploration in deep-marine carbonate reservoirs. Reservoir pressures are between 75 MPa (10,878 psi) and 85 MPa (12,328 psi), with pressure coefficients between 1.2 and 1.7 and temperatures ranging between 140°C (284°F) and 172°C (342°F). Charging and accumulation of petroleum occurred during the late Hercynian orogeny, followed by subsequent gradual deep burial, which took place before rapid subsidence beginning circa 5 Ma. Following subsidence, the thickness of overlying strata increased by more than 2000 m (〉6600 ft) before finally attaining current depth. Therefore, this oil accumulation represents a well-preserved ancient petroleum system. Based on the geochemical features of oils and gases, the crude oils can be classified as mature, sourced from mixed marine organofacies of shale, marl, and carbonate, whereas the gases were cogenerated with oils. Despite very high present-day reservoir temperatures, no oil cracking has occurred because of the relatively short exposure of oils to high temperatures in a low geothermal gradient regime. Thus, there is significant exploration potential under similar conditions for liquid petroleum in superdeep strata. Faults and reservoirs are major factors controlling petroleum accumulation. Interlayer karsts with excellent fracture-cavity connectivity developed adjacent to faults, generally resulting in the enrichment of oil and gas along fault zones. High-quality reservoirs in this area are easy to identify because they exhibit strong bead-like amplitude features in seismic sections. Wells located near faults produce relatively large amounts of oil and gas. Effective karst fracture-cave reservoirs with noncracked oil may exist below 8000 m (26,000 ft) in the Tarim Basin and represent a significant exploration target in China.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Substantial amounts of petroleum were recently discovered in the Carboniferous andesite, tuff, breccia, and basalt reservoirs of the Chepaizi uplift in the western Junggar Basin. However, the charging history of the Carboniferous petroleum reservoir is poorly understood. Oil–oil correlation studies indicate that all of the oils were mainly derived from the middle Permian Wuerhe Formation source rocks, possibly mixed with a small contribution from Carboniferous Baogutu Formation source rocks in the neighboring Changji sag. Based on the petrographic and microthermometry of fluid inclusions, two hydrocarbon charging episodes are defined; these episodes were characterized by a low-peak-range homogenization temperature (〈span〉Th〈/span〉) distribution (80°C–90°C) and high salinity (13.22–13.42 wt. % NaCl) and a high-peak-range 〈span〉Th〈/span〉 distribution (120°C–130°C) and low salinity (4.89–11.72 wt. % NaCl), respectively. Through one-dimensional basin modeling and pressure–volume–temperature–composition simulation, the burial-thermal histories for wells P61, P66, P668, and P663 were reconstructed, and their trapping temperatures of the hydrocarbon inclusions were calculated to be higher than their corresponding highest paleotemperature (i.e., 56.8°C, 53.7°C, 60.9°C, and 58.1°C, respectively), implying fast hydrocarbon charging processes promoted by deep hydrothermal fluids. Associated with the hydrocarbon generation history, sealing process of the Hongche fault, and regional tectonic evolution, these two hydrocarbon charging events were deduced as the adjustments of oils previously accumulated along the Hongche fault zone, because of the tectonic extension in the Paleogene and regional tilting in the Neogene, respectively. The general direction of oil charging was traced from south to north and from east to west, as indicated by the molecular parameters of nitrogen-bearing compounds and C〈sub〉20〈/sub〉 + C〈sub〉21〈/sub〉 triaromatic steroids/C〈sub〉20〈/sub〉 + C〈sub〉21〈/sub〉 + C〈sub〉26〈/sub〉–C〈sub〉28〈/sub〉 triaromatic steroids (TA(I)/TA(I+II)), which roughly coincided with the active fracturing.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Modern oil and gas seismic surveys commonly use areal arrays that record continuously, and thus routinely collect “excess” data that are not needed for the conventional common reflection point imaging that is the primary goal of exploration. These excess data have recently been recognized to have utility not only in resource exploration but also for addressing a diverse range of scientific issues.Here we report processing of such discarded data from recent exploration surveys carried out in southeastern New Mexico. These have been used to produce new three-dimensional (3-D) seismic reflection imagery of a layered complex within the crystalline basement as well as elements of the underlying crust. This enigmatic basement layering is similar to that found on industry and academic seismic reflection surveys at many sites in the central United States. Correlation of these reflectors with similar features encountered by drilling in northwestern Texas suggest that they may be part of an extensive, continental-scale network of tabular mafic intrusions linked to Keweenawan rifting of the igneous eastcentral Unites States during the late Proterozoic. More importantly, this analysis clearly demonstrates that the new generation of continuously recorded 3-D exploration datasets represent a valuable source of fresh information on basement structure that should be examined rather than discarded. Such basement information is not only important to understanding crustal evolution, it is directly relevant to assessing risks associated with fossil fuel extractions, such as induced seismicity related to waste water injection.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The detailed depositional systems and basin evolution of Lower Cretaceous coal-bearing strata in the Erlian Basin of northeastern China were analyzed based on extensive borehole and outcrop data. A total of 7 facies associations are interpreted and consist of 14 distinct lithofacies, with lithologies including conglomerates, sandstones, siltstones, mudstones, shales, and coals. Five third-order sequences were recognized, and their internal lowstand, transgressive, and highstand systems tracts were defined based on six key sequence stratigraphic boundaries. These boundaries were represented by regional unconformities, basal erosional surfaces of incised valley fills, interfluvial paleosols, and abrupt depositional facies-reversal surfaces. Sequences I–V correspond to the rift-initiation stage, the early-rift climax stage, the late-rift climax stage, the immediate postrift stage, and the late postrift stage of the basin, respectively. The preferred sites for coal accumulation were braided fluvial delta plain, meandering fluvial delta plain, and littoral–shallow lake environments. The major coal seams formed during the early and late transgressive systems tract of sequences III, IV, and V, which were well developed in the eastern, northeastern, and northeastern parts of the Erlian Basin, respectively. Three coal depositional models were summarized in the sequence stratigraphic framework, including types 1, 2, and 3, corresponding to the Newark type, Newark–Richmond type, and Richmond type, respectively. These coal depositional models were closely related to the basin evolution. These results could provide preferred depositional environments and favorable areas of coal and coalbed methane (CBM) for the exploration and development of coal and CBM in the Erlian Basin, with the Jiergalangtu, Huolinhe, Baiyinhua, and A’nan sags recommended as the key sags.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Different from most tight oil reservoirs, the tuffaceous tight oil of the Tiaohu Formation is not in situ oil, and no close contact exists between the source rock and reservoir in the Malang sag (Santanghu Basin, China). This study determined the mechanism of hydrocarbon accumulation of this tuffaceous tight oil reservoir through an integrated analysis of oil–source rock correlation, reservoir characteristics, and rock wettability combined with a comprehensive analysis of geological conditions. An oil–source rock correlation using biomarkers and stable carbon isotopes shows that the crude oil originated from underlying source rocks in the Lucaogou Formation. The oil in the tuffaceous tight reservoir is not indigenous but has migrated over a long distance to accumulate in these reservoirs. Faults and fractures that developed at the end of the Cretaceous are the oil migration pathways. Vitric and crystal-vitric tuffs constitute the main rock types of the tuffaceous tight reservoir. Matrix-related pores in the tuffs mainly comprise interparticle pores between minerals and dissolution intraparticle pores formed by devitrification. The adsorption of polar components of the oil generated from original organic matter in the tuff leads to wettability of lipophilicity, which is the main reason for hydrocarbon charging and accumulation. To our knowledge, this is the first comprehensive study reporting this finding.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Predicting vertical and lateral facies variations in various depositional environments remains a major challenge in the oil and gas industry because it impacts petroleum system assessments and the associated exploration-risking phase. The use of multidisciplinary constraints (geomorphology, geology, geophysics) in forward stratigraphic models sheds light on the complex interaction of local, regional, and global driving mechanisms that influence sediment transport and deposition along continuously evolving landscapes. In this paper, we develop an integrated statistical approach to examine the sensitivity of forward stratigraphic models in complex salt provinces to several parameters, including water discharge, sedimentary load, grain size and associated diffusion coefficients, and slope. This statistical analysis was applied to the Barremian–Albian sequence of the central Scotian Basin (Canada) and highlights the influence of complex salt kinematics on sediment pathway diversion and accumulation around salt domes and canopies. Forward stratigraphic modeling results point to regions of higher probability of Lower Cretaceous sandy reservoirs. Automating simulation runs significantly reduced the time required to achieve a statistically valid number of simulations and allowed the sensitivity of the model to be evaluated.〈/span〉

Description: 〈span〉In his comment on our paper about the hyperpycnites in the Triassic Yanchang Formation, G. Shanmugam puts forward that hyperpycnites do not exist. Consequently, he considers our interpretation that hyperpycnal flows are an important depositional process in the Yanchang Formation to be invalid. We unravel his arguments and demonstrate that evidence supports our assertion that hyperpycnal flows were an important sedimentary process in the lake in which the Yanchang Formation accumulated. Moreover, we provide proof from modern observations that hyperpycnal flows do exist in lakes.〈/span〉

Description: 〈span〉〈a href="https://pubs.geoscienceworld.org/aapgbull#b27"〉Yang et al. (2017b)〈/a〉 have advocated the importance of hyperpycnites by using a genetic facies model proposed for deposits of hyperpycnal flows by 〈a href="https://pubs.geoscienceworld.org/aapgbull#b14"〉Mulder et al. (2003)〈/a〉. The problem is that the authors have ignored experimental flume results and other empirical field data that discredited the model. This discussion is a rigorous evaluation of data, documentation, and the facies model.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉We analyze western Caribbean structural styles and depositional controls associated with Late Cretaceous–Cenozoic deformational events using a 1600-km (994-mi)-long, regional, northwest–southeast transect extending from the Cayman Trough in Honduras to northern Colombia. Different structural provinces defined along the transect include (1) the Cayman Trough and adjacent Honduran borderlands marking the North American–Caribbean transtensional plate boundary characterized by late Eocene–Holocene fault-controlled depocenters; (2) the Nicaraguan Rise that includes continental Paleocene–Eocene rocks deposited in sag basins, which are overlain by relatively undeformed Miocene–Holocene carbonate and clastic shelf deposits of the northern Nicaraguan Rise, following a Late Cretaceous convergent phase; (3) the Colombian Basin that includes thick Miocene clastic depocenters and the localized presence of Upper Cretaceous rocks overlying the basement and where much of the subsidence is likely isostatic and flexurally driven given its proximity to the subduction zone of northern Colombia; (4) the south Caribbean deformed belt, an active, accretionary prism produced by the subduction of the Caribbean large igneous province beneath the South American plate, which has deformed the Cenozoic prism and fore-arc section and produced thrust-fault–controlled accommodation space for upper Miocene–Holocene piggyback deposits; and (5) the onshore Cesar–Rancheria Basin in northern Colombia, which has recorded the uplift of its bounding mountain ranges, the Sierra de Santa Marta massif to the west and Perija Range to the east. Plate reconstructions place the various crustal provinces along the transect into the context of the Late Cretaceous–Cenozoic deformation events that can be partitioned into strike-slip, convergent, and extensional components.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The lacustrine shale of the Upper Cretaceous Qingshankou Formation is the principal prospective unconventional target lithology, acting as source, reservoir, and seal. Lithofacies and associated storage capacity are two significant factors in shale oil prospectivity. This paper describes an investigation of the lower Qingshankou Formation lacustrine shale based on detailed description and analysis of cores, shale lithofacies characteristics, depositional setting, and stacking patterns.Seven lithofacies are recognized based on organic matter content, sedimentary structure, and mineralogy, all exhibiting rapid vertical and lateral changes controlled by the depositional setting and basin evolution. An overall trend from shallow-water to deep-water depositional environments is interpreted from the characteristics of the infilling sequences, characterized by increasing total organic carbon (〈span〉TOC〈/span〉) and total clay content and decreasing layer thickness (i.e., from bedded to laminated then to massive sedimentary structures). Periods of deposition during shallowing cycles show a reverse trend in the sedimentary characteristics described above. The sedimentary rocks in the studied interval show three complete short-term cycles, each one containing progressive and regressive system tracts.Massive siliceous mudstones with both high and moderate 〈span〉TOC〈/span〉 are considered to have the best hydrocarbon generation potential. Laminated siliceous mudstones, bedded siltstones, and calcareous mudstones with moderate and low 〈span〉TOC〈/span〉 could have the same high hydrocarbon saturations as the high-〈span〉TOC〈/span〉 massive siliceous mudstones, but these lithologies contain more brittle minerals than the massive mudstones. Several siltstone samples show low or zero saturation of in situ hydrocarbons; this is considered to be related to a combination of fair to poor hydrocarbon generation potential and extremely low permeability, limiting migration. Moderate-〈span〉TOC〈/span〉 laminated siliceous mudstones were also observed to have connective pore-fracture networks. It can be demonstrated that successive thick sequences of moderate-〈span〉TOC〈/span〉 laminated siliceous mudstones, showing high volumes of hydrocarbon in situ, a high mineral brittleness index, and good permeability, combine to form shale oil exploration “sweet spots.”〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Helium and nitrogen variations in Panhandle–Hugoton field (PHF) gases are products of interaction between hydrocarbon gas from the Anadarko basin and at least two water masses with dissolved nitrogen and helium. The two most distinct water masses are from the Palo Duro basin (highest He/N〈sub〉2〈/sub〉) and the Hugoton embayment (lowest He/N〈sub〉2〈/sub〉). Geochemical data indicate several hundred million years of helium generation in porous rock. Helium migrated to the gas by diffusion through water-saturated rock and by west-to-east water flow.Sediment and basement helium generation and helium migration were modeled to validate timing and source of PHF helium. Models indicate a predominantly sedimentary helium source with some basement helium charge on the Amarillo uplift. Helium in the central and eastern PHF diffused from underlying rocks, whereas gases on the west and southwest sides were enriched in nitrogen and helium delivered by hydrodynamic water flow.Nitrogen in high-nitrogen gases was probably sourced as ammonium released from clays by cation exchange with brines derived from overlying salt units. The amount of mudrock (nitrogen and helium source) relative to other potential helium sources (arkose, radioactive dolomite) correlates to decreasing gas He/N〈sub〉2〈/sub〉.The high helium concentrations in PHF gases result from multiple favorable circumstances. Old pore water accumulated dissolved helium during hundreds of millions of years of helium generation in sediment. High water/gas and low pressure favored higher helium concentrations in gas. Hydrodynamic flow delivered helium-rich pore water from basins west of the PHF.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉As an important unconventional and alternative resource, shale gas has attracted worldwide attention. The breakthrough pressure is a major factor in the generation and migration of shale gas as well as in the evaluation of the caprock sealing capacity. Carboniferous shales are considered to have great potential for the exploitation of shale gas; thus, investigations of the breakthrough pressure and gas effective permeability are significant. Two shale samples taken from the Carboniferous Hurleg Formation in the eastern Qaidam Basin, China, were chosen to conduct breakthrough experiments to investigate the effects of water saturation and CO〈sub〉2〈/sub〉–CH〈sub〉4〈/sub〉 mixed mole fractions on the breakthrough pressure and gas effective permeability. Prior to the experiments, various relevant parameters (e.g., the porosity, mineral composition, and organic geochemistry; the total organic content, thermal maturity and kerogen type; and microstructure) of these samples were also measured.The results of our breakthrough experiments show that the breakthrough pressure increases with the water saturation and decreases with the CO〈sub〉2〈/sub〉 mole fraction in the gas mixture. The situation for the gas effective permeability is just the opposite. Pore-size distribution measurements indicate that there are many nanoscale micropores that can easily be blocked by water molecules. This results in the reduced connectivity of gas pathways; thus, the breakthrough pressure increases and the gas effective permeability decreases with increasing water saturation. The breakthrough pressure decreases with the CO〈sub〉2〈/sub〉 mole fraction because the interfacial tension of the CO〈sub〉2〈/sub〉–water system is smaller than that of the CH〈sub〉4〈/sub〉–water system. The viscosity of the CO〈sub〉2〈/sub〉–CH〈sub〉4〈/sub〉 mixture was found to increase with the CO〈sub〉2〈/sub〉 mole fraction by fitting a series of values under the same temperature and pressure conditions, leading to an increase in the gas effective permeability. Furthermore, CO〈sub〉2〈/sub〉 molecules are smaller than CH〈sub〉4〈/sub〉 molecules, making it easier for CO〈sub〉2〈/sub〉 to move across pathways. After each breakthrough experiment, the CO〈sub〉2〈/sub〉 mole fraction in the effluent was less than that in the injected gas, and it increased over time until reaching the initial injected gas composition. This is because the adsorption and solubility of CO〈sub〉2〈/sub〉 in water are greater than those of CH〈sub〉4〈/sub〉. This study provides practical information for further investigations of shale gas migration and extraction and the sealing capacities of caprocks.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉A comprehensive study on rift stratigraphy requires a solid understanding of sequence architecture along the steep margins of rift basins. This study analyzes an Eocene lacustrine sequence along the steep margin of the Dongying depression in eastern China through integrated core, well-log, and three-dimensional seismic analyses. The lacustrine sequence is bounded by unconformities and their correlative conformities at the base and top and consists of three systems tracts, namely an early expansion systems tract (EEST), late expansion–early contraction systems tract (LEECST), and late contraction systems tract (LCST), which record a lake expansion–contraction cycle. These systems tracts differ in thickness and development of depositional systems. The EEST is the thickest and contains well-developed marginal and basinal fan systems with an overall retrogradational stacking pattern. The well-developed fan systems are the most striking features within the sequence. The LEECST is the most widespread and contains dominantly profundal–sublittoral deposits. The LCST is the thinnest, with poorly developed fan systems, and is characterized by significant erosion by fluvial incision. The variable thickness and development of depositional systems in the three systems tracts are the responses to the interplay of sediment supply and accommodation space. Accommodation space establishes the framework for sedimentary infill, and sediment supply determines spatial distribution and temporal evolution of depositional systems within each systems tract. This study provides a lake expansion–contraction scheme to divide a lacustrine stratigraphic sequence into systems tracts and highlights the feasibility of applying this approach in studying sequence stratigraphy along the steep margin of a lacustrine rift basin. The results also provide understandings for the development, distribution, and evolution of depositional systems and their controlling factors along the steep margin of other rift basins in the world.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The purpose of this work is to identify genetic affinities among 48 crude oil samples from the onshore and offshore Santa Maria basins. A total of 21 source-related biomarker and stable carbon isotope ratios among the samples were assessed to assure that they were unaffected by secondary processes. Chemometric analysis of these data identifies six oil families with map and stratigraphic distributions that reflect organofacies variations within the Miocene Monterey Formation source rock. The data comprise a training set that was used to create a chemometric decision tree to classify newly collected oil samples. Three onshore families originated from two synclines, which may contain one or more pods of thermally mature source rock. Multiple biomarker parameters indicate that the six oil families achieved early oil window maturity in the range of 0.6%–0.7% equivalent vitrinite reflectance. The offshore oil samples consist of one family from Point Pedernales field and two families from the “B” prospect. Geochemical characteristics of these families indicate origins under differing water column and sediment oxicity and carbonate versus siliceous and detrital input in ‘carbonate,’ ‘marl,’ and ‘shale’ organofacies like those in the lower calcareous–siliceous, carbonaceous marl, and clayey–siliceous members of the Monterey Formation elsewhere in coastal California. The corresponding lithofacies and organofacies appear to be linked to the early–middle Miocene climate optimum and subsequent paleoclimatic cooling after circa 14 Ma, a systematic up-section increase in the stable carbon isotope composition of related oil samples, decreased preservation of calcium carbonate shells from planktic foraminifera and coccoliths, and increased preservation of clay-sized siliceous shells of diatoms and radiolarians. The results show that organofacies within the Monterey source rock are responsible for many of the geochemical differences between the oil families. This paleoclimate–organofacies model for crude oil from the Monterey Formation can be used to enhance future exploration efforts in many areas of coastal California.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Twenty-four oils produced from the Woodford Shale and overlying Mississippian strata in central Oklahoma were characterized geochemically to determine their possible source(s). The 168 core samples from the Woodford and Mississippian sections of 14 wells in central Oklahoma were initially characterized by total organic carbon (TOC), Rock-Eval, and vitrinite reflectance, and select samples (TOC 〉 1.0 wt. %) were subjected to biomarker analyses to characterize source input, depositional environment, maturity, and oil-to-source rock correlations. Thermal maturity parameters indicate the Woodford Shale is immature to marginally mature in Payne County, Oklahoma, and shows a progressive increase in maturity toward the southwest. Close to the Nemaha uplift, the Woodford is in the main stage of oil generation. It is proposed that the oils in this area have three possible origins: (1) Oils produced from the Woodford and overlying Mississippian strata have similar fingerprints, suggesting the Woodford Shale and overlying Mississippian strata are in communication; (2) oils produced near the Nemaha uplift (Logan and western Payne Counties) were sourced from the Woodford but had a significant Mississippian source contribution based on source-specific biomarkers; (3) oils east of the Cherokee platform (eastcentral Payne County) share strong Woodford source characteristics, and they were not generated in situ from the immature Woodford Shale but probably migrated from the Woodford Shale in the deeper part of the Anadarko Basin in southern Oklahoma. These results are consistent with the findings that indicate abundant marine coarse-grained biogenic silica (radiolarian-rich) chert facies found in eastcentral Payne County may contribute to good reservoir petrophysical properties, suggesting the Woodford Shale may not be a source in this area but simply a tight reservoir.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Primary depositional mineralogy has a major impact on sandstone reservoir quality. The spatial distribution of primary depositional mineralogy in sandstones is poorly understood, and consequently, empirical models typically fail to accurately predict reservoir quality. To address this challenge, we have determined the spatial distribution of detrital minerals (quartz, feldspar, carbonates, and clay minerals) in surface sediment throughout the Ravenglass Estuary, United Kingdom. We have produced, for the first time, high-resolution maps of detrital mineral quantities over an area that is similar to many oil and gas reservoirs. Spatial mineralogy patterns (based on x-ray diffraction data) and statistical analyses revealed that estuarine sediment composition is primarily controlled by provenance (i.e., the character of bedrock and sediment drift in the source area). The distributions of quartz, feldspar, carbonates, and clay minerals are controlled by a combination of the grain size of specific minerals (e.g., rigid vs. brittle grains) and estuarine hydrodynamics. The abundance of quartz, feldspar, carbonates, and clay minerals is predictable as a function of depositional environment and critical grain-size thresholds. This study may be used, by analogy, to better predict the spatial distribution of sandstone composition and thus reservoir quality in ancient and deeply buried estuarine sandstones.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Petroleum types in the Eagle Ford resource play span the range from black oil to dry gas and are produced along regional trends that are largely maturity controlled. A total of 61 shale samples covering all maturity zones were evaluated to document organic richness, organic matter type, and maturation characteristics using established geochemical parameters. Pyrolysis experiments were then performed to simulate the generation of petroleum fluids. Termed the “PhaseSnapShot” approach, one or more target wells with known fluid properties were used as reference; a match with that composition was made using next-formed fluids generated from the shale in a closely located well of slightly lower thermal maturity than the target well(s). Phase behavior predictions from the model were calibrated using a regional pressure–volume–temperature (PVT) database compiled from the public domain. The conceptual model that best matched the PVT data were comprised of two reactive components: (1) a mixture of kerogen and bitumen that generated petroleum within the low permeability shale matrix and (2) bitumen in zones of enhanced porosity within the matrix. The combined generation of gas from both of these components as well as the strong retention of C〈sub〉7+〈/sub〉 fluids in the matrix during production were required to match the calibration data. Retention of oil was needed over a broad thermal maturity range (Rock-Eval 〈span〉Tmax〈/span〉 release: 440°C –475°C). A key result of this forward model is that phase behavior and bulk compositional properties of hydrocarbons can be quickly and effectively predicted using mature shale samples as long as calibration data from PVT reports are available.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Mapping of seismic and lithological facies is a very complex process, especially in regions with low seismic resolution caused by extensive salt layers, even when only an exploratory view of the distribution of the reservoir facies is required. The aim of this study was to apply multi-attribute analysis using an unsupervised classification algorithm to map the carbonate facies of an exploratory presalt area located in the Outer high region of the Santos Basin. The interval of interest is the Barra Velha Formation, deposited during the Aptian, which represents an intercalation of travertines, stromatolites, grainstones and spherulitic packstones, mudstones, and authigenic shales, which were deposited under hypersaline lacustrine conditions during the sag phase. A set of seismic attributes, calculated from a poststack seismic amplitude volume, was used to characterize geological and structural features of the study area. We applied k-means clustering in an approach for unsupervised seismic facies classification. Our results show that at least three seismic facies can be differentiated, representing associations of buildup lithologies, aggradational or progradational carbonate platforms, and debris facies. We quantitatively evaluated the seismic facies against petrophysical properties (porosity and permeability) from available well logs. Seismic patterns associated with the lithologies helped identify new exploration targets.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Dibei gas field is a large tight gas field located in the Kuqa subbasin, Tarim Basin, northwestern China. The reservoir is within the Lower Jurassic Ahe Formation (J〈sub〉1〈/sub〉a) and has porosity and permeability ranges of 2%–8% and 0.01–1 md, respectively. Two episodes of hydrocarbon charge are identified based on a detailed study of fluid-inclusion petrography and microthermometry, fluorescence spectroscopy characteristics, and the thermal maturity of both gas and light oil. Low-maturity oil as represented by hydrocarbon inclusions with yellow-green fluorescence entered the reservoir circa 23–12 Ma, whereas high-maturity hydrocarbons, as indicated by hydrocarbon inclusions with blue-white fluorescence, have charged the reservoir since 5 Ma. The hydrocarbon charge process combined with porosity evolution determined the present gas–water distribution characteristics in the Dibei gas field. Porosity in the J〈sub〉1〈/sub〉a sandstone reservoir was relatively high during the first episode of hydrocarbon charge, which allowed oil to migrate upward and accumulate in structural highs under buoyancy. From 5 Ma to the present, the Dibei gas field experienced strong tectonic compression associated with intense thrust-fault reactivation, causing deformation and oil leakage from the reservoir. Continuous tight sand deposits along the slope areas, located far away from the active faults, became favorable accumulation sites for gas derived from the underlying Triassic source rocks. Hydrocarbon accumulation along the slope area in the Ahe Formation is dominantly controlled by equilibrium between hydrocarbon-generation pressure and capillary pressure.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉This study examines the influences on fluid flow within a shale outcrop where the networks of two distinct paleoflow episodes have been recorded by calcite-filled veins and green alteration halos. Such direct visualization of flow networks is relatively rare and provides valuable information of fluid-flow behavior between core and seismic scale.Detailed field mapping, fracture data, and sedimentary logging were used over a 270 m〈sup〉2〈/sup〉 (2910 ft〈sup〉2〈/sup〉) area to characterize the paleo–fluid-flow networks in the shale. Distal remnants of turbidite flow deposits are present within the shale as very thin (1–10 mm [0.04–0.4 in.]) fine-grained sandstone bands. The shale is cut by a series of conjugate faults and an associated fracture network, all at a scale smaller than seismic detection thresholds. The flow episodes used fluid-flow networks consisting of subgroups of both the fractures and the thin turbidites. The first fluid-flow episode network was mainly comprised of thin turbidites and shear fractures, whereas the network of the second fluid-flow episode was primarily small joints (opening mode fractures) connecting the turbidites.The distribution of turbidite thicknesses follows a negative exponential trend. which reflects the distribution of thicker turbidites recorded in previous studies. Fracture density varies on either side of faults and is highest in an area between closely spaced faults. Better predictions of hydraulic properties of sedimentary-structural networks for resource evaluation can be informed from such outcrop subseismic scale characterization. These relationships between the subseismic features could be applied when populating discrete fracture networks models, for example, to investigate such sedimentary-structural flow networks in exploration settings.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The middle Cambrian Maryville–Basal sands in the interval of 4600–4720 ft (1402.1–1438.7 m) in the Kentucky Geological Survey 1 Hanson Aggregates well (i.e., muddy sandstones separated by sandy mudstones) were evaluated to determine effective porosity (ϕ〈sub〉〈span〉e〈/span〉〈/sub〉), clay volume (〈span〉Vc〈/span〉), and supercritical CO〈sub〉2〈/sub〉 storage capacity. Average porosity and permeability measured in core plugs were 8.71% porosity and 2.17 md permeability in the Maryville sand and 10.61% porosity and 15.79 md permeability in the Basal sand. The ϕ〈sub〉〈span〉e〈/span〉〈/sub〉 and 〈span〉Vc〈/span〉 were calculated from the density log using a multiple-matrix shaly sand model to identify four formation lithologies: muddy sandstone, sandy mudstone, dolomitic mudstone, and dolomitic claystone. Average ϕ〈sub〉〈span〉e〈/span〉〈/sub〉 and 〈span〉Vc〈/span〉 calculated in the Maryville sand were 8.9% and 35.3%, respectively, and an average of 8.7% and 41.2% in the Basal sand, respectively. Calculated ϕ〈sub〉〈span〉e〈/span〉〈/sub〉 exhibits a good match with porosity measured in core plugs. Prior to step-rate testing, static reservoir pressure was 2020 psi (13.9 MPa), representing a 0.435 psi/ft (9.8 kPa/m) hydrostatic gradient, which is consistent with other underpressured reservoirs in Kentucky. The interval fractured at 2698 psi (18.0 MPa), yielding a fracture gradient of 0.581 psi/ft (12.7 kPa/m). Pressure falloff analysis suggests a dual-porosity/dual-permeability reservoir consistent with core data. Estimated 50th percentile supercritical CO〈sub〉2〈/sub〉 storage volume supercritical CO〈sub〉2〈/sub〉 storage volume, using 7% porosity cutoff for determining net reservoir volume, is 0.538 tons/ac (1.33 t/ha). Thin reservoir sands, low porosity and permeability, and low fracture gradient, however, preclude the Maryville–Basal sands as large-volume deep-saline CO〈sub〉2〈/sub〉 storage reservoirs in this area.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Porosity–permeability transforms were generated using an extensive data set covering two oil-bearing formations in Ohio: the Clinton Sandstone in eastern Ohio and the Copper Ridge Dolomite in central Ohio. The reservoirs were selected because of their historical importance as oil producers and their potential as targets for CO〈sub〉2〈/sub〉 use for enhanced oil recovery and associated geological storage. The porosity-permeability transforms generated in this study have coefficients of determination that are nearly double those in the published literature. Methods applying other information (e.g., lithofacies type and reservoir depth) to improve the transforms are also discussed. Ultimately, it was determined that although subdividing the Clinton Sandstone data by geologically similar areas constrained the porosity and permeability values, the data for most areas were too limited to yield robust correlations. Thus, the range of possible outcomes should be determined using the transform derived from all available data. The Copper Ridge values were largely not constrained when subdivided by depth.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The Pennsylvanian–Middle Triassic Cooper Basin is Australia’s premier conventional onshore hydrocarbon-producing province. The basin also hosts a range of unconventional gas play types, including basin-centered gas and tight gas accumulations, deep dry coal gas associated with the Patchawarra and Toolachee Formations, and the Murteree and Roseneath shale gas plays.This study used petroleum systems analysis to investigate the maturity and generation potential of 10 Permian source rocks in the Cooper Basin. A deterministic petroleum systems model was used to quantify the volume of expelled and retained hydrocarbons, estimated at 1272 billion BOE (512 billion bbl and 760 billion BOE) and 977 billion BOE (362 billion bbl and 615 billion BOE), respectively. Monte Carlo simulations were used to quantify the uncertainty in volumes generated and to demonstrate the sensitivity of these results to variations in source-rock characteristics.The large total generation potential of the Cooper Basin and the broad distribution of the Permian source kitchen highlight the basin’s significance as a world-class hydrocarbon province. The large disparity between the calculated volume of hydrocarbons generated and the volume so far found in reservoirs indicates the potential for large volumes to remain within the basin, despite significant losses from leakage and water washing. The hydrocarbons expelled have provided abundant charge to both conventional accumulations and to the tight and basin-centered gas plays, and the broad spatial distribution of hydrocarbons remaining within the source rocks, especially those within the Toolachee and Patchawarra Formations, suggests the potential for widespread shale and deep dry coal plays.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The strike-slip fault systems in the central Tarim Basin, China, afford an exceptional opportunity to document the structural characteristics and evolution process of small displacement intracratonic strike-slip faults using three-dimensional seismic reflection data. These strike-slip faults display subvertical segments at depth and en echelon normal fault zones where relatively shallow. Fault segmentation and flower structures can be commonly observed in plan view and cross-section view, respectively.Consistent with the notion that segment coalescence is the fundamental process for fault evolution, the mean segment length of representative strike-slip faults examined in this study is positively correlated to the measured fault offset. The width of the en echelon normal fault zone is positively correlated with the estimated maximum overburden thickness. The integrated data sets suggest that the evolution of the conjugate fault array followed a sequential evolution process instead of forming simultaneously. The switch in slip direction of the master fault of the conjugate fault array is attributed to the change of stress orientation. Regarding individual strike-slip faults, increase in displacement induces the formation of faults with lower fault-array angles linking initially formed en echelon normal faults. In cross sections, throughgoing fault surfaces can also form, connecting the lower subvertical fault segment and the upper en echelon normal faults.The presented data sets and evolution models established in this study can be used as tools to better predict the structural attributes of subsurface strike-slip fault systems with important consequences for reservoir formation and hydrocarbon accumulation in the Tarim Basin in particular, and in ancient marine basins in general.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Static formation temperature (〈span〉SFT〈/span〉) can be estimated from temperatures measured during wire-line logging (〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉). A large number of correction models for obtaining 〈span〉SFT〈/span〉 from 〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉 have been suggested. Several studies have shown that 〈span〉SFT〈/span〉s yielded by such models are off by an average of 6°C–10°C (43°F–50°F) at burial depths of 1.5–3.5 km (0.9–2.2 mi) and thus have the potential to cause serious issues in thermal and hydrocarbon generation models. This paper explores the causes for erroneous 〈span〉SFT〈/span〉 predictions generated from 〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉 measurements and identifies factors that should be addressed to generate a globally applicable correction model. We also present an improved empirical correction model for 〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉 data from eight oil and gas fields, located on the Norwegian continental shelf. The new empirical model was designed to give correct average 〈span〉SFT〈/span〉 predictions and is applicable to single 〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉 measurements. It has been validated against temperatures recorded during drill-stem testing, which closely represent local 〈span〉SFT〈/span〉s. The expression yields improved results compared with other correction models applied to the data set. However, the average error in computed 〈span〉SFT〈/span〉 values varies by up to 10°C (18°F) between the investigated hydrocarbon fields. We conclude that these variations result from differences in operational practices such as fluid circulation and drilling velocities. Therefore, current empirical and physical models for 〈span〉SFT〈/span〉 prediction from 〈span〉T〈/span〉〈sub〉〈span〉m〈/span〉〈/sub〉 require local calibration. It is also suggested that more accurate compilations and analyses of operational data could lead to improved and more globally applicable models.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The petroleum system concept spans the spatial and temporal extent of all elements and processes required for the generation and preservation of petroleum. The critical moment of a petroleum system is the moment with the highest probability for the generation–migration–accumulation of hydrocarbons. It is an important concept in petroleum exploration risk assessment because the stratigraphic and geographic extents of a petroleum system are determined at the critical moment. In petroleum systems, thermal history data, burial history data, and vitrinite reflectance data may be unavailable, unreliable, or incomplete; this introduces significant uncertainty in the choice of the critical moment. We present here a quantitative probabilistic framework for estimating the critical moment and quantifying the associated uncertainty in such cases. We define a probabilistic early bound and late bound for the critical moment (which, combined together, we term the critical range) and then estimate the moment with the highest numerical probability of generation–migration–accumulation. We define the uncertainty associated with the critical moment as half the absolute value of the critical range. In cases with little ambiguity or duplicity in the timing of petroleum system elements and processes, the critical range converges to one point, which is also the critical moment. The probabilistic framework introduces consistency to the critical moment estimation problem and quantifies the level of uncertainty in the estimation. This reduces the risk involved in petroleum exploration assessment.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Very limited literature is available relating to gas production from ultradeep (〉9000 ft [〉2700 m]) coal seams. This paper investigates permeability enhancement in ultradeep coal seams of the late Carboniferous and early Permian to Late Triassic Cooper Basin in central Australia, using a time-lapse pressure transient analysis (PTA) approach for a pilot well. The gas production history and three extended shut-in periods are used to construct the time-lapse PTA for the study well. A new approach is introduced to construct a permeability ratio function. This function allows the calculation of permeability change resulting from competition between the compaction and coal-matrix shrinkage effects.Pressure transient analysis indicates that gas flow is dominated by a bilinear flow regime in all extended pressure buildup tests. Hence, reservoir depletion is restricted to the stimulated area near the hydraulic fracture. This implies that well-completion practices that create a large contact area with reservoirs, such as multistage hydraulically fractured horizontal wells, may be required for achieving economic success in these extremely low-permeability reservoirs. The permeability ratio is constructed using the slope of the straight lines in bilinear flow analysis. Because of uncertainty in average reservoir pressure, probabilistic analysis is used and a Monte Carlo simulation is performed to generate a set of possible permeability ratio values. The permeability ratio values indicate that coal permeability has increased during the production life of the wellbore because of the coal-matrix shrinkage effect. Permeability enhancement in this ultradeep coal reservoir has offset the effect of permeability reduction caused by compaction, which is beneficial to gas production.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Nanometer to micrometer mica and illite separates of indurated Cambrian and Ordovician oil-bearing sandstones from the Hassi Messaoud field (Algeria) were extracted, x-rayed, observed by scanning and transmission electron microscopy, and K-Ar dated. Electron microscope observations revealed typical euhedral shapes for the mica to illite particles of most size fractions; almost no odd-shaped detrital crystals were detected. The combined results document several generations of mineralogical and morphological identical mica to illite crystals that could not be differentiated by the traditional identification methods. Illite and mica genesis was multiphased with crystallization episodes at 340 ± 10 (ca. Middle Mississippian), 280 ± 10 Ma (ca. early Permian), and 170 ± 10 Ma (ca. Middle Jurassic). Younger than the stratigraphic age of the host rocks, which is incompatible with a detrital origin, the two older mica ages confirm that the hydrocarbon generation and emplacement had to start after the Variscan tectonothermal event and before exhumation of the meta-sediments. The younger K-Ar ages at 135 to 110 Ma (ca. Early Cretaceous) relate to further crystallization episodes, whereas those at circa 295, 265, and 210 Ma probably correspond to variable mixtures of the older and younger mica to illite end-members. Three average K-Ar values are statistically significant: the oldest at 340 ± 10 Ma corresponds to the start of the Variscan tectonic activity, and the intermediate at 280 ± 10 Ma sets its end, both episodes probably modifying the reservoir capacities of the potential hydrocarbon host rocks. The ages at 170 ± 10 Ma identify a further diagenetic activity characterized by illitization of dickite-type precursors in local reservoirs. These younger ages could correspond to the hydrocarbon charge into reservoirs, which stopped diagenetic illitization at a present-day depth of approximately 4000 m (∼13,000 ft).〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉This paper analyzes regional hydrogeological conditions and divides the study area into three hydrogeological types and seven hydrogeological units, to investigate hydrogeology and its effect on coalbed methane (CBM) enrichment in the southern Junggar Basin, China. From this work, it is found that the groundwater flow paths in the study area are the joint effects of south-to-north and west-to-east flows. This study also shows that microbial gases are widely developed, although the depth limit of microbial gas occurrence is still unclear in the study area. Microbial CO〈sub〉2〈/sub〉 reduction is the leading formation path in the study area, except for the Houxia region, where fermentation is the formation mechanism. The abnormally high CO〈sub〉2〈/sub〉 in stagnant zones (i.e., water flow is slow and stagnant) is mainly associated with methanogenesis, whereas relatively low CO〈sub〉2〈/sub〉 (microbial or thermogenic) is present where water flow is active. The average CBM content within the Xishanyao Formation changes within various hydrogeological units; moreover, the average CBM content within the Badaowan Formation of the same hydrogeological unit (e.g., Fukang) suggests that the hydrogeological and CBM enrichment conditions are different within various structural types. Overall, the hydrogeological conditions exert control on the gas content in the study area; that is, the gas content is high in stagnant zones. Finally, influenced by supplemental microbial gases, changes in the CBM oxidation zone are relatively complex in the study area, the depth of which has no obvious correlation with hydrogeological conditions and changes significantly from west to east.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉How and when sediment moves from terrestrial sources to deep-water sinks is a significant area of research. We have used an array of seismic, borehole, and gravity core data sets to explore the timing and magnitude of sediment-routing to Pearl River slope over the last 478 k.y. As predicted by existing sequence stratigraphic models, most sediment dispersal to deep water is shown to have occurred during glacial sea-level falls; however, clastic detritus was still being transported into deep water during interglacial sea-level rises. We suggest that sediment routing to deep water during interglacial sea-level rise is caused by summer monsoon strengthening and resultant warmer and wetter climates, both of which have enhanced effective precipitation and sediment supply. Although some models for the delivery of sediment to deep-water basins stress the importance of proximity of canyon heads and coeval shorelines, we observed that sediment routing to deep water could occur regardless of the distance between channel head and coeval shorelines. In the present case, the success of delivery is related to the combined effects of (1) the short duration and high amplitude of sea-level oscillations during the past 478 k.y. and (2) the enhanced sediment supply caused by more humid climates and greater temperature difference between glacial and interglacial period. This hypothesis is supported by (1) observations that outer Pearl River deltas prograded as an apron over preexisting shelf edges for 10–15 km (6–9 mi) and (2) the occurrence of slope channels extending back to prodelta reaches of Pearl River shelf-edge deltas.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The discovery of carbonate gas fields in the Middle Triassic Leikoupo Formation of the Sichuan Basin has a complex history. In recent years, a series of structural fields have been discovered in the western Sichuan Basin. Their discovery confirms the immense exploration potential of the Leikoupo Formation. In this study, we analyze the characteristics of Leikoupo Formation exploration plays using exploration wells and test data, aiming to provide a reference for further discoveries. The Leikoupo Formation represents the uppermost unit in the Sichuan marine carbonate platform succession. During its deposition, the whole basin was characterized by a restricted and evaporitic platform. Two classes of reservoirs developed. One is pore–fracture reservoirs, in marginal platform and intraplatform shoals, and another is fracture–vug reservoirs in the karstic weathering crust of the formation-capping unconformity. Three hydrocarbon accumulation models were established for the Leikoupo Formation based on the spatial and temporal relationship among the source, reservoir, and cap rocks. Two types of exploration plays are present in the Leikoupo Formation, that is, shoal (including intraplatform shoal and marginal platform shoal) dolomite plays and karstic dolomite weathering crust plays (including intraplatform shoal karst and marginal platform shoal karst). The western Sichuan depression in the karstic slope belt presents immense exploration potential because of a proximal hydrocarbon supply, charging via an extensive fracture network, shoals and karstic reservoir, a good seal rock of terrestrial mudstone, and potential composite hydrocarbon accumulations in stratigraphic traps, making it a promising area for future exploration.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Carbonaceous debris (CD) within uranium-bearing strata has been studied in the Daying uranium deposit of the northern Ordos Basin, northern China. The influence of radiogenic heat from uranium on organic matter maturation was investigated through a series of tests including measurements of vitrinite reflectance (〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉), fission-track (FT) analysis in quartz grains, and the calculation of the radiogenic heat production rate of the samples. The results show that 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉 in uranium-bearing strata generally increases as the burial depth increases, indicating that CD experienced normal burial coalification. However, 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉 values of the samples rich in uranium are 0.062% 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉 higher than those without uranium mineralization. Vitrinite reflectance bears a positive relationship with uranium content, and an inverse relationship with distance to the closest sandstone rich in uranium, indicating that uranium enrichment enhances organic matter maturation. The production of uranium decay makes FT observable in quartz grains, and the intensity of decay increases with proximity to the uranium ore body. The calculated radioactive heat production rate from the uranium ore body is 6.857 × 10〈sup〉−5〈/sup〉 W/m〈sup〉3〈/sup〉. During the long-term stable decay, as the uranium ore body theoretically results in an abnormal increase in temperature of 52°C without consideration of the loss of heat conduction, heat convection, and thermal radiation, this would yield a theoretical 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉 increase of 0.209% 〈span〉R〈/span〉〈sub〉〈span〉o〈/span〉〈/sub〉, reasonably greater than the observed. Therefore, the long-term stable radiogenic heat produced by uranium ore body can slightly enhance organic matter maturation, which is instructive in uranium prospecting.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The purpose of this study is to deconstruct the relationship between the Leaf River anticline and the preglacial bedrock paleotopography at the eastern terminus of the Plum River Fault Zone in Ogle County, Illinois, using a geostatistical approach. The contour maps derived from the elevation models provided detailed depictions of the ancient bedrock landscape and subsurface structure in the study area. The Leaf River anticline is interpreted to be a component of hanging-wall anticline at the terminus of the Plum River Fault Zone. The topographic high created by the anticline controlled local drainage and led to the development of the Leaf River paleovalley prior to the Pleistocene. The catastrophic failure of an ice damn during the Illinois glacial episode carved a glacial spillway into the north flank of the Leaf River anticline that interfaced with a tributary of the Leaf River paleovalley. This rerouted the preglacial drainage network and permanently diverted the ancient Rock River to its modern-day position. Ultimately, the subsurface geometry of the Leaf River anticline and its relationship to the local bedrock paleotopography were revealed by the elevation models. The position and development of the Leaf River paleovalley and glacial spillway interpreted in this study aligned with the regional interpretations for the evolution of the ancient bedrock landscape established in prior works. However, this study revealed that the Leaf River anticline and, by association, the terminus of the Plum River Fault Zone extend farther east into the region than indicated by prior works.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉In the last decade, production of shale gas has tremendously increased, and the need for local pre-exploitation baseline data on dissolved natural gas in aquifers has been stressed. This study investigated the origin of hydrocarbons naturally present in shallow aquifers of the Saint-Édouard area (Québec, eastern Canada), where the underlying Utica Shale is known to contain important gas resources that have not yet been exploited. Groundwater and shallow bedrock gas samples were collected and analyzed for isotopic composition of alkanes (δ〈sup〉13〈/sup〉C and δ〈sup〉2〈/sup〉H〈sub〉C1–C3〈/sub〉), dissolved inorganic carbon (δ〈sup〉13〈/sup〉C〈sub〉DIC〈/sub〉), and radiocarbon in methane and DIC (〈sup〉14〈/sup〉C〈sub〉DIC〈/sub〉, 〈sup〉14〈/sup〉C〈sub〉CH4〈/sub〉). This multi-isotope approach proved enlightening, and results revealed that (1) most of the methane in the region is of microbial origin; (2) partial contribution of thermogenic gas occurs in 15% of the wells; (3) processes such as late-stage methanogenesis and methane oxidation are responsible for ambiguous methane isotopic compositions; and (4) both microbial and thermogenic gas originate from the shallow bedrock aquifer, with the exception of one sample likely coming from deeper units. The thick succession of shales overlying the Utica Shale thus appears to act as an effective migration barrier for the shallow aquifers. However, evidence of upward migration of old brines near major fault zones indicates that these may serve as a preferential migration pathway over a certain depth but most likely no more than approximately 200–500 m (∼650–1640 ft). The geochemical framework presented here will hopefully be useful in other research projects, especially when conventional indicators of natural gas origin provide ambiguous results.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Urbanization modifies the natural water cycle. In this study, a weighted-rating multicriteria analysis was adopted to quantify the runoff index and to assess the impact of urbanization on the water cycle. The considered parameters are (1) slope, (2) permeability of soil, and (3) rainfall. Using the land use map, a runoff risk map was established. The approach was applied to Manouba catchment. The main results revealed that between 2004 and 2014, the area with a high runoff index increased from 32% to 39%. The runoff risk increased; in 2004, the high class covered 18% of the watershed area. This value became 30% in 2014. Results demonstrate that urbanization affects hydrological processes. This method is appropriate in other similar watersheds to estimate the runoff index.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉To better understand controls on the origin and evolution of brackish groundwater, the hydrogeochemistry of brackish groundwaters was studied within the Triassic Dockum Group across the Midland Basin in Texas. The suitability of Dockum Aquifer water for use in hydraulic fracturing fluid was examined because the area overlies the largest and most productive tight oil province in the United States. Groundwater generally flows southward and eastward across the basin. Transmissivities indicate that water yield from the Dockum Aquifer is mixed. Higher salinity (up to ∼100 g/L), group I water is found mainly in the center and western parts of the basin; chemistry of these meteoric waters is controlled by water–rock interaction with salinity increasing along its flow path via dissolution of halite and anhydrite, followed by salinity-enhanced carbonate dissolution and/or cation release from clays. Along the down-gradient basin margins, lower salinity (〈7.5 g/L), group II waters of various ion compositions are more commonly found. Group II waters are also meteoric but from local recharge including downward flow from the Edwards–Trinity or other aquifers. Despite having lower salinity, the water in the down-gradient southern and eastern margins of the basin can exceed acceptable SO〈sub〉4〈/sub〉 limits for cross-linked gel fluids. Generally, the majority of the water in the basin is suitable for use with slick-water hydraulic fracturing. Findings from this research provide important information on the complex controls on the chemistry of brackish groundwater and their potential beneficial uses in the oil and gas industry.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉Groundwater is the major source of drinking water in both urban and rural India. Estimation of natural groundwater recharge is essential for the sustainable development of groundwater. Natural recharge was estimated by various methods, such as the water level fluctuation method, water balance method, linear regression model, and nonlinear regression model. The recharge estimates by the water balance method was compared with the recharge obtained from the water level fluctuation method for the study area and found to be in good agreement.Estimation of recharge by the water level fluctuation method is laborious, and envisaging the difficulties in the availability and reliability of data, the water balance method is taken as the standard for developing regression equations in the present study. Simpler linear and nonlinear regression models were developed for the study area to estimate natural recharge by correlating with the water balance model. The models were calibrated with 10-yr data and validated with 5-yr data. The statistical analysis showed that no significant difference exists between the recharge estimate by the water balance method and the two estimates of natural recharges, such as linear regression and nonlinear regression models. The average recharge percentages from the water level fluctuation method, water balance method, linear regression model, and nonlinear regression model are 15.09%, 14.92%, 14.62%, and 14.57%, respectively, for the watershed during the study period. The study proves that regression equations can be efficiently used in recharge computation with proper calibration for ungauged basins, and laborious data-intensive computation methods can be eliminated.〈/span〉

Description: 〈span〉〈div〉ABSTRACT〈/div〉The White River watershed encompasses four major tributaries within a basin area of 130 km〈sup〉2〈/sup〉 (1595 mi〈sup〉2〈/sup〉) in extreme northwestern Nebraska. An examination of the historical (1968–1975) aqueous geochemistry data (major cations and anions and total dissolved solids [TDS]) supplied by the Nebraska Department of Environmental Quality revealed that the TDS is relatively low (130–1200 mg/L), excluding Big Cottonwood Creek (BCC), with a basin-wide median of 340 mg/L. The median TDS for the BCC is 1880 mg/L (brackish); the median values for Na and SO〈sub〉4〈/sub〉 are 385 and 897 mg/L, respectively. Mineralization in the river increases steadily downstream. The scatter plots of meq/L concentrations for selected anions and cations reveal the impact of silicate mineral (e.g., feldspar) weathering on the aqueous geochemistry throughout the watershed. These ubiquitous feldspar minerals most likely originated along the eastern slope of the Front Range during the Late Cretaceous and Tertiary (Laramide orogeny). Twenty-nine samples for three White River stations and the BCC exceed the US Environmental Protection Agency secondary maximum contaminant levels for TDS and/or SO〈sub〉4〈/sub〉 in drinking water supplies at 500 and 250 mg/L, respectively. Uncontaminated streams that drain marine shales (typically containing S-bearing minerals) nationwide typically show an excess of Na and a deficiency of Ca and Mg. This is due in part to cation exchange of Ca in solution for Na on clay minerals. Consequently, the weathering of shale terrains commonly produces an Na-SO〈sub〉4〈/sub〉 brackish surface-water runoff as is the case with BCC, which drains the Pierre Shale hills.〈/span〉