Search NASA⌕ Search

SEARCH · Search NASA

Results for “Relative permeability”

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 451 records · Page 25

Conservative solute transport processes and associated transient storage mechanisms: Comparing streams with contrasting channel morphologies, land use and land cover

Abstract Land use within a watershed impacts stream channel morphology and hydrology and, therefore, in‐stream solute transport processes and associated transient storage mechanisms. This study evaluated transport processes in two contrasting stream sites where channel morphology was influenced by the surrounding land use, land cover, climate and geologic controls: Como Creek, CO, a relatively undisturbed, high gradient, forested stream with a gravel bed and complex channel morphology, and Clear Creek, IA, an incised, low‐gradient stream with low‐permeability substrate draining an agricultural landscape. We performed conservative stream tracer injections at these sites to address the following questions: (1) How does solute transport vary between streams with differing morphologies? and (2) How does solute transport at each stream site change as a function of discharge? We analysed in‐stream tracer time series data and compared results quantifying solute attenuation in surface and subsurface transient storage zones. Significant trends were observed in these metrics with varying discharge conditions at the forested site but not at the agricultural site. There was a broad range of transport mechanisms and evidence of substantial exchange with both surface and hyporheic transient storage in the relatively undisturbed, forested stream. Changing discharge conditions activated or deactivated different solute transport mechanisms in the forested site and greatly impacted advective travel time. Conversely in the simplified agricultural stream, there was a narrow range of solute transport behaviour across flows and predominantly surface transient storage at all measured discharge conditions. These results demonstrate how channel simplification inhibits available solute transport mechanisms across varying discharge conditions.

54 ENVIRONMENTAL SCIENCES↗

Effects of nonionic surfactants on pigment excretion and cell morphology in extractive fermentation of Monascus sp. NJ1

Abstract BACKGROUND Different nonionic surfactants in submerged fermentation of Monascus sp. demonstrate significant differences regarding increasing pigment yield. In this study, 15 surfactants from five series were analyzed to investigate the influence of nonionic surfactants on Monascus pigments, with the aim of simultaneously obtaining a novel nonionic surfactant. RESULTS Addition of the novel surfactant Brij 35 greatly enhanced pigment excretion and demonstrated good biocompatibility. Extracellular red, orange and yellow pigments increased by 1.47‐, 1.71‐ and 2.07‐fold respectively. Production of extracellular pigments was not only related to the hydrophile–lipophile balance value (HLB) but also affected by the cloud point temperature (CP) of the fermentation medium. It was found that nonionic surfactants can improve cell membrane permeability and cell storage capacity by modifying the cell walls of Monascus mycelium and by increasing lipid droplet levels, enhancing pigment excretion. Different nonionic surfactants modify Monascus mycelium to different degrees. CONCLUSION The novel surfactant Brij 35, which has good capacity for pigment extraction and biocompatibility, was identified in the analysis. The effects of nonionic surfactants on the secretion of pigments are related to not only the modification of the cell wall and internal structure but also the CP and HLB. © 2018 Society of Chemical Industry

Yang, Xuelian↗

Thermal regimes in the detachment fault environment as deduced from fluid inclusions

Extensional tectonism, which dominates middle- and late-Tertiary geology in western Arizona, southeastern California, and southern Nevada, is characterized by normal regionally extensive, low-angle detachment faults. The decollement movement of Fupper plate rocks relative to lower plate assemblages created extensive zones of dilatency, including synthetic and antithetic listric normal faults, tear faults, tectonic crush breccias, shatter breccias, and gash veins in lithologic units above and below the detachment. The tectonically enhanced permeability above and below the detachment fault permitted mass migration of large volumes of hydrothermal solutions along the fault zone during and following upper plate movement. Major quantities of MgO, CaO, K2O, FeO/Fe2O3, SiO2 and CO2 were added to rocks in and near the detachment and related structures. Also introduced were varying amounts of trace elements including Mn, Cu, S, Mo, Ba, Au, Pb, Zn, U and/or Ag. Minerals containing fluid incusions were collected from all of these loci at locations in detachment faulted terranes in western Arizona and southeastern California.

