Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Shale”

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.

At least 451 records · Page 25

Marcellus assessment shows continued expansion

Increasing uncertainty regarding the remaining volume of recoverable gas in Appalachian shale plays provided the original motivation for this study. Assessments of the Marcellus and Utica shale plays published by the US Geological Survey (USGS) in 2019 established resource volumes dramatically lower than ongoing well performance would indicate. Further, as production in the Marcellus increased and wells matured, operators observed that wells were producing more than the estimated gas-in-place (GIP).

03 NATURAL GAS↗

Deposition, Diagenesis, and Sequence Stratigraphy of the Pennsylvanian Morrowan and Atokan Intervals at Farnsworth Unit

Farnsworth Field Unit (FWU), a mature oilfield currently undergoing CO2-enhanced oil recovery (EOR) in the northeastern Texas panhandle, is the study area for an extensive project undertaken by the Southwest Regional Partnership on Carbon Sequestration (SWP). SWP is characterizing the field and monitoring and modeling injection and fluid flow processes with the intent of verifying storage of CO2 in a timeframe of 100–1000 years. Collection of a large set of data including logs, core, and 3D geophysical data has allowed us to build a detailed reservoir model that is well-grounded in observations from the field. This paper presents a geological description of the rocks comprising the reservoir that is a target for both oil production and CO2 storage, as well as the overlying units that make up the primary and secondary seals. Core descriptions and petrographic analyses were used to determine depositional setting, general lithofacies, and a diagenetic sequence for reservoir and caprock at FWU. The reservoir is in the Pennsylvanian-aged Morrow B sandstone, an incised valley fluvial deposit that is encased within marine shales. The Morrow B exhibits several lithofacies with distinct appearance as well as petrophysical characteristics. The lithofacies are typical of incised valley fluvial sequences and vary from a relatively coarse conglomerate base to an upper fine sandstone that grades into the overlying marine-dominated shales and mudstone/limestone cyclical sequences of the Thirteen Finger limestone. Observations ranging from field scale (seismic surveys, well logs) to microscopic (mercury porosimetry, petrographic microscopy, microprobe and isotope data) provide a rich set of data on which we have built our geological and reservoir models.

04 OIL SHALES AND TAR SANDS↗

Pyrolysis of Oils from Unconventional Resources

In this study, oils from various sources were subjected to pyrolysis conditions; that is, without oxidizer, as the samples were heated to 500 °C, and held at that temperature. The oils studied included: (1) heavy oil from Grassy Creek, Missouri; (2) oil from tar sands of Asphalt Ridge in Utah; (3) mid-continent oil shales of three formations (two of Chattanooga formation, Pennsylvanian (age) formation, and Woodford formation); and (4) a Colorado Piceance Basin shale. Differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA) with either gas chromatography (GC) or mass spectrometry (MS) were used to quantify the produced gases evolved in the tests. Purge gases of helium, argon, and humid carbon dioxide were utilized. Larger scale pyrolysis tests were conducted in a tube furnace coupled to a MS and a GC. The results consistently showed that pyrolysis occurred between 300 °C and 500 °C, with the majority of gases being mainly hydrogen and light alkanes. This behavior was essentially consistent, regardless of the oil source.

42 ENGINEERING↗

Temporal patterns of fine‐root dynamics have little influence on seasonal soil CO 2 efflux in a mixed, mesic forest

