Showing papers in "AAPG Bulletin in 2009"

Pore-throat sizes in sandstones, tight sandstones, and shales

Abstract: Pore-throat sizes in siliciclastic rocks form a continuum from the submillimeter to the nanometer scale. That continuum is documented in this article using previously published data on the pore and pore-throat sizes of conventional reservoir rocks, tight-gas sandstones, and shales. For measures of central tendency(mean,mode,median),pore-throatsizes(diameters) are generally greater than2 mm in conventionalreservoir rocks, range from about 2 to 0.03 mm in tight-gas sandstones, and rangefrom0.1to0.005 mminshales.Hydrocarbonmolecules, asphaltenes, ring structures, paraffins, and methane, form another continuum, ranging from 100 A ˚ (0.01 mm) for asphaltenes to 3.8 A ˚ (0.00038 mm) for methane. The pore-throat size continuum provides a useful perspective for considering (1) the emplacement of petroleum in consolidated siliciclastics and (2) fluid flow through fine-grained source rocks now being exploited as reservoirs.

Three common source rock evaluation errors made by geologists during prospect or play appraisals

Abstract: Geologists are frequently called on to evaluate the source rocks associated with their exploration prospects or plays. The three most common questions asked and answered about the source rock during project reviews are What's the total organic carbon (TOC)?, What kerogen type does Rock-Eval indicate?, and What maturity level does the vitrinite reflectance data point to? The answers to these seemingly innocuous questions may, in fact, be providing a false sense of security about the source rock in question. Understanding how this line of questioning can lead you astray and make you the victim of the TOC myth (“If I have high TOC, I have a good source rock.”), the Rock-Eval fallacy (“The Rock-Eval data tell me what kind of kerogen is in my source rock.”), and the vitrinite reflectance deficiency (“Vitrinite reflectance will tell me if my source rock is generating.”) is important. Some of the solutions to these problems include fully integrating TOC and Rock-Eval data, supplementing Rock-Eval data with pyrolysis-gas chromatography, and using burial history diagrams to help interpret vitrinite reflectance. Harry Dembicki Jr. is a senior geological advisor for geochemistry in the Geological Technology Group at Anadarko Petroleum Corporation, where he provides technical support to both exploration and development teams. He holds a B.S. degree (1973) in geology from the State University College of New York at New Paltz and a Ph.D. (1977) in geology, with emphasis in organic geochemistry, from Indiana University. He previously worked as an organic geochemist for Conoco and Marathon.

Mechanical and fracture stratigraphy

Abstract: Using examples from core studies, this article shows that separate identification of mechanical stratigraphy and fracture stratigraphy leads to a clearer understanding of fracture patterns and more accurate prediction of fracture attributes away from the wellbore. Mechanical stratigraphy subdivides stratified rock into discrete mechanical units defined by properties such as tensile strength, elastic stiffness, brittleness, and fracture mechanics properties. Fracture stratigraphy subdivides rock into fracture units according to extent, intensity, or some other observed fracture attribute. Mechanical stratigraphy is the by-product of depositional composition and structure, and chemical and mechanical changes superimposed on rock composition, texture, and interfaces after deposition. Fracture stratigraphy reflects a specific loading history and mechanical stratigraphy during failure. Because mechanical property changes reflect diagenesis and fractures evolve with loading history, mechanical stratigraphy and fracture stratigraphy need not coincide. In subsurface studies, current mechanical stratigraphy is generally measurable, but because of inherent limitations of sampling, fracture stratigraphy is commonly incompletely known. To accurately predict fractures in diagenetically and structurally complex settings, we need to use evidence of loading and mechanical property history as well as current mechanical states.

Natural fracture characterization in tight gas sandstones: Integrating mechanics and diagenesis

Abstract: Accurate predictions of natural fracture flow attributes in sandstones require an understanding of the underlying mechanisms responsible for fracture growth and aperture preservation. Poroelastic stress calculations combined with fracture mechanics criteria show that it is possible to sustain opening-mode fracture growth with sublithostatic pore pressure without associated or preemptive shear failure. Crack-seal textures and fracture aperture to length ratios suggest that preserved fracture apertures reflect the loading state that caused propagation. This implies that, for quartz-rich sandstones, the synkinematic cement in the fractures and in the rock mass props fracture apertures open and reduces the possibility of aperture loss on unloading and relaxation. Fracture pattern development caused by subcritical fracture growth for a limited range of strain histories is demonstrated to result in widely disparate fracture pattern geometries. Substantial opening-mode growth can be generated by very small extensional strains (on the order of 104); consequently, fracture arrays are likely to form in the absence of larger scale structures. The effective permeabilities calculated for these low-strain fracture patterns are considerable. To replicate the lower permeabilities that typify tight gas sandstones requires the superimposition of systematic cement filling that preferentially plugs fracture tips and other narrower parts of the fracture pattern.

Joint sets that enhance production from Middle and Upper Devonian gas shales of the Appalachian Basin

Abstract: The marine Middle and Upper Devonian section of the Appalachian Basin includes several black shale units that carry two regional joint sets (J1 and J2 sets) as observed in outcrop, core, and borehole images. These joints formed close to or at peak burial depth as natural hydraulic fractures induced by abnormal fluid pressures generated during thermal maturation of organic matter. When present together, earlier J1 joints are crosscut by later J2 joints. In outcrops of black shale on the foreland (northwest) side of the Appalachian Basin, the east-northeast–trending J1 set is more closely spaced than the northwest-striking J2 set. However, J2 joints are far more pervasive throughout the exposed Devonian marine clastic section on both sides of the basin. By geological coincidence, the J1 set is nearly parallel the maximum compressive normal stress of the contemporary tectonic stress field (SHmax). Because the contemporary tectonic stress field favors the propagation of hydraulic fracture completions to the east-northeast, fracture stimulation from vertical wells intersects and drains J2 joints. Horizontal drilling and subsequent stimulation benefit from both joint sets. By drilling in the north-northwest–south-southeast directions, horizontal wells cross and drain J1 joints, whenever present. Then, staged hydraulic fracture stimulations, if necessary, run east-northeast (i.e., parallel to the J1 set) under the influence of the contemporary tectonic stress field thereby crosscutting and draining J2 joints.