Beane, R. E.↗

Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO 2 in fractured rocks

Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheres and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO 2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for development of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.

54 ENVIRONMENTAL SCIENCES↗

Fundamental Structure/Property Studies of Gas Separation Membrane Polymers

Solution-processable polybenzimidazoles (PBIs) present opportunities for new gas separation membranes. We plan to develop systematic structure/property relations for this relatively new class of membrane materials. The fundamental questions we will address are: (1) How do chemical structure modifications of PBIs influence their gas solubility, diffusivity and permeability properties for a systematic series of gases? (2) Is the separation performance of PBIs affected by plasticization under operating conditions that are relevant to, for example, natural gas purification or olefin/paraffin separation and, if so, can this family of materials be tuned structurally to minimize deleterious effects of plasticization on separation properties? (3) What is the effect of water vapor on PBI separation performance? (4) What is the influence of physical aging on permeation properties of thin (i.e., <1 µm) PBI polymer films, and how does the rate of physical aging relate to polymer chemical structure? The primary focus of this research program is to develop a fundamental, molecular-based mechanistic understanding of the influence of chemical structure on the gas transport properties of PBIs and to establish structure/property relationships for this unique platform of materials.

08 HYDROGEN↗

The timescales of mineral redox reactions

Redox-reactive minerals can serve as electron donors or acceptors for abiotic reactants or microbial metabolic processes, and hence can play important roles in terrestrial and aquatic environments, particularly if their reaction rates are comparable to those of other biogeochemical processes. Under such circumstances, their reactions can control metal and contaminant bioavailability and change the permeability of soils and sediments. While the thermodynamic driving force for a mineral redox reaction is frequently a good predictor of relative rates of reaction, there are many examples in which kinetic factors limit reaction rate. Understanding the thermodynamic and kinetic aspects of mineral reaction rates, and their sensitivity to environmental conditions such as temperature or pH, is important for anticipating the biogeochemical evolution of natural environments subjected to change. Achieving this goal requires knowledge of the reaction pathway, the timescales of intermediate steps, and the lifetimes of metastable reaction states. Mineral redox reactions proceed through a combination of steps that can include electron and proton transfer, the breaking or formation of bonds, and mineral dissolution or phase transformation. The combination of conventional kinetics approaches, newly-developed ultrafast time-resolved methods and molecular simulation can provide elucidation of such complex reaction pathways. This chapter summarizes key concepts in mineral and interfacial redox reactions with illustrations from the rich geochemistry of iron and ironbearing minerals.

58 GEOSCIENCES↗

Creating Novel Gas Separation Constructs by Manipulating Free-Volume Distributions in Polymer-Grafted Nanoparticles

Polymer membranes are critical to sustainability applications, and in this work we focused on one key application, the efficient separation of gas mixtures. Despite their widespread use, important challenges remain in the use of polymer membranes in this context – i.e., finding membrane materials with selectively improved transport of the desired mixture component(s) while possessing enhanced strength and improved aging behavior relative to the best current benchmarks. The important separation figures of merit are the gas flux, which is proportional to the permeability, P i = D i × S i (Di is the gas diffusivity and S i its solubility coefficient) and selectivity (i.e., gas purity,α ij = P i /P j ). The goal is to simultaneously maximize P i and α ij . Most research to date has empirically targeted the development of new glassy polymers with improved separation performance. These include thermally rearranged (TR) polymers and polymers of intrinsic microporosity (PIMs).

36 MATERIALS SCIENCE↗

Hydrologic Connectivity and Patch‐To‐Hillslope Scale Relations in Dryland Ecosystems

Abstract In drylands, runoff during storms redistributes water and nutrients from bare soil areas to vegetated patches, subsidizing vegetation with additional resources. The extent of this redistribution depends on the interplay between surface roughness and permeability; greater permeability in vegetated patches promotes run‐on to vegetation, but greater surface roughness diverts runoff, producing tortuous flow paths that bypass vegetation. Here, this interplay is examined in virtual experiments using the 2D Saint Venant Equations to measure runoff connectivity. Flowpaths are delineated using tracers advected by the flow. Distances between tracer sources and sinks along flowpaths measure hydrologic connectivity at two lengthscales: connectivity to the hillslope outlet and within‐slope source‐sink connectivity. Differences between these connectivity lengthscales indicate how flow may “by‐pass” vegetated patches within hillslopes. At the hillslope scale, a derived power‐law relation between the runoff coefficient and outlet connectivity describes hillslope water losses, providing a foundation for identifying landscapes likely to shed water.

