Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SANDSTONES”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Challenges of Maintaining Institutional Controls: Maintenance of an Alternate Water Supply- A Case Study - 19583

The US Department of Energy (DOE) cofunded an alternate water-supply system (AWSS) with the Indian Health Service that was installed in the vicinity of the former uranium mill site near Riverton, Wyoming, as one of the institutional controls (ICs) at the site. The AWSS was constructed in 1998 to supply drinking water to residents within the IC boundary in lieu of drinking untreated groundwater that could potentially be impacted by the contaminated shallow aquifer. The AWSS derives water from three source wells that pump groundwater from the sandstones of the Wind River Formation to a 3.8 x 106 liter (1,000,000 gallon) storage tank. Water from the storage tank is distributed throughout the AWSS service area by approximately 18 kilometers (11 miles) of PVC pipeline. The AWSS currently serves 103 residential customers and 1 commercial customer. DOE and the Northern Arapaho tribe (owner of the AWSS) have partnered over the years to address numerous challenges with the AWSS that included engineering assessments, scientific investigations, sampling support, a hydrant flushing program, and funding. Challenges with the AWSS have included the following: (1) Transmission pipe installed during initial construction did not have a pressure rating high enough to account for system pressure surges. (2) AWSS appurtenances were initially installed without accounting for the effects of corrosive soils on metal components. (3) Natural radium from the Wind River Formation groundwater builds up in the AWSS transmission pipes to levels above the primary drinking water standard because of the lack of a flushing program, dead-end lines, and low system demand. A 2-year study showed that a unidirectional hydrant-flushing program will control the radium buildup. (4) Aging system components including hydrants, valves, and pumps are nearing their effective lifespan and are in need of replacement. (5) With the addition of large commercial customers, demand on the system has increased since the initial installation of the system, and future growth in the area will continue to increase demand. (6) Additional connections and integration of the AWSS with an areawide water system in the near future will require changes to the current flushing plan and add additional complexity regarding areas of responsibility. (7) A well-organized financial and administrative infrastructure is not in place, and additional revenue is needed to operate, maintain, and improve the AWSS. To be a viable IC, the AWSS will require management, maintenance, and upgrades to provide safe drinking water until groundwater meets drinking water standards. Responsibilities and funding for maintaining the AWSS as a viable IC are the focus of current discussions, and decisions regarding responsibilities and funding will be formalized in future cooperative agreements and memoranda of understanding with the owners. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

On the origin of nonlinear elasticity in disparate rocks

Dynamic acousto-elastic (DAE) studies are performed on a set of 6 rock samples (four sandstones, one soapstone, and one granite). From these studies, at 20 strain levels 10 -7 < ϵ < 10 -5 , four measures characterizing the nonlinear elastic response of each sample are found. Additionally, each sample is tested with nonlinear resonant ultrasonic spectroscopy (NRUS) and a fth measure of nonlinear elastic response is found. The ve measures of the nonlinear elastic response of the samples (approximately 3 x 6 x 20 x 5 numbers as each measurement is repeated 3 times) are subjected to careful analysis using model independent statistical methods, principal component analysis and fuzzy clustering. This analysis reveals di erences among the samples and di erences among the nonlinear measures. Four of the nonlinear measures are sensing much the same physical mechanism in the samples. The fth is seeing something di erent. This is the case for all samples. Although the same physical mechanisms (two) are operating in all samples there are distinctive features in the way the physical mechanisms present themselves from sample to sample. This suggests classi cation of the samples into two groups. The numbers in this study and the classi cation of the measures/samples constitute an empirical characterization of rock nonlinear elastic properties that can serve as a valuable testing ground for physically based theories that relate rock nonlinear elastic properties to microscopic elastic features.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Resonant ultrasound spectroscopy for materials with high damping and samples of arbitrary geometry

