Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “logging”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Brief Announcement: Communication Optimal Sparse LU Factorization for Planar Matrices

We introduce a new parallel algorithm for solving sparse LU factorization of planar matrices, which commonly arise in the finite element method for 2D PDEs. Existing scalable methods, such as the multifrontal approach with subtree-to-subcube mapping by Gupta et al. [1] and right-looking with 3D mapping by Sao et al. [2] fail to achieve optimal communication costs for these matrices. Our new algorithm combines 3D mapping and subtree-to-subcube mapping to minimize communication costs while allowing trade-offs between extra memory and reduced communication. We demonstrate that our proposed algorithm attains the communication lower bound up to a factor of O(log log n) in the memory-optimal case and up to a factor of O(log P) in the memory-independent case for an n-dimensional planar sparse matrix on P processors.

Sao, Piyush↗
DOE OSTI2022

Properly Learning Decision Trees in almost Polynomial Time

We give an n O (log log n ) -time membership query algorithm for properly and agnostically learning decision trees under the uniform distribution over { ± 1} n . Even in the realizable setting, the previous fastest runtime was n O (log n ) , a consequence of a classic algorithm of Ehrenfeucht and Haussler. Our algorithm shares similarities with practical heuristics for learning decision trees, which we augment with additional ideas to circumvent known lower bounds against these heuristics. To analyze our algorithm, we prove a new structural result for decision trees that strengthens a theorem of O’Donnell, Saks, Schramm, and Servedio. While the OSSS theorem says that every decision tree has an influential variable, we show how every decision tree can be “pruned” so that every variable in the resulting tree is influential.

Computer Science↗
DOE OSTI2024

Detecting fractures and monitoring hydraulic fracturing processes at the first EGS Collab testbed using borehole DAS ambient noise

Enhanced geothermal systems (EGS) require cost-effective monitoring of fracture networks. We validate the capability of using borehole distributed acoustic sensing (DAS) ambient noise for fracture monitoring using core photos and core logs. The EGS Collab project has conducted 10 m scale field experiments of hydraulic fracture stimulation using 50–60 m deep experimental wells at the Sanford Underground Research Facility (SURF) in Lead, South Dakota. The first EGS Collab testbed is located at 1616.67 m (4850 ft) depth at SURF and consists of one injection well, one production well, and six monitoring wells. All wells are drilled subhorizontally from an access tunnel called a drift. The project uses a single continuous fiber-optic cable installed sequentially in the six monitoring wells to record DAS data for monitoring hydraulic fracturing during stimulation. We analyze 60 s time records of the borehole DAS ambient noise data and compute the noise root-mean-square (rms) amplitude on each channel (points along the fiber cable) to obtain DAS ambient noise rms amplitude depth profiles along the monitoring wellbore. Our noise rms amplitude profiles indicate amplitude peaks at distinct depths. We compare the DAS noise rms amplitude profiles with borehole core photos and core logs and find that the DAS noise rms amplitude peaks correspond to the locations of fractures or lithologic changes indicated in the core photos or core logs. We then compute the hourly DAS noise rms amplitude profiles in two monitoring wells during three stimulation cycles in 72 h and find that the DAS noise rms amplitude profiles vary with time, indicating the fracture opening/growth or closing during the hydraulic stimulation. Our results demonstrate that borehole DAS passive ambient noise can be used to detect fractures and monitor fracturing processes in EGS reservoirs.

58 GEOSCIENCES↗
DOE OSTI2026

2014 Southern Nevada Household Travel Survey

The 2014 Southern Nevada Household Travel Survey collected information from residents in the Las Vegas area to update the regional travel demand model and assess travel behavior. The Regional Transportation Commission of Southern Nevada contracted with Westat to conduct the survey. The survey was conducted in two phases—from March to May 2014 and from August to October 2014. Participants provided demographic and travel data (via a travel log). A 10% subsample (1,694 participants) was randomly selected to take part in a wearable global positioning system (GPS) technology-based component of the study, the purpose of which was to assess the level of trip under-reporting in the self-reported travel logs. Participants in the GPS portion of the study were instructed to record trips in their travel logs on the first day only, while passively recording their travel for three full days.

1Hz data↗
DOE OSTI2021

Diurnal observations of basal stem CO2 efflux in three canopy dominant trees in the central Amazon

