Engineering topics
Laubach, Stephen E.
Publications and source records attributed to Laubach, Stephen E..
Multiscale spatial analysis of fracture nodes in two dimensions
Spatial arrangement of fractures as a function of scale is an important component of fracture quantification for inferential and predictive modeling. Available methods that analyze fracture spatial arrangement are based on one-dimensional spacing data; therefore, they are limited to semi-parallel fractures. Such methods cannot be applied to fracture networks in higher dimensions, particularly when fractures have different orientations. Here to characterize fracture arrangements in two dimensions, we propose using Ripley’s K-function, as a method of point pattern analysis, to quantify spatial arrangement of fracture nodes. Fracture nodes, such as barycenters, intersection points, and tips, are point-based representations of fracture locations and connectivity within the fracture network. We introduce formulations for isotropic as well as directional analyses of spatial arrangement. In addition, we derive formulations for edge correction in circular and rectangular study domains. Finally, we demonstrate applications of Ripley’s K-function on two natural fracture datasets. Our proposed method supports quantification and characterization of fracture spatial arrangements that allow practitioners to build representative models of fractures in the subsurface.
Opening-mode fracturing and cementation during hydrocarbon generation in shale: An example from the Barnett Shale, Delaware Basin, West Texas
Relative timing of fracturing is a key input for predictive fracture models, but timing information for fractures is commonly difficult to obtain. In this study, we used crosscutting relations and fluid inclusion assemblage temperatures from fracture cements from a few well-documented sampled fractures, combined with a one-dimensional burial history model, to establish timing for three generations of opening-mode fractures in a Barnett Shale core from the southern part of the Delaware Basin, Pecos County, West Texas. A burial history model is presented for the cored well and matched to measured vitrinite reflectance in samples from the core, and bottomhole temperature in the well. Here, the earliest fractures (group 1) likely formed due to early fluid-expulsion events (ca. 300 Ma) and were folded during host-rock compaction. Later group 2 fractures are sealed with fibrous barite containing primary, liquid hydrocarbon inclusions (mean homogenization temperature [T h ] = –9°C) and aqueous fluid inclusions (mean T h = 108.1°C). Group 2 fractures likely formed in response to fluid overpressure associated with cracking of type II kerogen to oil. Group 3 vertical fractures are up to 2 m in height with kinematic apertures ranging from less than 0.05 to 1.4 mm, partly open, and strike dominantly 010°–020°. Sequentially trapped aqueous fluid inclusions in fracture-spanning quartz cement bridges (mean T h = 110°C in crack-seal texture and 128°C in post-crack-seal fracture cement) record fracture opening under increasing temperature, inferred to reflect increasing burial, with continued overpressuring during the Triassic to Late Cretaceous. Some group 3 fractures may have continued to fill during Cenozoic uplift.
Mechanisms for the Generation of Complex Fracture Networks: Observations From Slant Core, Analog Models, and Outcrop
We use observations of hydraulic fractures in core, outcrop attributes of natural hydraulic fractures, and analogue models, to address how hydraulic fracture networks evolve. A slant core from the Wolfcamp Formation—an unconventional shale hydrocarbon reservoir in the Permian Basin of West Texas—collected within 18 m and 30 m of two hydraulically stimulated horizontal wells, provided an opportunity to examine hydraulic fractures directly. In approximately 183 m of core, 309 calcite-sealed natural opening-mode fractures and 375 hydraulic fractures were identified. Many hydraulic fractures in the core show complex morphology, including twist-hackle segmentation, diversion, and bifurcation; these structures most commonly develop at lithological bed boundaries and mechanical heterogeneities such as natural fractures and concretions. An outcrop of bed-parallel pavements in the Cretaceous Boquillas Formation in West Texas contains opening-mode fractures that likely formed by natural hydraulic fracturing. Fracture traces provide evidence of twist-hackle segmentation, and are typically associated with bed boundaries and preexisting bed-parallel stylolites. A laboratory study of hydraulic fracturing of 33 synthetic blocks of gypsum and hydrostone revealed fracture steps, diversions, twist hackles, and multiple overlapping fractures together with information on fracture growth directions. These complexities in the fracture network were dominantly nucleated at inclusions used to simulate pre-existing fractures, and as a result of mechanical heterogeneity introduced by the wellbore and perforations. Collectively, our results show that complex fracture networks are produced in hydraulic fracturing of self-sourced reservoir strata. Mechanical stratigraphic boundaries and other heterogeneities are likely to enhance fracture network complexity through the processes of segmentation, diversion, and bifurcation. These processes create multiple fracture strands, resulting in an increased number of hydraulic fractures over those initiated, thereby increasing total fracture surface area. Our study provides insight into hydraulic fracture network propagation, and has applications for evaluation, completion, production, and fracture modeling of unconventional reservoirs.
Value-of-Information Analysis of a Fracture Prediction Method
A core-based fracture prediction method is used to illustrate a Value of Information (VOI) decision-analysis protocol to inform completion decisions in tight-gas sandstones. The ratio of late host-rock cement to available pore volume or degradation index uses petrographic observations of cement distributions in core (including sidewall cores) to predict whether nearby but unsampled fractures (widths > 0.5 to 1 mm) are sealed (non-conductive) or open (conductive). Measurements from four sandstone plays suggest that the index correctly predicts open vs. sealed fractures with an accuracy in excess of 80%. The value added is calculated using Bayesian inference in which accuracy of the index serves as the likelihood of the prior distribution of open fractures to assess the posterior probability that data represent a useful predictor of producibility. VOI of the prediction method is more than three times the cost to acquire the data. VOI is most sensitive to play-specific geologic and cost parameters including cost to drill, expected revenue from a successful well, cost of completion, cost of acquiring data for the index, and fracture probability distributions. The approach provides a way to value acquiring fracture data and points to a need for zone-specific production data in tight-gas sandstones.
The nature and origins of decametre-scale porosity in Ordovician carbonate rocks, Halahatang oilfield, Tarim Basin, China
At >7 km depths in the Tarim Basin, hydrocarbon reservoirs in Ordovician rocks of the Yijianfang Formation contain large cavities ( c . 10 m or more), vugs, fractures and porous fault rocks. Although some Yijianfang Formation outcrops contain shallow (formed near surface) palaeokarst features, cores from the Halahatang oilfield lack penetrative palaeokarst evidence. Outcrop palaeokarst cavities and opening-mode fractures are mostly mineral filled but some show evidence of secondary dissolution and fault rocks are locally highly ( c . 30%) porous. Cores contain textural evidence of repeated formation of dissolution cavities and subsequent filling by cement. Calcite isotopic analyses indicate depths between c . 220 and 2000 m. Correlation of core and image logs shows abundant cement-filled vugs associated with decametre-scale fractured zones with open cavities that host hydrocarbons. A Sm–Nd isochron age of 400 ± 37 Ma for fracture-filling fluorite indicates that cavities in core formed and were partially cemented prior to the Carboniferous, predating Permian oil emplacement. Repeated creation and filling of vugs, timing constraints and the association of vugs with large cavities suggest dissolution related to fractures and faults. In the current high-strain-rate regime, corroborated by velocity gradient tensor analysis of global positioning system (GPS) data, rapid horizontal extension could promote connection of porous and/or solution-enlarged fault rock, fractures and cavities. Supplementary material: Stable isotopic analyses and the velocity gradient tensor and principal direction and magnitude calculation are available at https://doi.org/10.6084/m9.figshare.c.4946046 Thematic collection: This article is part of the The Geology of Fractured Reservoirs collection available at: https://www.lyellcollection.org/cc/the-geology-of-fractured-reservoirs
Fracture, Dissolution, and Cementation Events in Ordovician Carbonate Reservoirs, Tarim Basin, NW China
Ordovician carbonate rocks of the Yijianfang Formation in the Tabei Uplift, Tarim Basin, contain deeply buried (>6000 m), highly productive oil and gas reservoirs associated with large cavities (>10 m). Previous workers inferred that large cavities are paleocaves (paleokarst) formed near the surface and subsequently buried. Alternately, caves may have formed by dissolution at depth along faults. Using 227 samples from 16 cores, we document textures and cement compositions bearing on cavity histories with petrographic, high-resolution scanning electron microscopy (SEM), isotopic, and fluid inclusion microthermometric observations. Results show that dissolution occurred at depth and was caused by (1) acidic fluids derived from Middle-Late Silurian and/or Devonian-Permian hydrocarbon generation and maturation, (2) high-temperature fluids, of which some were associated with Late Permian igneous activity, and (3) Mg-rich fluids that accompanied Jurassic-Cretaceous deformation and the formation of partially open fractures and stylobreccias (fault breccias). The relative paragenetic sequence of the structure-related diagenesis suggests seven stages of fracturing, dissolution, and cementation. Mottle fabrics in the Yijianfang Formation contain argillaceous carbonate-rich silt and are bioturbation features formed within the marine environment. Those mottled fabrics differ from clearly karstic features in the overlying Lianglitage Formation, which formed by near-surface dissolution and subsequent infilling of cavities by allochthonous sediment. Mottle fabrics are crosscut by compacted fractures filled with phreatic-vadose marine cements and followed by subsequent generations of cement-filled fractures and vugs indicating that some fractures and vugs became cement filled prior to later dissolution events. Calcite cements in fractures and vugs show progressively depleted values of δ18O documenting cement precipitation within the shallow (~220 m), intermediate (~625 m), and deep (~2000 m) diagenetic environments. Deep (mesogenetic) dissolution associated with fractures is therefore the principal source of the high porosity-permeability in the reservoir, consistent with other pieces of evidence for cavities localized near faults.