Quantification of pore structure and its effect on sonic velocity and permeability in carbonates

Abstract: Carbonate rocks commonly contain a variety of pore types that can vary in size over several orders of magnitude. Traditional pore-type classifications describe these pore structures but are inadequate for correlations to the rock's physical properties. We introduce a digital image analysis (DIA) method that produces quantitative pore-space parameters, which can be linked to physical properties in carbonates, in particular sonic velocity and permeability. The DIA parameters, derived from thin sections, capture two-dimensional pore size (DomSize), roundness (), aspect ratio (AR), and pore network complexity (PoA). Comparing these DIA parameters to porosity, permeability, and P-wave velocity shows that, in addition to porosity, the combined effect of microporosity, the pore network complexity, and pore size of the macropores is most influential for the acoustic behavior. Combining these parameters with porosity improves the coefficient of determination (R2) velocity estimates from 0.542 to 0.840. The analysis shows that samples with large simple pores and a small amount of microporosity display higher acoustic velocity at a given porosity than samples with small, complicated pores. Estimates of permeability from porosity alone are very ineffective (R2 = 0.143) but can be improved when pore geometry information PoA (R2 = 0.415) and DomSize (R2 = 0.383) are incorporated. Furthermore, results from the correlation of DIA parameters to acoustic data reveal that (1) intergrain and/or intercrystalline and separate-vug porosity cannot always be separated using sonic logs, (2) P-wave velocity is not solely controlled by the percentage of spherical porosity, and (3) quantitative pore geometry characteristics can be estimated from acoustic data and used to improve permeability estimates.

Origin of overpressure and pore-pressure prediction in the Baram province, Brunei

Abstract: Accurate pore-pressure prediction is critical in hydrocarbon exploration and is especially important in the rapidly deposited Tertiary Baram Delta province where all economic fields exhibit overpressures, commonly of high magnitude and with narrow transition zones. A pore-pressure database was compiled using wireline formation interval tests, drillstem tests, and mud weights from 157 wells in 61 fields throughout Brunei. Overpressures are observed in 54 fields both in the inner-shelf deltaic sequences and in the underlying prodelta shales. Porosity vs. vertical effective stress plots from 31 fields reveal that overpressures are primarily generated by disequilibrium compaction in the prodelta shales but have been generated by fluid expansion in the inner-shelf deltaic sequences. However, the geology of Brunei precludes overpressures in the inner-shelf deltaics being generated by any conventional fluid expansion mechanism (e.g., kerogen-to-gas maturation), and we propose that these overpressures have been vertically transferred into reservoir units, via faults, from the prodelta shales. Sediments overpressured by disequilibrium compaction exhibit different physical properties to those overpressured by vertical transfer, and hence, different pore-pressure prediction strategies need to be applied in the prodelta shales and inner-shelf deltaic sequences. Sonic and density log data detect overpressures generated by disequilibrium compaction, and pore pressures are accurately predicted using an Eaton exponent of 3.0. Sonic log data detect vertically transferred overpressures even in the absence of a porosity anomaly, and pore pressures are reasonably predicted using an Eaton exponent of 6.5.

Fractures in sandstone reservoirs with ultra-low permeability: A case study of the Upper Triassic Yanchang Formation in the Ordos Basin, China

Abstract: The Upper Triassic Yanchang Formation in the Ordos Basin, central China, is a typical sandstone reservoir with an ultra-low permeability. High-angle tectonic fractures and diagenetic fractures, such as near-horizontal bedding fractures, intragranular fractures, and grain-boundary fractures, are abundant. Fractures are major pathways and enhance fluid flow in sandstone reservoirs with ultra-low permeability. Because of their weak lateral continuity and their small apertures under lithostatic pressure, bedding fractures make a relatively small contribution to the overall permeability of reservoirs. As they are both of small size and low permeability, intragranular fractures and grain-boundary fractures mainly improve the connectivity of reservoirs by connecting the matrix pores. High-angle tectonic fractures control the fluid movement in the ultra-low-permeability reservoirs. Under the effect of different successive tectonic stress fields, four assemblages of high-angle tectonic fractures developed in the sandstone reservoirs. Under the present-day stress, differently oriented fractures have different connectivities, apertures, and permeabilities. The northeast–southwest fractures parallel to the maximum principal stress have a good connectivity, large apertures, and a high permeability, forming the dominant flow paths. Knowledge of these paths can be used for optimizing well placement. Because of their sensitivity to pressure, fractures in different directions will show varying apertures as the formation pressure decreases. Therefore, the permeability of the fractures of different orientations and its impact on reservoirs vary at different developmental stages.

Porosity-preserving chlorite cements in shallow-marine volcaniclastic sandstones: Evidence from Cretaceous sandstones of the Sawan gas field, Pakistan