Among the contributors to soil CO 2 efflux, there remains uncertainty about the contribution of root activity to the overall soil efflux. Soil water and temperature frequently have been used to predict a large portion of the variation in soil CO 2 efflux. We hypothesized that fine-root dynamics explain most of the remaining variability in soil CO 2 efflux that cannot be explained by soil temperature and water content. We anticipated that seasonal increases in root production, mortality via decomposition, and standing crop would result in corresponding increases in soil CO 2 efflux. We tested our hypotheses by collecting and analyzing two years of minirhizotron and soil chamber CO 2 flux data from plots distributed throughout the Shale Hills Catchment of the Susquehanna-Shale Hills Critical Zone Observatory in Central Pennsylvania, USA. Here we showed that: (1) seasonal fluctuations in fine-root dynamics yielded only a very small increase in the predictability of soil CO 2 efflux; (2) fine-root mortality effects on soil CO 2 efflux were strongly tied to soil temperature; (3) fluctuations in fine-root presence or standing mass independent of temperature and moisture had little effect on soil CO 2 efflux; and (4) new fine-root length and root length mortality had limited impacts on soil CO 2 efflux rates. We conclude that, at least in temperate forests on rocky soils, characterizing fine-root dynamics may provide only limited improvement in the estimation of soil CO 2 efflux.

54 ENVIRONMENTAL SCIENCES↗

Manganese uptake by red maples in response to mineral dissolution rates in soil

Manganese (Mn) is an essential plant micronutrient that influences photosynthesis, ecosystem productivity, and soil carbon storage. Our objective was to quantify how Mn uptake by forest vegetation relates to Mn release into soil solution through mineral dissolution. In this work, a greenhouse pot experiment was conducted to quantify Mn uptake by red maple saplings as a function of mineral solubility to test whether Mn uptake was limited by the supply of Mn to soil solution. Differences in soil microbial community composition between treatments, particularly amongst Mn cycling bacteria and fungi, were also evaluated to assess potential microbial impacts on observed Mn fluxes. Plant Mn uptake was highest in systems supplied with dissolved Mn(II) because it was not kinetically limited by mineral weathering. Mn uptake was also higher in systems supplied with a fast-weathering substrate (shale containing Mn(II)-bearing pyrite) than a slow-weathering substrate (Mn(IV)-oxide). However, vegetation accumulated a decreasing proportion of available Mn with increasing solubility, indicating that uptake was tempered relative to solubility. The presence of bacterial phyla containing known Mn-oxidizing bacteria indicates the potential for Mn-oxidizing bacteria to influence Mn solubility within Mn-oxide and dissolved Mn treatments. A relatively low abundance of these bacteria points to their limited capacity to oxidize and retain Mn in the shale treatments, consistent with substantial Mn leaching. We conclude that Mn uptake was primarily controlled by dissolution rates of Mn-bearing minerals but modified by biological processes. The quantitative framework presented here can guide understanding of how biogeochemical processes control element cycling between plants and soils.

59 BASIC BIOLOGICAL SCIENCES↗

Incorporating geological structure into sensitivity analysis of subsurface contaminant transport

Simulating subsurface contaminant transport at the kilometer-scale often entails modeling reactive flow and transport within and through complex geologic structures. These structures are typically meshed by hand and as a result geologic structure is usually represented by one or a few deterministically generated geological models for uncertainty studies of flow and transport in the subsurface. Uncertainty in geologic structure can have a significant impact on contaminant transport. In this study, the impact of geologic structure on contaminant tracer transport in a shale formation is investigated for a simplified generic deep geologic repository for permanent disposal of spent nuclear fuel. An open-source modeling framework is used to perform a sensitivity analysis study on transport of two tracers from a generic spent nuclear fuel repository with uncertain location of the interfaces between the stratum of the geologic structure. The automated workflow uses sampled realizations of the geological structural model in addition to uncertain flow parameters in a nested sensitivity analysis. Concentration of the tracers at observation points within, in line with, and downstream of the repository are used as the quantities of interest for determining model sensitivity to input parameters and geological realization. Finally, the results of the study indicate that the location of strata interfaces in the geological structure has a first-order impact on tracer transport in the example shale formation, and that this impact may be greater than that of the uncertain flow parameters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A critical literature review of the advances in methane dehydroaromatization over multifunctional metal-promoted zeolite catalysts