54 ENVIRONMENTAL SCIENCES↗

Quantifying Groundwater Response and Uncertainty in Beaver‐Influenced Mountainous Floodplains Using Machine Learning‐Based Model Calibration

Abstract Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model‐data integration workflow to address uncertainty in groundwater response to beaver‐induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning‐based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down‐valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down‐valley underflow, from 2 during dry pond periods to 20 during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down‐valley underflow, minimizes the effect of beaver‐induced inundation on water quality. We emphasize the need for field‐scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response. Plain Language Summary Beavers change the flow of water in river corridors by creating ponds, expanding wetlands, and flooding floodplains. This increases surface water area, promotes plant growth, and enhances biodiversity. However, the impact of this flooding on groundwater flow is not well understood, especially in mountainous areas with gravel layers where water moves easily beneath soil. In this study, we used numerical modeling to investigate how beaver ponds influence groundwater in a mountainous floodplain of the Upper Colorado River Basin. We adapted a machine learning method to validate our numerical models using multiple field data sets. Our findings show that beaver ponds significantly increase vertical water flow from the soil to the gravel during wet periods, compared to when the ponds are fully drained. The study also highlights the importance of floodplain structure in controlling both water flow in gravel layers along the river direction and vertical flow from the soil to the gravel with the presence of beavers. To reduce uncertainty in groundwater response, we emphasize the need for more field‐scale measurements of floodplain structure, hydraulic properties, and evapotranspiration changes. Key Points Floodplain structures and hydraulic conductivities are important for groundwater response with beaver ponds in mountainous floodplains Large down‐valley underflow in permeability‐stratified floodplains reduces beaver‐induced impacts on groundwater storage and water quality Machine learning‐based model calibration methods are effective for estimating posterior distributions of groundwater model parameters

Wang, Lijing↗

An Integrated Feasibility Study of Reservoir Thermal Energy Storage (RTES) in Portland, OR, USA

In regions with long cold overcast winters and sunny summers, Deep Direct-Use (DDU) can be coupled with Reservoir Thermal Energy Storage (RTES) technology to take advantage of pre-existing subsurface permeability to save summer heat for later use during cold seasons. Many aquifers worldwide are underlain by permeable regions (reservoirs) containing brackish or saline groundwater that has limited beneficial use due to poor water quality. We investigate the utility of these relatively deep, slow flowing reservoirs for RTES by conducting an integrated feasibility study in the Portland Basin, Oregon, USA, developing methods and obtaining results that can be widely applied to RTES systems elsewhere. As a case study, we have conducted an economic and social cost-benefit analysis for the Oregon Health and Science University (OHSU), a teaching hospital that is recognized as critical infrastructure in the Portland Metropolitan Area. Our investigation covers key factors that influence feasibility including 1) the geologic framework, 2) heat and fluid flow modeling, 3) capital and maintenance costs, 4) the regulatory framework, and 5) operational risks. By pairing a model of building seasonal heat demand with an integrated model of RTES resource supply, we determine that the most important factors that influence RTES efficacy in the study area are operational schedule, well spacing, the amount of summer heat stored (in our model, a function of solar array size), and longevity of the system. Generally, heat recovery efficiency increases as the reservoir and surrounding rocks warm, making RTES more economical with time. Selecting a base-case scenario, we estimate a levelized cost of heat (LCOH) to compare with other sources of heating available to OHSU and find that it is comparable to unsubsidized solar and nuclear, but more expensive than natural gas. Additional benefits of RTES include energy resiliency in the event that conventional energy supplies are disrupted (e.g., natural disaster) and a reduction in fossil fuel consumption resulting in a smaller carbon footprint. Key risks include reservoir heterogeneity and a possible reduction in permeability through time due to scaling (mineral precipitation). Lastly, a map of thermal energy storage capacity for the Portland Basin yields a total of 43,400 GWh, suggesting tremendous potential for RTES in the Portland Metropolitan Area.

