Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SAND”

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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Corrosion Behavior of B and Ti Grain-Refined Sr-Modified A356

Corrosion behavior of A356 aluminum alloy with different grain refiners was evaluated. For this purpose, the Sr-modified A356 alloy was grain refined with Al5Ti1B and Al3B, and then, the corrosion characteristics were determined. For the mechanical properties, impact and Vickers hardness tests were carried out. To determine the corrosion behavior, both the electrochemical impedance spectroscopy technique and by potentiodynamic polarization plots were carried out in a 3.5 wt.% NaCl solution at 25 °C. The corrosion resistance of A356 alloy was dramatically improved in sequence with the following modifications compared to the unmodified castings: Sr, Ti and B additions. This situation applies to both different cooling rates, namely die and sand casting. The increase in corrosion rate was more than 80% for all alterations of the microstructure. Additionally, charge transfer resistance was increased ten times in die cast samples, while this value was 2.5 for sand cast samples compared to the unmodified alloy. It is important to note that sand casting of the unmodified alloy was improved four times with regard to the die cast unmodified A356.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Mobilization and Transport of Organic Compounds from Reservoir Rock and Caprock in Geological Carbon Sequestration Sites

Supercritical CO2 (scCO2) is an excellent solvent for organic compounds, including benzene, toluene, ethyl-benzene, and xylene (BTEX), phenols, and polycyclic aromatic hydrocarbons (PAHs). Monitoring results from geological carbon sequestration (GCS) field tests has shown that organic compounds are mobilized following CO2 injection. Such results have raised concerns regarding the potential for groundwater contamination by toxic organic compounds mobilized during GCS. Knowledge of the mobilization mechanism of organic compounds and their transport and fate in the subsurface is essential for assessing risks associated with GCS. Extraction tests using scCO2 and methylene chloride (CH2Cl2) were conducted to study the mobilization of volatile organic compounds (VOCs, including BTEX), the PAH naphthalene, and n-alkanes (n-C20 – n-C30) by scCO2 from representative reservoir rock and caprock obtained from depleted oil reservoirs and coal from an enhanced coal-bed methane recovery site. More VOCs and naphthalene were extractable by scCO2 compared to the CH2Cl2 extractions, while scCO2 extractable alkane concentrations were much lower than concentrations extractable by CH2Cl2. In addition, dry scCO2 was found to extract more VOCs than water saturated scCO2, but water saturated scCO2 mobilized more naphthalene than dry scCO2. In sand column experiments, moisture content was found to have an important influence on the transport of the organic compounds. In dry sand columns the majority of the compounds were retained in the column except benzene and toluene. In wet sand columns the mobility of the BTEX was much higher than that of naphthalene. Based upon results determined for the reservoir rock, caprock, and coal samples studied here, the risk to aquifers from contamination by organic compounds appears to be relatively low; however, further work is necessary to fully evaluate risks from depleted oil reservoirs.

Zhong, Lirong↗
DOE OSTI2019

Additively Manufactured Single-Use Molds and Reusable Patterns for Large Automotive and Hydroelectric Components

Additive manufacturing (AM) has the potential to reduce costs in the casting industry, particularly for large parts (i.e., components with any dimension > 12 in.). In this study, pilot-scale production case studies of hydroelectric and automotive parts were used to evaluate two AM-enabled casting methods: the direct printing of sand molds by binder jet and printing of reusable cope and drag patterns for sand casting. An additional benefit was found when additively manufactured patterns were used in combination with heat-treat free aluminum alloys, to produce castings of complex hydrodynamic surfaces. These parts previously required production via subtractive machining due to part distortion induced by the quench step of heat treatments. The machine types used for the three case studies are sand-binder jet, high-resolution polymer fused deposition modeling (FDM), and big area additive manufacturing (BAAM) FDM in combination with CNC machining. Each method demonstrated distinct advantages over traditional casting practices in particular use cases. Single-use molds show great reductions in start-up cost to produce one-off or legacy parts, and additively manufactured impression patterns show promise for innovating the tooling design process for complex geometry and large castings.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Numerical analysis of experimental studies of methane hydrate formation in a sandy porous medium