Abstract: Sandstones that have high porosity and permeability at great burial depth and high temperatures are of economic importance because a significant amount of hydrocarbons have been discovered in such reservoirs The Sawan gas field, with an expected ultimate recovery of more than 1 tcf, lies in the middle Indus Basin The reservoir rocks, Cretaceous volcaniclastic sandstones of the lower Goru Formation, show very high porosities at a reservoir temperature of 175C and depths of 3000 to 3500 m (9842 to 11,483 ft) The sandstones are mostly feldspathic litharenites Strongly altered volcanic rock fragments are the most important lithic component The clay fraction consists of Fe-rich chlorite (chamosite) and illite Diagenetic features such as compaction, quartz overgrowths, carbonate cements, and feldspar dissolution are observed The most distinguishing feature is a double layer of authigenic chlorite, lining the pores of the sandstones Chlorite additionally occurs as a pore-filling cement and as chloritized detrital components, all having similar chemical composition The pore-lining cement clearly developed in two stages: an earlier, poorly crystallized, and a later, better crystallized growth Missing rims at grain contacts show that precipitation occurred after an initial stage of compaction but early relative to other diagenetic phases Both chlorite rims grew by direct precipitation from pore waters, using products derived from volcanic rock fragments In areas with no, thin, or discontinuous chlorite rims, quartz cementation is common Well-developed chlorite rims inhibited quartz cementation, preserved porosities of up to 20%, and good permeabilities Porosity-preserving chlorite cementation in Sawan is restricted to sediments of a shallow-marine environment

Petroleum geology of the giant Elm Coulee field, Williston Basin

Abstract: The Elm Coulee field of the Williston Basin is a giant oil discovery in the middle Bakken Formation (Devonian–Mississippian) discovered in 2000. Horizontal drilling began in the field in 2000, and to date, more than 600 wells have been drilled. The estimated ultimate recovery for the field is more than 200 million bbl of oil (31.8 million m3). The Bakken Formation in the field area consists of three members: (1) upper shale, (2) middle silty dolostone, and (3) lower siltstone. The total Bakken interval ranges in thickness from 10 to 50 ft (3.1 to 15.3 m) over the field area. The upper shale is dark-gray to black, hard, siliceous, slightly calcareous, pyritic, and fissile. The shale consists of dark organic kerogen, minor clay, silt-size quartz, and some calcite and dolomite. The kerogen consists mainly of amorphous material, and the organic material is distributed evenly throughout the shale interval (not concentrated in laminations or lenses). The upper shale ranges in thickness from 6 to 10 ft (1.8 to 3.1 m) over the field area. The middle member consists of a silty dolostone and ranges in thickness from 10 to 40 ft (3.1 to 12.2 m). The lower member in the Elm Coulee field consists of brownish-gray, argillaceous, organic-rich siltstone. Burrowing and brachiopod fragments are common in the lower member. This facies is equivalent to the lower Bakken black shale facies on the northern side of the field and is interpreted to be an updip-landward equivalent to the deeper-water, black shale facies. The lower member ranges in thickness from 2 to 6 ft (0.61 to 1.8 m). Based on the abundance of fossil content and amount of burrowing, the members of the Bakken Formation are interpreted to have been deposited under aerobic (middle member, common burrows, and rare fossils), dsyaerobic (lower member, common fossils, and lesser amount of burrows), and anaerobic conditions (upper member, rare fossils, and burrows). The main reservoir in Elm Coulee field is the middle member, which has low matrix porosity and permeability and is found at depths of 8500 to 10,500 ft (2593 to 3203 m). The current field limits cover approximately 450 mi2 (1165 km2). The middle Bakken porosities range from 3 to 9%, and permeabilities average 0.04 md. Overall, the reservoir quality in the middle Bakken improves upward as the amount of clay matrix decreases. The middle Bakken is interpreted to be a dolomitized carbonate-shoal deposit based on subsurface mapping and dolomite lithology. The main production is interpreted to come from matrix permeability in the field area. Occasional vertical and horizontal fractures are noted in cores. The vertical pay ranges in thickness from 8 to 14 ft (2.4 to 4.3 m). The Bakken is slightly overpressured with a pressure gradient of 0.53 psi/ft (0.02 kPa/m). Horizontal wells are drilled on 640- to 1280-ac (259- to 518.4-ha) spacing units. Long single laterals, dual laterals, and trilaterals have all been drilled in the field. The horizontal intervals are sand-, gel-, and water-fracture stimulated. Initial production ranges from 200 to 1900 BOPD (31.8 to 302.1 m3/day). Initial potential rates for vertical wells are generally less than 100 BOPD (15.9 m3/day). The upper Bakken shale probably also contributes to the overall production in the field. The exact contribution is unknown but estimated to be less than 20% of the total production. The Elm Coulee field illustrates that the Bakken petroleum system has enormous potential for future oil discoveries in the Williston Basin.

Geochemistry, origin, and deep-water exploration potential of natural gases in the Pearl River Mouth and Qiongdongnan basins, South China Sea

Abstract: More than 10 gas pools have been discovered since 1983 in the shallow-water region of the Pearl River Mouth (PRM) Basin and the Qiongdongnan (QDN) Basin, offshore South China Sea. Gases produced from the QDN basin are characterized by high contents of benzene and toluene and relatively heavy 13C2 values (25 to 27). The associated condensates have a high abundance of bicadinanes and oleanane, providing a good correlation with the coal-bearing sequence of the Oligocene Yacheng Formation in the basin. In contrast, the gases from the PRM basin contain lower amounts of benzene and toluene and lighter 13C2 values (24 to 34). Widely variable concentrations of bicadinane and oleanane were identifiied from the associated condensates, which can be mostly correlated with the lower Oligocene Enping Formation source rocks formed in a swamp to shallow lake environment. Oil-cracked gases sourced from the Eocene oil-prone source rock may also provide some contribution to the PRM basin gases. The available geochemical data indicate that both the Yacheng and Enping formations contain mainly type III and II2 kerogens with dominant gas potential. Regional geological information indicates that the deep-water regions of the two basins share the same hydrocarbon source sags with the shallow-water areas, and they developed massive sandstone reservoirs during the Oligocene and Miocene. Fluid-flow modeling results show that the deep-water regions were on the pathway of lateral migrating gases, and the interpreted reservoirs in these zones have developed abundant seismic bright spots, which may reflect the presence of gas. The deep-water regions of the offshore South China Sea are believed to have great gas exploration potential.

