SEARCH · Engineering Papers
Results for “outcrop”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Diagenesis is key to unlocking outcrop fracture data suitable for quantitative extrapolation to geothermal targets
Exceptionally large, well-exposed sandstone outcrops in New York provide insights into folds, deformation bands, and fractures that could influence permeability, heat exchange, and stimulation outcomes of geothermal reservoir targets. Cambrian Potsdam Sandstone with <5% porosity contains decimeter-scale open, angular-limbed monoclines <0.5 km apart with associated low-porosity mm-wide cataclastic deformation bands. Crossing and abutting relationships among sub-vertical opening-mode fractures show four chronological Sets A–D, striking NNW, NE, NW, and ENE, respectively. Fracture lengths and heights range from millimeters to tens of meters. Sets A and C macro-fractures, and possibly B and D, contain quartz deposits. All sets have abundant associated quartz cemented microfractures that also record set orientations and crosscutting relations. Quartz cement deposits—evidence of diagenesis—are the key to identifying attributes of outcrop fractures suitable for extrapolation to geothermal targets in sandstones because they show which fractures formed in the subsurface. Set A fluid inclusion homogenization temperatures (120°C–129°C) are compatible with fracture at >3 km depth. Fractures are stiff and those ≥0.05 mm (Set C) and ≥0.1 mm (Set A) are open and potentially conducive to flow. Sets A and D are abundant in outcrops with close fracture spacing—0.18 m and 0.68 m, respectively—and define a rectangular connectivity network dominated by crossing and abutting X and Y nodes. Set A aperture distributions follow a power law with slope –0.8 up to 0.15 mm; other sets have lognormal distributions. Set A and D microfractures are weakly clustered, while macro-fractures commonly have 1D anticlustered (regular or periodic) arrangements at shorter length scales (<0.2 m). Sub-horizontal fractures are barren and may have formed near the surface. Fracture heights, lengths, and spatial arrangements show good trace connectivity but low open connectivity. For geothermal applications, outcrop results predict low initial well-test permeabilities owing to quartz disconnecting open fractures, but stimulation of closely spaced microfractures and partly open macro-fractures could yield high surface area for heat exchange. Quantitative extrapolation of key fracture attributes like abundance, orientation, spatial arrangement, length, and open fracture connectivity is possible from outcrops to fractured reservoirs if differing thermal histories and diagenesis are accounted for.
Establishing Outcrop Time Series at Various Depths in a Nonlinear Soil Column: An Iterative Approach
Seismic soil-structure interaction (SSI) analysis of nuclear facilities is an important consideration during design and retrofit. SSI tools used in the nuclear industry are currently based on an equivalent linear approach. Procedures for developing input ground motion for equivalent linear approaches are well established. However, the procedures for establishing input ground motion for Nonlinear Soil-Structure Interaction (NLSSI) analysis of nuclear facilities is not well established. A collaborative research group at Idaho National Laboratory (INL) has recently developed analytical methods and numerical tools for using NLSSI analysis for nuclear facility seismic calculations. NLSSI analysis for a nuclear facility allows for calculation of seismic wave motion through a near field soil domain using either, (1) vertically propagating shear and compressive waves, which is the current industry practice, or (2) a three-dimensional non-vertical wave field. This paper presents an iterative procedure for establishing outcrop motion at a depth in the soil column for NLSSI analysis that uses vertically propagating shear waves. The approach presented in this paper starts with a known ground motion at the surface, deconvolves to a depth, and then convolves the obtained motion up to a different desired location of input for the NLSSI model. To demonstrate the validity of the approach a finite element soil column, that is representative of a nuclear facility site in the US, is used to produce compatible outcrop seismic time series for reduced nonlinear soil mesh depths. The developed approach for reducing the nonlinear soil column model depth is a two-step iterative method; 1) the first step is establishing an outcrop time series at the lowest depth considered that produces the top-of-soil response spectrum of an actual recorded ground motion, and 2) the second step is providing compatible outcrop time series at a shallower depth based on the information from the first step. The comparison of the 5% damped response spectrum from the resulting acceleration time series based on the iterated outcrop motions and the original acceleration time series are conducted. The study showed that the proposed iterative approach produced comparable results within 1% range of the original recorded time series results when sufficient iterations were performed.
Enhanced Outcrop Methane Capture (Final Project Report)
This report provides the final status of the “Enhanced Outcrop Methane Capture (EOMC) Project” which was focused on drilling two new methane capture wells along the Fruitland Formation Outcrop for capture of methane and carbon dioxide which are naturally seeping from the Outcrop. The EOMC Project was focused on evaluating whether novel horizontal drilling technology will result in improved capture of methane and carbon dioxide which are naturally seeping from the Fruitland Formation Outcrop on Southern Ute tribal trust lands within the Southern Ute Indian Reservation (the “Reservation”).
Fallon FORGE: Analogue Outcrop Samples
Compilation of results for mechanical and fluid flow properties of analogue outcrop samples - experimental data for compressional and shear wave velocities, tensile strengths, and compressive strengths. Outcrop location and sample orientation data are documented in a separate csv file.
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.
Outcrops as guides to subsurface natural fractures: Example from the Nikanassin Formation tight-gas sandstone, Grande Cache, Alberta foothills, Canada
Not Available
Alabama Coalfields Project Maps, Outcrop Traces and Structural Features SoApp Basin FE0032045
Arc GIS files of coal maps, traces and structural features for the Alabama Coalfields Project, Southern Appalachian Basin.
Geomechanical properties of the Meletta sandstone - the high-temperature heat storage reservoir rock of DeepStor
The DeepStor project aims at storing excess heat at temperatures up to 140 °C in the depleted Leopoldshafen oil field at a depth of about 1300 m. In order to gain knowledge on the target horizons, the different layers of the Meletta sandstone, samples cored in a block retrieved in a quarry near Nussloch where the Meletta sandstone outcrops were studied. Several petrophysical properties were investigated including mineralogy, porosity, permeability, thermal conductivity, P- and S-wave velocities and quality factor. A mechanical study focused on compressive strength under uniaxial and triaxial stress conditions, tensile strength, critical pressure and stress-dependence of physical properties. Our results show that the Meletta sandstone is heterogeneous, anisotropic, mechanically weak, stress-sensitive and prone to water weakening. Petrophysical measurements on few Meletta sandstone cores retrieved in boreholes at about 1250 m depth showed that the outcrop samples are significantly less cohesive and weaker. Therefore, one should consider the Nussloch outcrop as a good reservoir analog only with great caution. Several models were applied to account for our experimental results on anisotropy, mechanical behavior and stress-dependence of the outcrop samples. Our study provides a valuable data set which can feed numerical models to simulate the behavior of the host formation during heat exchange operations in the DeepStor project.
Geothermal Exploration North of Mount St. Helens: Washington State Play-Fairway Project
Active seismicity within a broad zone in proximity to volcanism in the St. Helens Seismic Zone north of Washington State’s Mount Saint Helens provide key ingredients for hydrothermal circulation at depth. This seismicity extends well north of the volcanic edifice below where several faults and associated fractures in outcrop record repeated slip, dilation and alteration indicative of localized fluid flow. Candidate reservoir rocks include marine metasediments overlain by extrusive volcanics. This study characterizes the faults and fractures in the extrusives which are available in outcrop and sampled along well TGH17-24 and Plate Boundary Observatory strainmeter borehole, B201. Outcrop mapping and borehole image log analysis down to 250 m document abundant volcaniclastic deposits and basalt flows that are highly fractured, although intervening ash layers truncate the vertical extent of many of these structures. However, strike slip faults with well-developed gouge zones and adjacent localized high fracture density provide potential for vertically extensive fracture flow paths. Where these rocks are altered by intrusion, they develop high fracture density that can support heat transfer. Borehole image logs further record breakout that indicate the azimuth of SHmax in the shear zone is approximately 054/234°±23.6°, consistent with seismically active strike slip faults in the seismic zone, faults mapped in outcrop and evident in the borehole image logs, and with of plate convergence data. Furthermore, temperature profiles in two wells indicate isothermal conditions at average groundwater temperatures, consistent with rapidly flowing water localized within fractures.
Aperture size distribution, length, and preferential location of bed-parallel veins in shale
Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.
Rates of subcritical cracking and long-term rock erosion
Bedrock fracture is a key element of rock erosion and subsequent surface processes. We test the hypothesis that rock's susceptibility to subcritical cracking, a specific type of fracturing, significantly drives and limits rock erosion. We measured 10 Be-derived erosion rates, compressive strength, and crack characteristics on 20 outcrops of different rock units (quartzite, granite, and two metasandstones) in the northern Blue Ridge Mountains of Virginia (USA). We also measured the subcritical cracking index (n), Charles’s law velocity constant (A), and fracture toughness (K IC ) of samples from four of the same outcrops, representative of each rock type. Erosion rates range from 1.16 ± 0.67 to 32.3 ± 7.8 m/m.y. We find strong correlations—across the four rock units—between average erosion rates and the three subcritical cracking parameters (R 2 > 0.85, p < 0.05), but not compressive strength (R 2 = 0.6; p > 0.1). We also find a correlative relationship between n and outcrop fracture length (R 2 = 0.91; p < 0.05). The latter correlation is consistent with that of published model predictions, further indicating a mechanistic link between subcritical cracking and rock erosion. We infer that subcritical cracking parameters closely tie to erosion rates, because subcritical cracking is the dominant process of mechanical weathering, leading to positive feedbacks relating subcritical cracking rates, crack length, porosity, and water accessibility. Finally, these data are the first that directly test and support the hypothesis that subcritical cracking can set the pace of long-term rock erosion.
Fracture length data for geothermal applications
Fracture lengths govern permeability and are unknowns in geothermal assessment. Along their lengths, fracture widths vary due to growth by linkage. Under the influence of diagenesis, narrow widths seal, breaking porosity continuity and reducing open length. The largest range of widths and thus susceptibility to fill occurs where fractures are linked by narrow segments. Outcrops of a geothermal target, Cambrian Potsdam quartz arenite, contain opening-mode fractures having lengths spanning five orders of magnitude from 0.082 mm to 17.9 m. Combined lengths measured at a range of scales can be described by power laws, but at a given image resolution, lengths are best fit by exponential functions. Owing to preferential sealing of small fractures, open fractures follow exponential functions, but values depend on rules for designating fractures as continuous. En échelon segments, offset 10 mm, are connected by narrow fractures or microfractures (hard linked) not evident on outcrop 1 m-elevation LiDAR or 30 m-height drone images. A rule that identifies where narrow and likely connected segments are located can yield lengths meaningful for flow simulation. Depending on diagenesis, continuity rules can halve or double average and maximum length values. Length values from outcrop for geothermal applications should be adjusted based on wellsite-specific diagenesis information.
Surface Water Mass Transformation in the Southern Ocean: The Role of Eddies Revisited
The water mass transformation (WMT) framework describes how water of one class, such as a discrete interval of density, is converted into another class via air–sea fluxes or interior mixing processes. Here, this paper investigates how this process is modified at the surface when mesoscale ocean eddies are present, using a state-of-the-art high-resolution climate model with reasonable fidelity in the Southern Ocean. The method employed is to coarse-grain the high-resolution model fields to remove eddy signatures, and compare the results with those from the full model fields. This method shows that eddies reduced the WMT by 2–4 Sv (10%–20%; 1 Sv ≡ 10 6 m 3 s −1 ) over a wide range of densities, from typical values of 20 Sv in the smoothed case. The corresponding water mass formation was reduced by 40% at one particular density increment, namely, between 1026.4 and 1026.5 kg m −3 , which corresponds to the lighter end of the range of Indian Ocean Mode Water in the model. The effect of eddies on surface WMT is decomposed into three terms: direct modulation of the density outcrops, then indirectly, by modifying the air–sea density flux, and the combined effect of the two, akin to a covariance. It is found that the first and third terms dominate, i.e., smoothing the outcrops alone has a significant effect, as does the combination of smoothing both outcrops and density flux distributions, but smoothing density flux fields alone has little effect. Results from the coarse-graining method are compared to an alternative approach of temporally averaging the data. Implications for climate model resolution are also discussed.
Applicability of 1D site response analysis to shallow sedimentary basins: A critical evaluation through physics-based 3D ground motion simulations
One-dimensional site response analysis (1D SRA) remains the standard practice in considering the effect of local soil deposits and predicting site-specific ground motions, although its range of applicability to realistic seismic wavefields is still in question. In this 1D approach, horizontal and vertical ground shaking are assumed to be induced by vertically propagating shear and compressional waves, respectively. A recent study based on analytical two-dimensional (2D) plane waves and simple point source earthquake simulations has shown two mechanistic limitations in this 1D modelling technique for general inclined seismic waves, that is, systematic over-prediction of the vertical motion and wave trapping in the 1D soil column. In this article, we evaluate in detail the applicability of this 1D modelling approach to realistic three-dimensional (3D) simulated seismic wavefields in shallow sedimentary basins. Linear-viscoelastic 1D SRA predictions using two types of input motions that are commonly used in practice—rock outcrop and in-column motions, are compared with the reference true site response results from 3D earthquake simulations in terms of various measures in the frequency and time domain. It is shown that the horizontal motion in the 3D seismic wavefield exhibits dominant shear wave propagation phenomenon, while the vertical motion is a combined effect of compressional and shear waves and can be over-predicted by the 1D approach when the incident seismic waves are inclined. Direct evidence of the wave refraction process that leads to the vertical motion over-prediction is provided. 1D SRA with in-column inputs can yield motions that have significantly longer duration compared to the true 3D site response solution due to trapped waves, casting in doubt the frequent need for increased soil damping in existing site studies to compensate for wave attenuation due to scattering alone. Sensitivity investigation on the increase of soil profile damping by a multiplier D mul shows D mul values compatible with those found in the literature for both horizontal and vertical motions. It is shown that the level of D mul optimized for a best match of the spectral acceleration is dependent on the characteristic of the input motion and a larger D mul is typically required for the vertical component. In contrast, 1D SRA with outcrop motions predicts motions with shorter significant duration due to its inability to capture the basin-edge generated surface waves. A suite of ground motion simulations was performed to assess the sensitivity of the observations to the basin geologic structure including the velocity gradient, rock-basin impedance contrast and basin depth. The analysis results show that the accuracy of the simplified 1D procedure is dependent on the wavefield composition of both the input motions and the true 3D site response solution. While the horizontal motions in shallow sedimentary basins can, to the first order, be reasonably captured by the simplified 1D approach, 1D SRA for the vertical component is in general not reliable and contributions from inclined shear waves should be accounted for in site-specific evaluation of the vertical design ground motion.
The Role of Stratigraphy and Loading History in Generating Complex Compaction Bands in Idealized Field-Scale Settings
The Buckskin Gulch locality in Utah is a landmark example of compaction localization. The outcrop of this locality involves distinct stratigraphic heterogeneity and was exposed to complex loading history. It features multiple sets of deformation bands with different kinematics and orientation. Similar formations were seen in the Valley of Fire, Nevada, and the Orange quarry, France, among other localities. The formation of such complex structures, their propagation mechanisms, and frequency is affected by numerous local and ambient factors whose impacts are not yet fully understood. The simulation of the above-mentioned localities is not feasible because of the limited amount of available information. This work, instead, investigates from a geomechanics standpoint how the interplay among material nonlinearity, outcrop stratigraphy, and loading history interconnects with specific spatiotemporal patterns of compaction band propagation. Our study shows that the system stratigraphy can be responsible for the emergence of coexisting compaction bands with different inclination and kinematics. Specifically, we show that stiffness contrasts induce nonlocal stress changes which may favor the initiation of secondary structures with different compaction localization characteristics. Furthermore, systems of inclined compaction bands induced by burial increase display secondary, noncontemporaneous sets of vertical compaction bands under the effects of postburial shortening. Our results indicate that stages of intermediate burial decrease prior to tectonic shortening can promote the formation of such complex systems. Despite the simplifications involved in our analyses, these findings show how geomechanics computations complement field observations and could provide a mechanics-based validation of site-specific reconstruction hypothesis.
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