Here we analyse numerically an earlier experimental study that involved the formation of methane hydrates by the excess water method in a small reactor filled with a sandy porous medium, and seek to address questions about the type of the hydration reaction and the phase heterogeneity in the resulting hydrate-bearing sand. Using a fine discretization describing the reactor assembly, the experimental process is faithfully replicated numerically. The multi-stage process of hydrate formation is subdivided in 7 steps. The experimental data from the continuously-monitored pressure and temperature during each step are used for comparison against the numerical predictions, the identification of the dominant processes and the determination of the associated parameters through a history-matching process that minimizes deviations between observations and simulation results. The results of this first-ever study on this subject demonstrate unequivocally that the hydration reaction is a kinetic (as opposed to an equilibrium) process, and that the spatial distributions of the various phases (aqueous, gas and hydrate) at the end of the formation process are strongly heterogeneous. This has serious implications in simulation studies of hydrate dissociation that assume uniform initial phase saturation distributions. Finally, the history-matching process indicates that (a) the system behaviour is sensitive to some flow parameters (porosity and irreducible water saturation) only during the first water injection, (b) it is insensitive to the sand intrinsic permeability during all steps of the study, and (c) thermal processes dominate after the first water injection, yielding estimates of the thermal properties of the sand and of time-variable key parameters of the kinetic reaction.

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

Modeling the impact of bubbling bed hydrodynamics on tar yield and its fluctuations during biomass fast pyrolysis

The impact of bubbling bed hydrodynamics on temporal variations in the exit tar yield for biomass fast pyrolysis was investigated using computational simulations of an experimental laboratory-scale reactor. A multi-fluid computational fluid dynamics model was employed to simulate the differential conservation equations in the reactor, and this was combined with a multi-component, multi-step pyrolysis kinetics scheme for biomass to account for chemical reactions. The predicted mean tar yields at the reactor exit appear to match corresponding experimental observations. Parametric studies predicted that increasing the fluidization velocity should improve the mean tar yield but increase its temporal variations. Increases in the mean tar yield coincide with reducing the diameter of sand particles or increasing the initial sand bed height. However, trends in tar yield variability are more complex than the trends in mean yield. The standard deviation in tar yield reaches a maximum with changes in sand particle size. As a result, the standard deviation in tar yield increases with the increases in initial bed height in freely bubbling state, while reaches a maximum in slugging state.

09 BIOMASS FUELS↗
DOE OSTI2019

Hydromechanical modelling to evaluate impact of fault structure on CO 2 migration in stacked storage system

Faulting processes produce a complex fault damage zone (FDZ) having permeability variations over several orders of magnitude within a very short distance. In this study we use an example normal fault system from a depleted oil field in southern Louisiana with multiple stacked sand beds along the fault structure. Geostatistical techniques are used to introduce heterogeneities of several orders of magnitude over very short distances that honor the across- and along-fault architecture of the system. A total of 15 cases are simulated to account for the impact of a fragmented core, an intact core and a naturally fractured core on inter-sand leakage rates. The presented results show the complex nature of CO 2 migration within the fault structure, where an interplay between overcoming the fault core’s high capillary entry pressure and stress propagation in the vertical direction dictates the leakage direction and rates. It is observed that hydromechanical evolution of fracture flow properties may result in significant leakage from the storage zone. In one case with small caprock thickness, a loss of 18% of the total injected mass is observed. The results presented highlight some unique features of representative fault related leakage in stacked sand systems and thus will be beneficial for storage integrity analysis in these systems.

58 GEOSCIENCES↗
DOE OSTI2023

Hydrogel-assisted self-healing of biomineralized living building materials

