Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “core formation”

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

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗
DOE OSTI2025

LABQ3: Bayesian method for quantification of mineral compositions and nano-scale elemental mapping of 3D synchrotron XCT data

Quantitative analysis of mineral compositions is essential in understanding geochemical, mineralogical and environmental processes. Fine-resolution 3D imaging is widely done using synchrotron X-ray computed tomography (XCT), but existing analyses are limited to visualization and segmentation. This paper presents a new method, Linear Attenuation Bayesian Quantitative 3D-mapper (LABQ3), based on the linearity of X-ray attenuation with respect to elemental concentrations. To address the random variability in attenuation measurements, LABQ3 employs Bayesian decision theory to minimize classification error, using reference attenuation distributions from scans of pure mineral standards. To demonstrate LABQ3 and test its performance, we studied precipitated carbonate samples. XCT scans were done at multiple energies using the transmission X-ray microscope (TXM) at beamline 32-ID-C of the Advanced Photon Source at Argonne National Laboratory. The reconstructed 3D images have a voxel size of 20 nm. Analyses revealed rich nano-scale compositional heterogeneity within individual particles. A mixture of calcium and cadmium produced an overall stoichiometric composition of (Ca 0.78 ,Cd 0.22 )CO 3 , with some voxels containing nearly pure CdCO 3 . The addition of zinc led to an overall stoichiometric composition of 33% Ca, 28% Cd, 39% Zn, with a nearly pure CaCO 3 core and compositional zonation through the rim. These compositional gradients are related to temporal sequences of carbonate mineral formation where Cd precipitated at the beginning in (Ca,Cd)CO 3 , while Cd and Zn precipitated at the end in (Ca, Cd,Zn)CO 3 . Results differ from bulk analyses using Inductively Coupled Plasma-Mass Spectrometry (ICP-MS), showing that LABQ3 provides particle-specific insights. LABQ3 distinguishes itself by quantifying chemical compositions along a continuum, making it different from XCT analyses based on segmentation. LABQ3 allows simultaneous acquisition of morphology and chemical composition in 3D, facilitating the interpretation of chemical gradients of trace elements, quantification of solid solution compositions, inferences about temporal sequences of mineral precipitation, and addressing other concerns about solid-phase chemistry.

58 GEOSCIENCES↗
DOE OSTI2020

P SI4 1.4: Open-source software for high-throughput quantum chemistry

PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree–Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4 's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. We announce a rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Experiment of Sea Breeze Convection, Aerosols, Precipitation, and Environment (ESCAPE) C-band Radar Deployment Field Campaign Report

Convective clouds play a critical role in the Earth’s climate system, acting as sinks of total water in the atmospheric column through precipitation, thereby contributing to the atmospheric energy balance. They serve as a primary mechanism for the transfer of heat, moisture, and momentum through the troposphere, significantly impacting the large-scale atmospheric circulation and local environment, and affecting the probability of subsequent cloud formation (e.g., Hartmann et al. 1984, Su et al. 2014, Sherwood et al. 2014). Convective characteristics including convective core size (reflectivity and updraft), cloud lifetime, precipitation intensity, total precipitation amount, precipitation efficiency, and lightning flash rates can be strongly related to environmental factors (e.g., wind shear, humidity). Since cumulus clouds evolve rapidly, their microphysical and dynamical properties and life cycles are challenging to resolve in models, and even in observations (e.g., Fridlind et al. 2017, Ladino et al. 2017).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evolution of the Thermodynamic Properties of Clusters of Galaxies out to Redshift of 1.8

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. Here, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant (z > 1.2) SPT-selected clusters, and compare them with those of the nearby clusters (z < 0.1) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1983

The nature of R136a, the superluminous central object of the 30 Doradus nebula

UV data show R136a to be an extremely luminous object having spectral characteristics similar to that of O3 stars. It is suggested that the bright component of R136a is either a single, supermassive hydrogen-burning star, or a compact multiple system of supermassive stars. In the former case, the observational data imply a mass exceeding that of the most massive stars known by a factor of 20 and raise important theoretical problems concerning stellar formation and stability. Formation may occur through ordinary collapse in a region with peculiar dust properties, or through coalescence in the core of a compact cluster.

Cassinelli, J. P.↗
NASA NTRS1984

Earth