We report the drastic rise in shale gas production has encouraged the quest for alternative uses of methane as a chemical feedstock in the manufacturing industry. While two-step syngas routes for methane valorization are deployed commercially, direct one-step routes for methane conversion are attracting much attention. As steam cracking installations have shifted from using oil-based naphtha to shale-based natural gas liquids, production of aromatics has dropped. Methane dehydroaromatization (MDA) is a one-step reaction capable of valorizing methane to hydrogen and benzene. Challenges with the MDA reaction are two-fold: the reaction is thermodynamically limited with low one-pass methane conversion and even the best catalytic systems, Mo/zeolites, suffer rapid deactivation from coking. A catalyst design strategy to improve stability is the use of multifunctional Mo-X/zeolite systems where X is a dopant capable of modulating the stability. In this paper we provide a complete overview of the main Mo-X/zeolite systems used in MDA and critically draw connections among the different types of dopants (X) employed, as a function of the role they play in the reaction/deactivation pathway. We have also dedicated a section to emerging trends with non-Mo based catalysts. The goal of this review article is to establish a basis that will facilitate the identification of useful multifunctional catalytic systems, and recognize gaps in the knowledge of these systems that deserve more attention. Improving MDA systems to the point to which they can be commercially deployed requires a multifaceted approach that combines optimization of the designs of both the catalyst and the reactor configuration. We therefore also provide a brief overview of the most recent advances in process intensification strategies employed with different reactor configurations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Determination of transition metal ions in fossil fuel associated wastewaters using chelation ion chromatography

Here, this study outlines the development and subsequent validation of a method using chelation ion chromatography (CIC) pretreatment followed by traditional ion chromatography (IC) and post column UV/vis detection to measure transition metals in fossil fuel wastewaters, such as oil & gas (O&G) brines and coal mine drainage (CMD) waters. Measurement of transition metals is often an important characterization step in the research of environmental and energy systems. IC represents one way to measure these metals with the advantages of being versatile, simple and relatively low cost compared to other analytical methods. However, high concentrations of alkali and alkaline earth metals present in fossil fuel wastewaters will decrease IC detectability of transition metals in these waters. In this study, a CIC method was developed for the analysis of transition metal ions (Fe 3+ , Cu 2+ , Ni 2+ , Zn 2+ , Co 2+ , Mn 2+ , and Fe 2+ ) in fossil fuel associated wastewaters such as Appalachian CMD and O&G wastewaters from the Permian and Bakken shale basins in the United States. CIC system incorporated an on-line chelator column (e.g., the MetPac CC-1) with high selectivity for transition metals over alkali and alkaline earth metals for salt matrix removal prior to transition metal separation and detection. Additional method developments also included acidifying all samples to 2% v/v HCl and using gradient elution rather than isocratic. The recoverability of transition metals in simple salt solutions commonly found in CMD and brine samples (e.g. NaCl, Na 2 SO 4 , CaCl 2 ) using CIC was evaluated and compared to that using traditional IC. Our results found that the CIC system significantly improved transition metal recoveries for samples in 10,000 mg/L CaCl 2 matrix, reaching 87%-108% recovery for all analytes, as opposed to 2-323% recovery in traditional IC. The limits of detection in this study achieved 10.09 – 161.2 μg/L, comparable to reported values in similar IC studies. The developed method was also verified with certified water samples, resulting in 89% - 111% recoveries in samples with higher analyte concentrations (i.e. >4x the LoDs). The developed method achieved 87% -112% recoveries for most analytes in CMD samples and 72%-138% recoveries for Bakken shale samples, relative to ICP-MS values. Overall, the current IC method can be a very good screening tool for fast and cheap analysis for transition metals at mg/L level, to facilitate selection of samples for more detailed ICP-MS analysis.

01 COAL, LIGNITE, AND PEAT↗

How to avoid multiple scattering in strongly scattering SANS and USANS samples