We explored mechanisms responsible for disturbance-induced shifts in carbon metabolism in forests north of Manaus, Brazil (S2º 38' 17", W60º 09' 25"), using data from a selective logging experiment (BIONTE). BIONTE (the BIOmass and NuTrient Experiment) selective logging treatments were initiated in the mid-1980s and comprised variable removal of commercial tree (not total) species basal area (T1 =32%, T2 = 42%, T3 = 69%), and control plots with no logging (T0). Raw data including tree base diameter (Db), CO2 concentration versus time inside a static stem chamber, and wood density were used to derive stem respiration rates (Rw), wood production rates (Pw), and wood carbon use efficiency (CUE). Changes in tree base diameter (Db; measured at 1.3 m height, or above the buttresses) was used to calculate wood production for individual trees (Pw). Wood density enabled calculation of stem growth rate in the same units as stem respiration (mmol C m-2 s-1). Four trees were randomly selected from five tree growth rate classes for each treatment block, for a total of 20 trees per treatment block, or 80 trees total from the BIONTE plots. Each of the 80 trees selected for the Rw study were outfitted with dendrometer bands (da Silva et al. 2002). Stem respiration was measured using the method described in Chambers et al. (2004). Briefly, an infra-red gas analyzer (LiCor 820) was operated as a closed dynamic chamber with a flow rate of 1.0 L min-1. Polyvinyl chloride (PVC) semi-cylindrical chambers (250–400 mL) were secured to the tree stem near the dendrometer bands using nylon straps. The measurement interval spanned 1–2 min, and the CO2 flux from the stem of each tree was quantified using the enclosed stem area and slope of the stem CO2 concentration versus time.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Growth, mortality, wood density, biomass data from BIONTE inventories in Manaus, Brazil

BIONTE (BIOmass and NuTrient Experiment) is a selective logging experiment established at the Experimental Station of Tropical Forestry (EEST, aka “ZF2”) field research station in the mid 1980s in the central Amazon (Higuchi et al. 1997, Amaral et al. 2019). Led by the National Institute for Amazon Research (INPA) in Brazil, the project aimed at assessing the effects of logging intensity on forest dynamics and enabling the creation of a model of forest management for the Central Amazon. The experiment included three levels of increasing selective logging intensity and controls, with 1 hectare sample plots (12 total) located at the center of 4 hectare treatment plots. The site’s Köppen classification is tropical rainforest (Af), characterized by high temperatures and humidity, with mean annual temperatures around 27 ℃ and mean annual precipitation around 2200 mm of rain. The vegetation has a high floristic diversity, the soils of the region are poor in nutrients, and the topography is characterized by plateaus (where BIONTE is located), and also valley bottoms and slopes. The inventory (growth and mortality) and biomass data included here covers the 1990 to 2019 period, with wood density being averaged from existing datasets. This dataset includes a data file in .csv file format and a .txt file, BIONTE_mortality-rates_headers.txt, that provides descriptions for the data file headers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

AmeriFlux US-Me2 Metolius-intermediate aged ponderosa pine

This is the AmeriFlux version of the carbon flux data for the site US-Me2 Metolius-intermediate aged ponderosa pine. Site Description - Site Description before Fire (January, 2002 - August, 2020): The mean stand age is 71 years old and the stand age of the oldest 10% of trees is about 108 years old. This site is one of the Metolius core cluster sites with different age and disturbance classes and part of the AmeriFlux network. The overstory is almost exclusively composed of ponderosa pine trees (Pinus ponderosa Doug. Ex P. Laws) with a few scattered incense cedars (Calocedrus decurrens (Torr.) Florin) and has a peak leaf area index (LAI) of 2.1 m2 m-2. Tree height is relatively homogeneous at about 18 m, and the mean tree density is approximately 339 trees ha-1 (Irvine et al., 2008). The understory is sparse with an LAI of 0.2 m2 m-2 and primarily composed of bitterbrush (Purshia tridentata (Push) DC.) and greenleaf manzanita (Arctostaphylos patula Greene). Soils at the site are sandy (69%/24%/7% sand/silt/clay at 0–0.2 m depth and 66%/27%/7% at 0.2–0.5 m depth, and 54%/ 35%/11% at 0.5–1.0 m depth), freely draining with a soil depth of approximately 1.5 m (Irvine et al., 2008; Law et al., 2001b; Schwarz et al., 2004). Green Ridge Fire: On August 20, 2020, the Green Ridge Fire burned through Us-Me2. The fire was ignited by lightning on August 16th, and grew rapidly to the east over the first few days driven by strong, downslope (westerly) afternoon winds. Fire behavior and observed fire effects were highly heterogeneous due to the localized wind pattern carrying the flaming head of the fire forward, and the efforts being made by suppression resources to contain the fire. The site experienced the full range of fire effects, from <1 m high surface fire that charred litter and duff and only consumed shrubs and herbaceous material to full tree (>15 m) crown fire that consumed 100% of needles, small limbs, and surface fuels at high intensity, leaving only ash and bare soil post-fire. Salvage Logging: From late March to late April 2021, salvage logging by the landowner occurred at the site. Almost all trees within the flux footprint were logged except a small area with lower burn severity, where sap flow and automatic soil respiration measurements are continued since the fire. In August 2022 the primary flux system was moved from the top of the damaged tall tower to a nearby shorter tower centering on the salvage/regenerating footprint.

Hanson, Chad↗
DOE OSTI2022

AmeriFlux FLUXNET-1F US-Me2 Metolius mature ponderosa pine

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-Me2 Metolius mature ponderosa pine. This is the FLUXNET version of the carbon flux data for the site US-Me2 Metolius mature ponderosa pine produced by applying the standard ONEFlux (1F) software. Site Description - Site Description before Fire (January, 2002 - August, 2020): The mean stand age is 71 years old and the stand age of the oldest 10% of trees is about 108 years old. This site is one of the Metolius core cluster sites with different age and disturbance classes and part of the AmeriFlux network. The overstory is almost exclusively composed of ponderosa pine trees (Pinus ponderosa Doug. Ex P. Laws) with a few scattered incense cedars (Calocedrus decurrens (Torr.) Florin) and has a peak leaf area index (LAI) of 2.1 m2 m-2. Tree height is relatively homogeneous at about 18 m, and the mean tree density is approximately 339 trees ha-1 (Irvine et al., 2008). The understory is sparse with an LAI of 0.2 m2 m-2 and primarily composed of bitterbrush (Purshia tridentata (Push) DC.) and greenleaf manzanita (Arctostaphylos patula Greene). Soils at the site are sandy (69%/24%/7% sand/silt/clay at 0–0.2 m depth and 66%/27%/7% at 0.2–0.5 m depth, and 54%/ 35%/11% at 0.5–1.0 m depth), freely draining with a soil depth of approximately 1.5 m (Irvine et al., 2008; Law et al., 2001b; Schwarz et al., 2004). Green Ridge Fire: On August 20, 2020, the Green Ridge Fire burned through Us-Me2. The fire was ignited by lightning on August 16th, and grew rapidly to the east over the first few days driven by strong, downslope (westerly) afternoon winds. Fire behavior and observed fire effects were highly heterogeneous due to the localized wind pattern carrying the flaming head of the fire forward, and the efforts being made by suppression resources to contain the fire. The site experienced the full range of fire effects, from <1 m high surface fire that charred litter and duff and only consumed shrubs and herbaceous material to full tree (>15 m) crown fire that consumed 100% of needles, small limbs, and surface fuels at high intensity, leaving only ash and bare soil post-fire. Salvage Logging: From late March to late April 2021, salvage logging by the landowner occurred at the site. Almost all trees within the flux footprint were logged except a small area with lower burn severity, where sap flow and automatic soil respiration measurements are continued since the fire. In August 2022 the primary flux system was moved from the top of the damaged tall tower to a nearby shorter tower centering on the salvage/regenerating footprint.

Hanson, Chad↗
DOE OSTI2017

Stan : A Probabilistic Programming Language

Stan is a probabilistic programming language for specifying statistical models. A Stan program imperatively defines a log probability function over parameters conditioned on specified data and constants. As of version 2.14.0, Stan provides full Bayesian inference for continuous-variable models through Markov chain Monte Carlo methods such as the No-U-Turn sampler, an adaptive form of Hamiltonian Monte Carlo sampling. Penalized maximum likelihood estimates are calculated using optimization methods such as the limited memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Stan is also a platform for computing log densities and their gradients and Hessians, which can be used in alternative algorithms such as variational Bayes, expectation propagation, and marginal inference using approximate integration. To this end, Stan is set up so that the densities, gradients, and Hessians, along with intermediate quantities of the algorithm such as acceptance probabilities, are easily accessible. Stan can also be called from the command line using the cmdstan package, through R using the rstan package, and through Python using the pystan package. All three interfaces support sampling and optimization-based inference with diagnostics and posterior analysis. rstan and pystan also provide access to log probabilities, gradients, Hessians, parameter transforms, and specialized plotting.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Development of Time Lapse VSP Integration Workflow: A Case Study at Farnsworth CO2-EOR Project

Abstract This study aims to develop a 4D Vertical Seismic Profile (VSP) integration workflow to improve the prediction of subsurface stress changes. The selected study site is a 5-spot pattern within the ongoing CO2-EOR operations at the Farnsworth Field Unit FWU in Ochiltree County, Texas. The specific pattern has undergone extensive geological and geomechanical characterization through the acquisition of 3D seismic data, geophysical well logs, and core. This workflow constrains a numerical hydromechanical model by applying a penalty function formed between "modeled" versus "observed" time-lapse compressional and shear seismic velocity changes. Analyses of geophysical logs and ultra-sonic measurements on core exhibit measurable sensitivities to changes in both fluid saturation and mean effective stress. These data are used to develop a site-specific rock physics model and stress-velocity relationship, which inform the numerical models used to generate the "modeled" portion of the penalty function. The "observed" portion of the penalty function is provided by a novel elastic full-waveform inversion of the available 3D baseline and three monitor surveys to produce high-quality estimates of time-lapse compressional and shear seismic velocity changes. The modeling workflow accounts sequentially for fluid substitution and stress impacts. Hydrodynamic and geomechanical properties of the 3D coupled numerical model are estimated through geostatistical integration of well log and core data with 3D seismic inversion products. Changes in seismic velocities due to fluid substitution are computed using the Biot-Gassmann workflow and site-specific rock physics. Stress impacts on time-lapse seismic velocity changes are modeled from the effective stress output of the hydromechanical model and are initially based on the velocity versus effective stress relationship extracted from core mechanical testing. Based on the principle of superposition of seismic wavefields, seismic velocity changes attributed to fluid substitution and that due to changes in mean effective stress are treated as linearly additive. The modeled results are upscaled using Backus averaging to reconcile scale discrepancies between the modeled and measured datasets to formulate the penalty function. This manuscript presents the forward modeling process and concludes that for the base case, the seismic velocity changes due to mean effective stress dominates over the seismic velocity changes attributed to fluid substitution because of the extensive range of the pressure perturbations. Successful minimization of this penalty function calibrates the coupled hydrodynamic geomechanical numerical model and affirms the suitability of acoustic time-lapse measurements such as 4D-VSP for geomechanical calibration.

02 PETROLEUM↗
DOE OSTI2020

Equations of state, phase relations, and oxygen fugacity of the Ru-RuO 2 buffer at high pressures and temperatures

Experimental studies and measurements of inclusions in diamonds show that ferric iron components are increasingly stabilized with depth in the mantle. To determine the thermodynamic stability of such components, their concentration needs to be measured at known oxygen fugacities. The metal-oxide pair Ru and RuO 2 are ideal as an internal oxygen fugacity buffer in high-pressure experiments. Both phases remain solid to high temperatures and react minimally with silicates, only exchanging oxygen. To calculate oxygen fugacities at high pressure and temperature, however, requires information on the phase relations and equation of state properties of the solid phases. We have made in situ synchrotron X-ray diffraction measurements in a multi-anvil press on mixtures of Ru and RuO 2 to 19.4 GPa and 1473 K with which we have determined phase relations of the RuO 2 phases and derived thermal equations of state (EoS) parameters for both Ru and RuO 2 . Rutile-structured RuO 2 was found to undergo two phase transformations, first at ~7 GPa to an orthorhombic structure and then above 12 GPa to a cubic structure. The phase boundary of the cubic phase was constrained for the first time at high pressure and temperature. We have derived a continuous Gibbs free energy expression for the tetragonal and orthorhombic phases of RuO 2 by fitting the second-order phase transition boundary and P-V-T data for both phases, using a model based on Landau theory. The transition between the orthorhombic and cubic phases was then used along with EoS terms derived for both phases to determine a Gibbs free energy expression for the cubic phase. We have used these data to calculate the oxygen fugacity of the Ru + O 2 = RuO 2 equilibrium, which we have parameterized as a single polynomial across the stability fields of all three phases of RuO 2 . The expression is log 10 f O 2 (Ru – RuO 2 ) = (7.782 – 0.00996P + 0.001932P 2 – 3.76 × 10 –5 P 3 ) + (–13 763 + 592P – 3.955P 2 )/T + (–1.05 × 10 6 – 4622P)/T 2 , which should be valid from room pressure up to 25 GPa and 773–2500 K, with an estimated uncertainty of 0.2 log units. Our calculated fO 2 is shown to be up to 1 log unit lower than estimates that use previous expressions or ignore EoS terms.

Armstrong, Katherine↗
DOE OSTI2017

SECARB Commercial Scale CO 2 Injection and Optimization of Storage Capacity in the Southeastern United States

The Commercial Scale Project is a US DOE/NETL funded initiative aimed at enhancing the knowledge-base and industry’s ability to geologically store vast quantities of anthropogenic carbon. In support of this goal, a large-scale, stacked reservoir geologic model was developed for Gulf Coast sediments centered on the Citronelle Dome in southwest Alabama, the site of the SECARB Phase III Anthropogenic Test. Characterization of regional geology to construct the model consists of an assessment of the entire stratigraphic continuum at Citronelle Dome, from surface to the depth of the Donovan oil-bearing formation. This project utilizes all available geologic data available, which includes: modern geophysical well logs from three new wells drilled for SECARB’s Anthropogenic Test; vintage logs from the Citronelle oilfield wells; porosity and permeability data from whole core and sidewall cores obtained from the injection and observation wells drilled for the Anthropogenic Test; core data obtained from the SECARB Phase II saline aquifer injection test; regional core data for relevant formations from the Geological Survey of Alabama archives. Cross sections, isopach maps, and structure maps were developed to validate the geometry and architecture of the Citronelle Dome for building the model, and assuring that no major structural defects exist in the area. A synthetic neural network approach was used to predict porosity using the available SP and resistivity log data for the storage reservoir formations. These data are validated and applied to extrapolate porosity data over the study area wells, and to interpolate permeability amongst these data points. Geostatistical assessments were conducted over the study area. In addition to geologic characterization of the region, a suite of core analyses was conducted to construct a depositional model and constrain caprock integrity. Petrographic assessment of core was conducted by OSU and analyzed to build a depositional framework for the region and provide modern day analogues. Stability of the caprock over several test parameters was conducted by UAB to yield comprehensive measurements on long term stability of caprocks. The detailed geologic model of the full earth volume from surface thru the Donovan oil reservoir is incorporated into a state-of-the-art reservoir simulation conducted by the University of Alabama at Birmingham (UAB) to explore optimization of CO 2 injection and storage under different characterizations of reservoir flow properties. The application of a scaled up geologic modeling and reservoir simulation provides a proof of concept for the large scale volumetric modeling of CO 2 injection and storage the subsurface.

58 GEOSCIENCES↗
DOE OSTI2017

Small Scale Field Test Demonstrating CO 2 Sequestration In Arbuckle Saline Aquifer And By CO 2 -Eor At Wellington Field, Sumner County, Kansas

The objectives of this project are to understand the processes that occur when a maximum of 70,000 metric tonnes of CO 2 are injected into two different formations to evaluate the response in different lithofacies and depositional environments. The evaluation will be accomplished through the use of both in situ and indirect MVA (monitoring, verification, and accounting) technologies. The project will optimize for carbon storage accounting for 99% of the CO 2 using lab and field testing and comprehensive characterization and modeling techniques. Site characterization and CO 2 injection should demonstrate state-of-the-art MVA tools and techniques to monitor and visualize the injected CO 2 plume and to refine geomodels developed using nearly continuous core, exhaustive wireline logs, and well tests and a multi-component 3-D seismic survey. Reservoir simulation studies will map the injected CO 2 plume and estimate tonnage of CO 2 stored in solution, as residual gas, and by mineralization and integrate MVA results and reservoir models shall be used to evaluate CO 2 leakage. A rapid-response mitigation plan was developed to minimize CO 2 leakage and provide a comprehensive risk management strategy. The CO 2 was intended to be supplied from a reliable facility and have an adequate delivery and quality of CO 2 . However, several unforeseen circumstances complicated this plan: (1) the initially negotiated CO 2 supply facility went offline and contracts associated with CO 2 supply had to be renegotiated, (2) a UIC Class VI permit proved to be difficult to obtain due to the experimental nature of the project. Both subjects are detailed in separate deliverables attached to this report. The CO 2 enhanced oil recovery (EOR) and geologic storage in Mississippian carbonate reservoir was sucessully deployed. Approximately 20,000 metric tons of CO 2 was injected in the upper part of the Mississippian reservoir to verify CO 2 EOR viability in carbonate reservoirs and evaluate a potential of transitioning to geologic CO 2 storage through EOR. A total of 1,101 truckloads, 19,803 metric tons—an average of 120 tonnes per day—were delivered over the course of injection that lasted from January 9 to June 21, 2016. After cessation of CO 2 injection, the KGS 2-32 well was converted to water injector and continues to operate. CO 2 EOR progression in the field was monitored weekly with fluid level, temperature, and production recording and formation fluid composition sampling. It is important to note that normally, CO 2 EOR pilots are less efficient than commercial operations due to lack of directional and precise well control, lack of surface facilities for CO 2 recycling, and other factors. As a result of this pilot CO 2 injection, the observed incremental average oil production increase was ~68% with only ~18% of injected CO 2 produced back. Decline curve analysis forecasts of additional cumulative oil produced were 32.44M STB to the end of 2027. Wellington Mississippian pilot efficiency by the end of forecast calculations is 11 MCF per barrel of produced oil. Using 32M STB oil production and $1,964,063 cost of CO 2 , CO 2 EOR cost per barrel of oil production is ~$60. Simple but robust monitoring technologies proved to be very efficient in detecting and locating CO 2 . High CO 2 reservoir retentions with low yields within an actively producing field could help to estimate real-world risks of CO 2 geological storage for future projects. The Wellington Field CO 2 EOR was executed in a controlled environment with high efficiency. This case study proves that CO 2 EOR could be successfully applied in Kansas carbonate reservoirs if CO 2 sources and associated infrastructure are available. Recent developments in unconventional resources development in Mid-Continent USA and associated large volume disposal of backflow water and the resulting seismic activity have brought more focus and attention to the Arbuckle Group in southern Kansas. Despite the commercial interest, limited essential information about reservoir properties and structural elements has impeded the management and regulation of disposal, an issue brought to the forefront by recent seismicity in and near areas of large volumes and rates of brine disposal. The Kansas Geological Survey (KGS) collected, compiled, and analyzed available data, including well logs, core data, step rate tests, drill stem tests, 2-D and 3-D seismic data, water level measurements, and others types of data. Several exploratory wells were drilled and core was collected and modern suites of logs were analyzed. Reservoir properties were populated into several site-specific geological models. The geological models illustrate the highly heterogeneous nature of the Arbuckle Group. Vertical and horizontal variability results in several distinct hydro-stratigraphic units that are the result of both depositional and diagenetic processes. During the course of this project, it has been demonstrated that advanced seismic interpretation methods can be used successfully for characterization of the Mississippian reservoir and Arbuckle saline aquifer. Analysis of post-stack 3-D seismic data at the Mississippian reservoir showed the response of a gradational velocity transition. Pre-stack gather analysis showed that porosity zones of the Mississippian and Arbuckle reservoirs exhibit characteristic amplitude versus offset (AVO) response. Simultaneous AVO inversion estimated P- and S-impedances. The 3-D survey gather azimuthal anisotropy analysis (AVAZ) provided information about the fault and fracture network and showed good agreement to the regional stress field and well data. Mississippian reservoir porosity and fracture predictions agreed well with the observed mobility of injected CO 2 in KGS well 2-32. Fluid substitution modeling predicted acoustic impedance reduction in the Mississippian carbonate reservoir introduced by the presence of CO 2 . Seismicity in the United States midcontinent has increased by orders of magnitude over the past decade. Spatiotemporal correlations of seismicity to wastewater injection operations have suggested that injection-related pore fluid pressure increases are inducing the earthquakes. In this investigation, we examine earthquake occurrence in southern Kansas and northern Oklahoma and its relation to the change in pore pressure. The main source of data comes from the Wellington Array in the Wellington oil field, in Sumner County, Kansas, which has monitored for earthquakes in central Sumner County, Kansas, since early 2015. The seismometer array was established to monitor CO 2 injection operations at Wellington Field. Although no seismicity was detected in association with the spring 2016 Mississippian CO 2 injection, the array has recorded more than 2,500 earthquakes in the region and is providing valuable understanding to induced seismicity. A catalog of earthquakes was built from this data and was analyzed for spatial and temporal changes, stress information, and anisotropy information. The region of seismic concern has been shown to be expanding through use of the Wellington earthquake catalog, which has revealed a northward progression of earthquake activity reaching the metropolitan area of Wichita. The stress orientation was also calculated from this earthquake catalog through focal mechanism inversion. The calculated stress orientation was confirmed through comparison to other stress measurements from well data and previous earthquake studies in the region. With this knowledge of the stress orientation, the anisotropy in the basement could be understood. This allowed for the anisotropy measurements to be correlated to pore pressure increases. The increase in pore pressure was monitored through time-lapse shear-wave anisotropy analysis. Since the onset of the observation period in 2010, the orientation of the fast shear wave has rotated 90°, indicating a change associated with critical pore pressure build up. The time delay between fast and slow shear wave arrivals has increased, indicating a corresponding increase in anisotropy induced by pore pressure rise. In-situ near-basement fluid pressure measurements corroborate the continuous pore pressure increase revealed by the shear-wave anisotropy analysis over the earthquake monitoring period. This research is the first to identify a change in pore fluid pressure in the basement using seismological data and it was recently published in the AAAS journal Science Advances (Nolte et al., 2017). The shear-wave splitting analysis is a novel application of the technique, which can be used in other regions to identify an increase in pore pressure. This increasing pore fluid pressure has become more regionally extensive as earthquakes are occurring in southern Kansas, where they previously were absent. These monitoring techniques and analyses provide new insight into mitigating induced seismicity’s impact on society.

02 PETROLEUM↗
DOE OSTI2019

Mid-Atlantic U.S. Offshore Carbon Storage Resource Assessment Project (Final Technical Report)

The greatest potential for carbon storage in the northeastern United States lies in the offshore geologic formations comprising the outer continental shelf (Monteverde et al., 2011). Offshore storage can be linked to large point-sources of carbon dioxide (CO 2 ) while avoiding many of the logistical difficulties and potential risks encountered when siting onshore projects, especially in densely populated areas of the East Coast. Offshore CO 2 storage assessment and research in the United States is underway to address uncertainty in potential storage resources, particularly in the mid- and north-Atlantic offshore area. Given the current knowledge base and access to publicly available data, the objectives of the Mid-Atlantic U.S. Offshore Carbon Storage Resource Assessment Project (MAOCSRAP) were fourfold: (1) complete a systematic carbon storage resource assessment of the mid-Atlantic offshore coastal region from the Georges Bank Basin (GBB) through the Long Island Platform (LIP) to the southern Baltimore Canyon Trough (BCT); (2) define key input parameters to reduce uncertainty for offshore storage resource and efficiency estimates; (3) perform a preliminary assessment of risk factors, uncertainties, and data gaps; and (4) engage industry and regulatory stakeholders through development of a Road Map to assist future project planning and implementation. The MAOCSRAP project has compiled, inventoried, and assimilated various publicly available data sets to provide a strong technical basis on which future carbon storage studies and applications can be built. The knowledge infrastructure necessary to support the development of full-scale offshore carbon storage include: Defining the geologic characteristics of candidate storage sites (Chapter 2). Using existing seismic data to better define the continuity of the storage zones and seals (Chapter 3). Cataloguing the hydrogeologic properties of mid-Atlantic offshore storage sites (Chapter 4). Calculating prospective CO 2 storage resources using net effective pore volumes and fluid displacement properties specific to offshore lithologies (Chapter 5). Examining risk factors related to offshore storage (Chapter 6). Communicating with industry and other stakeholders about the future prospects for offshore storage (Chapter 7). Ensuring technology transfer to industry and other stakeholders (Chapter 8). Led by Battelle, this project was conducted by public and private entities with expertise in offshore geology and resources for the study region, including state geological surveys of Delaware, Maryland, and Pennsylvania; United States Geological Survey-Woods Hole Coastal and Marine Science Center and Haifa University; Rutgers University; Harvard University; and Lamont-Doherty Earth Observatory (LDEO) at Columbia University. The storage resource assessment was completed for a broad region offshore of the U.S. East Coast, from Massachusetts to Virginia. The results include high-level storage resource estimates for areas not previously characterized and improved storage resource estimates for geographically expansive portions of offshore geologic units. Results of the project provide a foundation for developing CO 2 storage along offshore areas of the mid-Atlantic United States. A combination of information was used to provide confidence in project results. Key findings of the major project tasks are summarized as follows. Subsurface Data Analysis: Legacy seismic, well-log, core, and biostratigraphic data were digitized, reprocessed, and analyzed using modern techniques, augmenting previous characterization efforts. The sequence stratigraphic framework for the Cretaceous and Jurassic strata was developed by correlating well-log depositional sequences to seismic reflectors identified in the BCT, LIP, and GBB. The integration of well-log stratigraphy with seismic stratigraphy and petrophysical analysis helped to identify and map potential storage zones within sequences of the Logan Canyon (LC), Missisauga, and Mohawk formations, as well as the regionally extensive caprock comprised of the Dawson Canyon shale. Approximately 4,000 kilometers (km) of seismic data was reprocessed to create a structural framework that was used to tie together all available well data, constrain sequence stratigraphic interpretations, and make regional assessments of offshore carbon storage resources. Prospective storage resource estimates were 150 to 1136 megatons (Mt) for the combined storage zones. This suggests mid-Atlantic U.S. offshore formations can store decades of CO 2 from industrial sources in the region. Preliminary local reservoir simulations suggest that injection rates of 1 Mt CO 2 /year may be sustained for 30 years in single injection wells in the BCT. Advanced geologic modeling and new data acquisition are needed to address data gaps and advance carbon capture and storage (CCS) in key offshore areas selected for further investigation. Offshore Risk Factors: Offshore geologic risk factors include soft sediment deformation, unit continuity, sedimentological and structural features, seismicity and hydrates. CO 2 storage risks include inadequate seals, migration/leakage, and chemical interactions leading to decreased storage. Sensitive habitats, environmental impacts, disturbance to seafloor, and other risks need to be identified in advance of project activities and integrated into detailed mitigation plans for all project phases. Stakeholder Engagement: Input and participation from government, industry, and environmental groups provided input into the Road Map and address next steps needed for project deployment. Early engagement and ongoing communication, as well as policy framework development, is key to large-scale deployment of CCS. Next Steps: This project represents an important first step by completing a high-level CO 2 storage resource assessment and building the knowledge infrastructure necessary to improve quantitative storage resource estimates. The data sets that have been curated under this project provide an opportunity to conduct research and development (R&D) needed to address data gaps and reduce risk and uncertainty. Offshore characterization and validation strategies that are systematically designed to provide data and infrastructure that can be upscaled to meet commercial requirements could be developed off of this foundation.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2020

In situ and laboratory evidence for high electrical anisotropy in marine gas hydrate

Marine gas hydrate is electrically resistive compared to the surrounding hydrate-free formation, a property that is exploited in borehole logging and also by marine controlled-source electromagnetic (CSEM) surveys. Well logging shows that hydrate can be electrically anisotropic, with vertical resistivities up to 10 times higher than horizontal resistivities, although the use of logging tools capable of measuring anisotropy is the exception rather than the rule. Marine CSEM methods are well known to be sensitive to anisotropy in seafloor sediments, and, for example, we have observed anisotropies of more than an order of magnitude using CSEM to map offshore relic permafrost.

02 PETROLEUM↗
DOE OSTI2022

Tuscaloosa Marine Shale Laboratory

The Tuscaloosa Marine Shale (TMS) in Louisiana and Mississippi is an Upper Cretaceous source rock formation sandwiched between the sands of the upper and lower Tuscaloosa sections. The TMS is believed to be the source rock for underlying prolific Tuscaloosa sand formation. The TMS has an unproven estimate of 7,000,000,000 bbls of recoverable oil while its current total average production is about 3,000 bbls of oil per day in 2017. In 2013 and 2014, more than 80 wells were drilled horizontally into the TMS that were fractured using multi-stage fracturing technology. The results from this have been mixed, but recent production for several wells show an appealing initial oil production rate of more than 1000 bbl/day. The preliminary core analysis by industry partners and a few literature studies shows that the TMS is one of the most clay-rich and sensitive shales to water. Due to these and other technical problems, there is high risk for the economic development of TMS compared to other shale plays. The experiences of major industrial players in the TMS show the necessity of open and collaborative efforts to better understand the critical gaps in the development of this challenging and potentially highly economic shale play to enable more cost-efficient and environmentally-sound recovery from this unconventional liquid-rich shale play. The overall objective of this project is to form a consortium of science and industry partners to address the following six major objectives using scientific and technical approaches: 1. To improve wellbore integrity by better understanding the sources of the wellbore instability issues, proposing innovative mud and cement design for the TMS. 2. To improve formation evaluation using laboratory techniques for the evaluation of mineralogical composition, organic content, and produced-water chemistry as well as well log and geophysical analysis. 3. To determine the role of geologic discontinuities on fracture growth and shale creep behavior using digital image correlation technique. 4. To investigate the application of stable CO 2 foam and super-hydrophobic proppants for improved reservoir stimulation. 5. To better understand the nature of water/hydrocarbon/CO 2 flow in clay and organic-rich formation and the role of water/fluid interaction on recovery. 6. To prepare better socio-economic environment for TMS development by community engagement. Subsequently, the TMS virtual laboratory conducted testing and analysis of various properties of rock and formation fluids from the TMS, including but not limited to the following: Analyzing reports and logs to better understand the source of wellbore instability in TMS wells; Experiments to design a customized cement based on TMS requirements; Experiments to obtain the mineralogical and geochemical composition of TMS samples; Seismic analysis of TMS geophysical data to better predict total organic carbon (TOC) content and brittleness in TMS; Well log analysis to better estimate the TOC and geo-mechanical properties of TMS; Experiments on formation water to understand the chemistry of produced water; Experiments to determine the role of lamination and natural fractures on fracture propagation or rock deformation using digital image correlation technique in in-direct tensile tests, semi-circular bend test and creep tests Experiments to determine the stability and rheological properties of nanoparticle-stabilized CO 2 foam in TMS rock samples; Experiments to determine fluid dynamics in un-propped TMS fractures and the role of nano-coating of proppants on fluid dynamics in fractures with proppants; Micro-fluidics experiments to enhance the understanding of fluid dynamics in tight liquidrich pores with high clay content; Socio-economic studies to better engage communities in TMS development.

58 GEOSCIENCES↗
DOE OSTI2022

Wabash CarbonSAFE. Geologic Analysis of the Potosi Dolomite Reservoir Interval and Potential Confining Units, at the Wabash #1 Well, Vigo Co., Indiana, Subtask 7.2, Technical Report

The Wabash CarbonSAFE project drilled the Wabash #1 stratigraphic test well (ID# 168045) at the Wabash Valley Resources (WVR) IGCC facility in Vigo County, Indiana, to evaluate the feasibility of commercial-scale CO 2 storage near the site. The Mt. Simon Sandstone was the initial target for storage evaluation and was found to have generally poor reservoir qualities in the Wabash #1 well. Focus was then placed on the secondary target, the Potosi Dolomite. For storage and containment evaluation an extensive suite of geophysical logs was collected throughout the borehole, and whole cores were obtained from the Mt. Simon Sandstone, Eau Claire Formation, and Maquoketa Group. The Maquoketa Group is identified as a regionally extensive secondary sealing interval above the Mt. Simon Sandstone and as a primary seal for the Potosi Dolomite storage complex. Due to concerns over potential lost circulation zones in the Potosi Dolomite during drilling, no core samples or Formation Micro Imager (FMI) logs were acquired in the Potosi Dolomite. However, cased hole well testing was performed over a 20 feet (6 m) interval within the Potosi Dolomite during which a fluid swab sample was obtained. Approximately 35 miles of 2D seismic information was acquired in the project area near Wabash #1 to aid in reservoir and caprock characterization. The highly porous and permeable zones in the Potosi Dolomite and lower Oneota Dolomite observed within the Wabash #1 well are, for the purpose of this report, considered the same reservoir package and referred to as the Potosi reservoir interval. This report presents a regional depositional context for the Potosi reservoir interval, the overlying confining zones in the Knox Group, and Maquoketa Group primary seal. This context is developed through detailed information obtained from the Wabash #1 well such as lithologic data collected from cuttings and core, geophysical logging, geomechanical analysis of core samples, and well testing and fluid sampling within the Potosi reservoir interval. Information discussed in this report was compiled as a basis for the development of static and dynamic geologic models of the Potosi Dolomite storage complex. The static and dynamic modeling of CO 2 injection in the Potosi reservoir interval are discussed in a separate report (Dessenberger et al., 2022) under the Wabash CarbonSAFE project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