Search NASA⌕ Search

SEARCH · Search NASA

Results for “Geological Modeling”

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 559 records · Page 31

Ishtar deformed belts: Evidence for deformation from below?

The mountain belts of Ishtar Terra are unique on Venus. Models for their formation include mantle upwelling, mantle downwelling, and horizontal convergence. The present forms of these models are too simple to predict surface strain, topography, or gravity. More detailed models will require specific constraints as imposed by geologic relations. In order to develop specific constraints for geodynamic models, we examine the geology of Ishtar Terra as viewed in Magellan SAR imagery in an attempt to interpret regional surface strain patterns. In this paper, we present geologic and structural relations that leads us to postulate that Ishtar deformed belts result from shear forces within the mantle acting on the lithosphere, and not by horizontal forces from colliding plates. We propose that the surface strains result from differential strain and displacement of domains within the upper mantle, and that further analysis of Ishtar deformation may allow us to identify individual domains within the mantle, and to constrain displacement trajectories between domains.

Hansen, V. L.↗

Thermal inertia imaging - A new geologic mapping tool

A thermal model of the earth's surface has been developed and used to determine the thermal inertia of a test site in the Mojave Desert, California. The model, which includes meteorological heating terms as well as radiation and conduction heating terms, is used with remotely sensed surface temperature and reflectance data to determine the thermal inertia of the surface materials at the test site. The thermal inertia is displayed in image form, and can aid in the differentiation of the various lithologic materials in the test site. Since this thermal property is representative of the upper several cm of the surface, it complements visible and reflected near-IR image data.

Kahle, A. B.↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

Building an EPA Class VI Permit Application

Summary To accelerate the commercialization of carbon capture and storage (CCS), the US Department of Energy (US DOE) is building on decades of characterization efforts and pilot-scale projects through their CarbonSAFE program. Administered through their National Energy Technology Laboratory, this program seeks to bring fully integrated projects to the sector that can store more than 50 million tonnes of CO2 over a 30-year period. The program, which was enacted before the enhancement of Internal Revenue Code Section 45Q, is in the capture assessment, characterization, and permitting phase. The objectives of this paper are to discuss (a) the injection permitting requirements of the CarbonSAFE projects; (b) information gathering in support of the permit; (c) the timelines of field development and permit-related activities; (d) the major technical components of the field development plan; and (e) early feedback from the regulators toward acceptance of the permit. In Mississippi, more than 30,000 acres have been characterized by six deep characterization wells, a deep groundwater well, and 92 line miles of 2D seismic as part of the CarbonSAFE Project ECO2S. During the acquisition of seismic data, all receiver lines were live, which resulted in the generation of a pseudo-3D seismic design. The incorporation of a 3D seismic survey was not included as part of this project due to logistical difficulties presented by the undulating, wooded surface terrain. A suite of openhole geophysical logs was taken from each well, allowing for a detailed interpretation of prospective storage reservoirs and confining intervals to complement the analysis carried out on the 290 ft of a whole core that was cut through the prospective confining zone and storage reservoir. The detailed geologic and reservoir data were assembled and entered into a 3D model to assess the injection capacity and the area of review (AoR). This information fed into the detailed corrective action, monitoring, testing, and postinjection site care (PISC) modeling. The results have been exceptional. The geologic assessment has revealed three primary storage targets, ranging in depth from 3,500 ft to 6,000 ft. These storage reservoirs net 1,300 ft of sandstone, with mean porosity and permeability of 29% and 3.6 darcies, respectively. Together, these reservoirs have storage capacities that may exceed 20 million tonnes per square mile, making this a gigatonne prospect. Forward modeling of the project resulted in an AoR of 16 sq miles, injecting about 8000 t/d, for 30 years, via two deep injection wells. The excellent confining characteristics of the caprock, relatively simple geologic structure, and lack of historical well drilling activity in this area provide excellent containment of the injected CO2. Based on this work, the project has proposed 20 years of PISC. To date, only two US CO2 injection permits have been granted. These projects relied on a singular capture point feeding a singular sequestration point (source to sink), and considerations have not been made to garner CO2 emissions from other industrial sources. The Kemper County Storage Complex is a first-of-its-kind storage hub concept that looks to develop an area capable of storing significant quantities of CO2 from the region. Also, this work will show how characterization efforts, geological and numerical modeling efforts, and plan development were constructed in support of permit and incentives acceptance.

Energy & Fuels↗