Structural and diagenetic control of fluid migration and cementation along the Moab fault, Utah

Abstract: The Moab fault, a basin-scale normal fault that juxtaposes Jurassic eolian sandstone units against Upper Jurassic and Cretaceous shale and sandstone, is locally associated with extensive calcite and lesser quartz cement. We mapped the distribution of fault-related diagenetic alteration products relative to the fault structure to identify sealing and conductive fault segments for fluid flow and to relate fault–fluid-flow behavior to the internal architecture of the fault zone. Calcite cement occurs as vein and breccia cement along slip surfaces and as discontinuous vein cement and concretions in fault damage zones. The cement predominates along fault segments that are composed of joints, sheared joints, and breccias that overprint earlier deformation bands. Using the distribution of fault-related calcite cement as an indicator of paleofluid migration, we infer that fault-parallel fluid flow was focused along fault segments that were overprinted by joints and sheared joints. Joint density, and thus fault-parallel permeability, is highest at locations of structural complexity such as fault intersections, extensional steps, and fault-segment terminations. The association of calcite with remnant hydrocarbons suggests that calcite precipitation was mediated by the degradation and microbial oxidation of hydrocarbons. We propose that the discontinuous occurrence of microbially mediated calcite cement may impede, but not completely seal, fault-parallel fluid flow. Fault-perpendicular flow, however, is mostly impeded by the juxtaposition of the sandstone units against shale and by shale entrainment. The Moab fault thus exemplifies the complex interaction of fault architecture and diagenetic sealing processes in controlling the hydraulic properties of faults in clastic sequences.

Three-dimensional seismic geomorphology of a deep-water slope-channel system: The Sequoia field, offshore west Nile Delta, Egypt

Abstract: Within the Nile Delta gas province, reservoirs are dominated by Pliocene slope-channel systems, which are spectacularly imaged on high-quality three-dimensional seismic data. This article deals with the detailed seismic geomorphology of the Sequoia channel system, focusing on the geometry and distribution of its component sand bodies and the impact they have on reservoir heterogeneity. The Sequoia reservoir serves as a potential analog for similar but less well-imaged, deep-water slope systems. The reservoir consists of a succession of sandstones and mudstones organized into a composite upward-fining profile. Sand bodies include laterally amalgamated channels, sinuous channels, channels with frontal splays, and leveed channels and are interpreted to be the products of deep-water gravity-flow processes. Above a major basal incision surface, the reservoir is highly sand prone and made up of laterally amalgamated channels. The medial section of the reservoir is more aggradational and exhibits laterally isolated and sinuous channels. Within the upper part of the reservoir, channels are smaller, straighter, and built of individual channels with associated frontal splay elements and less common leveed channels. The main channel system is buried by a prograding slope succession that includes lobate sand-sheet elements. The stacking of facies within the Sequoia channel system implies a punctuated waning of sediment supply prior to eventual abandonment. The Sequoia channel is interpreted to be the late lowstand to transgressive infilling of a third-order early lowstand slope incision. The channel fill is overlain by a mudstone unit, which delineates a major correlatable hot gamma-ray event, and on seismic data, is a prominent downlap surface and therefore a possible maximum flooding surface. The Sequoia channel system shows evidence for synsedimentary faulting, including a large-scale downdip widening of the channel and small-scale channel diversions and intraslope ponding of flows. Understanding reservoir architecture in terms of sand-body geometries and connectivity is vital within Sequoia because the gas column occupies the most complex and heterogeneous upper part of the reservoir. Correspondingly, the basal sand-rich part of the reservoir will significantly influence aquifer behavior during production.

Present-day stress and neotectonics of Brunei: Implications for petroleum exploration and production

Abstract: The present-day state of stress in Tertiary deltas is poorly understood but vital for a range of applications such as wellbore stability and fracture stimulation. The Tertiary Baram Delta province, Brunei, exhibits a range of contemporary stress values that reflect the competing influence of the northwest Borneo active margin (situated underneath the basin) and local stresses generated within the delta. Vertical stress (v) gradients at 1500-m (4921-ft) depth range from 18.3 MPa/km (0.81 psi/ft) at the shelf edge to 24.3 MPa/km (1.07 psi/ft) in the hinterland, indicating a range in the shallow bulk density across the delta of 2.07–2.48 g/cm3. The maximum horizontal stress (Hmax) orientation rotates from margin parallel (northeast–southwest; deltaic) in the outer shelf to margin normal (northwest–southeast; basement associated) in the inner shelf. Minimum horizontal stress (hmin) gradients in normally pressured sequences range from 13.8 to 17.0 MPa/km (0.61–0.75 psi/ft) with higher gradients observed in older parts of the basin. The variation in contemporary stress across the basin reveals a delta system that is inverting and self-cannibalizing as the delta system rapidly progrades across the margin. The present-day stress in the delta system has implications for a range of exploration and production issues affecting Brunei. Underbalanced wells are more stable if deviated toward the hmin direction, whereas fracture stimulation in mature fields and tight reservoirs can be more easily conducted in wells deviated toward Hmax. Finally, faults near the shelf edge are optimally oriented for reactivation, and hence exploration targets in this region are at a high risk of fault seal breach.