This paper describes resonant ultrasound spectroscopy (RUS) as a powerful and established technique for measuring elastic constants of a material with general anisotropy. The first step of this technique consists of extracting resonance frequencies and damping from the vibrational frequency spectrum measured on a sample with free boundary conditions. An inversion technique is then used to retrieve the elastic tensor from the measured resonance frequencies. As originally developed, RUS has been mostly applicable to (i) materials with small damping such that the resonances of the sample are well separated and (ii) samples with simple geometries for which analytical solutions exist. In this paper, these limitations are addressed with a new RUS approach adapted to materials with high damping and samples of arbitrary geometry. Resonances are extracted by fitting a sum of exponentially damped sinusoids to the measured frequency spectrum. The inversion of the elastic tensor is achieved with a genetic algorithm, which allows searching for a global minimum within a discrete and relatively wide solution space. First, the accuracy of the proposed approach is evaluated against numerical data simulated for samples with isotropic symmetry and transversely isotropic symmetry. Subsequently, the applicability of the approach is demonstrated using experimental data collected on a composite structure consisting of a cylindrical sample of Berea sandstone glued to a large piezoelectric disk. In the proposed experiments, RUS is further enhanced by the use of a 3-D laser vibrometer allowing the visualization of most of the modes in the frequency band studied.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Effect of pore pressure buildup on slowness of rupture propagation

Pore fluid pressure is known to play an important role in brittle fracture initiation and propagation, yet the underlying mechanisms remain unclear. We conducted triaxial experiments on saturated porous sandstones to investigate effects of pore pressure buildup on the slowness of shear rupture propagation at different confining pressures. At low to intermediate confinements, rocks fail by brittle faulting, and pore pressure buildup causes a reduction in rock’s shear strength but does not induce measurable differences in slip behavior. When the confinement is high enough to prohibit dynamic faulting, rocks fail in the brittle-ductile transitional regime. In the transitional regime, pore pressure buildup promotes slip instability on an otherwise stably sliding fracture. Compared to those observed in the brittle regime, the slip rate, stress drop, and energy dissipated during rupture propagation with concurrent pore pressure buildup in the transitional regime are distinctively different. When decreasing confining pressure instead, the slip behavior resembles the ones of the brittle regime, emphasizing h ow the observed slowness is related to excess pore pressure beyond the effective pressure phenomenon. Analysis of the mechanical data using existing theoretical models confirms these observations. Quantitative microstructural analyses reveal that increasing pore pressure lessens the dilatancy hardening during failure, thus enhances slip along the localized zone in the transitional regime. Lastly, our experimental results suggest that pore pressure buildup induces slow slip in the transitional regime, and slip rates along a shear fracture may vary considerably depending on effective stress states.

58 GEOSCIENCES↗
DOE OSTI2016

Observation of > 5 wt % zinc at the Kimberley outcrop, Gale crater, Mars

Zinc-enriched targets have been detected at the Kimberley formation, Gale crater, Mars, using the Chemistry Camera (ChemCam) instrument. The Zn content is analyzed with a univariate calibration based on the 481.2 nm emission line. The limit of quantification for ZnO is 3 wt % (at 95% confidence level) and 1 wt % (at 68% confidence level). The limit of detection is shown to be around 0.5 wt %. As of sol 950, 12 targets on Mars present high ZnO content ranging from 1.0 wt % to 8.4 wt % (Yarrada, sol 628). Those Zn-enriched targets are almost entirely located at the Dillinger member of the Kimberley formation, where high Mn and alkali contents were also detected, probably in different phases. Zn enrichment does not depend on the textures of the rocks (coarse-grained sandstones, pebbly conglomerates, and resistant fins). The lack of sulfur enhancement suggests that Zn is not present in the sphalerite phase. Zn appears somewhat correlated with Na 2 O and the ChemCam hydration index, suggesting that it could be in an amorphous clay phase (such as sauconite). On Earth, such an enrichment would be consistent with a supergene alteration of a sphalerite gossan cap in a primary siliciclastic bedrock or a possible hypogene nonsulfide zinc deposition where Zn, Fe, Mn would have been transported in a reduced sulfur-poor fluid and precipitated rapidly in the form of oxides.

58 GEOSCIENCES↗
DOE OSTI2016

The potassic sedimentary rocks in Gale Crater, Mars, as seen by ChemCam on board Curiosity : Potassic Sedimentary Rocks, Gale Crater