Subsurface Characterization of Hydraulic Fracture Test Site-2 (HFTS-2), Delaware Basin

Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp Formation of the Permian (Delaware) Basin. This paper focuses on integration, advanced geological characterization, and 3D subsurface modeling of the comprehensive HFTS-2 dataset. The study showcases a multidisciplinary reservoir characterization approach that incorporates geology, petrophysics, geochemistry, geomechanics, microseismic, and subsurface engineering analysis. Subsurface characterization of organic-rich mudstone formations requires understanding complex hydraulic fracture network growth in relation to inherent lithology, geomechanical properties, and interaction with pre-existing natural fractures. This paper presents a characterization workflow incorporating pre- and post-stimulation subsurface data, unique to the HFTS-2 dataset. The study integrated: (1) rock properties from logs, cores, and thin sections; (2) natural and hydraulic fracture descriptions from cores and image logs; (3) local and regional stresses; (4) geomechanics; (5) microseismic; (6) fiber optic (FO) and bottomhole pressure gauge (BHPG) response; and (7) produced fluids analysis. During a stimulation treatment, creation of the stimulated rock volume (SRV) is influenced by several subsurface factors. Key contributing factors include structural context, stress conditions, lithology, facies architecture, pre-existing natural fractures, and geomechanical properties. The HFTS-2 subsurface data integration indicates that the SRV is comprised of a complex juxtaposition of hydraulic fracture swarms, as evidenced by image logs analysis, core description, and microseismic monitoring. The HFTS-2 microseismic event density was used to generate 3D heat maps that serve as a representative SRV footprint, corroborated by secondary datasets. These maps were further integrated with petrophysical and geomechanical characteristics, as well as responses from FO and BHPG, to estimate the lateral and vertical dimensions of the effective fractures. The geological characterization for the HFTS-2 dataset combined with 3D modeling for petrophysical and geomechanical properties provides a strong foundation for subsurface simulation and optimization studies. Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S005R001/2477501/urtec-2021-5243-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5243 The workflow improved our understanding of HFTS-2 hydraulic fracture propagation and characteristics in relation to offset pressure depletion and interaction with pre-existing natural fractures. Analysis showed that fracture geometry varies by stage and by well, and a complex fracture network is generated with varying fracture density. The multidisciplinary workflow presented herein for integration and characterization serves as a foundation to evaluate completion efficiency and estimate areal and vertical stimulation and depletion extent for the project. Furthermore, the workflow and learnings can also be transferred to other unconventional plays.

58 GEOSCIENCES↗

Ion-Exchange Interdiffusion Model with Potential Application to Long-Term Nuclear Waste Glass Performance

Ion exchange is an integral mechanism influencing the corrosion of glasses. Due to the formation of alteration layers in aqueous conditions, it is difficult to conclusively deconvolute the process of ion exchange from other processes, principally dissolution of the glass matrix. Therefore, we have developed a method to isolate alkali diffusion that involves contacting glass coupons with a solution of 6 LiCl dissolved in functionally inert dimethyl sulfoxide. Here, we employ the method at temperatures ranging from 25 to 150 °C with various glass formulations. Glass compositions include simulant nuclear waste glasses, such as SON68 and the international simple glass (ISG), glasses in which the nature of the alkali element was varied, and glasses that contained more than one alkali element. An interdiffusion model based on Fick’s second law was developed and applied to all experiments to extract diffusion coefficients. The model expands established models of interdiffusion to the case where multiple types of alkali sites are present in the glass. Activation energies for alkali ion exchange were calculated and the results are in agreement with those obtained in glass strengthening experiments but are nearly five times higher than values reported for diffusion-controlled processes in nuclear waste glass corrosion experiments. A discussion of the root causes for this apparent discrepancy is provided. The interdiffusion model derived from laboratory experiments is expected to be useful for modeling glass corrosion in a geological repository when the silicon concentration is high.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Abstracts of the Annual Meeting of Planetary Geologic Mappers, Flagstaff, AZ, 2008