Living building materials (LBMs) are an emergent class of structural materials that leverage the biomineralization capability of microorganisms within sand-hydrogel scaffolds to produce living, load-bearing structures. Here we produced LBMs using a physically crosslinkable sand-hydrogel scaffold and two microorganisms with different biomineralization pathways - Synechococcus sp. PCC 7002 (photosynthetic) and S. pasteurii (ureolytic) - and investigated their self-healing capacity. Our results reveal that both Synechococcus sp. PCC 7002 and S. pasteurii demonstrated exceptional viability within all LBMs for more than 20 days. Damaged LBMs containing Synechococcus sp. PCC 7002 exhibited 103% and 112% of their original compressive and flexural strengths, respectively, after three days of healing at 50% relative humidity (RH) (i.e., ambient conditions). In contrast, LBMs containing S. pasteurii exhibited 71% and 66% of their original compressive and flexural strengths, respectively, after three days of healing at 50% RH. The compressive and flexural strengths of all LBMs rebounded to 93-100% after seven days of healing at 50% RH. Data substantiate that the self-healing ability of the hydrogel plays a critical role in facilitating healing of LBMs, as evidenced by the 82-118% rebounds in compressive and flexural strengths by the sand-hydrogel scaffold alone after three or seven days of healing at 50% RH. Healing was less effective at 100% RH for all LBMs investigated herein.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Estimating biofuel contaminant concentration from 4D ERT with mixing models

Here, we present the results of a lab-scaled feasibility study to assess the performance of electrical resistivity tomography for detection, characterization, and monitoring of fuel grade ethanol releases to the subsurface. Further, we attempt to determine the concentration distribution of the ethanol from the electrical resistivity tomography data using mixing-models. Ethanol is a renewable fuel source as well as an oxygenate fuel additive currently used to replace the known carcinogen methyl tert-butyl ether; however, ethanol is preferentially biodegraded and a cosolvent. When introduced to areas previously impacted by nonethanol-based fuels, it will facilitate the persistence of carcinogenic fuel compounds like benzene and ethylbenzene, as well as remobilize them to the ground water. These compounds would otherwise be retained in the soil column undergoing active or passive remediation processes such as soil vapor extraction or natural attenuation. Here, we introduce ethanol to a saturated Ottawa sand in a tank instrumented for four-dimensional geoelectrical measurements. Forward model results suggest pure phase ethanol released into a water saturated silica sand should present a detectable target for electrical resistivity tomography relative to a saturated silica sand only. We observe the introduction of ethanol to the closed hydraulic system and subsequent migration over the duration of the experiment. One-dimensional and three–dimensional temporal data are assessed for the detection, characterization, and monitoring of the ethanol release. Results suggest one-dimensional geoelectrical measurements may be useful for monitoring a release, while three-dimensional geoelectrical field imaging would be useful to characterize, monitor, and design effective remediation approaches for an ethanol release, assuming field conditions do not preclude the application of geoelectrical methods. We then attempt to use predictive mixing models to calculate the distribution of ethanol concentration within the measurement domain. For this study we examine four different models: a nested parallel mixing model, a nested cubic mixing model, the complex refractive index model (CRIM), and the Lichtenecker-Rother (L-R) model. The L-R model, modified to include an electrical formation factor geometry term, provided the best agreement with expected EtOH concentrations.

09 BIOMASS FUELS↗
DOE OSTI2010

Structural resistance of chemically modified 1-D nanostructured titanates in inorganic acid environment