14 SOLAR ENERGY↗

Preface

High-temperature geothermal resources are often the manifestations of cooling magmatic intrusions, and are vertically zoned. Beneath the low-permeability cap and the hydrothermal circulation volume, there exists a zone of limited permeability where water does not circulate freely. At present, almost all of the commercial production of geothermal energy is derived from the relatively shallow (usually < 2–3 km) hydrothermal zone (i.e. zone where water can flow freely). Reliance on hydrothermal resources is dictated by currently available technology and by cost considerations.

15 GEOTHERMAL ENERGY↗

Evolution and Transport of Water in the Upper Regolith of Mars

Long standing theoretical predictions [1-3], as well as recent spacecraft observations [4] indicate that large quantities of ice is present in the high latitudes upper decimeters to meters of the Martian regolith. At shallower depths and warmer locations small amounts of H2O, either adsorbed or free, may be present transiently. An understanding of the evolution of water based on theoretical and experimental considerations of the processes operating at the Martian environment is required. In particular, the porosity, diffusivity, and permeability of soils and their effect on water vapor transport under Mars-like conditions have been estimated, but experimental validation of such models is lacking. Goal: Three related mechanisms may affect water transport in the upper Martian regolith. 1) diffusion along a concentration gradient under isobaric conditions, 2) diffusion along a thermal gradient, which may give rise to a concentration gradient as ice sublimes or molecules desorb from the regolith, and 3) hydraulic flow, or mass motion in response to a pressure gradient. Our combined theoretical and experimental investigation seeks to disentangle these mechanisms and determine which process(es) are dominant in the upper regolith over various timescales. A detailed one-dimensional model of the upper regolith is being created which incorporates water adsorption/ desorption, condensation, porosity, diffusivity, and permeability effects. Certain factors such as diffusivity are difficult to determine theoretically due to the wide range of intrinsic grain properties such as particle sizes, shapes, packing densities, and emergent properties such as tortuosity. An experiment is being designed which will allow us to more accurately determine diffusivity, permeability, and water desorption isotherms for regolith simulants.

Hudson, T. L.↗

Highly CO 2 -permeable membranes derived from a midblock-sulfonated multiblock polymer after submersion in water

To mitigate the effect of atmospheric CO 2 on global climate change, gas separation materials that simultaneously exhibit high CO 2 permeability and selectivity in gas mixtures must be developed. In this study, CO 2 transport through midblock-sulfonated block polymer membranes prepared from four different solvents is investigated. The results presented here establish that membrane morphology and accompanying gas transport properties are sensitive to casting solvent and relative humidity. We likewise report an intriguing observation: submersion of these thermoplastic elastomeric membranes in liquid water, followed by drying prior to analysis, promotes not only a substantial change in membrane morphology, but also a significant improvement in both CO 2 permeability and CO 2 /N 2 selectivity. Measured CO 2 permeability and CO 2 /N 2 selectivity values of 482 Barrer and 57, respectively, surpass the Robeson upper bound, indicating that these nanostructured membranes constitute promising candidates for gas separation technologies aimed at CO 2 capture.

36 MATERIALS SCIENCE↗

Numerical Simulations in Support of a Long-Term Test of Gas Production From Hydrate Accumulations on the Alaska North Slope: Water Production and Associated Design and Management Issues