Topics discussed include: Merging of the USGS Atlas of Mercury 1:5,000,000 Geologic Series; Geologic Mapping of the V-36 Thetis Regio Quadrangle: 2008 Progress Report; Structural Maps of the V-17 Beta Regio Quadrangle, Venus; Geologic Mapping of Isabella Quadrangle (V-50) and Helen Planitia, Venus; Renewed Mapping of the Nepthys Mons Quadrangle (V-54), Venus; Mapping the Sedna-Lavinia Region of Venus; Geologic Mapping of the Guinevere Planitia Quadrangle of Venus; Geological Mapping of Fortuna Tessera (V-2): Venus and Earth's Archean Process Comparisons; Geological Mapping of the North Polar Region of Venus (V-1 Snegurochka Planitia): Significant Problems and Comparisons to the Earth's Archean; Venus Quadrangle Geological Mapping: Use of Geoscience Data Visualization Systems in Mapping and Training; Geologic Map of the V-1 Snegurochka Planitia Quadrangle: Progress Report; The Fredegonde (V-57) Quadrangle, Venus: Characterization of the Venus Midlands; Formation and Evolution of Lakshmi Planum (V-7), Venus: Assessment of Models using Observations from Geological Mapping; Geologic Map of the Meskhent Tessera Quadrangle (V-3), Venus: Evidence for Early Formation and Preservation of Regional Topography; Geological Mapping of the Lada Terra (V-56) Quadrangle, Venus: A Progress Report; Geology of the Lachesis Tessera Quadrangle (V-18), Venus; Geologic Mapping of the Juno Chasma Quadrangle, Venus: Establishing the Relation Between Rifting and Volcanism; Geologic Mapping of V-19, V-28, and V-53; Lunar Geologic Mapping Program: 2008 Update; Geologic Mapping of the Marius Quadrangle, the Moon; Geologic Mapping along the Arabia Terra Dichotomy Boundary: Mawrth Vallis and Nili Fossae, Mars: Introductory Report; New Geologic Map of the Argyre Region of Mars; Geologic Evolution of the Martian Highlands: MTMs -20002, -20007, -25002, and -25007; Mapping Hesperia Planum, Mars; Geologic Mapping of the Meridiani Region, Mars; Geology of Holden Crater and the Holden and Ladon Multi-Ring Impact Basins, Margaritifer Terra, Mars; Geologic Mapping of Athabasca Valles; Geologic Mapping of MTM -30247, -35247 and -40247 Quadrangles, Reull Vallis Region of Mars; Geologic Mapping of the Martian Impact Crater Tooting; Geology of the Southern Utopia Planitia Highland-Lowland Boundary Plain: First Year Results and Second Year Plan; Mars Global Geologic Mapping: Amazonian Results; Recent Geologic Mapping Results for the Polar Regions of Mars; Geologic Mapping of the Medusae Fossae Formation on Mars (MC-8 SE and MC-23 NW) and the Northern Lowlands of Venus (V-16 and V-15); Geologic Mapping of the Zal, Hi'iaka, and Shamshu Regions of Io; Global Geologic Map of Europa; Material Units, Structures/Landforms, and Stratigraphy for the Global Geologic Map of Ganymede (1:15M); and Global Geologic Mapping of Io: Preliminary Results.

Bleamaster, Leslie F., III↗

Science Plan for Improving 3-D Seismic Velocity Models in the San Francisco Bay Region, 2019-24

This five-year science plan outlines short-term and long-term goals for improving 3-D seismic velocity models in the greater San Francisco Bay region as well as how to foster a community effort in reaching those goals. The short-term goals focus on improving the current USGS San Francisco Bay region geologic and seismic velocity model using existing data. The long-term goals focus on acquiring new data and leveraging better analysis tools to improve the model and characterize the uncertainty. The plan describes opportunities for contributions by members of the community to development of these seismic velocity models, provides current and potential users with general information on where efforts will likely be focused to improve these models and how new versions will be released, and outlines funding needs and obstacles for improving and maintaining such models. Several aspects of this plan, including how to foster a community effort, are independent of the geographic region and apply to other, similar efforts.

58 GEOSCIENCES↗

Photogrammetry using Apollo 16 orbital photography, part B

Discussion is made of the Apollo 15 and 16 metric and panoramic cameras which provided photographs for accurate topographic portrayal of the lunar surface using photogrammetric methods. Nine stereoscopic models of Apollo 16 metric photographs and three models of panoramic photographs were evaluated photogrammetrically in support of the Apollo 16 geologic investigations. Four of the models were used to collect profile data for crater morphology studies; three models were used to collect evaluation data for the frequency distributions of lunar slopes; one model was used to prepare a map of the Apollo 16 traverse area; and one model was used to determine elevations of the Cayley Formation. The remaining three models were used to test photogrammetric techniques using oblique metric and panoramic camera photographs. Two preliminary contour maps were compiled and a high-oblique metric photograph was rectified.

Wu, S. S. C.↗

A new process sensitivity index to identify important system processes under process model and parametric uncertainty