The following aspects of the planet Earth are discussed: plate tectonics, the interior of the planet, the formation of the Earth, and the evolution of the atmosphere and hydrosphere. The Earth's crust, mantle, and core are examined along with the bulk composition of the planet.

Michael H Carr↗
DOE OSTI2020

Evolution of the Thermodynamic Properties of Clusters of Galaxies out to Redshift of 1.8

The thermodynamic properties of the hot plasma in galaxy clusters retains information on the processes leading to the formation and evolution of the gas in their deep, dark matter potential wells. These processes are dictated not only by gravity but also by gas physics, e.g. AGN feedback and turbulence. In this work, we study the thermodynamic properties, e.g. density, temperature, pressure, and entropy, of the most massive and the most distant ($z > 1.2$) SPT-selected clusters, and compare them with those of the nearby clusters ($z<0.1$) to constrain their evolution as a function of time and radius. We find that thermodynamic properties in the outskirts of high redshift clusters are remarkably similar to the low redshift clusters, and their evolution follows the prediction of the self-similar model. Their intrinsic scatter is larger, indicating that the physical properties that lead to the formation and virialization of cluster outskirts show evolving variance. On the other hand, thermodynamic properties in the cluster cores deviates significantly from self-similarity indicating that the processes that regulate the core are already in place in these very high redshift clusters. This result is supported by the unevolving physical scatter of all thermodynamic quantities in cluster cores.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Galaxy Size and Rotation Curve Diversity in ΛCDM with Baryons

The observed rotation curves of dwarf galaxies exhibit significant diversity at fixed halo mass, challenging galaxy formation within the cold dark matter (CDM) model. Previous cosmological galaxy formation simulations with baryonic physics fail to reproduce the full diversity of rotation curves, suggesting that there is a flaw in baryonic feedback models, observational bias, or that an alternative to CDM must be invoked. In this work, we use the Marvelous Massive Dwarf zoom-in simulations, a suite of high-resolution dwarf simulations with M 200 ∼ 10 10 –10 11 M ⊙ and M * ∼ 10 7 –10 9 M ⊙ , designed to target the mass range where the galaxy rotation curve diversity is maximized, i.e., between and 100 km s −1 . We add to this a set of low-mass galaxies from the Marvel Dwarf Zoom Volumes to extend the galaxy mass range to lower values. Our fiducial star formation and feedback models produce simulated dwarfs with a broader range of rotation curve shapes, similar to observations. These simulations both create dark matter cores via baryonic feedback, reproducing the slower-rising rotation curves, while also allowing for compact galaxies and steeply rising rotation curves. Our simulated dwarfs also reproduce the observed size–M * relation, including scatter, producing both extended and compact dwarfs for the first time in simulated field dwarfs. However, the slowly rising and high baryon mass fraction, as well as the steeply rising and low baryon mass fraction, remain missing. We explore star formation and feedback models and conclude that previous simulations may have had feedback that was too strong to produce compact dwarfs.