The Mars Science Laboratory rover Curiosity encountered potassium-rich clastic sedimentary rocks at two sites in Gale Crater, the waypoints Cooperstown and Kimberley. These rocks include several distinct meters thick sedimentary outcrops ranging from fine sandstone to conglomerate, interpreted to record an ancient fluvial or fluvio-deltaic depositional system. Furthermore, from ChemCam Laser-Induced Breakdown Spectroscopy (LIBS) chemical analyses, this suite of sedimentary rocks has an overall mean K 2 O abundance that is more than 5 times higher than that of the average Martian crust. The combined analysis of ChemCam data with stratigraphic and geographic locations then reveals that the mean K 2 O abundance increases upward through the stratigraphic section. Chemical analyses across each unit can be represented as mixtures of several distinct chemical components, i.e., mineral phases, including K-bearing minerals, mafic silicates, Fe-oxides, and Fe-hydroxide/oxyhydroxides. Possible K-bearing minerals include alkali feldspar (including anorthoclase and sanidine) and K-bearing phyllosilicate such as illite. Mixtures of different source rocks, including a potassium-rich rock located on the rim and walls of Gale Crater, are the likely origin of observed chemical variations within each unit. Physical sorting may have also played a role in the enrichment in K in the Kimberley formation. The occurrence of these potassic sedimentary rocks provides additional evidence for the chemical diversity of the crust exposed at Gale Crater.

58 GEOSCIENCES↗
DOE OSTI2016

Oxidation of manganese in an ancient aquifer, Kimberley formation, Gale crater, Mars

We report that the Curiosity rover observed high Mn abundances (>25 wt % MnO) in fracture-filling materials that crosscut sandstones in the Kimberley region of Gale crater, Mars. The correlation between Mn and trace metal abundances plus the lack of correlation between Mn and elements such as S, Cl, and C, reveals that these deposits are Mn oxides rather than evaporites or other salts. On Earth, environments that concentrate Mn and deposit Mn minerals require water and highly oxidizing conditions; hence, these findings suggest that similar processes occurred on Mars. In conclusion, based on the strong association between Mn-oxide deposition and evolving atmospheric dioxygen levels on Earth, the presence of these Mn phases on Mars suggests that there was more abundant molecular oxygen within the atmosphere and some groundwaters of ancient Mars than in the present day.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Stress distribution during cold compression of a quartz aggregate using synchrotron X-ray diffraction: Observed yielding, damage, and grain crushing

Herein, we report new experimental results that quantify the stress distribution within a quartz aggregate during pore collapse and grain crushing. The samples were probed with synchrotron X-ray diffraction as they were compressed in a multianvil deformation apparatus at room temperature from low pressure (tens of megapascal) to pressures of a few gigapascal. In such a material, stress is likely to concentrate at grain-to-grain contacts and vanish where grains are bounded by open porosity. Therefore, internal stress is likely to vary significantly from point to point in such an aggregate, and hence, it is important to understand both the heterogeneity and anisotropy of such variation with respect to the externally applied stress. In our quartz aggregate (grain size of ~4 μm), the measured diffraction peaks broaden asymmetrically at low pressure (tens of megapascal), suggesting that open pores are still a dominant characteristic of grain boundaries. In contrast, a reference sample of novaculite (a highly dense quartz polycrystal, grain size of ~6–9 μm) showed virtually no peak broadening with increasing pressure. In the quartz aggregate, we observed significant deviation in the pressure-volume curves in the range of P = 400–600 MPa. We suggest that this marks the onset of grain crushing (generally denoted as P* in the rock mechanic literature), which is commonly reported to occur in sandstones at pressures of this order, in general agreement with a Hertzian analysis of fracturing at grain contacts.

58 GEOSCIENCES↗
DOE OSTI2018

Monitoring and Modeling Caprock Integrity at the In Salah Carbon Dioxide Storage Site, Algeria

The In Salah gas development project, initiated in 2004, involves the long-term storage of waste carbon dioxide associated with natural gas production at several central Algerian fields. Rather then vent the extracted CO 2 into the atmosphere the operators, BP, Sonatrach, and Statoil chose to have the gas compressed, dehydrated, transported, and injected into a deep saline aquifer leg of the Krechba gas reservoir (Matheison et al. 2011). For this effort three horizontal state-of-the-art injection wells were drilled to a depth of 1500-1800m, using geo-steering technology to follow the roughly 20m thick reservoir. To maximize the injectivity, the wells were oriented perpendicular to the direction of maximum stress, the dominant fracture direction. Long-term plans were to store up to 17 million tonnes of carbon dioxide during the life of the program (Matheison et al. 2010). The project at In Salah was a pioneering effort in several respects and involved the investment of over 100 million dollars by the member companies. One feature that distinguishes In Salah from several other industrialscale CO 2 sequestration projects, such as Sleipner (Chadwick et al. 2012), Weyburn (White et al. 2011), Decatur (Kaven et al. 2014), and Aquastore (Worth et al. 2014) is the generally smaller storage volume in terms of overall reservoir thickness and lower average porosities of around 15%. With the exception of Weyburn, these other sites are characterized by generally thick sequences of high porosity sandstones. Therefore, one might expect that the pressure buildup near the wells could be larger at In Salah, with consequences for reservoir deformation and caprock strain. The role of fractures in providing storage and increasing reservoir permeability is one factor that would likely impact the pressure distribution around the well. However, the detailed effect of existing fractures on flow was not established prior to the injection at In Salah. The unique features of the In Salah storage site point to the importance of wellbore, geophysical, and geochemical monitoring to ensure that the integrity of the storage complex is maintained both during and after injection. In 2005 a Joint Industry Project (JIP) was set up, in an effort to use various geophysical, geochemical, and production techniques to monitor the fate of the injected carbon dioxide (Matheison et al. 2011). Here in this paper we describe the efforts made to monitoring the fate of the injected carbon dioxide, the state of the caprock, the estimated outcome of the storage effort, and some of the uncertainties that remain.