Small Angle Neutron Scattering (SANS) and Ultra Small Angle Neutron Scattering (USANS) are the only available experimental techniques to provide seamless non-destructive measurements of the geometry of the accessible and inaccessible pore structure of rocks from sub-nanopore size to the scale of macropores. They have therefore become the measurement of choice for tight reservoir rocks such as organic rich shales. A simplifying assumption in the analysis is, however, that during the path of neutrons through the sample each neutron is only scattered once. Shales are samples with a high scattering power and Multiple Scattering (MS) may occur which requires special modelling for deconvolution of the results. The approach to avoid MS is to simply reduce the sample thickness to <0.15–0.5 mm. Here, in this work, we present a systematic method on wavelength selection and preparation of samples to optimise extraction of microstructural data and minimise parasitic errors. Experimentally measured SAS transmission (TSAS) values are used as a practical criterion for estimation of the extent of MS. Generous beamtime allocations allowed robust testing revealing that sample thicknesses can be twice as thick as predicted using the standard protocol. Analysing thicker samples is particularly beneficial for statistically relevant characterisation of heterogeneous samples making the new protocol the method of choice for such samples.

(U)SANS↗

A dynamic solvent chamber propagation estimation framework using RNN for warm solvent injection in heterogeneous reservoirs

Warm solvent injection (WSI), injecting low-temperature solvent into formations to reduce the viscosity of heavy oil, is a clean technology for heavy oil production through reducing greenhouse gas emissions and water usage. The success of WSI operation depends on the uniform development and propagation of solvent chambers in reservoirs. However, reservoir heterogeneity stemming from shale barriers plays a detrimental role in the conformance of solvent chamber development and oil production rate. In this work, we developed a novel recurrent neural network (RNN)-based framework with the capability of efficiently tracking and estimating the solvent chamber positions in heterogeneous reservoirs based on only production time-series data. The developed estimation model utilizes the “sequence-to-sequence" mapping methodology to correlate observed production time-series sequence and solvent chamber edge sequence via a long short-term memory (LSTM) algorithm. The trained RNN models exhibit high accuracy, evidenced by the predicted dynamic solvent chamber locations match the corresponding true locations from numerical simulation, with a high coefficient of determination (R 2 ) and a low mean squared error. Specifically, the achieved R 2 values exceed 0.98 on both the training and testing data. The developed RNN-based workflow was tested via several cases from both regularly- and irregularly-shaped shale barriers, and the results were promising. The predicted solvent chambers showed strong agreement with those obtained from numerical simulations. The major benefits of this workflow include reducing computational time and saving overall monitoring and tracking costs for conventional techniques. In conclusion, the present work would provide a good demonstration of the capability of practical integration of machine learning methods in solving engineering problems.

58 GEOSCIENCES↗

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗

Role of interaction between hydraulic and natural fractures on production

One of the main purposes of hydraulic fracturing in unconventional resources such as shale is to improve the connectivity between existing natural fractures and the production well. Since the inherent permeability of shale is extremely low, this newly formed fracture connectivity provides faster paths for the hydrocarbon stored in preexisting natural fractures. In this study, we performed numerical simulations of the free gas production at large scale fractured reservoir with different pre-existing natural fracture network intensities and with variable hydraulic fracture sizes and staging. We explored production characteristics as a function of the properties of both, natural fractures (intensity) and hydraulic fractures (size and staging). Our hypothesis is that once maximum fracture connectivity is achieved, additional stimulation does not increase production significantly. In this numerical study we use an advanced modeling tool to study fractures connectivity, Discrete Fracture Network (DFN) approach, which can represent fracture networks similar to those observed on individual reservoir sites. In our model the horizontal production well and vertical hydraulic fractures were deterministically defined and surrounded by stochastically generated natural fracture networks. The performed numerical experiments show how fracture connectivity depends on hydraulic fracture settings in fractured reservoir with different preexisting fracture densities. We numerically showed that cumulative free gas production is not affected significantly by hydraulic fractures size (25 m–150 m) or spacing (15 m–45 m) for dense natural fracture networks. However, the hydraulic fracture size is important in reservoirs with sparse natural fractures, although spacing (number of stages) does not play a significant role.