Cruz, Akaxia [Flatiron Institute, New York, NY (Un↗
DOE OSTI2021

On the Generation of Weddell Sea Polynyas in a High-Resolution Earth System Model

Larger Weddell Sea polynyas (WSPs), differentiated in this study from the smaller Maud Rise Polynyas (MRPs), forming to the east of the prime meridian in the proximity of the Maud Rise seamount, have last been observed in the 1970s. We investigate WSPs that grow realistically out of MRPs in a high-resolution (HR) preindustrial simulation with the Energy Exascale Earth System Model version 0.1. The formation of MRPs requires HR to simulate the detailed flow around Maud Rise, while the realistic formation of WSPs requires a model to produce MRPs. Furthermore, WSPs tend to follow periods of a prolonged build-up of a heat reservoir at depth and weakly negative wind-stress curl in association with the core of the southern hemisphere westerlies at an anomalously northern position. While this scenario also leads to drier conditions over the central Weddell Sea, which some literature claims to be a necessary condition for the formation of WSPs, our model results indicate that open-ocean polynyas do not occur during periods of weakly negative wind-stress curl despite drier atmospheric conditions. Our study supports the hypothesis noted in earlier studies that a shift from a weakly negative to a strongly negative wind-stress curl over the Weddell Sea is a prerequisite for WSPs to form, together with a large heat reservoir at depth. However, the ultimate trigger is a pronounced MRP; whose associated convection creates high surface salinity anomalies that propagate westward with the flow of the Weddell Gyre. If large enough, these anomalies trigger the formation of a WSP and a pulse of newly formed Antarctic Bottom Water.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Computed Tomography Scanning and Geophysical Measurements of the CarbonSAFE Seal Integrity Wells in the Illinois Basin

The U.S. Department of Energy’s (DOE) National Energy Technology Laboratory (NETL) researchers in Morgantown, West Virginia, utilized computed tomography (CT) facilities and the Multi-Sensor Core Logger (MSCL) to evaluate the integrity of three cores (Geophysical Monitor #2 (GM2), Verification Well #1 (VW #1), and Wabash #1 for the Illinois Basin Carbon Storage Assurance Facility Enterprise (CarbonSAFE) efforts. This work was carried out in collaboration with the Illinois State Geological Survey (ISGS) as part of their efforts to analyze cores from two field locations, including the Illinois Basin-Decatur Project (IBDP) and second stage Illinois Industrial Carbon Capture and Storage (IL-ICCS) sites at Decatur, Illinois, and the Wabash CarbonSAFE project at the Wabash Valley Resources Integrated Gasification Combined Cycle plant in Terre Haut, Indiana. The results of this study are presented in several formats and are available online on the Energy Data eXchange (EDX). The rock characterization was conducted using non-destructive techniques, allowing for future analysis of the cores. While the equipment used did not provide direct visualization of shale pore space, it allowed for detection of fractures and discontinuities. Low resolution CT imagery with the NETL medical CT scanner was performed on the entire core. Qualitative analysis of the medical CT images, coupled with x-ray fluorescence (XRF), P-wave, and magnetic susceptibility measurements from the MSCL were utilized to identify areas of interest for further study as well as fractured zones. The ability to quickly identify key areas for more detailed study with higher resolution will save time and resources in future studies. The combination of methods used provided a multi-scale analysis of this core and provided both a macro and micro description of the core that is relevant for many subsurface energy-related examinations traditionally performed at NETL.

58 GEOSCIENCES↗
NASA NTRS2003

Terrestrial Planet Formation in Binary Star Systems

Most stars reside in multiple star systems; however, virtually all models of planetary growth have assumed an isolated single star. Numerical simulations of the collapse of molecular cloud cores to form binary stars suggest that disks will form within such systems. Observations indirectly suggest disk material around one or both components within young binary star systems. If planets form at the right places within such circumstellar disks, they can remain in stable orbits within the binary star systems for eons. We are simulating the late stages of growth of terrestrial planets within binary star systems, using a new, ultrafast, symplectic integrator that we have developed for this purpose. We show that the late stages of terrestrial planet formation can indeed take place in a wide variety of binary systems and we have begun to delineate the range of parameter space for which this statement is true. Results of our initial simulations of planetary growth around each star in the alpha Centauri system and other 'wide' binary systems, as well as around both stars in very close binary systems, will be presented.

Lissauer, Jack J.↗
DOE OSTI2024

Application of performance portability solutions for GPUs and many-core CPUs to track reconstruction kernels

Next generation High-Energy Physics (HEP) experiments are presented with significant computational challenges, both in terms of data volume and processing power. Using compute accelerators, such as GPUs, is one of the promising ways to provide the necessary computational power to meet the challenge. The current programming models for compute accelerators often involve using architecture-specific programming languages promoted by the hardware vendors and hence limit the set of platforms that the code can run on. Developing software with platform restrictions is especially unfeasible for HEP communities as it takes significant effort to convert typical HEP algorithms into ones that are efficient for compute accelerators. Multiple performance portability solutions have recently emerged and provide an alternative path for using compute accelerators, which allow the code to be executed on hardware from different vendors. We apply several portability solutions, such as Kokkos, SYCL, C++17 std::execution::par and Alpaka, on two mini-apps extracted from the mkFit project: p2z and p2r. These apps include basic kernels for a Kalman filter track fit, such as propagation and update of track parameters, for detectors at a fixed z or fixed r position, respectively. The two mini-apps explore different memory layout formats. We report on the development experience with different portability solutions, as well as their performance on GPUs and many-core CPUs, measured as the throughput of the kernels from different GPU and CPU vendors such as NVIDIA, AMD and Intel.

Kwok, Ka Martin↗
DOE OSTI2024

Application of performance portability solutions for GPUs and many-core CPUs to track reconstruction kernels

Next generation High-Energy Physics (HEP) experiments are presented with significant computational challenges, both in terms of data volume and processing power. Using compute accelerators, such as GPUs, is one of the promising ways to provide the necessary computational power to meet the challenge. The current programming models for compute accelerators often involve using architecture-specific programming languages promoted by the hardware vendors and hence limit the set of platforms that the code can run on. Developing software with platform restrictions is especially unfeasible for HEP communities as it takes significant effort to convert typical HEP algorithms into ones that are efficient for compute accelerators. Multiple performance portability solutions have recently emerged and provide an alternative path for using compute accelerators, which allow the code to be executed on hardware from different vendors. We apply several portability solutions, such as Kokkos, SYCL, C++17 std::execution::par, Alpaka, and OpenMP/OpenACC, on two mini-apps extracted from the mkFit project: p2z and p2r. These apps include basic kernels for a Kalman filter track fit, such as propagation and update of track parameters, for detectors at a fixed z or fixed r position, respectively. The two mini-apps explore different memory layout formats. We report on the development experience with different portability solutions, as well as their performance on GPUs and many-core CPUs, measured as the throughput of the kernels from different GPU and CPU vendors such as NVIDIA, AMD and Intel.

Kwok, Ka Hei Martin↗
NASA NTRSDate not supplied

SURFACE EXOSURE TIMESCALES OF APOLLO CORE SAMPLE 73002 SPACE WEATHERED GRAINS

Introduction: Space weathering causes the surface regolith on airless bodies like the Moon to be morphologically, microstructurally, and chemically altered due to micrometeoroid bombardment and solar wind irradiation[1]. These processes produce a multitude of microstructural and chemical changes in individual soil grains that accumulate as grains are exposed on the surface over time. One microstructural feature produced by exposure to the solar wind is the development of ion-damaged rims on regolith grains via H+ and He+ ion irradiation. Also present are solar energetic particle (SEP) tracks, which are nanoscale lineations of ionization damage within grain interiors formed by high energy solar flare ions(primarily Fe group nuclei) which penetrate millimeters below the surface [2]. Recent work has confirmed that the thickness of solar wind-damaged amorphous rims on anorthite grains, nanocrystalline rims on olivine grains, and their respective SEP track densities are correlated with their surface exposure ages [3].Core sample 73002, recently released under the Apollo Next Generation Sample Analysis (ANGSA) Program, has provided an opportunity to study material collected from the light mantle formation during Apollo 17. The light mantle is thought to have been deposited via a landslide originating from the neighboring South Massif [4]. Spectral profiles and ferromagnetic resonance measurements of bulk soils sampled at cm-intervals from 73002 indicate the existence of mature regolith extending ~8 cm below the surface, suggesting the presence of an in-situ reworking zone in the core[5,6].Here we present the distribution of SEP track densities, solar wind damaged rim widths, and the corresponding surface exposure ages of grains residing in the proposed in-situ reworking zone of core 73002. Methods: Bulk samples of regolith from the first eight intervals and every following fourth interval down the core (dissection Pass 2) were delivered to Purdue University as <45 μm size fractions. The first eight0.5 cm intervals, representing the top 4cmof regolith, were individually dry sieved to a<20 μm size fraction, and grains were prepared by ultramicrotomy for analysis in the scanning transmission electron microscope (STEM). Bright field (BF) and dark field (DF) STEM images of solar flare tracks and solar wind damaged rims were acquired on a JEOL 2500SE TEM, equipped with a 60 mm2ultra-thin window silicon drift energy dispersive X-ray (EDX) spectrometer at NASA John-son Space Center. Grain compositions were determined by EDX compositional spectrum imaging. SEP track densities and amorphous rim thicknesses were measured on BF and DF images. Results: The fraction of grains with discernable space weathering features decreased with depth, with 100% in the first interval exhibiting space weathering features as opposed to ~60% in Interval 8.BF and DF STEM images show splash-melt and vapor deposited rims on outermost grain margins with embedded Fe-bearing nanoparticles ranging in size up to ~10 nm in diameter. Images also show solar-wind damaged rims below the vapor deposits, and SEP tracks present in grain interiors. The rim thicknesses and SEP track densities were determined using the methods of [3] for 54grains spanning the top 4 cm. ~80% of the grains analyzed were anorthite and ~20% were olivine. The track production rate from[3] was used to estimate the surface exposure times of the grains in this study. The majority of exposure times are in the 1 -5 MY range for all intervals. The lowest value was 3.8x 10¬5years in Interval 8and the highest value was1.1x 107years in Interval 6. Grains with the highest solar wind damage rim thicknesses and SEP track densities are found in Interval 6.The highest rim thickness and track density observed for anorthiteis126.5 nm and4.7 x 1011tracks/cm2, while for olivine they are 98.6 nm and 1.83 x 1011tracks/cm2, respectively. Discussion: Intervals from the top 4 cm of 73002 show similar SEP track density and solar wind damaged rim thickness distributions. Spectral and ferromagnetic resonance measurements of 73002 indicate an in-situ reworking zone up to~8 cm below the surface[5,6].The depth analyzed in this work is within the core’s uppermost in-situ reworking zone and may serve as an explanation to the homogenous distribution of exposure ages. The decreasing abundance of space weathered features with depth is consistent with regolith mixing models showing an exponential decay in regolith maturity with depth over time[7].A maximum surface exposure age determined via SEP track density is approximately 10million years. This timescale is consistent with the lower estimate of surface exposure ages of the light mantle determined in previous studies, which can range from 10s to over 100 Ma [4].It also falls within the core regolith reworking development timescale of ~17 Ma [5,6]. Additional analysis will be performed for a further seven intervals in 73002.

J A McFadden↗
DOE OSTI2026

A Review and Outlook on Experimental Advances and Innovations in Geological CO 2 Storage: Insights from Depleted Gas Reservoirs and Saline Aquifers

Geological storage of carbon dioxide (CO 2 ) in depleted gas reservoirs and deep saline aquifers is a key part of global decarbonization efforts. As carbon capture and storage advances toward commercial-scale deployment, the credibility and scalability of laboratory experiments are increasingly vital for guiding safe and effective field implementation. This review offers a comprehensive, cross-scale evaluation of experimental methodologies, including core flooding, high-pressure, high-temperature systems, microfluidic visualization, and emerging systems such as multilayer commingled/compartmentalized core flooding, 3D-printed micromodels, and AI-powered digital twins. These innovations are demonstrated to enhance representativeness, reproducibility, and real-time insight, thereby addressing the limitations of conventional workflows. A critical analysis of methodological gaps, such as inconsistent pressure–temperature conditions, oversimplified brine chemistry, and a lack of standardization, reveals experimental sources of scale translation errors and performance uncertainty. By comparing the unique challenges of depleted gas reservoirs (such as low water saturation and legacy well leakage) to those of saline aquifers (including pressure buildup and caprock integrity), this review identifies formation-specific priorities for experimental design. Novel contributions include a synthesis of best practices, integration strategies for model calibration, and recommendations for standardizing core handling, saturation procedures, and reporting protocols. Furthermore, this work serves as a guide for developing robust, field-relevant experimental strategies that can increase the deployment and regulatory acceptance of CO 2 storage technologies at scale.

58 GEOSCIENCES↗
DOE OSTI2020

Potential Formation Damage: An Integrated Reservoir Characterization Study of the Naturally Fractured Carbonate Middle Duperow Formation at the Kevin Dome, Montana

In this study, we integrate geologic and engineering data of a naturally fractured carbonate reservoir at the Kevin Dome, Montana. Well test data are correlated with core description, geochemical and lithology study to determine the flow behavior and communication within the injection test interval and to the surrounding area. Based on a dual-continuum geologic model, numerical brine injection simulations are carried out to validate the interpretation results from our well test analytical models and forecast the probability of CO 2 injection success using current reservoir properties. As a result, our well test analytical models as well as lithology/core description suggest that fluid flow may be mainly restricted to the injection interval and the assumption of radial (horizontal) flow may be appropriate. The well test models also indicate that there is potentially formation damage with a positive skin factor although prior to brine injection well tests, well stimulation through acid treatment was performed. Our numerical simulation results appear to confirm this formation damage by showing additional pressure buildup in the injection data during later test periods. To explain this, acid may have dissolved dolomite then dolomite or calcite may have been formed again further into the matrix/fracture system. Another possible explanation is mechanical clogging of the fractures due to acid dissolving dolomite and dislodging fine grains. Our work also predicts that if no additional well stimulation is performed, the project will have a lower probability of successfully injecting 1 million tons of CO 2 into the Middle Duperow formation over 4 years.

42 ENGINEERING↗