58 GEOSCIENCES↗
DOE OSTI2020

Discretization limits of lattice‐Boltzmann methods for studying immiscible two‐phase flow in porous media

Summary Digital images of porous media often include features approaching the image resolution length scale. The behavior of numerical methods at low resolution is therefore important even for well‐resolved systems. We study the behavior of the Shan‐Chen (SC) and Rothman‐Keller (RK) multicomponent lattice‐Boltzmann models in situations where the fluid‐fluid interfacial radius of curvature and/or the feature size of the medium approaches the discrete unit size of the computational grid. Various simple, small‐scale test geometries are considered, and a drainage test is also performed in a Bentheimer sandstone sample. We find that both RK and SC models show very high ultimate limits: in ideal conditions the models can simulate static fluid configuration with acceptable accuracy in tubes as small as three lattice units across for RK model (six lattice units for SC model) and with an interfacial radius of curvature of two lattice units for RK and SC models. However, the stability of the models is affected when operating in these extreme discrete limits: in certain circumstances the models exhibit behaviors ranging from loss of accuracy to numerical instability. We discuss the circumstances where these behaviors occur and the ramifications for larger‐scale fluid displacement simulations in porous media, along with strategies to mitigate the most severe effects. Overall we find that the RK model, with modern enhancements, exhibits fewer instabilities and is more suitable for systems of low fluid‐fluid miscibility. The shortcomings of the SC model seem to arise predominantly from the high, strongly pressure‐dependent miscibility of the two fluid components.

Li, Zhe↗
DOE OSTI2023

The Kimberlina synthetic multiphysics dataset for CO 2 monitoring investigations

Abstract We present a synthetic multi‐scale, multi‐physics dataset constructed from the Kimberlina 1.2 CO 2 reservoir model based on a potential CO 2 storage site in the Southern San Joaquin Basin of California. Among 300 models, one selected reservoir‐simulation scenario produces hydrologic‐state models at the onset and after 20 years of CO 2 injection. Subsequently, these models were transformed into geophysical properties, including P‐ and S‐wave seismic velocities, saturated density where the saturating fluid can be a combination of brine and supercritical CO 2 , and electrical resistivity using established empirical petrophysical relationships. From these 3D distributions of geophysical properties, we have generated synthetic time‐lapse seismic, gravity and electromagnetic responses with acquisition geometries that mimic realistic monitoring surveys and are achievable in actual field situations. We have also created a series of synthetic well logs of CO 2 saturation, acoustic velocity, density and induction resistivity in the injection well and three monitoring wells. These were constructed by combining the low‐frequency trend of the geophysical models with the high‐frequency variations of actual well logs collected at the potential storage site. In addition, to better calibrate our datasets, measurements of permeability and pore connectivity have been made on cores of Vedder Sandstone, which forms the primary reservoir unit. These measurements provide the range of scales in the otherwise synthetic dataset to be as close to a real‐world situation as possible. This dataset consisting of the reservoir models, geophysical models, simulated time‐lapse geophysical responses and well logs forms a multi‐scale, multi‐physics testbed for designing and testing geophysical CO 2 monitoring systems as well as for imaging and characterization algorithms. The suite of numerical models and data have been made publicly available for downloading on the National Energy Technology Laboratory's (NETL) Energy Data Exchange (EDX) website.

58 GEOSCIENCES↗
DOE OSTI2016

A multiscale hydro-geochemical-mechanical approach to analyze faulted CO 2 reservoirs: Original Research Article: A multiscale hydro-geochemical-mechanical approach to analyze faulted CO 2 reservoirs