Mechanisms for oil depletion and enrichment on the Shijiutuo uplift, Bohai Bay Basin, China

Abstract: The Shijiutuo uplift is a major uplift to the north of the Bozhong depression, the largest generative kitchen in the Bozhong subbasin, Bohai Bay Basin. Although the N35-2 trap on this uplift contains a medium-size oil accumulation and the Q32-6 trap contains China's third largest offshore oil accumulation, the Q31-1 trap between the N35-2 and Q32-6 traps with very similar evolution history was confirmed to be dry. Biomarker associations of crude oil and source rock samples were analyzed, and three-dimensional migration pathway modeling was conducted to investigate the origin of oils and mechanisms for oil enrichment and depletion on the uplift. Multiple-parameter oil-source correlation and hierarchical cluster analysis using 10 selected biomarker parameters allowed the identification of four source-related oil classes. Almost all oils from the Shijiutuo uplift are derived from the Eocene Shahejie Formation, whereas oils found between the Shijiutuo uplift and the Bozhong depression either are derived from or have important contributions from the Oligocene Dongying Formation. Variations in oil classes and biomarker parameters suggest sequential migration of oil generated from the Shahejie and then Dongying formations in the Bozhong depression, which is reasonably supported by petroleum migration pathway modeling. Oil charge from two oil-prone source rock intervals and, more importantly, focusing of oil originating from a large area of the Bozhong generative kitchen into the same trap accounted for oil enrichment and formation of China's third largest offshore oil field in the Q32-6 structure. The complexity and primary control of the sealing surface (top surface of the carrier bed) morphology on the positions of migration pathways caused the Q31-1 trap to be shielded from migration of oil originating from the Bozhong depression, resulting in oil depletion in this trap. Shadows to petroleum migration may occur because of the three-dimensional behavior of petroleum migration, and two-dimensional migration modeling may be misleading in predicting petroleum occurrences.

Fault facies and its application to sandstone reservoirs

Abstract: The concept of fault facies is a novel approach to fault description adapted to three-dimensional reservoir modeling purposes. Faults are considered strained volumes of rock, defining a three-dimensional fault envelope in which host-rock structures and petrophysical properties are altered by tectonic deformation. The fault envelope consists of a varying number of discrete fault facies originating from the host rock and organized spatially according to strain distribution and displacement gradients. Fault facies are related to field data on dimensions, geometry, internal structure, petrophysical properties, and spatial distribution of fault elements, facilitating pattern recognition and statistical analysis for generic modeling purposes. Fault facies can be organized hierarchically and scale independent as architectural elements, facies associations, and individual facies. Adding volumetric fault-zone grids populated with fault facies to reservoir models allows realistic fault-zone structures and properties to be included. To show the strength of the fault-facies concept, we present analyses of 26 fault cores in sandstone reservoirs of western Sinai (Egypt). These faults all consist of discrete structures, membranes, and lenses. Measured core widths show a close correlation to fault displacement; however, no link to the distribution of fault facies exists. The fault cores are bound by slip surfaces on the hanging-wall side, in some cases paired with slip surfaces on the footwall side. The slip surfaces tend to be continuous and parallel to the fault core at the scale of the exposure. Membranes are continuous to semicontinuous, long and thin layers of fault rock, such as sand gouge, shale gouge, and breccia, with a length/thickness ratio that exceeds 100:1. Most observed lenses are four sided (Riedel classification of marginal structures) and show open to dense networks of internal structures, many of which have an extensional shear (R) orientation. The average lens long axis/short axis aspect ratio is about 9:1.

Spatial variation of microstructure and petrophysical properties along deformation bands in reservoir sandstones

Abstract: A series of deformation bands from various reservoir sandstones deformed at different burial depths have been studied with respect to microstructural and petrophysical variations. In many of the examples explored, the internal microstructure, porosity, and permeability vary along the bands at the centimeter or even millimeter scale, changing and in most cases reducing the ability of the bands to act as barriers to fluid flow. Porosity varies by up to 18% and permeability by up to two orders of magnitude. Such petrophysical variations are found along different types of deformation bands, but the range depends upon the deformation mechanisms, in particular on the degree of cataclasis and dissolution in cataclastic and dissolution bands, and on the phyllosilicate content in disaggregation bands. For cataclastic bands, the grain-size distribution changes along the bands with regard to the degree of cataclasis. Furthermore, the increased specific surface area of the pore-grain interface as a result of cataclasis causes higher permeability reduction in cataclastic bands than in other types of deformation bands. Phyllosilicate content can influence the thickness of phyllosilicate bands. However, no apparent correlation between thickness and intensity of cataclasis in the studied cataclastic deformation bands is observed.

Sequence stratigraphy of experimental strata under known conditions of differential subsidence and variable base level