Here, we investigated numerical simulation strategies for a long-term test of depressurization-induced gas production from the B1 Sand of Unit B at the Hydrate-01 Stratigraphic Test Well. The main objective of this study was to estimate fluid production rates (with emphasis on water production) under a variety of conditions and production scenarios and contribute new insights to the design and management of the field test. In the first part of the study, we investigated the system response to a three-step depressurization process using two limiting sets of flow properties─the expected maximum and minimum intrinsic and effective permeabilities─for the very heterogeneous reservoir. In the second part, we investigated the effect of the production interval length and placement within the formation relative to the boundaries of the hydrate-bearing unit. The best performing well configuration was used in the third part of the study, which used the most representative subsurface flow properties to investigate the effect of the depressurization strategy on the production performance. The best overall performance (largest gas production with modest water production and a strong response at the observation wells) was obtained with a 10 m-long well situated 3 m below the top of the formation and a three-step depressurization scheme at 15-day intervals to a terminal bottomhole pressure of 2.8 MPa. The overall production performance was enhanced by a faster rate of depressurization. Estimated water production rates in all cases were limited and easily manageable. None of the tested well configurations or depressurization strategies significantly reduced water production without also severely reducing gas production. In all the investigated cases, 95% of the long-term fraction of produced water was replenished by inflows from the boundaries and could not be reduced. These substantial water inflows are an unavoidable feature of HU-B and cannot be easily mitigated by a hydraulic control.

02 PETROLEUM↗

Core-flood Effluent and Shale Surface Chemistries in Predicting Interaction between Shale, Brine, and Reactive Fluid

We report field and laboratory observations to date indicate that the efficiency of hydraulic fracturing, as it relates to hydrocarbon recovery, depends significantly on geochemical alterations to rock surfaces that diminish accessibility by partial or total plugging of the pore and fracture networks. This is caused by mineral scale deposition such as coating of fracture surfaces with precipitates, particle migration, and/or crack closure due to dissolution under stress. In reactive flow-through experiments, mineral reactions in response to acidic fluid injection significantly reduced system porosity and core permeability. The present study focuses on changes to fluid chemistry and shale surfaces (inlet and fracture walls) resulting from shale-fluid interactions and integrating these findings for an improved estimate of transport-related consequences. The reacted shale surfaces were examined by spatially-resolved scanning electron microscopy - energy dispersive spectroscopy (SEM-EDS) analysis. Importantly, inductively coupled plasma - mass spectrometry/optical emission spectroscopy (ICP-MS/OES) was utilized to probe the chemical evolution of the core-flood effluents. The three study cores selected from the Marcellus formation represent different mineralogies and structural features. In flow-through experiments, lab-generated brine and HCl-based fracture fluid (pH=2) were injected sequentially under effective stress (up to 500 psi) at reservoir temperature (80°C). SEM-EDS results confirmed by the ICP concentration trends showed significant Fe hydroxide precipitates in clay- and pyrite-rich outcrop samples due to partial oxidation of Fe-bearing phases in the case of intrusion of low salinity water-based fluids. Porosity reduction in the MSEEL (Marcellus Shale Energy and Environmental Laboratory) carbonate-rich sample is related to compaction of cores under stress due to matrix softening with substantial dissolution, and pore-filling by hydroxides, as well as barite and salts. Despite the same fluid compositions and experimental conditions used for both MSEEL samples, barite precipitation was much more intense in the MSEEL clay-rich sample due to its greater sorption capacity and additional sulfate source as well as fissile nature with multiple lengthwise cracks. ICP tests revealed time-resolved concentration trends in produced brine and reactive fluids that in turn complemented the pre-/post-reaction SEM-EDS observations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SYNTHESIS OF NOVEL ALL-DIELECTRIC GRATING FILTERS USING GENETIC ALGORITHMS