Sodium containing one-dimensional nanostructured layered titanates (1-D NSLT) were produced both from commercial anatase powder and Brazilian natural rutile mineral sands by alkali hydrothermal process. The 1-D NSLT were chemically modified with proton, cobalt or iron via ionic exchange and all products were additionally submitted to intensive inorganic acid aging (pH = 0.5) for 28 days. The morphology and crystal structure transformations of chemically modified 1-D NSLT were followed by transmission electron microscopy, powder X-ray diffraction, selected area electron diffraction and energy dispersive spectroscopy. It was found that the original sodium rich 1-D NSLT and cobalt substituted 1-D NSLT were completely converted to rutile nanoparticles, while the protonated form was transformed in a 70%-30% (by weight) anatase-rutile nanoparticles mixture, very similar to that of the well-known TiO{sub 2}-photocatalyst P25 (Degussa). The iron substituted 1-D NSLT presented better acid resistance as 13% of the original structure and morphology remained, the rest being converted in rutile. A significant amount of remaining 1-D NSLT was also observed after the acid treatment of the product obtained from rutile sand. The results showed that phase transformation of NSLT into titanium dioxide polymorph in inorganic acid conditions were controllable by varying the exchanged cations. Finally, the possibility to transform, through acid aging, 1-D NSLT obtained from Brazilian natural rutile sand into TiO{sub 2}-polymorphs was demonstrated for the first time to the best of authors' knowledge, opening path for producing TiO{sub 2}-nanoproducts with different morphologies through a simple process and from a low cost precursor.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Influence of living grass Roots and endophytic fungal hyphae on soil hydraulic properties

Soil hydraulic properties are often estimated based on laboratory data or pedotransfer functions dependent on soil physical properties, which often do not consider potential impacts of soil roots or fungal hyphae. Here, we first review current knowledge of how these soil biotic components affect hydraulic properties, then we conducted laboratory experiments to specifically test if the presence of roots and mycorrhizal fungi had a significant effect on the hydraulic properties of two soils with contrasting textures: Flint sand and Hamblen silt loam. Soil cores were seeded with (Panicum virgatum) and grown in a greenhouse over three separate growth periods. The endophytic fungus Serendipita indica was injected as liquid inoculant into designated mycorrhizal cores. Saturated hydraulic conductivity (K sat ) measurements were made with a constant head permeameter, and soil water retention curves were obtained by the evaporation method, supplemented at the dry end for Hamblen silt loam with water activity meter data. Retention curve parameters were obtained by fitting the van Genuchten equation to the resulting measurements. Mean root volume ratios were higher in the mycorrhizal inoculated treatment than in the uninoculated treatment for both soils. For Flint sand, analysis of variance revealed that K sat was reduced by the presence of roots as compared to bare soil. This was likely due to roots clogging soil pores. Results also indicated the presence of roots changed the shape of the water retention curve for Flint sand by increasing water content at saturation and by reducing the slope of the curve. These changes suggested roots created additional porosity and broadened the pore-size distribution. The presence of mycorrhizal fungi accentuated the root effects. The influence of roots and mycorrhizal fungi on hydraulic properties was less obvious for the Hamblen silt loam, as none of the treatments differed from each other at p < 0.05. The results highlight the necessity to consider the impact of root and fungal structures on models of soil hydraulic properties.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

In situ treatment of arsenic contaminated groundwater by aquifer iron coating: Experimental Study

In situ arsenic removal from groundwater by an iron coating method has great potential to be a cost effective and simple groundwater remediation technique, especially in rural and remote areas where groundwater is used as the main source of drinking water. The in situ arsenic removal technique was first optimized by simulating arsenic removal in various quartz sand columns under anoxic conditions., Its effectiveness was then evaluated in an actual high-arsenic groundwater environment. The mechanism of arsenic removal by the iron coating was investigated under different conditions using scanning electron microscopy (SEM)/X-ray absorption spectroscopy, an electron microprobe, and Fourier transformation infrared spectroscopy. A 4-step alternative cycle aquifer iron coating method was developed. A continuous injection of 5 mmol/L FeSO4 and 2.5 mmol/L NaClO for 96 hours can create a uniform coating of crystalline goethite on the surface of quartz sand in the columns without causing clogging. At a flow rate of 0.45 cm/min of the injection reagents (vi), the time for arsenic (as Na2HAsO4) to pass through the iron-coated quartz sand column was approximately 35 hours, which was much longer than that for tracer fluorescein sodium (approximately 2 hours). The retardation factor of arsenic was 23, and its adsorption capacity was 0.11 mol As per mol Fe, leading to an excellent arsenic removal. In situ arsenic removal from groundwater in an aquifer was achieved by simultaneous injections of As (V) and Fe (II) reagents. When the arsenic content in the groundwater was 233 μg/L, the aqueous phase arsenic was completely removed with an arsenic adsorption of 0.05 mol As per mol Fe. Arsenic fixation resulted from a process of adsorption/co-precipitation, in which arsenic and iron likely formed the arsenic-bearing iron mineral phases with poor crystallinity by way of bidentate binuclear complexes. Thus, the high arsenic removal efficiency of the technique likely resulted from the expanded specific iron oxide/hydroxide surface area with poor crystallinity and from coprecipitation.