Hydrological models are always composed of multiple components that represent processes key to intended model applications. When a process can be simulated by multiple conceptual-mathematical models (process models), model uncertainty in representing the process arises. While global sensitivity analysis methods have been widely used for identifying important processes in hydrologic modeling, the existing methods consider only parametric uncertainty but ignore the model uncertainty for process representation. To address this problem, this study develops a new method to probe multimodel process sensitivity by integrating the model averaging methods into the framework of variance-based global sensitivity analysis, given that the model averaging methods quantify both parametric and model uncertainty. A new process sensitivity index is derived as a metric of relative process importance, and the index includes variance in model outputs caused by uncertainty in both process models and model parameters. For demonstration, the new index is used to evaluate the processes of recharge and geology in a synthetic study of groundwater reactive transport modeling. The recharge process is simulated by two models that converting precipitation to recharge, and the geology process is also simulated by two models of different parameterizations of hydraulic conductivity; each process model has its own random parameters. The new process sensitivity index is mathematically general, and can be applied to a wide range of problems in hydrology and beyond.

Dai, Heng↗

Polar Heat Flow on Io

Recently, Galileo spacecraft data have revealed Io's polar regions to be much warmer than previously expected. This unexpected development came from Photo-Polarimeter Radiometer (PPR) data which show that the minimum night temperatures are in the range of 90-95 K virtually everywhere on Io. The minimum night temperatures show no dependence upon latitude and, when away from the sunset terminator, they show no dependence upon time of night. This is indeed bizarre behavior for surface units which generally had been assumed to be passive with respect to Io's pervasive volcanism. Night temperatures of 90-95 K at high, polar latitudes are particularly hard to explain. Even assuming infinite thermal inertia, at these latitudes there is insufficient sunlight to support these warm night temperatures. Thus, through the process of elimination of other possibilities, we come to the conclusion that these surfaces are volcanically heated. Taking previously passive units and turning them into new sources of heat flow is a radical departure from previous thermophysical model paradigms. However, the geological interpretation is straight forward. We are simply seeing the effect of old, cool lava flows which cover most of the surface of Io but yet have some heat to radiate. Under these new constraints, we have taken on the challenge of formulating a physical model which quantitatively reproduces all of the observations of Io's thermal emission. In the following we introduce a new parametric model which suffices to identify a previously unrecognized polar component of Io's heat flow.

Veeder, G. J.↗

Simulation of Spent Nuclear Fuel loading into a Final Waste Repository

Previous work towards studying repository loading have used: - spent nuclear fuel (SNF) with an average burnup composition; - a lumped capacitance thermal model for calculating temperature in a Cyclus repository model. The goal of this work is to improve on the repository models and use U.S. historical SNF inventory data in simulations to more accurately study the loading of a waste repository. These goals will be achieved by: - using UNF-ST-and-DARDS Unified Database (UDB) that has historic assembly-specific data (e.g, isotopic composition, heat) in Cyclus simulations; - implementing a more accurate thermal model within a Cyclus repository model. Objectives: - Create a Cyclus spent fuel conditioning model that packages spent fuel bundles into packages which have user-defined properties. - Create a Cyclus medium-fidelity repository model that gives accurate time and spatial dependent temperature values and loads the repository based on a user-selected loading strategy. Cyclus is an agent-based extensible framework for modeling flow of material through user-defined nuclear fuel cycles. In Cyclus, each facility in the fuel cycle is modeled individually and the facilities interact with one another as independent agents. The spent fuel conditioning model accepts spent fuel bundles and puts them into a cylindrical waste package. In the spent fuel conditioning model, the user can define variables: For each layer, radius, thermal conductivity, thermal diffusivity; For each package, Number of spent fuel bundles, Radius and height. The waste repository model accepts waste packages and emplaces them into available positions within the waste repository based on a thermal criteria. The thermal criteria is a temperature limit at the interface between the waste package surface and the host geology, that is set based on the repository's host geology. In the waste repository model, the user can define the variables: Capacity, Distance between waste packages, Distance between drifts, Repository host geology, Loading Strategy. After the addition of new waste packages at each time step, the waste repository model recalculates the temperature at each location in the repository. If the addition of this new package causes its temperature to exceed the thermal limit, it will be placed back into the buffer. A thermal model that relies on a transient 'outside' model and quasi-steady-state 'inside' model is used to accurately determine the temperature in the repository. The 'outside' model assumes a homogenous medium with the Engineered Barrier System (EBS) replaced by the geologic medium. Temperature solutions for the central waste package, adjacent point and line sources are superimposed to calculate the temperature at specific points in the repository. The central drift consists of one finite line source which represents the central waste package. The central drift also consists of point sources that represent neighboring waste packages in the central drift. The neighboring drifts are represented by infinite line sources. The 'inside' model is considered to be at a quasi-steady-state condition because EBS has a relatively low thermal mass compared to the infinite geologic medium. The steady state calculation is performed at each time step with the heat source and interface temperature as boundary conditions. Future Work: Run Cyclus simulations with U.S. historical SNF inventory data, the spent fuel conditioning and repository models to study how waste package acceptance strategies impact repository loading.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Bayesian seismic inversion based on rock-physics prior modeling for the joint estimation of acoustic impedance, porosity and lithofacies