We are concerned with the design of inhomogeneous, all dielectric (lossless) periodic structures which act as filters. Dielectric filters made as stacks of inhomogeneous gratings and layers of materials are being used in optical technology, but are not common at microwave frequencies. The problem is then finding the periodic cell's geometric configuration and permittivity values which correspond to a specified reflectivity/transmittivity response as a function of frequency/illumination angle. This type of design can be thought of as an inverse-source problem, since it entails finding a distribution of sources which produce fields (or quantities derived from them) of given characteristics. Electromagnetic sources (electric and magnetic current densities) in a volume are related to the outside fields by a well known linear integral equation. Additionally, the sources are related to the fields inside the volume by a constitutive equation, involving the material properties. Then, the relationship linking the fields outside the source region to those inside is non-linear, in terms of material properties such as permittivity, permeability and conductivity. The solution of the non-linear inverse problem is cast here as a combination of two linear steps, by explicitly introducing the electromagnetic sources in the computational volume as a set of unknowns in addition to the material unknowns. This allows to solve for material parameters and related electric fields in the source volume which are consistent with Maxwell's equations. Solutions are obtained iteratively by decoupling the two steps. First, we invert for the permittivity only in the minimization of a cost function and second, given the materials, we find the corresponding electric fields through direct solution of the integral equation in the source volume. The sources thus computed are used to generate the far fields and the synthesized triter response. The cost function is obtained by calculating the deviation between the synthesized value of reflectivity/transmittivity and the desired one. Solution geometries for the periodic cell are sought as gratings (ensembles of columns of different heights and widths), or combinations of homogeneous layers of different dielectric materials and gratings. Hence the explicit unknowns of the inversion step are the material permittivities and the relative boundaries separating homogeneous parcels of the periodic cell.

Zuffada, Cinzia↗

Superhydrophobic composite asymmetric electrospun membrane for sustainable vacuum assisted air gap membrane distillation

Membrane distillation (MD) relies on superhydrophobic membranes that must achieve the competing goals of high liquid entry pressure (LEP) and high permeability. Here, in this study, we developed the first superhydrophobic composite electrospun asymmetric membrane using fluorinated graphite (FG). We also provided among the first comparisons of the relative benefits of membrane asymmetry and composite coatings via 3D electrospinning and spraying. FG particles increased the roughness and reduced the zeta potential of the membranes. Superhydrophobic composite membranes had a contact angle of 157° (sliding angle of 7°) and an LEP of 327 kPa. Long-term MD testing revealed that FG particles could increase MD's sustainability while decreasing membrane wetting. Vacuum-assisted air gap membrane distillation (V-AGMD) showed that the composite asymmetric membrane outperformed the non-asymmetric version (composite nanofiber) membrane by 21 % in average water flux and 25 % in average thermal efficiency because of the asymmetric structure of the membrane. The numerical model developed for the V-AGMD system also successfully predicted the water flux generated experimentally from the composite nanofiber and composite asymmetric membranes. Our results confirmed that electrospun asymmetric FG membranes with high contact angle, low sliding angle, high permeability, and high gained output ratio (GOR) would be promising candidates for sustainable MD.

36 MATERIALS SCIENCE↗

India National Gas Hydrate Program Expedition 02 Summary of Scientific Results: Gas Hydrate Systems Along the Eastern Continental Margin of India

The India National Gas Hydrate Program Expedition 02 (NGHP-02) discovered gas hydrate at high saturation in sand reservoirs at several sites in the deepwater Bay of Bengal. To assess the potential response of those deposits to scientific depressurization experiments, comprehensive geologic models were constructed to enable numerical simulation for two sites. Both sites (NGHP-02-09 and NGHP-02-16) feature thick sequences of thinly-interbedded reservoir and non-reservoir facies at sub-seafloor depths less than 300m and sub-sea depths of 2400m or more. These settings pose significant challenges to current modeling capabilities. First, the thinly-interbedded reservoir architecture complicates the determination of basic reservoir parameters from both log and core data due to measurement resolution issues. Secondly, the fine scale variation in sediment properties imparts great contrasts in key parameters over very short distances, creating high gradients at multiple scales and varying orientations that necessitate careful design of high-definition simulation grids. Thirdly, the deposits include internal sources of water, as well as a range of complex boundary conditions, including variable permeability within the overlying mud-rich “seals,” that complicate reservoir depressurization. Lastly, because of the unique combination of great water depth and relatively shallow sub-seafloor depth: models designed to maximize the dissociation rate impose large pressure drawdowns on relatively low-strength sediments. This condition renders the proper evaluation and integration of the geomechanical response to hydrate dissociation critical. In this report, we review the history of gas hydrate reservoir simulation, discuss methods for creating geologic input models, and summarize the key findings and implications of the collaborative NGHP-02 numerical simulation effort. Together, the studies confirm the viability of the modeled accumulations for scientific testing and identify key challenges related to the selection of specific test sites and the design of test wells.

Collett, Tim↗