This paper applies a multiscale hydro-geochemical-mechanical approach to analyze faulted CO 2 reservoirs using the STOMP-CO 2 -R code that is coupled to the ABAQUS® finite element package. STOMP-CO 2 -R models the reactive transport of CO 2 causing mineral volume fraction changes that are captured by an Eshelby-Mori-Tanka model implemented in ABAQUS®. A three-dimensional (3D) STOMP-CO 2 -R model for a reservoir containing an inclined fault was built to analyze a formation containing a reaction network with 5 minerals: albite, anorthite, calcite, kaolinite and quartz. A 3D finite element mesh that exactly maps the STOMP-CO 2 -R grid is developed for coupled hydro-geochemical-mechanical analyses. The model contains alternating sandstone and shale layers. The impact of reactive transport of CO 2 on the geomechanical properties of reservoir rocks and seals are studied in terms of mineral composition changes that affect their geomechanical responses. Simulations assuming extensional and compressional stress regimes with and without coupled geochemistry are performed to study the stress regime effect on the risk of hydraulic fracture. The tendency for the fault to slip is examined in terms of stress regime, geomechanical and geochemical-mechanical effects as well as fault inclination. The results show that mineralogical changes due to long-term injection of CO 2 reduce the permeability and elastic modulus of the reservoir, leading to increased risk of hydraulic fracture in the injection location and at the caprock seal immediately above the injection zone. Fault slip is not predicted to occur. However, fault inclination and stress regime have an important impact on the slip tendency factor.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

CO 2 -induced chemo-mechanical alteration in reservoir rocks assessed via batch reaction experiments and scratch testing

Here, the injection of carbon dioxide (CO 2 ) into geological formations results in a chemical re-equilibration between the mineral assemblage and the pore fluid, with ensuing mineral dissolution and re-precipitation. Hence, target rock formations may exhibit changes of mechanical and petrophysical properties due to CO 2 exposure. We conducted batch reaction experiments with Entrada Sandstone and Summerville Siltstone exposed to de-ionized water and synthetic brine under reservoir pressure (9–10 MPa) and temperature (80°C) for up to four weeks. Samples originate from the Crystal Geyser field site, where a naturally occurring CO 2 seepage alters portions of these geologic formations. We conducted micro-scratch tests on rock samples without alteration, altered under laboratory conditions, and naturally altered over geologic time. Scratch toughness and hardness decrease as a function of exposure time and water salinity up to 52% in the case of Entrada and 87% in the case of Summerville after CO 2 -induced alteration in the laboratory. Imaging of altered cores with SEM-EDS and X-ray microCT methods show dissolution of carbonate and silica cements and matrix accompanied by minor dissolution of Fe-oxides, clays, and other silicates. Parallel experiments using powdered samples confirm that dissolution of carbonate and silica are the primary reactions. The batch reaction experiments in the autoclave utilize a high fluid to rock volume ratio and represent an end member of possible alteration associated with CO 2 storage systems. These types of tests serve as a pre-screening tool to identify the susceptibility of rock facies to CO 2 -related chemical-mechanical alteration during long-term CO 2 storage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Data-driven nonlocal model for fragmentation in the crushing of solids

A technique is proposed for reproducing particle size distributions in three-dimensional simulations of the crushing and comminution of solid materials. The method is designed to produce realistic distributions over a wide range of loading conditions, especially for small fragments. In contrast to most existing methods, the new model does not explicitly treat the small-scale process of fracture. Instead, it uses measured fragment distributions from laboratory tests as the basic material property that is incorporated into the algorithm, providing a data-driven approach. The algorithm is implemented within a nonlocal peridynamic solver, which simulates the underlying continuum mechanics and contact interactions between fragments after they are formed. Finally, the technique is illustrated in reproducing fragmentation data from drop weight testing on sandstone samples.

58 GEOSCIENCES↗
DOE OSTI2016

Data-Driven, In Situ, Relative Sensor Calibration Based on Waveform Fitting Moment Tensor Inversion

