SEARCH · Engineering Papers
Results for “geoscience”
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.
Vertical effective stress and temperature as controls of quartz cementation in sandstones: Evidence from North Sea Fulmar and Gulf of Mexico Wilcox sandstones
Here we present quantitative petrographic data, high spatial resolution oxygen isotope analyses of quartz cement, basin modelling and a kinetic model for quartz precipitation for two Paleocene-Eocene Wilcox Group sandstones from Texas and two Jurassic Fulmar Formation sandstones from the Central North Sea. At each location, one sandstone has been buried to ca. 145 °C and one to ca. 185 °C. A key difference between the Wilcox and Fulmar burial histories is that the Wilcox sandstones are currently at higher vertical effective stresses and, from basin modelling studies, have been subjected to generally higher vertical effective stresses through their burial history. The amounts of quartz cement in the Wilcox sandstones are between 12 and 18%, and between 2 and 6% in the Fulmar sandstones. High-spatial-resolution oxygen isotope data obtained from the quartz cements suggest temperature ranges for quartz precipitation from 60 to 80 °C to values approaching maximum burial temperature. Factors such as grain coatings or the timing of petroleum emplacement cannot explain the differences in the amounts of quartz cement. Petrographic data show that most of the silica for quartz cement can be derived from intergranular pressure dissolution. Although the sample set is small, we interpret the results to suggest that the differences in quartz cementation in Fulmar and Wilcox sandstones can be explained better by differences in their vertical effective stress history than their temperature history; in this case, the supply of silica rather than the precipitation of quartz becomes an important control on the rate and extent of cementation.
A pseudo-kinetic model to simulate phase changes in gas hydrate bearing sediments
Herein, modeling of the phase transitions anticipated in gas hydrate bearing sediments (GHBS) is critical for a proper understanding of time-dependent changes in states and volumes (e.g. the production of methane from this type of soils). We propose a new pseudo-kinetic approach to simulate the typical phase changes anticipated in GHBS, using published experimental results involving gas hydrate dissociation that are the basis of a widely used kinetic model. The proposed pseudo-kinetic model is formulated in the pressure-temperature (P-T) plane and assumes a rate of gas hydrate dissociation (or formation) proportional to the distance between the current state and the phase boundary. The model consists of only one parameter and is simple to implement in numerical simulators. A similar concept is used to model ice formation/thawing phenomena, but based on the ice/liquid-water phase boundary. We implemented the pseudo-kinetic model in a fully coupled thermo-hydro-chemo-mechanical (THCM) finite element code and validated it against experimental results performed on the dissociation of synthetic gas hydrate. We also evaluated the pseudo-kinetic model using synthetic cases covering several scenarios associated with gas hydrate formation/dissociation and ice formation/thawing. The model successfully reproduced the gas production test from a natural GHBS core from Korea (scaled gas venting experiment over 14 h), and also the formation of gas hydrate and ice in permafrost in Alaska (over 2 × 10 6 years). -The analyses show the versatility of the proposed pseudo-kinetic approach by applying it to model the different types of phase transitions typically encounter in GHBS. The simple formulation, easy implementation in numerical simulator, and reduced number of parameters (only one per phase change) make this model an attractive option for simulating phase transformations in problems involving GHBS.
Bulk and position-specific isotope geochemistry of natural gases from the Late Cretaceous Eagle Ford Shale, south Texas
Here, a dataset of bulk and position-specific isotope compositions of shale gases from the Late Cretaceous Eagle Ford Shale, south Texas is reported. The chemical and bulk isotopic compositions of the seven samples in this study and those available from the literature show that across the play, produced natural gas from the Eagle Ford Shale preserves a wide range of gas compositions (wetness, ~5 to >90%) and δ 13 C values of C 1 –C 3 . The depth profiles of the isotopic compositions suggest that the western region of the Eagle Ford Shale deposit had experienced as much as 700–800 m of uplift-erosion or gas migration. The timing of and lost-gas fractions by gas expulsion events very likely affected their gas compositions and bulk isotope compositions. Some deep (>3000–3500 m), matured (%Ro > 1.5) gases with heavy δ 13 C(C 1 –C 3 ) values indicate significant loss (>50%) of the early-stage gases. For the position-specific isotope deviations of propane, the ΔC 2-1 values of the five samples show small decrease with well depths, while the ΔH 2-1 values have a general increasing trend. Although some samples fall very close to the equilibrium model trajectories with reasonable calculated temperatures (138–148 °C), the position-specific isotope compositions of propane from the Eagle Ford Shale gases are likely results of thermal cracking of various organic molecules within the source rocks, which have different activation energy for cracking and non-statistical distributions of C/H isotope within them.
Episodic hydrocarbon charge in tight Mississippian reservoirs of Central Oklahoma, USA: Insights from oil inclusion geochemistry
The unforeseen discovery of substantial hydrocarbons in ultra-tight Mississippian reservoirs in Oklahoma testifies to the complexity of the hydrocarbon charge across the flanks of the Nemaha Uplift. Based on crude oil and fluid-inclusion chemistry, the hydrocarbons at the western side of the Nemaha Uplift consist of highly-mature, thermally-cracked and less-mature uncracked hydrocarbons (0.50–1.4% calculated vitrinite reflectance, R c ). The oils are less mature on the eastern flank of the structure (0.5–0.7% R c ). Based on distributions of biomarkers and higher diamondoid (composed of three adamantane cages or more) distributions, the hydrocarbons entrapped within reservoir fluid-inclusions are sourced from the Woodford Shale and also from Mississippian carbonates. To our knowledge, this is the first work to describe oil correlations using diamondoids inside fluid inclusions. Petroleum compositions in the ultra-tight Mississippian reservoirs are highly variable, which highlights the need for a thorough assessment of mechanisms controlling charge. By integrating oil and fluid inclusion geochemistry and fluid-inclusion microthermometry with basin modeling, we show that the primary controls are: (1) proximity to the source-rock kitchen(s), (2) episodic charging from Permian through Cretaceous time, and (3) mixing.
Geochemical advances in Arctic Alaska oil typing – North Slope oil correlation and charge history
The Arctic Alaska petroleum province is geologically and geochemically complex. Mixed hydrocarbon charge from multiple source rocks and/or levels of thermal maturity is common within an individual oil pool. Biomarker and chemometric statistical analyses were used to correlate twenty-nine oils to five oil families derived from: (1) Triassic Shublik Formation (calcareous organofacies), (2) Triassic Shublik Formation (shaly organofacies), (3) Jurassic Kingak Shale, (4) Cretaceous shale (pebble shale unit and Hue Shale), and (5) Paleogene shale (Canning Formation). Age-diagnostic and source-related oil biomarker parameters establish clear genetic relationships between the normal oil-window components and their putative source designations. However, application of diamondoid analyses reveals mixed-oil accumulations with postmature charge contributions (diamondoid-rich and biomarker-poor) in many oils. Most sampled reservoirs contain a predominant charge derived from a single oil-window source plus a minor contribution from one or more higher maturity source(s). Variations in source organofacies also are recognized in the Shublik, Kingak, and Cretaceous oil families. In some cases, oils from multiple pools within a single field display relatively homogeneous geochemical profiles, suggesting a common source and migration pathway. For example, oil from the significant Pikka discovery is inferred to originate mainly from the calcareous Shublik Formation. In other cases, variability among oils from multiple pools within a single field (e.g., Milne Point, Colville River, and Northstar) likely indicates a more complex source, migration, and charge history. The results of this study may be useful for anticipating the composition of oil charge in stratigraphic traps with low-permeability sandstone reservoirs, where oil gravity and other chemical parameters may influence economic viability.
Spatial heterogeneity analyses of pore structure and mineral composition of Barnett Shale using X-ray scattering techniques
The spatial heterogeneity of pore structure and mineral composition at mu m to mm scales in mm to cm-sized shale samples complicates the laboratory measurements of porosity & pore size distribution and consequently reservoir quality evaluation. Traditional fluid invasion characterization methods typically take hours to days to get data and can only test relatively small samples (usually 2.54 cm diameter core plugs and mm-sized rock chips). Furthermore, the Xray scattering techniques at the Advanced Photon Source of the Argonne National Laboratory have the capability to non-destructively determine the pore structure and mineral composition of mm x mm area in a short time of seconds and to hold several to hundreds of mm 2 sized sample to conduct large area heterogeneity analysis in the short duration of hours. With the use of ultra-small angle X-ray scattering (USAXS), small-angle X-ray scattering (SAXS), and wide-angle X-ray scattering (WAXS), the porosity, surface area and pore size distribution of two Barnett Shale samples have been mapped out and the mineral compositions are determined qualitatively. Overall, the porosity and pore surface area of the Barnett Shale Sample A over a 36 cm2 area ranges from 1.70 % to 8.13 % and 8.87 m 2 /g to 40.2 m 2 /g respectively. The porosity and surface area of Sample B over a 24 cm 2 area varies from 3.15 % to 11.7 % and 18.3 m 2 /g to 60.2 m 2 /g. In Sample A, the mineral composition shows a shift from carbonate-rich to siliciclastic-rich over a distance of several millimeters. For Sample B, the mineral compositions are mainly siliciclastic and do not show an obvious lithological change. The direct observations of high spatial heterogeneity of pore structure and mineral composition in shale illustrate the utility of integrated X-ray scattering techniques to provide valuable insights into shale studies and reservoir evaluation.
Diatom influence on the production characteristics of hydrate-bearing sediments: Examples from Ulleung Basin, offshore South Korea
The Ulleung Basin Gas Hydrate field expeditions in 2007 (UBGH1) and 2010 (UBGH2) sought to assess the Basin's gas hydrate resource potential. Coring operations in both expeditions recovered evidence of gas hydrate, primarily as fracture-filling (or vein type) morphologies in mainly silt-sized, fine-grained sediment, but also as pore-occupying hydrate in the coarser-grained layers of interbedded sand and fine-grained systems. A commonality across many of these occurrences is the presence of diatoms in the fine-grained sediment. Here we tested fine-grained sediment (median grain size <12.5 μm) associated with hydrate occurrences at four UBGH2 sites (UBGH2-2-2, UBGH2-3, UBGH2-6 and UBGH2-11) to investigate potential impacts of diatoms on efforts to extract methane from hydrate, or to tap hydrocarbon reservoirs beneath hydrate-bearing sediment. Two key considerations are: the extent to which diatoms control sediment mechanical properties, and the extent to which pore-water freshening, which occurs as gas hydrate breaks down during resource extraction, alters the diatom control on sediment mechanical properties. We conducted experiments to measure sediment index properties, sedimentation behavior and compressibility to address these considerations. We relied on scanning electron microscope (SEM) imagery and X-ray powder diffraction (XRD) to characterize the sediment mineralogy. Our high-level findings are that at the ~20–45% (by volume) diatom concentrations observed at these UBGH2 sites, sediment compressibility increases with diatom content, but diatoms only appear to increase porosity and permeability at the highest diatom concentration (~45%). Our measurements suggest in situ compression indices of 0.35–0.55 and permeabilities on the order of 0.01milliDarcies (1 × 10 -17 m 2 ) can be anticipated at these sites. Importantly, these properties are not expected to vary significantly upon pore water freshening that accompanies gas hydrate dissociation during production.
Illitization in the Mt. Simon Sandstone, Illinois Basin, USA: Implications for carbon dioxide storage
Clay mineral cementation is one of the most important controls on sandstone reservoir properties. Here, in this paper, the diagenetic history of the Mt. Simon reservoir complex is studied to reveal the origin, timing, and controls of illitization. Samples of sandstone and shale from Mt. Simon reservoir complex were acquired from the Illinois Basin–Decatur Project (IBDP), a CO 2 storage demonstration project in the central Illinois Basin. Petrographic, SEM, TEM, XRD analyses and K/Ar age dating were completed to identify the major detrital and diagenetic components of the samples and reveal illite to be the major clay component in all samples. Illitic clay coatings in the lower Mt. Simon reservoir are identified as a major control on reservoir properties by inhibiting major precipitation of authigenic quartz during illitization, resulting in highly permeable sandstone, essential for CO 2 storage. The coating box-work morphology and mineralogy are indicative of a detrital smectite origin with subsequent illite growth associated with feldspar dissolution and kaolinite alteration. The mineralogy of bulk material and separate grain size fractions (2–0.6 μm; 0.6–0.2 μm; < 0.2 μm) with illite polytypes 2M 1 , 1M, and 1M d were quantified and age dated via 40 K- 40 Ar methods. The shales or reservoir seals contain the highest proportions of detrital illite with illite in the lower Mt. Simon Sandstone reservoir identified as solely diagenetic. Two major events of illitization are identified throughout the Mt. Simon with more porous reservoir rock exhibiting the older event from approximately 360 to 315 Ma and tighter sandstone units with reservoir properties too low to be considered reservoir exhibiting illite dates from 250 to 220 Ma. This partitioning of illitization is attributed to varied mineralogy controlled by depositional changes and evolution of the greater basin.
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.
Hydrate-bearing sands in the Terrebonne Basin record the transition from ponded deposition to bypass in the deep-water Gulf of Mexico
Here, we integrated well and seismic data to interpret the stratigraphic architecture and the paleogeographic evolution of two hydrate-bearing sand reservoirs in the Terrebonne Basin in the deep-water Gulf of Mexico. These sands were deposited during the final stage of ponded fill within the Terrebonne Basin. The Green sand filled the final accommodation in the basin as a perched apron; it was formed coincidentally with downstream erosion and sediment bypass which cut through a canyon formed on the southeastern edge of the basin. Subsequently, a low energy turbidite channel continued to bypass sediment through the basin and aggrade. The Orange sand was then deposited as a second perched apron and it was once again incised by the low energy turbidite channel. The Green and Orange sands are regional sheet sands that extend and thicken deep into the Terrebonne Basin. They provide a connected regional aquifer that captured and focused gas migration to the crest of the structure where hydrate was formed. These sands were formed by large, unconfined, turbidity flows. We envision that they will be composed of coarser grained deposits relative to previously explored hydrate reservoirs within levee deposits. As such, they have the potential to be relatively coarse grained and hence have higher permeability: a condition favorable to production of hydrate deposits.
Hydrocarbon production reduces natural methane seeps in the Santa Barbara channel
Abundant natural hydrocarbon seepage in the Santa Barbara Channel includes large amounts of gas, mainly methane. Numerous seeps are focused along a faulted anticlinal axis where in 1982 ARCO placed two large (900 m 2 each) steel gas capture tents on the seabed adjacent to Platform Holly. The seep rate into the tent at high tide is about 5 percent less compared with low tide (2m tidal range). An earthquake in 2013 (M = 4.8) did not affect seepage rates into the steel tents, indicating that gas was not released by fault dilation. The CO 2 and H 2 S content of the captured gas are anomalously low compared to typical Holly gas due to bacterial and sea water-gas interaction at the seep tent. Platform Holly production had an overall effect of reducing gas seepage into the tents over time. Wells beneath the tents at 1-km depth had a much greater influence on the seepage rate than more distal wells. The nearby wells had a direct effect on seepage rate when they were put on production (seepage decreased) versus shut in (seepage increased). The response time was within hours, indicating a direct hydraulic connection between the reservoir at 1 km subsea and the seabed. Permeability of the approximate 1-km fault pathway between the reservoir and the seep tent is calculated to be about 3 × 10 -14 m 2 (30 millidarcys). Gradually, the seep rate declined from peak flow in 1990 (4.5×10 4 m 3 dy-1) to less than 7.1×10 3 m 3 dy -1 in 2013. In June of 2013, a horizontal well was completed high on the structural anticline between Platform Holly and directly beneath the seep tents. Within three weeks of turning on this new well, seepage into the tents completely ceased, ending more than 30 years of natural methane seepage captured by the tents.
A solid-state extrusion pathway extensible to upcycling aluminum scrap and end-of-life permanent magnets into ductile Al-6061/SmCo 5 structural magnetic composites
Permanent magnets offer strong potential for multifunctional structural components, but their inherent brittleness limits direct integration. Here, SolidStir® Extrusion (SSE) is used to embed SmCo5 powder into AA6061 in the solid state, and the resulting microstructure, interfacial chemistry, mechanical properties, and magnetic response are systematically characterized. SSE fragments the initial SmCo 5 powder into predominantly fine particles dispersed within an ultrafine Al matrix, with a volume fraction of ∼ 4–6 % in this proof-of-concept study. Interfacial reactions form a continuous Al–Co reaction layer, with local atomic-resolution indexing consistent with Al 13 Co 4 , together with Sm–Si-rich reaction products, while the matrix contains Sm–Co–Si-rich precipitates and Mg 2 Si, producing features from nanometers to micrometers. Compared with an Al extrudate, this architecture increases yield strength by ∼ 42 % and tensile strength by ∼ 31 %, while preserving uniform elongation, and it sustains work hardening in compression from 150 to 400 °C. The composite retained hard-magnetic behavior, with axial magnetization of 12.7 A·m 2 kg −1 at 6 T and intrinsic coercivity of 9.6 kOe. These findings demonstrate that SSE can produce Al-based load-bearing magnetic composites while preserving useful SmCo5 functionality.
Pd modified prussian blue frameworks: Multiple electron transfer pathways for improving catalytic activity toward hydrogenation of nitroaromatics
Prussian blue analogs (PBAs) exhibit potential as low-cost and eco-friendly nanocatalysts that can be fabricated with ease. However, the PBA framework structure suffers from poor electronic conductivity, which limits the catalytic efficiency for this class of materials. Noble metals represent an alternative class of materials that display inherent catalytic activity but suffer from aggregation, ultimately reducing the amount of accessible catalytic sites. Herein, we demonstrate a combinatory approach that circumvents the known disadvantages with these classes of catalytic materials in which PBA-supported nanocatalysts were synthesized. These composite materials exhibit excellent catalytic activity for the reduction of nitroaromatics to aminoaromatics while displaying long-term cycling stability, which is attributed to the availability of multiple electron transfer pathways. Overall, this work opens the study on the assembly of PBA-supported heterogeneous nanocatalysts and potentially paves the way toward future applications.
Fabric controls on fracture surface roughness of an architected rock material
Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.
Development of a Flotation Simulator That Can Predict Grade vs. Recovery Curves from Mineral Liberation Data
Flotation is a kinetic process designed to separate mineral fines by control of hydrophobicity. In general, both the recovery and selectivity of a flotation process increase with contact angle. It is difficult, however, to model flotation kinetics using the thermodynamic parameter. In the present work, we consider flotation as heterocoagulation so that the bubble-particle interactions occurring in the pulp phase of a flotation cell can be modeled using the hydrophobic force as a kinetic parameter that can represent the role of contact angles in flotation. Here, the pulp-phase model developed in this manner has been combined with a froth-phase model developed earlier by considering the effect of particles in film thinning kinetics and bubble coarsening (Park et al, 2018). The combined model has been used to develop a flotation simulator that can be used to construct grade vs. recovery curves from feed characteristics (e.g., mineral liberation and particle size) and various operating parameters. The simulator has been tested successfully against a set of bench-scale continuous flotation test results reported in the literature (dos Santos and Galery, 2018).