Abstract: Sequence stratigraphy has been applied from reservoir to continental scales, providing a scale-independent model for predicting the spatial arrangement of depositional elements. We examine experimental strata deposited in the Experimental EarthScape facility at St. Anthony Falls Laboratory, focusing on stratigraphic surfaces defined by discordant contact geometries, surfaces analogous to those delineated in the original work on seismic sequence stratigraphy. In this controlled setting, we directly evaluate critical sequence-stratigraphic issues, such as stratigraphic horizon development and time significance, as well as the internal geometry and migration of the bounded strata against the known boundary conditions and depositional history. Four key stratigraphic disconformities defined by marine downlap, marine onlap, fluvial erosion, and fluvial onlap are mapped and vary greatly in their relative degree of time transgression. Marine onlap and downlap contacts closely parallel topographic surfaces (time surfaces) and, prior to burial, approximate the instantaneous offshore topography. These stratal-bounding surfaces are also robust stratigraphic signals of relative base-level fall and rise, respectively. Marine onlap surfaces are of special interest. They tend to be the best preserved discordance, where widespread, allogenic-based onlap surfaces subdivide otherwise amalgamated depositional cycles amidst cryptic stacks of marine foresets; however, local, autogenic-based marine onlap discordances are present throughout the fill. A critical distinguishing feature of allogenic onlap is the greater lateral persistence of the discordance. Surfaces defined by subaerial erosional truncation and fluvial onlap do not have geomorphic equivalence because channel processes continually modify the surface as the stratigraphic horizons are forming. Hence, they are strongly time transgressive. Last, the stacking arrangement of the preserved bounded strata is found to be a good time-averaged representation of the mass-balance history.

Three-dimensional modeling of a shoreface-shelf parasequence reservoir analog: Part 1. Surface-based modeling to capture high-resolution facies architecture

Abstract: Conventional reservoir modeling approaches are developed to account for uncertainty associated with sparse subsurface data but are not equipped for detailed reconstruction of high-resolution geologic data sets. We present a surface-based modeling procedure that enables explicit representation of heterogeneity across a hierarchy of length scales. Numerous surfaces are used to construct complex facies-body geometries and distributions prior to generating a grid, allowing sampled and conceptual data to be fully incorporated within field-scale models. Our approach is driven by the improved efficiency that surfaces introduce to reservoir modeling through their geologically intuitive design, rapid construction, and ease of manipulation. Cornerpoint gridding of the architecture defined by the surfaces reduces the number of cells required to represent complex geometries, thus preserving geologic detail and rendering upscaling unnecessary for fluid-flow simulations. The application of surface-based modeling is demonstrated by reconstructing the detailed three-dimensional facies architecture of a wave-dominated shoreface-shelf parasequence from a rich outcrop data set. The studied outcrop data set describes reservoir architecture in a generic analog for many shallow-marine reservoirs. The process of model construction has demonstrated the function of (1) shoreface-shelf clinoforms, (2) paleogeographic changes in shoreline orientation, and (3) storm-event-bed amalgamation in controlling facies architecture. These subtle geometric features cannot be accurately represented using conventional stochastic reservoir modeling algorithms, which results in poor estimation of facies proportions and associated hydrocarbon volumes in place. In contrast, the surface-based modeling approach honors all data and captures subtle geometric facies relationships, thus allowing detailed and robust reservoir characterization.

Charging of the Neogene Penglai 19-3 field, Bohai Bay Basin, China: Oil accumulation in a young trap in an active fault zone

Abstract: The Penglai 19-3 (PL19-3) oil field, the largest offshore oil field in China, was found in shallow reservoirs (700–1700 m, 2297–5577 ft) within the most active fault zone in east China. The PL19-3 anticline was not finally formed until about 2.0 Ma and is cut by densely distributed faults. Source rock and crude oil samples from the PL19-3 field were analyzed to determine the origin and formation mechanisms of this large oil field. Three organic-rich, oil-prone source rock intervals exist in the Bozhong subbasin, each of which has a distinct biomarker assemblage. Oil samples from different wells have different biomarker associations, and three source-related oil classes were identified within the PL19-3 field based on biomarker compositions and multivariate analysis of the data. The PL19-3 field displays considerable compositional heterogeneity. The compositional heterogeneity within the field and comparison between oil samples from the PL19-3 field and those from nearby structures suggest three field-filling directions, which is consistent with the results of migration pathway modeling. The PL19-3 field was charged in the north by oil generated from Dongying Formation source rocks in the eastern Bozhong depression and Bodong depression, in the southeast by oil generated from Shahejie Formation source rocks in the Miaoxi depression, and in the northwest by oil generated from Shahejie Formation source rocks in the central Bozhong depression. Oil charge from multiple source rock intervals and multiple generative kitchens and focusing of oil originating from a large area of the Bozhong depression into the same trap resulted in rapid oil accumulation in the PL19-3 structure and the formation of this large oil field in a very young trap within an active fault zone.

Pore-throat characterization in highly porous and permeable sandstones

Abstract: Several methods have been developed to characterize the pore spaces in sandstone reservoirs using data on the pore-throat-size distribution obtained from mercury injection tests. The Winland equation, the threshold pressure, the displacement pressure, and Pittman's equation are mostly used for this purpose to delineate the stratigraphic traps and seals. This study examines the reliability of these methods applied to the highly permeable Nubia sandstones in their type section in southern Egypt. These sandstones are composed mainly of siliceous sandstones and constitute the main Paleozoic–Cretaceous aquifers and reservoirs in Egypt. Routine core analysis and mercury injection tests were conducted to delineate the pore network characteristics for these rocks. The relationships between helium porosity and the uncorrected air permeability from the routine core analysis, and the various parameters derived from mercury injection–capillary pressure curves were established using multiple regressions. This study indicates the high reliability of the displacement pressure at 10% mercury saturation and also reveals the apex of Pittman's hyperbole at 45% mercury saturation as a complexity apex at which the pore network becomes highly chaotic. Despite the great benefits of such types of measurements, they are not commonly used because of their high cost. This study introduces a series of empirical equations for constructing a partial pore-aperture-size distribution curve from routine core analysis for the highly permeable rocks.

Petroleum reservoir porosity versus depth: Influence of geological age