This article presents a deconvolution-based methodology for the relative calibration of acoustic emission sensors. The method uses the observations from multiple acoustic emission arrivals to estimate correction functions that account for instrument response and site effects. The corrections act on the measurements reshaping the wave arrivals, in many cases allowing visible discrimination between compressional and shear phases. The methodology is illustrated with the application on a dataset obtained during the hydraulic fracturing at triaxial stress conditions of a Colton sandstone block. A workflow for the estimation of full moment tensor solutions, including the generation of quality control attributes complements the presentation. The moment tensor solutions display fracture geometries mostly aligned with the main hydraulic fracture’s plane. The style of activation can be associated with an extensive stress regime before breakdown and to a compressional stress regime after breakdown, when the pumps were reversed to withdraw fluid from the block. In conclusion, the source types commonly deviate from the pure double-couple type; however, the deviations are under 30° in the case of opening fractures, and under 20° in the case of closing fractures. Therefore, the acoustic emissions in this experiment are the result of predominantly shearing dislocations.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Laboratory Evaluation of the Thermal Breakout Method for Maximum Horizontal Stress Measurement

Measuring in situ stress is essential for many problems in geomechanics, and the maximum horizontal stress is the most difficult to constrain. We are developing an extension of the breakout method to measure maximum horizontal stress in regions where natural breakouts do not occur. In the novel thermal breakout method, additional compression which leads to breakout development is induced by heating the borehole wall. In the present study, we validated the method experimentally in a true-triaxial apparatus on samples with predrilled boreholes. Two rocks were selected for laboratory testing: high-porosity Berea sandstone and low-porosity Niagaran dolomite. Prior to main true-triaxial tests, we carried out standard testing to characterize the strength, elasticity and thermal properties. The true-triaxial experiments consisted of: (1) room-temperature tests where samples were first loaded mechanically until the breakout formed, and (2) elevated-temperature tests where samples were loaded mechanically within the elastic range with additional compression induced thermally. Breakout initiation was monitored by acoustic emission sensors mounted on the pistons that applied horizontal stresses. Furthermore, the magnitude of induced thermal stress was calculated from temperature measurements around the borehole wall. In both rock types, we created thermally induced breakouts and examined analytical expressions to constrain maximum horizontal stress based on strength, elastic and thermal properties of the rocks.

58 GEOSCIENCES↗
DOE OSTI2024

Review: Pre-Darcy flows in low-permeability porous media

The widely used Darcy’s law specifies a linear relation between the Darcy velocity of fluid flow and the pressure gradient that drives the flow. However, studies have shown that Darcy velocity can exhibit a nonlinear dependence on the pressure gradient in low-permeability porous media such as clay and shale when the pressure gradient is adequately low. This phenomenon is referred to as low-velocity non-Darcian flow or pre-Darcy flow. This paper provides a comprehensive review of the theories, experimental data, and modeling methods for pre-Darcy flow in low-permeability porous media. The review begins by outlining the fundamental mechanisms underlying pre-Darcy flow that regulate the unique characteristics such as nonlinear dependence of the Darcy velocity on the pressure gradient and its relevance to fluid–rock interactions. The review then proceeds to present a thorough compilation of experimental investigations performed in various low-permeability geomaterials including tight sandstones, shales, and clays. Next, empirical and theoretical models and simulation methods that have been developed to fit and interpret experimental data are reviewed. Finally, the review underscores the challenges encountered in conducting and interpreting pre-Darcy flow experiments and suggests future research directions. By analyzing previous experimental investigations, this review aims to offer a valuable resource for researchers and practitioners seeking to enhance their understanding of fluid dynamics in low-permeability geomaterials. This provides insights into the application of pre-Darcy flow in numerous natural and engineered processes such as shale oil and gas recovery, contaminant transport in low-permeability aquifers, and geological disposal of nuclear waste.

58 GEOSCIENCES↗
DOE OSTI2024

An integrated approach to derive relative permeability from capillary pressure

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability curve can be determined. Relative permeability is a numerical description of the interaction between two or more fluids and the porous media. It is a critical parameter for various tools that characterize subsurface multiphase flow systems, such as numerical simulation for carbon sequestration, oil and gas development, and groundwater contamination remediation. Therefore, it is critical to get a good statistical distribution of relative permeability in the porous media under study. Empirical formula for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility that is required to match laboratory-derived relative permeability curves. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability curves. Mercury intrusion capillary pressure (MICP) data from the Morrow 'B' Sandstone coupled with interfacial tension and contact angle measurements were used to create a suite of relative permeability curves. Furthermore, these curves were then calibrated to a small sample of existing laboratory curves to elucidate common fitting parameters for the formation that were then used to create relative permeability curves from MICP data that does not have an associated laboratory-measured relative permeability curve.

58 GEOSCIENCES↗