arsenic, remediation, iron, coating↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2019

Wettability Effects on Primary Drainage Mechanisms and NAPL Distribution: A Pore-Scale Study

The pore-scale processes governing water drainage behavior in porous media have implications for geoscience multiphase scenarios including carbon capture and storage, contaminant site remediation, oil recovery, and vadose zone processes. However, few studies report directly observed pore-scale water drainage phenomena in 3-D soils. This knowledge gap limits our ability to verify assumptions underlying existing models and develop optimal solutions. As such, this paper utilizes synchrotron X-ray microtomography to present an experimental pore-scale examination of nonaqueous phase liquid (NAPL)/water distribution along a primary drainage front as dense NAPL was injected upward into water wetting (WW) and intermediate wetting (IW) sand-packed columns. Pore-network structures were extracted from imaged data sets and mapped onto segmented NAPL/water data sets which allowed quantitative examinations of wettability impacts on (a) the extent to which NAPL fills individual pore bodies and (b) relationships between pore size and the phase occupying the pore, with both considered as a function of distance (and capillary pressure) relative to the NAPL front. These results revealed that several hypotheses treating IW sand similarly to WW sands are simplistic. IW systems exhibited a sequence of pore filling that deviated from traditional capillary pressure-based model predictions: NAPL invades smaller pores, while larger, adjacent pores are bypassed leaving multipore residual water ganglia. NAPL pore saturations were close to 1 and did not change with capillary pressure in IW systems. Overall, the results illustrate how a relatively small change in operative contact angle alters NAPL distribution during water drainage, with important implications for geoscience multiphase flow scenarios.

58 GEOSCIENCES↗
DOE OSTI2020

Laboratory Electrical Conductivity of Marine Gas Hydrate

Methane hydrate was synthesized from pure water ice and flash frozen seawater, with varying amounts of sand or silt added. Electrical conductivity was determined by impedance spectroscopy, using equivalent circuit modeling to separate the effects of electrodes and to gain insight into conduction mechanisms. Silt and sand increase the conductivity of pure hydrate; we infer by contaminant NaCl contributing to conduction in hydrate, to values in agreement with resistivities observed in well logs through hydrate–saturated sediment. The addition of silt and sand lowers the conductivity of hydrate synthesized from seawater by an amount consistent with Archie's law. Here, all samples were characterized using cryogenic scanning electron microscopy and energy dispersive spectroscopy, which show good connectivity of salt and brine phases. Electrical conductivity measurements of pure hydrate and hydrate mixed with silt during pressure–induced dissociation support previous conclusions that sediment increases dissociation rate.

58 GEOSCIENCES↗
DOE OSTI2024

Microtomographic Measurements of Total Air‐Water Interfacial Areas for Soils

Abstract Synchrotron X‐ray microtomography (XMT) was used to measure total air‐water interfacial areas (A aw ) as a function of water saturation (S w ) for several soils that comprise a range of physical and geochemical properties. Measurements were also conducted for glass beads and quartz sands for comparison. Apparent near‐linearA aw ‐S w relationships are observed for the three sands and the three sandy soils. In contrast, the measured interfacial areas for two soils that contain greater proportions of silt and clay are strongly nonlinear functions of water saturation. The greater degree of nonlinearity observed for these two soils is due to their much greater particle‐size distributions (i.e., uniformity coefficients) and their concomitant greater range in pore sizes. Interfacial areas determined with the thermodynamic method were used to benchmark the XMT measurements. XMT‐measured interfacial areas compare well to the thermodynamic‐determined values for the sands and sandy soils. In contrast, the XMT‐measured interfacial areas for the two soils with larger particle‐size distributions are not fully congruent with the thermodynamic‐determined values. Both of these soils have large fractions of pore space comprising nominal pore diameters smaller than the resolution of the XMT imaging. These results suggest that air‐water interfacial area may not always be fully characterized by standard XMT for soils with large particle‐size distributions.