Abstract: Average porosity values for the producing zones of oil and gas fields worldwide are examined as a function of the present depth for sandstone and carbonate lithologies divided into 10 groupings by reservoir depositional age (Precambrian–Silurian to Pliocene–Pleistocene). The wide variations in average reservoir porosity within each depth range reflect the extreme ranges in porosity-controlling factors such as depositional facies, early

Impact of interlayer slip on fracture prediction from geomechanical models of fault-related folds

Abstract: Understanding and interpreting the timing, location, orientation, and intensity of natural fractures within a geologic structure are commonly important to both exploration and production planning activities. Here we explore the application of finite-element-based geomechanical models to fracture prediction. Our approach is based on the idea that natural fractures can be interpreted or inferred from the geomechanical-model-derived permanent strains. For this analysis, we model an extensional fault-tip monocline developed in a mechanically stratified limestone and shale sequence because field data exist that can be directly compared with model results. The approach and our conclusions, however, are independent of the specific structural geometry. The presence or absence of interlayer slip is shown to strongly control the distribution and evolution of strain, and this control has important implications for interpreting fractures from geomechanical models.

A new approach for outcrop characterization and geostatistical analysis of a low-sinuosity fluvial-dominated succession using digital outcrop models: Upper Triassic Oukaimeden Sandstone Formation, central High Atlas, Morocco

Abstract: The use of digital outcrop models in combination with traditional sedimentological field data improves the accuracy and efficiency of qualitative and quantitative characterization of outcrop analogs for reservoir modeling purposes. In this article, we apply an innovative methodology of outcrop characterization to an Upper Triassic fluvial-dominated system, exposed in extensive outcrops with limited three-dimensional (3-D) exposure. Qualitative analysis of the study outcrop allows the subdivision of the formation into three architectural intervals. Each interval can be further subdivided into subintervals on the basis of architectural style. This subdivision provides information on reservoir compartmentalization, which is used for zonation of the geocellular model. Qualitative analysis also provides valuable information on reservoir facies distribution. A new technique termed "perpendicular projection plane" is presented as a tool for quantitative analysis of outcrops with reduced 3-D exposure. This technique improves the accuracy of apparent width measurements of geobodies exposed in outcrops, which are subparallel to paleoflow. The quantitative analysis provides a detailed data set of geobody dimensions to use as conditioning data for analog reservoir models. Statistical analysis of the dimensions provides empirical relationships to apply in subsurface analog systems to reduce uncertainty related to stochastic modeling approaches.

Gas geochemistry of the Mobile Bay Jurassic Norphlet Formation: Thermal controls and implications for reservoir connectivity

Abstract: The Mobile Bay gas field is located offshore Alabama in the northern Gulf of Mexico. Production is from eolian dunes of the Jurassic Norphlet sandstone at depths exceeding 6100 m (20,000 ft) and temperatures greater than 200C. Reservoir connectivity and compositional variation, including the distribution of nonhydrocarbon gases (H2S and CO2), are critical factors in production strategy. To evaluate the controls on compositional variation and connectivity, detailed molecular and isotopic analyses were conducted for 29 wells. Analysis of volatiles in fluid inclusions suggests that the field was originally filled with oil that subsequently cracked to gas. In addition to the thermal destruction (cracking) of oil, the process of thermochemical sulfate reduction (TSR) continues to destroy the remaining hydrocarbons through oxidation of gas and reduction of sulfate to form H2S and CO2. The variable extent of the TSR process at Mobile Bay results in a wide range of hydrocarbon and H2S compositions. Condensates are almost exclusively composed of diamondoids whose composition appears controlled by H2S concentrations. In contrast to hydrocarbon and H2S contents, CO2 concentrations are relatively constant throughout the field. Carbon isotopic ratios for CO2 correlate positively with those for wet-gas hydrocarbons but are heavier than expected for CO2 originating from hydrocarbon oxidation via TSR. The narrow range of CO2 contents and heavy isotope ratios suggests that CO2 is regulated by water-rock equilibration and carbonate precipitation. The destruction of the hydrocarbon gas and mineralization of the carbon dioxide product create a volume reduction and an associated drop in reservoir pressure. This process creates several internal sinks (or exits) that may control the spill direction for gas in the field.

Quantifying and predicting naturally fractured reservoir behavior with continuous fracture models

Abstract: This article describes the workflow used in continuous fracture modeling (CFM) and its successful application to several projects. Our CFM workflow consists of four basic steps: (1) interpreting key seismic horizons and generating prestack and poststack seismic attributes; (2) using these attributes along with log and core data to build seismically constrained geocellular models of lithology, porosity, water saturation, etc.; (3) combining the derived geocellular models with prestack and poststack seismic attributes and additional geomechanical models to derive high-resolution three-dimensional (3-D) fracture models; and (4) validating the 3-D fracture models in a dynamic reservoir simulator by testing their ability to match well performance. Our CFM workflow uses a neural network approach to integrate all of the available static and dynamic data. This results in a model that is better able to identify fractured areas and quantify their impact on well and reservoir flow behavior. This technique has been successfully applied in numerous sandstone and carbonate reservoirs to both understand reservoir behavior and determine where to drill additional wells. Three field case studies are used to illustrate the capabilities of the CFM approach.

Overlapping faults and their effect on fluid flow in different reservoir types: A LIDAR-based outcrop modeling and flow simulation study