We propose a Bayesian approach for seismic inversion to estimate acoustic impedance, porosity and lithofacies within the reservoir conditioned to post-stack seismic and well data. The link between elastic and petrophysical properties is given by a joint prior distribution for the logarithm of impedance and porosity, based on a rock-physics model. The well conditioning is performed through a background model obtained by well log interpolation. Two different approaches are presented: in the first approach, the prior is defined by a single Gaussian distribution, whereas in the second approach it is defined by a Gaussian mixture to represent the well data multimodal distribution and link the Gaussian components to different geological lithofacies. The forward model is based on a linearized convolutional model. For the single Gaussian case, we obtain an analytical expression for the posterior distribution, resulting in a fast algorithm to compute the solution of the inverse problem, i.e. the posterior distribution of acoustic impedance and porosity as well as the facies probability given the observed data. For the Gaussian mixture prior, it is not possible to obtain the distributions analytically, hence we propose a Gibbs algorithm to perform the posterior sampling and obtain several reservoir model realizations, allowing an uncertainty analysis of the estimated properties and lithofacies. Both methodologies are applied to a real seismic dataset with three wells to obtain 3D models of acoustic impedance, porosity and lithofacies. The methodologies are validated through a blind well test and compared to a standard Bayesian inversion approach. Using the probability of the reservoir lithofacies, we also compute a 3D isosurface probability model of the main oil reservoir in the studied field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Noachian/Hesperian Hiatus and Erosive Reactivation of Martian Valley Networks

Despite new evidence for persistent flow and sedimentation on early Mars, it remains unclear whether valley networks were active over long geologic timescales (10(exp 5)-10(exp 8) yr), or if flows were persistent only during multiple discrete episodes of moderate (approx. 10(exp 4) yr) to short (<10 yr) duration. Understanding the long-term stability/variability of valley network hydrology would provide an important control on paleoclimate and groundwater models. Here we describe geologic evidence for a hiatus in highland valley network activity while the fretted terrain formed, followed by a discrete reactivation of persistent (but possibly variable) erosive flows. Additional information is included in the original extended abstract.

Irwin, R. P., III.↗

Optimization of radar imaging system parameters for geological analysis

The use of radar image simulation to model terrain variation and determine optimum sensor parameters for geological analysis is described. Optimum incidence angle is determined by the simulation, which evaluates separately the discrimination of surface features possible due to terrain geometry and that due to terrain scattering. Depending on the relative relief, slope, and scattering cross section, optimum incidence angle may vary from 20 to 80 degrees. Large incident angle imagery (more than 60 deg) is best for the widest range of geological applications, but in many cases these large angles cannot be achieved by satellite systems. Low relief regions require low incidence angles (less than 30 deg), so a satellite system serving a broad range of applications should have at least two selectable angles of incidence.

Waite, W. P.↗

Transport of perfluorocarbon tracers in the Cranfield Geological Carbon Sequestration Project

A field–scale carbon dioxide (CO 2 ) injection pilot project was conducted by the Southeast Regional Carbon Sequestration Partnership (SECARB) at Cranfield, Mississippi. Two associated campaigns in 2009 and 2010 were carried out to co–inject perfluorocarbon tracers (PFTs) and sulfur hexafluoride (SF 6 ) with CO 2 . Tracers in gas samples from two observation wells were analyzed to construct breakthrough curves. In this work, we present the field data and numerical modeling of the flow and transport of CO 2 , brine, and tracers. A high–resolution static model of the formation geology in the detailed area study (DAS) was used to capture the impact of connected flow pathways created by fluvial channels on breakthrough curves and breakthrough times of PFTs and SF6 tracers. We use the cubic–plus–association (CPA) equation of state, which takes into account the polar nature of water molecules, to describe the phase behavior of CO 2 –brine–tracer mixtures. Our simulated results show good agreement for the 2009 tracer campaign in Cranfield, while a larger discrepancy emerges by 2010. Here, the combination of multiple tracer injection pulses with detailed numerical simulations proves to be a powerful tool in constraining both formation properties and how complex flow paths develop over time.

54 ENVIRONMENTAL SCIENCES↗