Environmental Sciences & Ecology↗
DOE OSTI2022

Trimethyloxonium-mediated methylation strategies for the rapid and simultaneous analysis of chlorinated phenols in various soils by electron impact gas chromatography–mass spectrometry

Abstract The efficient methylation of a panel of five industrial and environmentally-relevant chlorophenols (CPs) employing trimethyloxonium tetrafluoroborate (TMO) for their qualitative detection and identification by electron impact gas chromatography–mass spectrometry (EI-GC–MS) is presented. The protocol’s execution is simple and smoothly converts the phenols into their O-methylated counterparts conveniently at ambient temperature. The efficiency of two versions of the protocol was successfully tested in their ability to simultaneously derivatize five CPs (2-chlorophenol, 2,4-dichlorophenol, 2,4,6-trichlorophenol, pentachlorophenol and triclosan) in six distinct, separate soil matrices (Nebraska EPA standard soil, Virginia Type A soil, Ottawa sand, Baker sand, Silt and Georgia EPA standard soil) when present at low levels (~ 10 μgg −1 ). The first version involves the direct derivatization of the spiked soils with the methylating salt while the second one involves an initial soil extraction step of the CPs followed by methylation. The MDL values for each methylated CP were determined and lower values were found (4.1–13.2 ng . mL −1 ) for both sand matrices (Ottawa and Baker) as well as for the Georgia EPA standard soil, while larger values (8.2–21.8 ng . mL −1 ) were found for the Virginia Type soil, Nebraska EPA standard soil and Silt. The presented protocol offers a safer and more practical alternative to the universally employed diazomethane method and can be readily applicable to matrices other than soils. Furthermore, the protocols described herein may find applicability to the methylation of other analytes bearing acidic protons.

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

Assessing the potential of composite confining systems for secure and long-term CO 2 retention in geosequestration

A potential geologic target for CO 2 storage should ensure secure containment of injected CO 2 . Traditionally, this objective has been achieved by targeting reservoirs with overlying seals-regionally extensive, low permeability units that have been proven capable of retaining buoyant fluid accumulations over geologic time. However, considering that the amount of CO 2 is limited by a decadal injection period, vertical migration of CO 2 can be effectively halted by a composite system of discontinuous shale/silt/mudstone barriers in bedded sedimentary rocks. Here, we studied the impact of depositional architectures in a composite confining system on CO 2 migration and confinement at reservoir scale. We stochastically generated lithologically heterogeneous reservoir models containing discontinuous barriers consistent with statistical distributions of net-sand-to-gross-shale ratio (NTG) and horizontal correlation lengths derived from well log data and observations of producing hydrocarbon fields in Southern Louisiana. We then performed an extensive suite of reservoir simulations of CO 2 injection and post-injection to evaluate the sensitivity of CO 2 plume migration and pressure response of the composite system to a series of geologic and fluid parameters including the lateral continuity of barriers, NTG, permeability anisotropy within the sand body, and capillary pressure contrast between the sand and shale facies. The results indicate that lateral continuity of barriers and NTG are the dominant controls on CO 2 plume geometry and pressure build-up in the reservoir, while the impact of NTG is particularly pronounced. The significance of intraformational barriers becomes apparent as they facilitate the local capillary trapping of CO 2 . Those barriers improve the pore space occupancy by promoting a more dispersed shape of the plume and ultimately retard the buoyancy-driven upward migration of the plume post injection.

58 GEOSCIENCES↗