Abstract: In this article we focus on the potential of fault-overlap zones as conduits for fluid flow in a variety of reservoir types. Light detection and ranging (LIDAR) technology were applied to collect a three-dimensional, spatially constrained data set from a well-exposed fault-overlap zone that crops out in the Devil's Lane area of the Canyonlands National Park in Utah. A virtual outcrop was generated and used to extract structural and stratigraphic data that were taken into a reservoir modeling software and reconstructed. The outcrop-based model was flow simulated and used to test fluid flow through a real-world fault-overlap zone. A structural framework was built based on collected outcrop data and combined with a series of nine different facies models. The different facies models included an eolian model based on the outcrop and a range of synthetic fluvial and shallow marine systems. Results show that, for certain depositional models, cross-fault reservoir communication may be poor despite the geometric connectivity of the relay beds. This was the case for low net/gross fluvial models and shoreface models. Conversely, high net/gross fluvial systems and eolian systems show good communication through the same relay zone. Overall, the results show that, in the presence of a fault-overlap zone, pressure communication across a relay ramp may still be poor depending on the scale of the faults and relay ramp as well as the geometry and volume of the sands.

Physical properties of Cenozoic mudstones from the northern North Sea: Impact of clay mineralogy on compaction trends

Abstract: Vertical and lateral changes in physical properties in Cenozoic mudstones from the northern North Sea Basin reflect differences in the primary mineralogical composition and burial history, which provides information about sedimentary facies and provenance. Integration of well-log data with mineralogical information shows the effect of varying clay mineralogy on compaction curves in mudstones. The main controlling factor for the compaction of Eocene to early Miocene mudstones within the North Sea is the smectite content, which is derived from volcanic sources located northwest of the North Sea. Mudstones with high smectite content are characterized by low P-wave velocities and bulk densities compared to mudstones with other clay mineral assemblages at the same burial depths. Smectitic clays are important during mechanical compaction because they are less compressible than other types of clay minerals. A comparison between well-log data and experimental work also shows that smectite may be a controlling factor for overpressure generation in the smectite-rich Eocene and Oligocene sediments. At greater burial depths and temperatures (70–80C), the dissolution of smectite and precipitation of illite and quartz significantly increases velocities and densities. Miocene and younger mudstones from the northern North Sea have generally low smectite contents and as a result have higher velocities and densities than Eocene and Oligocene mudstones. Lateral differences in the compaction trends between the north and south for these sediments also exist, which may be related to two different source areas in the Pliocene. The log-derived petrophysical data from the northern North Sea Basin show that mudstone lithologies have very different compaction trends depending on the primary composition. Simplified compaction curves may therefore affect the outcomes from basin modeling. The amplitude-versus-offset response of hydrocarbon sands and the seismic signature on seismic sections are also dependent on the petrophysical properties of mudstones and will vary depending on the mineralogical composition.

Three-dimensional modeling of a shoreface-shelf parasequence reservoir analog: Part 2. Geologic controls on fluid flow and hydrocarbon recovery

Abstract: Wave-dominated, shoreface-shelf parasequences are commonly modeled as simple layer-cake reservoirs, yet analysis of modern and ancient analogs demonstrates that these intervals contain a more complex physical stratigraphy. We investigate the impact of depositional and diagenetic heterogeneity associated with gently dipping clinoform surfaces on fluid flow and recovery during water flooding, using a three-dimensional model reconstructed from a well-exposed outcrop analog. We demonstrate that the volume of oil in place is affected by variations in facies thickness associated with interfingering along clinoforms, whereas waterflood sweep efficiency is affected by barriers to flow along clinoform surfaces, such as calcite-cemented layers, mudstones, and siltstones. Sweep efficiency is low when water flooding is down depositional dip because oil is bypassed at the toe of each clinothem as water flows preferentially through high-quality sandstone facies in the upper part of the parasequence. Sweep efficiency is higher when water flooding is up depositional dip because the gravity-driven, downward flow of water sweeps poorer-quality sandstone facies in the lower part of the parasequence. In both cases, injectors may offer limited pressure support to producers. Water flooding along depositional strike yields good pressure support but poor sweep because the gravity-driven flow of water into the lower part of the parasequence is significantly reduced. This yields highly variable fluid saturations but a uniform pressure gradient, which is consistent with pressure and fluid saturation data from the mature Rannoch Formation reservoir, Brent field, United Kingdom North Sea. Simple layer-cake models fail to capture the range of flow behaviors described above and overpredict recovery by up to 20% as a result.

Fracture characterization at multiple scales using borehole images, sonic logs, and walkaround vertical seismic profile

Abstract: We present a quantitative forward-modeling methodology to link and interpret several measurements relevant to mechanical properties of fractures such as borehole images, sonic anisotropy logs, and borehole seismic anisotropy. The analysis is applied to a case study from a north African tight gas field using data from a vertical well. Two studies are conducted independently using the same geological fracture data to model fracture-induced anisotropy. In the first study, we use the orientation of the natural and drilling-induced fractures interpreted on the image log to model the azimuthal fracture-induced anisotropy at the sonic scale. The mechanical effects of natural and drilling-induced fractures are treated using different compliance parameters for each fracture type. We show that modeled sonic fast shear azimuths could be biased by the presence of noncompliant fractures in each fracture type, and we propose an empirical selection criterion to reject noncompliant fractures prior to compliance estimation. Then, we estimate the fracture compliances and confirm that natural open fractures have larger compliances than drilling-induced fractures. In the second study, we apply interpreted borehole images toward modeling of the azimuthal vertical seismic profile (VSP) attributes as a function of source azimuthal position. Natural fractures inside a window of height, h, and located at depth, d, are included, and several volume sizes and positions (i.e., h and d) are considered. We find a good agreement between modeled and observed transverse-over-radial displacement trends using natural fractures within windows located at the depth of the VSP receiver, and having window heights on the order of one to two VSP shear wavelengths.