58 GEOSCIENCES↗

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

Characterization of the mud displacement in an enlarged wellbore: An integrated rock-fluid model

Cement-mud displacement plays a crucial role in the sealability of cement sheaths. Irregular geometric features of a wellbore due to washout can have a negative impact on mud and cement mobilization. An unstable interface between two fluids always leads to mud channeling, interfluid mixing, and cement contamination, degrading the cement quality. Many factors, such as mechanical and rheological properties of fluids, annulus geometry, flow pattern, and flow rate, significantly influence the displacement efficiency. This study investigates the characterization of the mud displacement in an irregular horizontal well using a 3D computational fluid dynamics (CFD) model. Mud is displaced in an enlarged wellbore by geopolymer and neat class G cement. The wellbore geometry is developed based on the caliper log data from an unconventional shale well in the Tuscaloosa Marine Shale (TMS) lithology. The effects of pump rate, density difference, and mud contamination are evaluated by numerical simulations. The results present those residual muds mainly exist in the upper annulus of the enlarged section. Geopolymer has a better sealing performance and can resist more water-based mud (WBM) contaminations than neat class G cement. The scenario with a low mud-cement density difference and high cement injection rate results in a high cement volume fraction, mitigating the gas migration.

58 GEOSCIENCES↗

Can CO 2 -assisted alkane dehydrogenation lead to negative CO 2 emissions?

The utilization of light alkanes to generate value-added chemicals has gained significant attention due to the U.S. shale gas boom in recent decades. For example, after methane, ethane (3-16 vol.%) and propane (up to 5 vol.%) are the most abundant hydrocarbon species in shale gas deposits. Further, these inexpensive, widely available, but often underutilized, light alkanes can be used as a feedstock to produce olefins, which are among the most important petrochemical building blocks in chemical industry. Ethylene and propylene are used in the production of polymers (e.g., polyethylene, polypropylene) and epoxides (e.g., ethylene oxide, propylene oxide), as well as intermediates in the manufacture of numerous other chemicals (e.g., ethylene glycol, ethylbenzene, acetone, isopropanol, acrylonitrile).

03 NATURAL GAS↗

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.

58 GEOSCIENCES↗

Pressure-stable supported ionic liquid membranes using isoporous supports for evaluating pure- and mixed-gas light paraffin fractionation

Advances in horizontal drilling and hydraulic fracturing have spurred the growth of domestic U.S. energy production. Membranes, typically silicone rubbers, have found utility in shale gas treatment as fuel gas conditioning units to selectively remove C 2 + hydrocarbons at pressures up to 30 bar, producing clean CH 4 for gas engines. However, more selective materials could be beneficial for broader shale gas treatment applications, such as dew point control units. Supported ionic liquid membranes (SILMs) offer a potential opportunity for improving C 3 H 8 /CH 4 selectivity, but they lack pressure-stability. Here, we report C 3 H 8 /CH 4 selective and pressure-stable SILMs using isoporous supports. SILMs with supports that had minimal defects remained stable up to 15 bar of transmembrane pressure. This stability allowed for pure- and mixed-gas testing of the resulting membranes at elevated pressures. Furthermore, these tests revealed that low viscosity ILs (<100 cp) may display mixed-gas permeances and permselectivity nearly identical to pure-gas results. On the other hand, higher viscosity ILs may display increasing permeances and permselectivity with increasing C 3 H 8 fugacity, similar to rubbery polymers. Ultimately, the SILMs demonstrated relatively high pressure-stability due to the isoporous supports and competitive mixed-gas C 3 H 8 /CH 4 permselectivity compared to silicone rubber.

Rosenthal, Justin J. [University of Texas at Austi↗