Understanding Potassium Isotope Fractionation During Authigenic Clay Formation in Pore-fluid Systems: Implications for the δ41K of Seawater

NASA Astrophysics Data System (ADS)

Santiago Ramos, D. P.; Higgins, J. A.

2015-12-01

Improvements in analytical precision on the latest generation multi-collector inductively coupled plasma mass spectrometers (MC-ICP-MS) have revealed a ~2‰ range in the ratios of stable potassium isotopes (41K/39K) in terrestrial materials (Morgan et al., in prep). Preliminary measurements of δ41K values indicate that seawater and silicate rocks are isotopically distinct reservoirs, with seawater having a δ41K value that is ~0.5‰ heavier than the silicate average (-0.5‰; Morgan et al., in prep). The heavy δ41K character of seawater might be related to 1) an isotopically enriched input flux (rivers and high-temperature hydrothermal reactions); or 2) a 41K-depleted sink associated with authigenic clay formation during low-temperature alteration of volcanic rocks. Here we present measurements of the δ41K values of pore-fluids from ODP site 1052 in order to constrain potassium isotope fractionation during secondary clay formation. We find that δ41K values and K concentrations both decline systematically with depth. Results from 1-D diffusion-advection-reaction modeling of potassium concentrations and isotopic compositions indicate that fractionation of K isotopes during diffusion (Bourg et al., 2010) can explain all of the change in δ41K values of the pore-fluid with depth. Although the size of the K sink at site 1052 is a trivial fraction of the global K sink in clay minerals, our results suggest that diffusive fractionation of K isotopes in shallow pore-fluids may be, in part, responsible for the elevated δ41K value of seawater.

Scaling of postinjection-induced seismicity: An approach to assess hydraulic fracturing related processes

NASA Astrophysics Data System (ADS)

Johann, Lisa; Dinske, Carsten; Shapiro, Serge

2017-04-01

Fluid injections into unconventional reservoirs have become a standard for the enhancement of fluid-mobility parameters. Microseismic activity during and after the injection can be frequently directly associated with subsurface fluid injections. Previous studies demonstrate that postinjection-induced seismicity has two important characteristics: On the one hand, the triggering front, which corresponds to early and distant events and envelops farthest induced events. On the other hand, the back front, which describes the lower boundary of the seismic cloud and envelops the aseismic domain evolving around the source after the injection stop. A lot of research has been conducted in recent years to understand seismicity-related processes. For this work, we follow the assumption that the diffusion of pore-fluid pressure is the dominant triggering mechanism. Based on Terzaghi's concept of an effective normal stress, the injection of fluids leads to increasing pressures which in turn reduce the effective normal stress and lead to sliding along pre-existing critically stressed and favourably oriented fractures and cracks. However, in many situations, spatio-temporal signatures of induced events are captured by a rather non-linear process of pore-fluid pressure diffusion, where the hydraulic diffusivity becomes pressure-dependent. This is for example the case during hydraulic fracturing where hydraulic transport properties are significantly enhanced. For a better understanding of processes related to postinjection-induced seismicity, we analytically describe the temporal behaviour of triggering and back fronts. We introduce a scaling law which shows that postinjection-induced events are sensitive to the degree of non-linearity and to the Euclidean dimension of the seismic cloud (see Johann et al., 2016, JGR). To validate the theory, we implement comprehensive modelling of non-linear pore-fluid pressure diffusion in 3D. We solve numerically for the non-linear equation of diffusion with a power-law dependent hydraulic diffusivity on pressure and generate catalogues of synthetic seismicity. We study spatio-temporal features of the seismic clouds and compare the results to theoretical values predicted by the novel scaling law. Subsequently, we apply the scaling relation to real hydraulic fracturing and Enhanced Geothermal System data. Our results show that the derived scaling relations well describe synthetic and real data. Thus, the methodology can be used to obtain hydraulic reservoir properties and can contribute significantly to a general understanding of injection related processes as well as to hazard assessment.

Susceptibility of experimental faults to pore pressure increase: insights from load-controlled experiments on calcite-bearing rocks

NASA Astrophysics Data System (ADS)

Spagnuolo, Elena; Violay, Marie; Nielsen, Stefan; Cornelio, Chiara; Di Toro, Giulio

2017-04-01

Fluid pressure has been indicated as a major factor controlling natural (e.g., L'Aquila, Italy, 2009 Mw 6.3) and induced seismicity (e.g., Wilzetta, Oklahoma, 2011 Mw 5.7). Terzaghi's principle states that the effective normal stress is linearly reduced by a pore pressure (Pf) increase σeff=σn(1 - αPf), where the effective stress parameter α, may be related to the fraction of the fault area that is flooded. A value of α =1 is often used by default, with Pf shifting the Mohr circle towards lower normal effective stresses and anticipating failure on pre-existing faults. However, within a complex fault core of inhomogeneous permeability, α may vary in a yet poorly understood way. To shed light on this problem, we conducted experiments on calcite-bearing rock samples (Carrara marble) at room humidity conditions and in the presence of pore fluids (drained conditions) using a rotary apparatus (SHIVA). A pre-cut fault is loaded by constant shear stress τ under constant normal stress σn=15 MPa until a target value corresponding roughly to the 80 % of the frictional fault strength. The pore pressure Pf is then raised with regular pressure and time steps to induce fault instability. Assuming α=1 and a threshold for instability τp_eff=μp σeff, the experiments reveal that an increase of Pf does not necessarily induce an instability even when the effective strength threshold is largely surpassed (e.g., τp_eff=1.3 μpσeff). This result may indicate that the Pf increase did not instantly diffuse throughout the slip zone, but took a finite time to equilibrate with the external imposed pressure increase due to finite permeability. Under our experimental conditions, a significant departure from α=1 is observed provided that the Pf step is shorter than about < 20s. We interpret this delay as indicative of the diffusion time (td), which is related to fluid penetration length l by l = √ κtd-, where κ is the hydraulic diffusivity on the fault plane. We show that a simple cubic law relates td to hydraulic aperture, pore pressure gradient and injection rate. We redefine α as the ratio between the fluid penetration length and sample dimension L resulting in α = min(√ktd,L) L. Under several pore pressure loading rates this relation yields an approximate hydraulic diffusivity κ ˜10-8 m2 s-1 which is compatible, for example, with a low porosity shale. Our results highlight that a high injection flow rate in fault plane do not necessarily induce seismogenic fault slip: a critical pore penetration length or fluid patch size is necessary to trigger fault instability.

Pore-scale modeling of phase change in porous media

NASA Astrophysics Data System (ADS)

Juanes, Ruben; Cueto-Felgueroso, Luis; Fu, Xiaojing

2017-11-01

One of the main open challenges in pore-scale modeling is the direct simulation of flows involving multicomponent mixtures with complex phase behavior. Reservoir fluid mixtures are often described through cubic equations of state, which makes diffuse interface, or phase field theories, particularly appealing as a modeling framework. What is still unclear is whether equation-of-state-driven diffuse-interface models can adequately describe processes where surface tension and wetting phenomena play an important role. Here we present a diffuse interface model of single-component, two-phase flow (a van der Waals fluid) in a porous medium under different wetting conditions. We propose a simplified Darcy-Korteweg model that is appropriate to describe flow in a Hele-Shaw cell or a micromodel, with a gap-averaged velocity. We study the ability of the diffuse-interface model to capture capillary pressure and the dynamics of vaporization/condensation fronts, and show that the model reproduces pressure fluctuations that emerge from abrupt interface displacements (Haines jumps) and from the break-up of wetting films.

A kinetic Monte Carlo approach to study fluid transport in pore networks

NASA Astrophysics Data System (ADS)

Apostolopoulou, M.; Day, R.; Hull, R.; Stamatakis, M.; Striolo, A.

2017-10-01

The mechanism of fluid migration in porous networks continues to attract great interest. Darcy's law (phenomenological continuum theory), which is often used to describe macroscopically fluid flow through a porous material, is thought to fail in nano-channels. Transport through heterogeneous and anisotropic systems, characterized by a broad distribution of pores, occurs via a contribution of different transport mechanisms, all of which need to be accounted for. The situation is likely more complicated when immiscible fluid mixtures are present. To generalize the study of fluid transport through a porous network, we developed a stochastic kinetic Monte Carlo (KMC) model. In our lattice model, the pore network is represented as a set of connected finite volumes (voxels), and transport is simulated as a random walk of molecules, which "hop" from voxel to voxel. We simulated fluid transport along an effectively 1D pore and we compared the results to those expected by solving analytically the diffusion equation. The KMC model was then implemented to quantify the transport of methane through hydrated micropores, in which case atomistic molecular dynamic simulation results were reproduced. The model was then used to study flow through pore networks, where it was able to quantify the effect of the pore length and the effect of the network's connectivity. The results are consistent with experiments but also provide additional physical insights. Extension of the model will be useful to better understand fluid transport in shale rocks.

A model of oscillatory transport in granular soils, with application to barometric pumping and earth tides.

PubMed

Neeper, D A

2001-04-01

A simple algebraic model is proposed to estimate the transport of a volatile or soluble chemical caused by oscillatory flow of fluid in a porous medium. The model is applied to the barometric pumping of vapors in the vadose zone, and to the transport of dissolved species by earth tides in an aquifer. In the model, the fluid moves sinusoidally with time in the porosity of the soil. The chemical concentration in the mobile fluid is considered to equilibrate with the concentration in the surrounding matrix according to a characteristic time governed by diffusion, sorption, or other rate processes. The model provides a closed form solution, to which barometric pressure data are applied in an example of pore gas motion in the vadose zone. The model predicts that the additional diffusivity due barometric pumping in an unfractured vadose zone would be comparable to the diffusivity in stagnant pore gas if the equilibration time is 1 day or longer. Water motion due to the M2 lunar tide is examined as an example of oscillatory transport in an aquifer. It is shown that the tidal motion of the water in an aquifer might significantly increase the vertical diffusivity of dissolved species when compared to diffusion in an absolutely stagnant aquifer, but the hydrodynamic dispersivity due to tidal motion or gravitational flow would probably exceed the diffusivity due to oscillatory advection.

Nucleation speed limit on remote fluid induced earthquakes

USGS Publications Warehouse

Parsons, Thomas E.; Akinci, Aybige; Malignini, Luca

2017-01-01

Earthquakes triggered by other remote seismic events are explained as a response to long-traveling seismic waves that temporarily stress the crust. However, delays of hours or days after seismic waves pass through are reported by several studies, which are difficult to reconcile with the transient stresses imparted by seismic waves. We show that these delays are proportional to magnitude and that nucleation times are best fit to a fluid diffusion process if the governing rupture process involves unlocking a magnitude-dependent critical nucleation zone. It is well established that distant earthquakes can strongly affect the pressure and distribution of crustal pore fluids. Earth’s crust contains hydraulically isolated, pressurized compartments in which fluids are contained within low-permeability walls. We know that strong shaking induced by seismic waves from large earthquakes can change the permeability of rocks. Thus, the boundary of a pressurized compartment may see its permeability rise. Previously confined, overpressurized pore fluids may then diffuse away, infiltrate faults, decrease their strength, and induce earthquakes. Magnitude-dependent delays and critical nucleation zone conclusions can also be applied to human-induced earthquakes.

Nucleation speed limit on remote fluid-induced earthquakes.

PubMed

Parsons, Tom; Malagnini, Luca; Akinci, Aybige

2017-08-01

Earthquakes triggered by other remote seismic events are explained as a response to long-traveling seismic waves that temporarily stress the crust. However, delays of hours or days after seismic waves pass through are reported by several studies, which are difficult to reconcile with the transient stresses imparted by seismic waves. We show that these delays are proportional to magnitude and that nucleation times are best fit to a fluid diffusion process if the governing rupture process involves unlocking a magnitude-dependent critical nucleation zone. It is well established that distant earthquakes can strongly affect the pressure and distribution of crustal pore fluids. Earth's crust contains hydraulically isolated, pressurized compartments in which fluids are contained within low-permeability walls. We know that strong shaking induced by seismic waves from large earthquakes can change the permeability of rocks. Thus, the boundary of a pressurized compartment may see its permeability rise. Previously confined, overpressurized pore fluids may then diffuse away, infiltrate faults, decrease their strength, and induce earthquakes. Magnitude-dependent delays and critical nucleation zone conclusions can also be applied to human-induced earthquakes.

Nucleation speed limit on remote fluid-induced earthquakes

PubMed Central

Parsons, Tom; Malagnini, Luca; Akinci, Aybige

2017-01-01

Earthquakes triggered by other remote seismic events are explained as a response to long-traveling seismic waves that temporarily stress the crust. However, delays of hours or days after seismic waves pass through are reported by several studies, which are difficult to reconcile with the transient stresses imparted by seismic waves. We show that these delays are proportional to magnitude and that nucleation times are best fit to a fluid diffusion process if the governing rupture process involves unlocking a magnitude-dependent critical nucleation zone. It is well established that distant earthquakes can strongly affect the pressure and distribution of crustal pore fluids. Earth’s crust contains hydraulically isolated, pressurized compartments in which fluids are contained within low-permeability walls. We know that strong shaking induced by seismic waves from large earthquakes can change the permeability of rocks. Thus, the boundary of a pressurized compartment may see its permeability rise. Previously confined, overpressurized pore fluids may then diffuse away, infiltrate faults, decrease their strength, and induce earthquakes. Magnitude-dependent delays and critical nucleation zone conclusions can also be applied to human-induced earthquakes. PMID:28845448

One-Dimension Diffusion Preparation of Concentration-Gradient Fe₂O₃/SiO₂ Aerogel.

PubMed

Zhang, Ting; Wang, Haoran; Zhou, Bin; Ji, Xiujie; Wang, Hongqiang; Du, Ai

2018-06-21

Concentration-gradient Fe₂O₃/SiO₂ aerogels were prepared by placing an MTMS (methyltrimethoxysilane)-derived SiO₂ aerogel on an iron gauze with an HCl atmosphere via one-dimensional diffusion, ammonia-atmosphere fixing, supercritical fluid drying and thermal treatment. The energy dispersive spectra show that the Fe/Si molar ratios change gradually from 2.14% to 18.48% with a height of 40 mm. Pore-size distribution results show that the average pore size of the sample decreases from 15.8 nm to 3.1 nm after diffusion. This corresponds well with TEM results, indicating a pore-filling effect of the Fe compound. In order to precisely control the gradient, diffusion kinetics are further studied by analyzing the influence of time and position on the concentration of the wet gel. At last, it is found that the diffusion process could be fitted well with the one-dimensional model of Fick’s second law, demonstrating the feasibility of the precise design and control of the concentration gradient.

Catalytic conversion in nanoporous materials: Concentration oscillations and spatial correlations due to inhibited transport and intermolecular interactions

DOE Office of Scientific and Technical Information (OSTI.GOV)

Garcia, Andres; Evans, James W.

2016-11-03

We show that steady-state catalytic conversion in nanoporous materials can occur in a quasi-counter-diffusion mode with the reactant (product) concentration strongly decaying (growing) into the pore, but also with oscillations in the total concentration. These oscillations reflect the response of the fluid to the transition from an extended to a confined environment near the pore opening. We focus on the regime of strongly inhibited transport in narrow pores corresponding to single-file diffusion. Here, limited penetration of the reactant into the pores and the associated low reaction yield is impacted by strong spatial correlations induced by both reaction (non-equilibrium correlations) andmore » also by intermolecular interactions (thermodynamic correlations). We develop a generalized hydrodynamic formulation to effectively describe inhibited transport accounting for the effect of these correlations, and incorporate this description of transport into appropriate reaction-diffusion equations. These equations accurately describe both shorter-range concentration oscillations near the pore opening and the longer-range mesoscale variation of concentration profiles in the pore (and thus also describe reaction yield). Success of the analytic theory is validated by comparison with a precise kinetic Monte Carlo simulation of an appropriate molecular-level stochastic reaction-diffusion model. As a result, this work elucidates unconventional chemical kinetics in interacting confined systems.« less

Laboratory triggering of stick-slip events by oscillatory loading in the presence of pore fluid with implications for physics of tectonic tremor

USGS Publications Warehouse

Bartlow, Noel M.; Lockner, David A.; Beeler, Nicholas M.

2012-01-01

The physical mechanism by which the low-frequency earthquakes (LFEs) that make up portions of tectonic (also called non-volcanic) tremor are created is poorly understood. In many areas of the world, tectonic tremor and LFEs appear to be strongly tidally modulated, whereas ordinary earthquakes are not. Anomalous seismic wave speeds, interpreted as high pore fluid pressure, have been observed in regions that generate tremor. Here we build upon previous laboratory studies that investigated the response of stick-slip on artificial faults to oscillatory, tide-like loading. These previous experiments were carried out using room-dry samples of Westerly granite, at one effective stress. Here we augment these results with new experiments on Westerly granite, with the addition of varying effective stress using pore fluid at two pressures. We find that raising pore pressure, thereby lowering effective stress can significantly increase the degree of correlation of stick-slip to oscillatory loading. We also find other pore fluid effects that become important at higher frequencies, when the period of oscillation is comparable to the diffusion time of pore fluid into the fault. These results help constrain the conditions at depth that give rise to tidally modulated LFEs, providing confirmation of the effective pressure law for triggering and insights into why tremor is tidally modulated while earthquakes are at best only weakly modulated.

Dual pore-connectivity and flow-paths affect shale hydrocarbon production

NASA Astrophysics Data System (ADS)

Hayman, N. W.; Daigle, H.; Kelly, E. D.; Milliken, K. L.; Jiang, H.

2016-12-01

Aided with integrated characterization approaches of droplet contact angle measurement, mercury intrusion capillary pressure, low-pressure gas physisorption, scanning electron microscopy, and small angle neutron scattering, we have systematically studied how pore connectivity and wettability are associated with mineral and organic matter phases of shales (Barnett, Bakken, Eagle Ford), as well as their influence on macroscopic fluid flow and hydrocarbon movement, from the following complementary tests: vacuum saturation with vacuum-pulling on dry shale followed with tracer introduction and high-pressure intrusion, tracer diffusion into fluid-saturated shale, fluid and tracer imbibition into partially-saturated shale, and Wood's metal intrusion followed with imaging and elemental mapping. The first three tests use tracer-bearing fluids (hydrophilic API brine and hydrophobic n-decane) fluids with a suite of wettability tracers of different sizes and reactivities developed in our laboratory. These innovative and integrated approaches indicate a Dalmatian wettability behavior at a scale of microns, limited connectivity (<500 microns from shale sample edge) shale pores, and disparity of well-connected hydrophobic pore network ( 10 nm) and sparsely connected hydrophilic pore systems (>50-100 nm), which is linked to the steep initial decline and low overall recovery because of the limited connection of hydrocarbon molecules in the shale matrix to the stimulated fracture network.

Dual pore-connectivity and flow-paths affect shale hydrocarbon production

NASA Astrophysics Data System (ADS)

Hu, Q.; Barber, T.; Zhang, Y.; Md Golam, K.

2017-12-01

Aided with integrated characterization approaches of droplet contact angle measurement, mercury intrusion capillary pressure, low-pressure gas physisorption, scanning electron microscopy, and small angle neutron scattering, we have systematically studied how pore connectivity and wettability are associated with mineral and organic matter phases of shales (Barnett, Bakken, Eagle Ford), as well as their influence on macroscopic fluid flow and hydrocarbon movement, from the following complementary tests: vacuum saturation with vacuum-pulling on dry shale followed with tracer introduction and high-pressure intrusion, tracer diffusion into fluid-saturated shale, fluid and tracer imbibition into partially-saturated shale, and Wood's metal intrusion followed with imaging and elemental mapping. The first three tests use tracer-bearing fluids (hydrophilic API brine and hydrophobic n-decane) fluids with a suite of wettability tracers of different sizes and reactivities developed in our laboratory. These innovative and integrated approaches indicate a Dalmatian wettability behavior at a scale of microns, limited connectivity (<500 microns from shale sample edge) shale pores, and disparity of well-connected hydrophobic pore network ( 10 nm) and sparsely connected hydrophilic pore systems (>50-100 nm), which is linked to the steep initial decline and low overall recovery because of the limited connection of hydrocarbon molecules in the shale matrix to the stimulated fracture network.

A fault constitutive relation accounting for thermal pressurization of pore fluid

USGS Publications Warehouse

Andrews, D.J.

2002-01-01

The heat generated in a slip zone during an earthquake can raise fluid pressure and thereby reduce frictional resistance to slip. The amount of fluid pressure rise depends on the associated fluid flow. The heat generated at a given time produces fluid pressure that decreases inversely with the square root of hydraulic diffusivity times the elapsed time. If the slip velocity function is crack-like, there is a prompt fluid pressure rise at the onset of slip, followed by a slower increase. The stress drop associated with the prompt fluid pressure rise increases with rupture propagation distance. The threshold propagation distance at which thermally induced stress drop starts to dominate over frictionally induced stress drop is proportional to hydraulic diffusivity. If hydraulic diffusivity is 0.02 m2/s, estimated from borehole samples of fault zone material, the threshold propagation distance is 300 m. The stress wave in an earthquake will induce an unknown amount of dilatancy and will increase hydraulic diffusivity, both of which will lessen the fluid pressure effect. Nevertheless, if hydraulic diffusivity is no more than two orders of magnitude larger than the laboratory value, then stress drop is complete in large earthquakes.

Multiscale modeling of fluid flow and mass transport

NASA Astrophysics Data System (ADS)

Masuoka, K.; Yamamoto, H.; Bijeljic, B.; Lin, Q.; Blunt, M. J.

2017-12-01

In recent years, there are some reports on a simulation of fluid flow in pore spaces of rocks using Navier-Stokes equations. These studies mostly adopt a X-ray CT to create 3-D numerical grids of the pores in micro-scale. However, results may be of low accuracy when the rock has a large pore size distribution, because pores, whose size is smaller than resolution of the X-ray CT may be neglected. We recently found out by tracer tests in a laboratory using a brine saturated Ryukyu limestone and inject fresh water that a decrease of chloride concentration took longer time. This phenomenon can be explained due to weak connectivity of the porous networks. Therefore, it is important to simulate entire pore spaces even those of very small sizes in which diffusion is dominant. We have developed a new methodology for multi-level modeling for pore scale fluid flow in porous media. The approach is to combine pore-scale analysis with Darcy-flow analysis using two types of X-ray CT images in different resolutions. Results of the numerical simulations showed a close match with the experimental results. The proposed methodology is an enhancement for analyzing mass transport and flow phenomena in rocks with complicated pore structure.

Computational approach to integrate 3D X-ray microtomography and NMR data

NASA Astrophysics Data System (ADS)

Lucas-Oliveira, Everton; Araujo-Ferreira, Arthur G.; Trevizan, Willian A.; Fortulan, Carlos A.; Bonagamba, Tito J.

2018-07-01

Nowadays, most of the efforts in NMR applied to porous media are dedicated to studying the molecular fluid dynamics within and among the pores. These analyses have a higher complexity due to morphology and chemical composition of rocks, besides dynamic effects as restricted diffusion, diffusional coupling, and exchange processes. Since the translational nuclear spin diffusion in a confined geometry (e.g. pores and fractures) requires specific boundary conditions, the theoretical solutions are restricted to some special problems and, in many cases, computational methods are required. The Random Walk Method is a classic way to simulate self-diffusion along a Digital Porous Medium. Bergman model considers the magnetic relaxation process of the fluid molecules by including a probability rate of magnetization survival under surface interactions. Here we propose a statistical approach to correlate surface magnetic relaxivity with the computational method applied to the NMR relaxation in order to elucidate the relationship between simulated relaxation time and pore size of the Digital Porous Medium. The proposed computational method simulates one- and two-dimensional NMR techniques reproducing, for example, longitudinal and transverse relaxation times (T1 and T2, respectively), diffusion coefficients (D), as well as their correlations. For a good approximation between the numerical and experimental results, it is necessary to preserve the complexity of translational diffusion through the microstructures in the digital rocks. Therefore, we use Digital Porous Media obtained by 3D X-ray microtomography. To validate the method, relaxation times of ideal spherical pores were obtained and compared with the previous determinations by the Brownstein-Tarr model, as well as the computational approach proposed by Bergman. Furthermore, simulated and experimental results of synthetic porous media are compared. These results make evident the potential of computational physics in the analysis of the NMR data for complex porous materials.

The influence of different diffusion pattern to the sub- and super-critical fluid flow in brown coal

NASA Astrophysics Data System (ADS)

Peng, Peihuo

2018-03-01

Sub- and super-critical CO2 flowing in nanoscale pores are recently becoming of great interest due to that it is closely related to many engineering applications, such as geological burial and sequestration of carbon dioxide, Enhanced Coal Bed Methane recovery ( ECBM), super-critical CO2 fracturing and so on. Gas flow in nanopores cannot be described simply by the Darcy equation. Different diffusion pattern such as Fick diffusion, Knudsen diffusion, transitional diffusion and slip flow at the solid matrix separate the seepage behaviour from Darcy-type flow. According to the principle of different diffusion pattern, the flow of sub- and super-critical CO2 in brown coal was simulated by numerical method, and the results were compared with the experimental results to explore the contribution of different diffusion pattern and swelling effect in sub- and super-critical CO2 flow in nanoscale pores.

Rock deformation models and fluid leak-off in hydraulic fracturing

NASA Astrophysics Data System (ADS)

Yarushina, Viktoriya M.; Bercovici, David; Oristaglio, Michael L.

2013-09-01

Fluid loss into reservoir rocks during hydraulic fracturing is modelled via a poro-elastoplastic pressure diffusion equation in which the total compressibility is a sum of fluid, rock and pore space compressibilities. Inclusion of pore compressibility and porosity-dependent permeability in the model leads to a strong pressure dependence of leak-off (i.e. drainage rate). Dilation of the matrix due to fluid invasion causes higher rates of fluid leak-off. The present model is appropriate for naturally fractured and tight gas reservoirs as well as for soft and poorly consolidated formations whose mechanical behaviour departs from simple elastic laws. Enhancement of the leak-off coefficient by dilation, predicted by the new model, may help explain the low percentage recovery of fracturing fluid (usually between 5 and 50 per cent) in shale gas stimulation by hydraulic fracturing.

Separating attoliter-sized compartments using fluid pore-spanning lipid bilayers.

PubMed

Lazzara, Thomas D; Carnarius, Christian; Kocun, Marta; Janshoff, Andreas; Steinem, Claudia

2011-09-27

Anodic aluminum oxide (AAO) is a porous material having aligned cylindrical compartments with 55-60 nm diameter pores, and being several micrometers deep. A protocol was developed to generate pore-spanning fluid lipid bilayers separating the attoliter-sized compartments of the nanoporous material from the bulk solution, while preserving the optical transparency of the AAO. The AAO was selectively functionalized by silane chemistry to spread giant unilamellar vesicles (GUVs) resulting in large continuous membrane patches covering the pores. Formation of fluid single lipid bilayers through GUV rupture could be readily observed by fluorescence microscopy and further supported by conservation of membrane surface area, before and after GUV rupture. Fluorescence recovery after photobleaching gave low immobile fractions (5-15%) and lipid diffusion coefficients similar to those found for bilayers on silica. The entrapment of molecules within the porous underlying cylindrical compartments, as well as the exclusion of macromolecules from the nanopores, demonstrate the barrier function of the pore-spanning membranes and could be investigated in three-dimensions using confocal laser scanning fluorescence imaging. © 2011 American Chemical Society

Effect of water phase transition on dynamic ruptures with thermal pressurization: Numerical simulations with changes in physical properties of water

NASA Astrophysics Data System (ADS)

Urata, Yumi; Kuge, Keiko; Kase, Yuko

2015-02-01

Phase transitions of pore water have never been considered in dynamic rupture simulations with thermal pressurization (TP), although they may control TP. From numerical simulations of dynamic rupture propagation including TP, in the absence of any water phase transition process, we predict that frictional heating and TP are likely to change liquid pore water into supercritical water for a strike-slip fault under depth-dependent stress. This phase transition causes changes of a few orders of magnitude in viscosity, compressibility, and thermal expansion among physical properties of water, thus affecting the diffusion of pore pressure. Accordingly, we perform numerical simulations of dynamic ruptures with TP, considering physical properties that vary with the pressure and temperature of pore water on a fault. To observe the effects of the phase transition, we assume uniform initial stress and no fault-normal variations in fluid density and viscosity. The results suggest that the varying physical properties decrease the total slip in cases with high stress at depth and small shear zone thickness. When fault-normal variations in fluid density and viscosity are included in the diffusion equation, they activate TP much earlier than the phase transition. As a consequence, the total slip becomes greater than that in the case with constant physical properties, eradicating the phase transition effect. Varying physical properties do not affect the rupture velocity, irrespective of the fault-normal variations. Thus, the phase transition of pore water has little effect on dynamic ruptures. Fault-normal variations in fluid density and viscosity may play a more significant role.

Degradation of porous poly(D,L-lactic-co-glycolic acid) films based on water diffusion.

PubMed

Huang, Ying-Ying; Qi, Min; Liu, Hong-Ze; Zhao, Hong; Yang, Da-Zhi

2007-03-15

Poly(D,L-lactic-co-glycolic acid) has been extensively used as a controlled release carrier for drug delivery due to its good biocompatibility, biodegradability, and mechanical strength. Effects of dense and porous film's degradation behavior have been systematically investigated up to 17 weeks in Hank's Simulated Body Fluid at 37 degrees C. The degradation of the films was studied by measuring changes in weight, molecular weight and its distribution, morphology, composition etc.. A special thing was that the differences in water diffusion in dense and porous structure films caused the different degradation behavior. According to the characteristic changes of various properties of films, the degradation process is suggested to be roughly divided into four stages, tentatively named as water absorption stage, dramatic loss of molecular weight or micro-pores formed stage, loss of weight or enlarged-pores formed stage, pores diminished or pores collapse stage.

Diffusion of Supercritical Fluids through Single-Layer Nanoporous Solids: Theory and Molecular Simulations.

PubMed

Oulebsir, Fouad; Vermorel, Romain; Galliero, Guillaume

2018-01-16

With the advent of graphene material, membranes based on single-layer nanoporous solids appear as promising devices for fluid separation, be it liquid or gaseous mixtures. The design of such architectured porous materials would greatly benefit from accurate models that can predict their transport and separation properties. More specifically, there is no universal understanding of how parameters such as temperature, fluid loading conditions, or the ratio of the pore size to the fluid molecular diameter influence the permeation process. In this study, we address the problem of pure supercritical fluids diffusing through simplified models of single-layer porous materials. Basically, we investigate a toy model that consists of a single-layer lattice of Lennard-Jones interaction sites with a slit gap of controllable width. We performed extensive equilibrium and biased molecular dynamics simulations to document the physical mechanisms involved at the molecular scale. We propose a general constitutive equation for the diffusional transport coefficient derived from classical statistical mechanics and kinetic theory, which can be further simplified in the ideal gas limit. This transport coefficient relates the molecular flux to the fluid density jump across the single-layer membrane. It is found to be proportional to the accessible surface porosity of the single-layer porous solid and to a thermodynamic factor accounting for the inhomogeneity of the fluid close to the pore entrance. Both quantities directly depend on the potential of mean force that results from molecular interactions between solid and fluid atoms. Comparisons with the simulations data show that the kinetic model captures how narrowing the pore size below the fluid molecular diameter lowers dramatically the value of the transport coefficient. Furthermore, we demonstrate that our general constitutive equation allows for a consistent interpretation of the intricate effects of temperature and fluid loading conditions on the permeation process.

Multi-Scale Multi-Physics Modeling of Matrix Transport Properties in Fractured Shale Reservoirs

NASA Astrophysics Data System (ADS)

Mehmani, A.; Prodanovic, M.

2014-12-01

Understanding the shale matrix flow behavior is imperative in successful reservoir development for hydrocarbon production and carbon storage. Without a predictive model, significant uncertainties in flowback from the formation, the communication between the fracture and matrix as well as proper fracturing practice will ensue. Informed by SEM images, we develop deterministic network models that couple pores from multiple scales and their respective fluid physics. The models are used to investigate sorption hysteresis as an affordable way of inferring the nanoscale pore structure in core scale. In addition, restricted diffusion as a function of pore shape, pore-throat size ratios and network connectivity is computed to make correct interpretation of the 2D NMR maps possible. Our novel pore network models have the ability to match sorption hysteresis measurements without any tuning parameters. The results clarify a common misconception of linking type 3 nitrogen hysteresis curves to only the shale pore shape and show promising sensitivty for nanopore structre inference in core scale. The results on restricted diffusion shed light on the importance of including shape factors in 2D NMR interpretations. A priori "weighting factors" as a function of pore-throat and throat-length ratio are presented and the effect of network connectivity on diffusion is quantitatively assessed. We are currently working on verifying our models with experimental data gathered from the Eagleford formation.

Controlling the Release of Indomethacin from Glass Solutions Layered with a Rate Controlling Membrane Using Fluid-Bed Processing. Part 2: The Influence of Formulation Parameters on Drug Release.

PubMed

Dereymaker, Aswin; Pelgrims, Jirka; Engelen, Frederik; Adriaensens, Peter; Van den Mooter, Guy

2017-04-03

This study aimed to investigate the pharmaceutical performance of an indomethacin-polyvinylpyrrolidone (PVP) glass solution applied using fluid bed processing as a layer on inert sucrose spheres and subsequently top-coated with a release rate controlling membrane consisting of either ethyl cellulose or Eudragit RL. The implications of the addition of a pore former (PVP) and the coating medium (ethanol or water) on the diffusion and release behavior were also considered. In addition, the role of a charge interaction between drug and controlled release polymer on the release was investigated. Diffusion experiments pointed to the influence of pore former concentration, rate controlling polymer type, and coating solvent on the permeability of the controlled release membranes. This can be translated to drug release tests, which show the potential of diffusion tests as a preliminary screening test and that diffusion is the main factor influencing release. Drug release tests also showed the effect of coating layer thickness. A charge interaction between INDO and ERL was demonstrated, but this had no negative effect on drug release. The higher diffusion and release observed in ERL-based rate controlling membranes was explained by a higher hydrophilicity, compared to EC.

Equation of state and some structural and dynamical properties of the confined Lennard-Jones fluid into carbon nanotube: A molecular dynamics study

NASA Astrophysics Data System (ADS)

Abbaspour, Mohsen; Akbarzadeh, Hamed; Salemi, Sirous; Abroodi, Mousarreza

2016-11-01

By considering the anisotropic pressure tensor, two separate equations of state (EoS) as functions of the density, temperature, and carbon nanotube (CNT) diameter have been proposed for the radial and axial directions for the confined Lennard-Jones (LJ) fluid into (11,11), (12,10), and (19,0) CNTs from 120 to 600 K using molecular dynamics (MD) simulations. We have also investigated the effects of the pore size, pore loading, chirality, and temperature on some of the structural and dynamical properties of the confined LJ fluid into (11,11), (12,10), (19,0), and (19,19) CNTs such as the radial density profile and self-diffusion coefficient. We have also determined the EoS for the confined LJ fluid into double and triple walled CNTs.

Computational approach to integrate 3D X-ray microtomography and NMR data.

PubMed

Lucas-Oliveira, Everton; Araujo-Ferreira, Arthur G; Trevizan, Willian A; Fortulan, Carlos A; Bonagamba, Tito J

2018-05-04

Nowadays, most of the efforts in NMR applied to porous media are dedicated to studying the molecular fluid dynamics within and among the pores. These analyses have a higher complexity due to morphology and chemical composition of rocks, besides dynamic effects as restricted diffusion, diffusional coupling, and exchange processes. Since the translational nuclear spin diffusion in a confined geometry (e.g. pores and fractures) requires specific boundary conditions, the theoretical solutions are restricted to some special problems and, in many cases, computational methods are required. The Random Walk Method is a classic way to simulate self-diffusion along a Digital Porous Medium. Bergman model considers the magnetic relaxation process of the fluid molecules by including a probability rate of magnetization survival under surface interactions. Here we propose a statistical approach to correlate surface magnetic relaxivity with the computational method applied to the NMR relaxation in order to elucidate the relationship between simulated relaxation time and pore size of the Digital Porous Medium. The proposed computational method simulates one- and two-dimensional NMR techniques reproducing, for example, longitudinal and transverse relaxation times (T 1 and T 2 , respectively), diffusion coefficients (D), as well as their correlations. For a good approximation between the numerical and experimental results, it is necessary to preserve the complexity of translational diffusion through the microstructures in the digital rocks. Therefore, we use Digital Porous Media obtained by 3D X-ray microtomography. To validate the method, relaxation times of ideal spherical pores were obtained and compared with the previous determinations by the Brownstein-Tarr model, as well as the computational approach proposed by Bergman. Furthermore, simulated and experimental results of synthetic porous media are compared. These results make evident the potential of computational physics in the analysis of the NMR data for complex porous materials. Copyright © 2018 Elsevier Inc. All rights reserved.

Data-driven fault mechanics: Inferring fault hydro-mechanical properties from in situ observations of injection-induced aseismic slip

NASA Astrophysics Data System (ADS)

Bhattacharya, P.; Viesca, R. C.

2017-12-01

In the absence of in situ field-scale observations of quantities such as fault slip, shear stress and pore pressure, observational constraints on models of fault slip have mostly been limited to laboratory and/or remote observations. Recent controlled fluid-injection experiments on well-instrumented faults fill this gap by simultaneously monitoring fault slip and pore pressure evolution in situ [Gugleilmi et al., 2015]. Such experiments can reveal interesting fault behavior, e.g., Gugleilmi et al. report fluid-activated aseismic slip followed only subsequently by the onset of micro-seismicity. We show that the Gugleilmi et al. dataset can be used to constrain the hydro-mechanical model parameters of a fluid-activated expanding shear rupture within a Bayesian framework. We assume that (1) pore-pressure diffuses radially outward (from the injection well) within a permeable pathway along the fault bounded by a narrow damage zone about the principal slip surface; (2) pore-pressure increase ativates slip on a pre-stressed planar fault due to reduction in frictional strength (expressed as a constant friction coefficient times the effective normal stress). Owing to efficient, parallel, numerical solutions to the axisymmetric fluid-diffusion and crack problems (under the imposed history of injection), we are able to jointly fit the observed history of pore-pressure and slip using an adaptive Monte Carlo technique. Our hydrological model provides an excellent fit to the pore-pressure data without requiring any statistically significant permeability enhancement due to the onset of slip. Further, for realistic elastic properties of the fault, the crack model fits both the onset of slip and its early time evolution reasonably well. However, our model requires unrealistic fault properties to fit the marked acceleration of slip observed later in the experiment (coinciding with the triggering of microseismicity). Therefore, besides producing meaningful and internally consistent bounds on in-situ fault properties like permeability, storage coefficient, resolved stresses, friction and the shear modulus, our results also show that fitting the complete observed time history of slip requires alternative model considerations, such as variations in fault mechanical properties or friction coefficient with slip.

A multi-scale Lattice Boltzmann model for simulating solute transport in 3D X-ray micro-tomography images of aggregated porous materials

NASA Astrophysics Data System (ADS)

Zhang, Xiaoxian; Crawford, John W.; Flavel, Richard J.; Young, Iain M.

2016-10-01

The Lattice Boltzmann (LB) model and X-ray computed tomography (CT) have been increasingly used in combination over the past decade to simulate water flow and chemical transport at pore scale in porous materials. Because of its limitation in resolution and the hierarchical structure of most natural soils, the X-ray CT tomography can only identify pores that are greater than its resolution and treats other pores as solid. As a result, the so-called solid phase in X-ray images may in reality be a grey phase, containing substantial connected pores capable of conducing fluids and solute. Although modified LB models have been developed to simulate fluid flow in such media, models for solute transport are relatively limited. In this paper, we propose a LB model for simulating solute transport in binary soil images containing permeable solid phase. The model is based on the single-relaxation time approach and uses a modified partial bounce-back method to describe the resistance caused by the permeable solid phase to chemical transport. We derive the relationship between the diffusion coefficient and the parameter introduced in the partial bounce-back method, and test the model against analytical solution for movement of a pulse of tracer. We also validate it against classical finite volume method for solute diffusion in a simple 2D image, and then apply the model to a soil image acquired using X-ray tomography at resolution of 30 μm in attempts to analyse how the ability of the solid phase to diffuse solute at micron-scale affects the behaviour of the solute at macro-scale after a volumetric average. Based on the simulated results, we discuss briefly the danger in interpreting experimental results using the continuum model without fully understanding the pore-scale processes, as well as the potential of using pore-scale modelling and tomography to help improve the continuum models.

Reverse weathering in marine sediments and the geochemical cycle of potassium in seawater: Insights from the K isotopic composition (41K/39K) of deep-sea pore-fluids

USGS Publications Warehouse

Santiago Ramos, Danielle P.; Morgan, Leah; Lloyd, Nicholas S.; Higgins, John A.

2018-01-01

In situ Al-silicate formation, also known as “reverse weathering,” is an important sink of many of the major and minor cations in seawater (e.g. Mg, K, and Li). However, the importance of this sink in global geochemical cycles and isotopic budgets of these elements remains poorly constrained. Here, we report on the potassium isotopic composition (41">41K/39">39K) of deep-sea sediment pore-fluids from four (Integrated) Ocean Drilling Program sites (1052, U1378, U1395 and U1403) to characterize potassium isotopic fractionation associated with the formation of authigenic Al-silicate minerals in marine sediments and its role in elevating the 41">41K/39">39K of seawater relative to bulk silicate Earth. Isotopic ratios are obtained by high-resolution multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) in cold plasma conditions with a long-term external reproducibility of ca. 0.17‰. We find that, although all sites are characterized by pore-fluid K concentrations that decline with increasing depth, their K isotopic profiles vary systematically from site-to-site; at sites characterized by rapid sedimentation rates, pore-fluid profiles of 41">41K/39">39K are relatively invariant whereas at sites characterized by slow sedimentation rates, 41">41K/39">39K declines with depth by up to 1.8‰. Results from 1-D diffusion-advection-reaction models suggest that these differences may result from a complex interplay between sedimentation rate and fractionation of K isotopes during diffusion, Al-silicate authigenesis, and ion exchange. Model simulations suggest fractionation factors between 0.9980 and 1.0000 for reverse weathering reactions in deep-sea sediments. Although deep-sea sites do not constitute major sinks of K in seawater, some of the processes responsible for K isotopic fractionation at these sites (diffusion and Al-silicate authigenesis) likely play a role in determining the 41">41K/39">39K of seawater.

Assessing the role of clay authigenesis in the seawater potassium cycle: A paired K and Mg isotope study of deep-sea pore fluids

NASA Astrophysics Data System (ADS)

Santiago Ramos, D. P.; Higgins, J. A.

2017-12-01

In situ formation of clays (clay authigenesis) in marine sediments and altered oceanic crust is an important sink of a number of seawater cations. In particular, clay authigenesis is a major, and yet unconstrained, flux in the global seawater potassium cycle. Potassium is the fourth most abundant cation in the ocean, which constitutes an isotopically enriched K reservoir (δ41K 0‰) compared to the solid Earth (δ41K -0.5‰). Understanding what processes control this isotopic offset is the main goal of this study. Here we use a multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) to measure the K and Mg isotope ratios (41K/39K and 26Mg/24Mg) of deep-sea pore fluids in order to assess the effects of clay formation in the K isotope composition of seawater. Mg isotopes are used as an independent proxy for clay formation, since marine authigenic clays are isotopically distinct from their detrital counterpart, an abundant component of marine sediments. Our study sites (ODP/IODP 1052, U1378, U1395, U1403) vary in location, lithology, age and sedimentation rates; however, pore-fluids from all sites show decreasing K concentrations with depth, suggesting potassium uptake into the sediments. We find that although K concentration trends are similar across all sites, measured δ41K values vary significantly. Results from 1-D diffusion-advection-reaction models suggest that these differences in isotopic profiles arise from a complex interplay between sedimentation rate and K isotopic fractionation during clay formation, aqueous K diffusion and ion exchange reactions. Further, model simulations yield fractionation factors between 0.9980 and 1.0000 for clay formation in deep-sea sediments. Despite the minor contribution of these deep-sea pore-fluids as sinks of seawater K, the processes responsible for K isotope fractionation in our study sites (clay formation and aqueous K diffusion) are also observed at shallow marine systems (major K sinks) and are thus likely responsible for setting the K isotopic composition of seawater.

Experimental Study of the Roles of Mechanical and Hydrologic Properties in the Initiation of Natural Hydraulic Fractures

NASA Astrophysics Data System (ADS)

French, M. E.; Goodwin, L. B.; Boutt, D. F.; Lilydahl, H.

2008-12-01

Natural hydraulic fractures (NHFs) are inferred to form where pore fluid pressure exceeds the least compressive stress; i.e., where the hydraulic fracture criterion is met. Although it has been shown that mechanical heterogeneities serve as nuclei for NHFs, the relative roles of mechanical anisotropy and hydrologic properties in initiating NHFs in porous granular media have not been fully explored. We designed an experimental protocol that produces a pore fluid pressure high enough to exceed the hydraulic fracture criterion, allowing us to initiate NHFs in the laboratory. Initially, cylindrical samples 13 cm long and 5 cm in diameter are saturated, σ1 is radial, and σ3 is axial. By dropping the end load (σ3) and pore fluid pressure simultaneously at the end caps, we produce a large pore fluid pressure gradient parallel to the long axis of the sample. This allows us to meet the hydraulic fracture criterion without separating the sample from its end caps. The time over which the pore fluid remains elevated is a function of hydraulic diffusivity. An initial test with a low diffusivity sandstone produced NHFs parallel to bedding laminae that were optimally oriented for failure. To evaluate the relative importance of mechanical heterogeneities such as bedding versus hydraulic properties, we are currently investigating variably cemented St. Peter sandstone. This quartz arenite exhibits a wide range of primary structures, from well developed bedding laminae to locally massive sandstone. Diagenesis has locally accentuated these structures, causing degree of cementation to vary with bedding, and the sandstone locally exhibits concretions that form elliptical rather than tabular heterogeneities. Bulk permeability varies from k=10-12 m2 to k=10-15 m2 and porosity varies from 5% to 28% in this suite of samples. Variations in a single sample are smaller, with permeability varying no more than an order of magnitude within a single core. Air minipermeameter and tracer tests document this variability at the cm scale. Experiments will be performed with σ3 and the pore pressure gradient both perpendicular and parallel to sub-cm scale bedding. The results of these tests will be compared to those of structurally homogeneous samples and samples with elliptical heterogeneities.

Peritoneal Fluid Transport rather than Peritoneal Solute Transport Associates with Dialysis Vintage and Age of Peritoneal Dialysis Patients.

PubMed

Waniewski, Jacek; Antosiewicz, Stefan; Baczynski, Daniel; Poleszczuk, Jan; Pietribiasi, Mauro; Lindholm, Bengt; Wankowicz, Zofia

2016-01-01

During peritoneal dialysis (PD), the peritoneal membrane undergoes ageing processes that affect its function. Here we analyzed associations of patient age and dialysis vintage with parameters of peritoneal transport of fluid and solutes, directly measured and estimated based on the pore model, for individual patients. Thirty-three patients (15 females; age 60 (21-87) years; median time on PD 19 (3-100) months) underwent sequential peritoneal equilibration test. Dialysis vintage and patient age did not correlate. Estimation of parameters of the two-pore model of peritoneal transport was performed. The estimated fluid transport parameters, including hydraulic permeability (LpS), fraction of ultrasmall pores (α u), osmotic conductance for glucose (OCG), and peritoneal absorption, were generally independent of solute transport parameters (diffusive mass transport parameters). Fluid transport parameters correlated whereas transport parameters for small solutes and proteins did not correlate with dialysis vintage and patient age. Although LpS and OCG were lower for older patients and those with long dialysis vintage, αu was higher. Thus, fluid transport parameters--rather than solute transport parameters--are linked to dialysis vintage and patient age and should therefore be included when monitoring processes linked to ageing of the peritoneal membrane.

Anisotropic changes in P-wave velocity and attenuation during deformation and fluid infiltration of granite

USGS Publications Warehouse

Stanchits, S.A.; Lockner, D.A.; Ponomarev, A.V.

2003-01-01

Fluid infiltration and pore fluid pressure changes are known to have a significant effect on the occurrence of earthquakes. Yet, for most damaging earthquakes, with nucleation zones below a few kilometers depth, direct measurements of fluid pressure variations are not available. Instead, pore fluid pressures are inferred primarily from seismic-wave propagation characteristics such as Vp/Vs ratio, attenuation, and reflectivity contacts. We present laboratory measurements of changes in P-wave velocity and attenuation during the injection of water into a granite sample as it was loaded to failure. A cylindrical sample of Westerly granite was deformed at constant confining and pore pressures of 50 and 1 MPa, respectively. Axial load was increased in discrete steps by controlling axial displacement. Anisotropic P-wave velocity and attenuation fields were determined during the experiment using an array of 13 piezoelectric transducers. At the final loading steps (86% and 95% of peak stress), both spatial and temporal changes in P-wave velocity and peak-to-peak amplitudes of P and S waves were observed. P-wave velocity anisotropy reached a maximum of 26%. Transient increases in attenuation of up to 483 dB/m were also observed and were associated with diffusion of water into the sample. We show that velocity and attenuation of P waves are sensitive to the process of opening of microcracks and the subsequent resaturation of these cracks as water diffuses in from the surrounding region. Symmetry of the orientation of newly formed microcracks results in anisotropic velocity and attenuation fields that systematically evolve in response to changes in stress and influx of water. With proper scaling, these measurements provide constraints on the magnitude and duration of velocity and attenuation transients that can be expected to accompany the nucleation of earthquakes in the Earth's crust.

Paired measurements of K and Mg isotopes and clay authigenesis in marine sediments

NASA Astrophysics Data System (ADS)

Santiago Ramos, D. P.; Dunlea, A. G.; Higgins, J. A.

2016-12-01

Despite its importance as a major sink for seawater K and Mg, estimates of clay authigenesis in marine sediments remain poorly constrained. Previous work on Mg isotope fractionation during clay formation has revealed a preferential uptake of 26Mg, yielding authigenic clay products with potentially distinct δ26Mg compared to the detrital component. In a similar manner, we aim to quantify the K isotope fractionation during authigenic clay formation and to use paired δ26Mg and δ41K measurements as proxies for the identification and quantification of authigenic clays in shallow and deep marine sedimentary systems. To better understand the behavior of paired Mg and K isotopes during authigenic clay formation in marine sediments, we measured δ26Mg and δ41K values of pore-fluids and sediments from ODP/IODP sites 1052, U1395, U1403 and U1366. We find that while pore-fluid K concentrations at sites 1052, U1395 and U1403 all decline with depth, δ41K profiles differ significantly. These differences might be a result of a complex interplay between clay authigenesis, sedimentation rate, and fractionation of K isotopes during diffusion. Results from 1-D diffusion-advection-reaction models suggest that, in contrast to Mg, diffusion may play an important role in determining the overall K isotope fractionation during clay authigenesis in sites with low-sedimentation rates. Sites with high sedimentation rates may act as close systems where diffusion is negligible. In such cases, K uptake can be modeled as a Rayleigh distillation process and K isotope fractionation can be estimated. Measurements of δ26Mg and δ41K of pore-fluids from site U1395 and bulk sediments from U1366 suggest that paired measurements of these isotopic systems in siliciclastic marine sediments can provide new insights into rates of marine clay authigenesis, a globally important but understudied component of many geochemical cycles.

A Systematic Procedure to Describe Shale Gas Permeability Evolution during the Production Process

NASA Astrophysics Data System (ADS)

Jia, B.; Tsau, J. S.; Barati, R.

2017-12-01

Gas flow behavior in shales is complex due to the multi-physics nature of the process. Pore size reduces as the in-situ stress increases during the production process, which will reduce intrinsic permeability of the porous media. Slip flow/pore diffusion enhances gas apparent permeability, especially under low reservoir pressures. Adsorption not only increases original gas in place but also influences gas flow behavior because of the adsorption layer. Surface diffusion between free gas and adsorption phase enhances gas permeability. Pore size reduction and the adsorption layer both have complex impacts on gas apparent permeability and non-Darcy flow might be a major component in nanopores. Previously published literature is generally incomplete in terms of coupling of all these four physics with fluid flow during gas production. This work proposes a methodology to simultaneously take them into account to describe a permeability evolution process. Our results show that to fully describe shale gas permeability evolution during gas production, three sets of experimental data are needed initially: 1) intrinsic permeability under different in-situ stress, 2) adsorption isotherm under reservoir conditions and 3) surface diffusivity measurement by the pulse-decay method. Geomechanical effects, slip flow/pore diffusion, adsorption layer and surface diffusion all play roles affecting gas permeability. Neglecting any of them might lead to misleading results. The increasing in-situ stress during shale gas production is unfavorable to shale gas flow process. Slip flow/pore diffusion is important for gas permeability under low pressures in the tight porous media. They might overwhelm the geomechanical effect and enhance gas permeability at low pressures. Adsorption layer reduces the gas permeability by reducing the effective pore size, but the effect is limited. Surface diffusion increases gas permeability more under lower pressures. The total gas apparent permeability might keep increasing during the gas production process when the surface diffusivity is larger than a critical value. We believe that our workflow proposed in this study will help describe shale gas permeability evolution considering all the underlying physics altogether.

Helium isotopes in matrix pore fluids from SAFOD drill core samples suggest mantle fluids cannot be responsible for fault weakening

NASA Astrophysics Data System (ADS)

Ali, S.; Stute, M.; Torgersen, T.; Winckler, G.

2008-12-01

To quantify fluid flow in the San Andreas Fault (SAF) (and since direct fracture fluid sampling of the fault zone was not available), we have adapted a method to extract rare gases from matrix fluids of whole rocks by diffusion. Helium was measured on drill core samples obtained from 3054 m (Pacific Plate) to 3990 m (North American Plate) through the San Andreas Fault Zone (SAFZ) ~3300 m during SAFOD Phases I (2004), II (2005), III (2007). Samples were typically collected as 2.54 cm diameter subcores drilled into the ends of the cores, or from the core catcher and drillcore fragments within <2hr after core recovery. The samples were placed into ultra high vacuum stainless steel containers, flushed with ultra high purity nitrogen and immediately evacuated. Helium isotopes of the extracted matrix pore fluids and the solid matrix were determined by mass spectrometery at LDEO. Matrix porefluid 3He/4He ratios are ~0.4 - 0.5xRa (Ra: atmospheric 3He/4He = 1.384 x 10-6) in the Pacific Plate, increasing toward the SAFZ, while pore fluids in the North American Plate have a 3He/4He range of 0.7-0.9Ra, increasing away from the SAFZ (consistent with results from mud gas samples (Wiersberg and Erzinger, 2007) and direct fluid samples (Kennedy et al., 2007)). Helium isotope ratios of the solid matrix are less than 0.06Ra across the SAF in samples from both the North American and the Pacific plates, thereby excluding the host matrix as source for the enhanced isotopic signature. If the system is assumed to be in steady state, then the flux of mantle helium must be from the North American Plate to the Pacific plate. The steeper gradient in the Pacific Plate relative to the North American plate is consistent with a porosity corrected effective diffusivity. The source for this mantle helium in the North American Plate is likely related to a low crustal conductivity zone identified by magnetotelluric signals (Becken et al., 2008) that provides a channel for transport of mantle helium within brittle crust under high strain rates (Kennedy et al., 2007). The helium isotope gradients suggest that fault weakening by mantle-derived fluid pressure is unlikely. More likely, mantle fluids "bleed" into the North American plate below seismogenic depths and are transported across the fault by nonseismic, diffusive processes.

A thermodynamically consistent model for granular-fluid mixtures considering pore pressure evolution and hypoplastic behavior

NASA Astrophysics Data System (ADS)

Hess, Julian; Wang, Yongqi

2016-11-01

A new mixture model for granular-fluid flows, which is thermodynamically consistent with the entropy principle, is presented. The extra pore pressure described by a pressure diffusion equation and the hypoplastic material behavior obeying a transport equation are taken into account. The model is applied to granular-fluid flows, using a closing assumption in conjunction with the dynamic fluid pressure to describe the pressure-like residual unknowns, hereby overcoming previous uncertainties in the modeling process. Besides the thermodynamically consistent modeling, numerical simulations are carried out and demonstrate physically reasonable results, including simple shear flow in order to investigate the vertical distribution of the physical quantities, and a mixture flow down an inclined plane by means of the depth-integrated model. Results presented give insight in the ability of the deduced model to capture the key characteristics of granular-fluid flows. We acknowledge the support of the Deutsche Forschungsgemeinschaft (DFG) for this work within the Project Number WA 2610/3-1.

Flux rates and sulfur isotopic composition of pore fluids from three mud volcanoes in the northern Gulf of Mexico

NASA Astrophysics Data System (ADS)

Gilhooly, W. P.; Ruppel, C. D.; Dickens, G. R.; Berg, P.; Macko, S. A.

2010-12-01

Chloride and sulfate pore water analyses were performed on a total of 29 piston and gravity cores collected along center to flank transects across 3 mud volcanoes, which were located on the Louisiana continental slope in Garden Banks (GB425), Green Canyon (CG185), and Mississippi Canyon (MC852). All three sites are known areas of oil and gas discharge. In addition, seepage at GC185 and GB425 supports highly developed chemosynthetic communities, whereas no known communities have been observed at MC852. Comparison of pore water chemistry (sulfur concentrations and sulfur isotope compositions) among these 3 sites provides initial insight about fluid migration processes and advection rates and about the connection between fluid flux and the establishment of chemosynthetic communities. Pore water advection velocities were calculated from chloride profiles using a steady-state one dimensional advection-diffusion model. In general, chloride concentrations increased with depth to more than four times seawater concentrations. Incidences of pore water freshening are likely associated with hydrate dissociation. Chloride profiles show characteristic concave-up shapes at the center of each mud volcano and concave-down shapes along the flanks, a pattern that we previously interpreted and modeled (doi:10.1029/2004GL021909; doi:10.1111/j.1468-8123.2007.00191.x) in terms of seawater recharge-discharge. The depth of the sulfate-methane interface (SMI) shoals toward the center of the mud volcanoes, indicating potentially rapid anaerobic methane oxidation in these areas. Where the SMI is shallow, pore water sulfide S-isotope values are correspondingly elevated (~ +10 ‰) relative to seawater sulfate (δ34S = +21‰) and presumably represent near-quantitative reduction of pore water sulfate at GB425 and MC852. There is no such pattern at GC185. Such differences potentially reflect advection rates, the ages of the fluids, timing of fluid efflux, and differences in their chemistry.

Fluid Absorption and Release of Nonwovens and their Response to Compression

NASA Astrophysics Data System (ADS)

Bateny, Fatemeh

Fluid handling is a key property in one of the major nonwoven applications in absorbent product such as wipes, hygiene products, and baby diapers. These products are subjected to various levels of compression in real-use. The aim of this study was to investigate the liquid absorption and release properties of nonwovens to establish the absorption structure-property relationship at various compression levels. A comprehensive methodology, considering various flow directions, was employed to establish the relationship by decoupling the effect of structural parameters and material properties in two phases of this study respectively. In the first phase, the mechanism of absorption by pore structure was investigated through considering various fiber cross-sectional size and shape, as well as heterogeneous layered structures having a pore size reduction and expansion. In the second phase, the mechanism of absorption by fiber and consequent swelling was evaluated in view of fluid diffusion into the rayon fibers in samples having different percentages of PET fiber (non-absorbent) and rayon fiber (absorbent). The analysis of absorption and release properties through the entire dissertation was based on the pore characteristics of the nonwovens by measuring the average pore sizes, pore size distribution, and solidity. The investigation revealed that the absorption and release properties of nonwovens are governed by their pore characteristics. In homogeneous non-layered nonwoven fabrics, maximum absorption is mainly governed by the available pore volume. Absorbency rate is determined according to pore size and the maximum rate of absorption is achieved at a specific range of pore sizes. This indicates that an in-depth understanding of the absorption and release properties brings about valuable information for the absorbent product engineering.

Driving forces for metamorphic vein filling during bauxite dehydration: insights from Li and Al transfer illustrated by LIBS compositional profiles (Western Alps)

NASA Astrophysics Data System (ADS)

Verlaguet, Anne; Brunet, Fabrice; Goffé, Bruno; Menut, Denis; Findling, Nathaniel; Poinssot, Christophe

2015-04-01

In subduction zones, the significant amounts of aqueous fluid released in the course of the successive dehydration reactions occurring during prograde metamorphism are expected to strongly influence the rock rheology, as well as kinetics of metamorphic reactions and mass transfer efficiency. Mineralized veins, ubiquitous in metamorphic rocks, can be seen as preserved witnesses of fluid and mass redistribution that partly accommodate the rock deformation (lateral segregation). However, the driving forces and mechanisms of mass transfer towards fluid-filled open spaces remain somewhat unclear. The aim of this study is to investigate the vein-forming processes and the modalities of mass transfer during local fluid-rock interactions, and their links with fluid production and rock deformation, with new insights from Laser Induced Breakdown Spectroscopy (LIBS) profiles. This study focuses on karstic pockets (metre scale) of Triassic metabauxites embedded in thick carbonate units, that have been isolated from large-scale fluid flow during HP-LT Alpine metamorphism (W. Vanoise, French Alps). These rocks display several generations of metamorphic veins containing various Al-bearing minerals, which give particular insights into mass transfer processes. It is proposed that the internally-derived fluid (~13 vol% produced by successive dehydration reactions) has promoted the opening of fluid-filled open spaces (euhedral habits of vein minerals) and served as medium for diffusive mass transfer from rock to vein. Based on mineralogical and textural features, two vein types can be distinguished: (1) some veins are filled with newly formed products of either prograde (chloritoid) or retrograde (chlorite) metamorphic reactions; in this case, fluid-filled open spaces seem to offer energetically favourable nucleation/growth sites; (2) the second vein type is filled with cookeite (Li-Al-rich chlorite) or pyrophyllite, which were present in the host-rock prior to the vein formation. In this closed chemical system, mass transfer from rock to vein was achieved through the fluid, in a dissolution-transport-precipitation process. To investigate the modalities of mass transfer towards this second vein type, LIBS profiles were performed in the host-rock, taking Li concentration as a proxy for cookeite distribution. Cookeite is highly concentrated (45-65 vol%) in regularly spaced veins, and the LIBS profiles show that cookeite is evenly distributed in the host-rock comprised between two veins. The absence of diffusion profiles suggests that the characteristic diffusion distance for Li, Al and Si is greater than or equal to the distance separating two cookeite veins (2-4 cm). This is in agreement with characteristic diffusion lengths calculated from both grain boundary and pore fluid diffusion coefficients, for the estimated duration of the peak of metamorphism. Which driving forces are responsible for cookeite selective transfer towards veins? Chemical potential gradients between host-rock pores and veins may have developed in response to either (1) a stress difference: thermochemical calculations show that pressure-solution processes may affect preferentially cookeite and pyrophyllite; (2) a difference in interfacial energy, phyllosilicates showing very different morphologies in host-rocks (fibers) compared to veins (euhedral crystals); fluid-mineral interfacial energy may be maximal in the small host-rock pores, which can maintain higher cookeite solubility than large fluid-filled open spaces (i.e., veins).

Wave propagation through elastic porous media containing two immiscible fluids

NASA Astrophysics Data System (ADS)

Lo, Wei-Cheng; Sposito, Garrison; Majer, Ernest

2005-02-01

Acoustic wave phenomena in porous media containing multiphase fluids have received considerable attention in recent years because of an increasing scientific awareness of poroelastic behavior in groundwater aquifers. To improve quantitative understanding of these phenomena, a general set of coupled partial differential equations was derived to describe dilatational wave propagation through an elastic porous medium permeated by two immiscible fluids. These equations, from which previous models of dilatational wave propagation can be recovered as special cases, incorporate both inertial coupling and viscous drag in an Eulerian frame of reference. Two important poroelasticity concepts, the linearized increment of fluid content and the closure relation for porosity change, originally defined for an elastic porous medium containing a single fluid, also are generalized for a two-fluid system. To examine the impact of relative fluid saturation and wave excitation frequency (50, 100, 150, and 200 Hz) on free dilatational wave behavior in unconsolidated porous media, numerical simulations of the three possible modes of wave motion were conducted for Columbia fine sandy loam containing either an air-water or oil-water mixture. The results showed that the propagating (P1) mode, which results from in-phase motions of the solid framework and the two pore fluids, moves with a speed equal to the square root of the ratio of an effective bulk modulus to an effective density of the fluid-containing porous medium, regardless of fluid saturation and for both fluid mixtures. The nature of the pore fluids exerts a significant influence on the attenuation of the P1 wave. In the air-water system, attenuation was controlled by material density differences and the relative mobilities of the pore fluids, whereas in the oil-water system an effective kinematic shear viscosity of the pore fluids was the controlling parameter. On the other hand, the speed and attenuation of the two diffusive modes (P2, resulting from out-of-phase motions of the solid framework and the fluids, and P3, the result of capillary pressure fluctuations) were closely associated with an effective dynamic shear viscosity of the pore fluids. The P2 and P3 waves also had the same constant value of the quality factor, and by comparison of our results with previous research on these two dilatational wave modes in sandstones, both were found to be sensitive to the state of consolidation of the porous medium.

Peritoneal Fluid Transport rather than Peritoneal Solute Transport Associates with Dialysis Vintage and Age of Peritoneal Dialysis Patients

PubMed Central

Waniewski, Jacek; Antosiewicz, Stefan; Baczynski, Daniel; Poleszczuk, Jan; Pietribiasi, Mauro; Lindholm, Bengt; Wankowicz, Zofia

2016-01-01

During peritoneal dialysis (PD), the peritoneal membrane undergoes ageing processes that affect its function. Here we analyzed associations of patient age and dialysis vintage with parameters of peritoneal transport of fluid and solutes, directly measured and estimated based on the pore model, for individual patients. Thirty-three patients (15 females; age 60 (21–87) years; median time on PD 19 (3–100) months) underwent sequential peritoneal equilibration test. Dialysis vintage and patient age did not correlate. Estimation of parameters of the two-pore model of peritoneal transport was performed. The estimated fluid transport parameters, including hydraulic permeability (LpS), fraction of ultrasmall pores (α u), osmotic conductance for glucose (OCG), and peritoneal absorption, were generally independent of solute transport parameters (diffusive mass transport parameters). Fluid transport parameters correlated whereas transport parameters for small solutes and proteins did not correlate with dialysis vintage and patient age. Although LpS and OCG were lower for older patients and those with long dialysis vintage, αu was higher. Thus, fluid transport parameters—rather than solute transport parameters—are linked to dialysis vintage and patient age and should therefore be included when monitoring processes linked to ageing of the peritoneal membrane. PMID:26989432

Fluctuation of Ultrafiltration Coefficient of Hemodialysis Membrane During Reuse

NASA Astrophysics Data System (ADS)

Arif, Idam; Christin

2010-12-01

Hemodialysis treatment for patient with kidney failure is to regulate body fluid and to excrete waste products of metabolism. The patient blood and the dialyzing solution (dialysate) are flowed counter currently in a dialyzer to allow volume flux of fluid and diffusion of solutes from the blood to the dialysate through a semipermiable membrane. The volume flux of fluid depends on the hydrostatic and the osmotic pressure difference between the blood and the dialysate. It also depends on the membrane parameter that represents how the membrane allows the fluid and the solutes to move across as a result of the pressure difference, known as the ultrafiltration coefficient Kuf. The coefficient depends on the number and the radius of membrane pores for the movement of the fluids and the solutes across the membrane. The measured membrane ultrafiltration coefficient of reused dialyzer shows fluctuation between one uses to another without any significant trend of change. This indicates that the cleaning process carried out before reuse does not cause perfect removal of clots that happen in the previous use. Therefore the unblocked pores are forced to work hardly to obtain targeted volume flux in a certain time of treatment. This may increase the unblocked pore radius. Reuse is stopped when there is indication of blood leakage during the hemodialysis treatment.

Three-Dimensional Modeling of Fluid and Heat Transport in an Accretionary Complex

NASA Astrophysics Data System (ADS)

Paula, C. A.; Ge, S.; Screaton, E. J.

2001-12-01

As sediments are scraped off of the subducting oceanic crust and accreted to the overriding plate, the rapid loading causes pore pressures in the underthrust sediments to increase. The change in pore pressure drives fluid flow and heat transport within the accretionary complex. Fluid is channeled along higher permeability faults and fractures and expelled at the seafloor. In this investigation, we examined the effects of sediment loading on fluid flow and thermal transport in the decollement at the Barbados Ridge subduction zone. Both the width and thickness of the Barbados Ridge accretionary complex increase from north to south. The presence of mud diapers south of the Tiburon Rise and an observed southward decrease in heat flow measurements indicate that the increased thickness of the southern Barbados accretionary prism affects the transport of chemicals and heat by fluids. The three-dimensional geometry and physical properties of the accretionary complex were utilized to construct a three-dimensional fluid flow/heat transport model. We calculated the pore pressure change due to a period of sediment loading and added this to steady-state pressure conditions to generate initial conditions for transient simulations. We then examined the diffusion of pore pressure and possible perturbation of the thermal regime over time due to loading of the underthrust sediments. The model results show that the sediment-loading event was sufficient to create small temperature fluctuations in the decollement zone. The magnitude of temperature fluctuation in the decollement was greatest at the deformation front but did not vary significantly from north to south of the Tiburon Rise.

The Capacity for Compaction Weakening in Fault Gouge in Nature and Experiment

NASA Astrophysics Data System (ADS)

Faulkner, D.; Boulton, C. J.; Sanchez Roa, C.; Den Hartog, S. A. M.; Bedford, J. D.

2017-12-01

As faults form in low permeability rocks, the compaction of fault gouge can lead to significant pore-fluid pressure increases. The pore pressure increase results from the collapse of the porosity through shear-enhanced compaction and the low hydraulic diffusivity of the gouge that inhibits fluid flow. In experiments, the frictional properties of clay-bearing fault gouges are significantly affected by the development of locally high pore-fluid pressures when compaction rates are high due to fast displacement rates or slip in underconsolidated materials. We show how the coefficient of friction of fault gouges sheared at different slip velocities can be explained with a numerical model that is constrained by laboratory measurements of contemporaneous changes in permeability and porosity. In nature, for compaction weakening to play an important role in earthquake nucleation (and rupture propagation), a mechanism is required to reset the porosity, i.e., maintain underconsolidated gouge along the fault plane. We use the observations of structures within the principal slip zone of the Alpine Fault in New Zealand to suggest that cyclic fluidization of the gouge occurs during coseismic slip, thereby resetting the gouge porosity prior to the next seismic event. Results from confined laboratory rotary shear measurements at elevated slip rates appear to support the hypothesis that fluidization leads to underconsolidation and, thus, to potential weakening by shear-enhanced compaction-induced pore-fluid pressurization.

FEM modeling of postseismic deformation of poroelastic material

NASA Astrophysics Data System (ADS)

Kawamoto, S.; Ito, T.; Hirahara, K.

2004-12-01

Following a large earthquake, postseismic deformation in the focal region has been observed by GPS, leveling measurements and the other geodetic measurements. To explain the postseismic deformation, researchers have proposed and well investigated two physical mechanisms of afterslip and viscoelastic relaxation. In some cases, however, there have been observed postseismic deformation which can not be explained by these mechanisms. Therefore, another mechanism has been proposed, where the crust is treated as "poroelastic material". This concept is called "poroelasticity". In this concept, postseismic deformation is caused by pore fluid flow due to the coseismic stress redistribution. We explored, therefore, the postseismic deformation due to pore fluid flow in a poroelastic material using finite element method (FEM), which can easily handle lateral variations of hydraulic diffusivity and elastic or plastic property. We used the FEM program 'CAMBIOT3D' originally developed by Geotech. Lab. Gunma University, Japan (2003). Because this program was developed for soil mechanics, we must have modified so as to calculate deformation due to earthquake faulting. We implemented the 'split node technique' (Melosh and Refsky, 1981) to calculate the coseismic deformation. In addition to this, we modified the program to calculate the deformation taking into account the Skempton's B. This coefficient B determines what fraction of the coseismic stress due to an earthquake is allotted to pore pressure. Without Skempton's B, coseismic pore pressure becomes too large and hence postseismic deformation is calculated too large. We evaluated the postseismic deformation in a poroelastic material to show that the poroelastic deformation is quite different from that of afterslip and viscoelastic relaxation models. In this presentation, we show the postseismic deformation due to pore fluids flow in a poroelastic material and the effect of Skempton's B. Especially, we discuss what different pattern of postseismic deformation is produced depending on the lateral variation of hydraulic diffusivity structures in and around the fault zone, which structures have been differently inferred from fault zone core sampling researches and so on.

Combined Neutron and X-Ray Radiographic/Tomographic Analysis of Dissolution Limestones under Acidic Conditions

NASA Astrophysics Data System (ADS)

Anovitz, L. M.; Cole, D. R.; Hussey, D. S.; LaManna, J.; Swift, A.; Jacobson, D. L.

2016-12-01

Carbon dioxide capture and sequestration in deep geological formations is an important option for reducing greenhouse gas emissions. While the importance of porosity and pore-evolution has long been recognized, the evolution of porosity and permeability in reactive carbonates exposed to CO2-loaded brines is not well constrained. A typical pH range for CO2-acidified brine is 3 to 4.5 depending on alkalinity. This represents a substantial perturbation of typical brines that range from pH 6 to 8. The key questions include how accessible are the pores to fluid transport and how does the pore network evolve as the matrix reacts with the acidic solution? Limestones and dolostones contain nano- to macroscale porosity comprised of cracks, grain boundaries, fluid inclusions, single pores, vugs and networks of pores of random shapes and orientations. Accessible, interconnected pores may act as pore throats, constraining overall flow and are the most likely locations for extensive rock alteration. Neutron imaging is well suited to interrogation of fluid flow in porous media. Because of the large scattering cross section of hydrogen it can be used to directly image water or hydrocarbons without an added contrast medium that might modify interfacial tension and fluid/fluid interactions. In order to understand the reaction of acidified fluids we used simultaneous neutron and X-ray tomography to study the uptake and reaction of water and an acidic fluid (pH 1 HCl) with two types of Indiana limestone, one with a permeability of 2-4 mD, and the other 70 mD. One set of experiments explored capillary uptake in a dry core. These documented rapid uptake and CO2 bubble formation at the inlet. A second set introduced at a constant forced flow rate of 10 ml/min. Both core types exhibited wormhole formation, but the low perm limestone wormhole consisted of one well-delineated channel with a few side "tributaries," whereas the high perm core exhibited a more diffuse array of channels. Post-flow neutron and X-ray tomography showed that grain boundaries and other initial parts of the porous network play an important role in controlling the dissolution process. Neutron radiography and tomography have the potential to significantly advance our understanding of fluid flow and reactive behavior relevant to a wide variety of subsurface applications.

Dehydration reactions, mass transfer and rock deformation relationships during subduction of Alpine metabauxites: insights from LIBS compositional profiles between metamorphic veins

NASA Astrophysics Data System (ADS)

Verlaguet, A.; Brunet, F.; Goffe, B.; Menut, D.; Findling, N.; Poinssot, C.

2011-12-01

In subduction zones, the significant amounts of aqueous fluid released in the course of the successive dehydration reactions occurring during prograde metamorphism are expected to strongly influence the rock rheology, as well as kinetics of metamorphic reactions and mass transfer efficiency. Mineralized veins, ubiquitous in metamorphic rocks, can be seen as preserved witnesses of fluid and mass redistribution that partly accommodate the rock deformation (lateral segregation). However, the driving forces and mechanisms of mass transfer towards fluid-filled open spaces remain somewhat unclear. The aim of this study is to investigate the modalities of mass transfer during local fluid-rock interactions, and their links with fluid production and rock deformation. This study focuses on karstic pockets (metre scale) of Triassic metabauxites embedded in thick carbonate units, that have been isolated from large-scale fluid flow during HP-LT Alpine metamorphism (W. Vanoise, French Alps). These rocks display several generations of metamorphic veins containing various Al-bearing minerals, which give particular insights into mass transfer processes. It is proposed that the internally-derived fluid (~13 vol% produced by successive dehydration reactions) has promoted the opening of fluid-filled open spaces (euhedral habits of vein minerals) and served as medium for diffusive mass transfer from rock to vein. Based on mineralogical and textural features, two vein types can be distinguished: (1) some veins are filled with newly formed products of either prograde (chloritoid) or retrograde (chlorite) metamorphic reactions; in this case, fluid-filled open spaces seem to offer energetically favourable nucleation/growth sites; (2) the second vein type is filled with cookeite (Li-Al-rich chlorite) or pyrophyllite, that were present in the host rock prior to the vein formation. In this closed chemical system, mass transfer from rock to vein was achieved through the fluid, in a dissolution-transport-precipitation process, possibly stress-assisted. Cookeite is highly concentrated (40-70 vol%) in regularly spaced veins. Laser Induced Breakdown Spectroscopy profiles show that cookeite is evenly distributed in the rock matrix comprised between two veins. The absence of diffusion profiles suggests that the characteristic diffusion length for Li, Al and Si is greater than or equal to the distance separating two cookeite veins (3-6 cm). This is in agreement with characteristic diffusion lengths calculated from both grain boundary and pore fluid diffusion coefficients, for the estimated duration of the peak of metamorphism. Phyllosilicates have very different morphologies in the rock matrix (fibers) compared to veins (euhedral crystals): fluid-mineral interfacial energy may be maximal in the small matrix pores, which can maintain higher cookeite solubility than in fluid-filled open spaces. Therefore, as soon as veins open, chemical potential gradients may develop and drive cookeite transfer from rock matrix to veins.

Harvesting liquid from unsaturated vapor - nanoflows induced by capillary condensation

NASA Astrophysics Data System (ADS)

Vincent, Olivier; Marguet, Bastien; Stroock, Abraham

2016-11-01

A vapor, even subsaturated, can spontaneously form liquid in nanoscale spaces. This process, known as capillary condensation, plays a fundamental role in various contexts, such as the formation of clouds or the dynamics of hydrocarbons in the geological subsurface. However, large uncertainties remain on the thermodynamics and fluid mechanics of the phenomenon, due to experimental challenges as well as outstanding questions about the validity of macroscale physics at the nanometer scale. We studied experimentally the spatio-temporal dynamics of water condensation in a model nanoporous medium (pore radius 2 nm), taking advantage of the color change of the material upon hydration. We found that at low relative humidities (< 60 % RH), capillary condensation progressed in a diffusive fashion, while it occurred through a well-defined capillary-viscous imbibition front at > 60 % RH, driven by a balance between the pore capillary pressure and the condensation stress given by Kelvin equation. Further analyzing the imbibition dynamics as a function of saturation allowed us to extract detailed information about the physics of nano-confined fluids. Our results suggest excellent extension of macroscale fluid dynamics and thermodynamics even in pores 10 molecules in diameter.

Evaluating the hydraulic and transport properties of peat soil using pore network modeling and X-ray micro computed tomography

NASA Astrophysics Data System (ADS)

Gharedaghloo, Behrad; Price, Jonathan S.; Rezanezhad, Fereidoun; Quinton, William L.

2018-06-01

Micro-scale properties of peat pore space and their influence on hydraulic and transport properties of peat soils have been given little attention so far. Characterizing the variation of these properties in a peat profile can increase our knowledge on the processes controlling contaminant transport through peatlands. As opposed to the common macro-scale (or bulk) representation of groundwater flow and transport processes, a pore network model (PNM) simulates flow and transport processes within individual pores. Here, a pore network modeling code capable of simulating advective and diffusive transport processes through a 3D unstructured pore network was developed; its predictive performance was evaluated by comparing its results to empirical values and to the results of computational fluid dynamics (CFD) simulations. This is the first time that peat pore networks have been extracted from X-ray micro-computed tomography (μCT) images of peat deposits and peat pore characteristics evaluated in a 3D approach. Water flow and solute transport were modeled in the unstructured pore networks mapped directly from μCT images. The modeling results were processed to determine the bulk properties of peat deposits. Results portray the commonly observed decrease in hydraulic conductivity with depth, which was attributed to the reduction of pore radius and increase in pore tortuosity. The increase in pore tortuosity with depth was associated with more decomposed peat soil and decreasing pore coordination number with depth, which extended the flow path of fluid particles. Results also revealed that hydraulic conductivity is isotropic locally, but becomes anisotropic after upscaling to core-scale; this suggests the anisotropy of peat hydraulic conductivity observed in core-scale and field-scale is due to the strong heterogeneity in the vertical dimension that is imposed by the layered structure of peat soils. Transport simulations revealed that for a given solute, the effective diffusion coefficient decreases with depth due to the corresponding increase of diffusional tortuosity. Longitudinal dispersivity of peat also was computed by analyzing advective-dominant transport simulations that showed peat dispersivity is similar to the empirical values reported in the same peat soil; it is not sensitive to soil depth and does not vary much along the soil profile.

Three-dimensional numerical simulation of water droplet emerging from a gas diffusion layer surface in micro-channels

NASA Astrophysics Data System (ADS)

Ding, Y.; Bi, H. T.; Wilkinson, D. P.

The dynamic formation of water droplets emerging from a gas diffusion layer (GDL) surface in micro-channels was simulated using the volume of fluid (VOF) method. The influence of GDL surface microstructure was investigated by changing the pore diameter and the number of pore openings on the GDL surface. Simulation results show that the microstructure of the GDL surface has a significant impact on the two-phase flow patterns in gas flow channels. For a non-uniform GDL surface, three stages were identified, namely emergence and merging on the GDL surface, accumulation on the channel sidewalls and detachment from the top wall. It was also found that if the pore size is small enough, the flow pattern in the channel does not change with further reduction in the pore diameter. However, the two-phase flow patterns change significantly with the wettability of the GDL surface and sidewalls, but remain the same when the liquid flow rate is reduced by two orders of magnitude from the reference case.

Investigating the relationship between seismicity and fluid injection in the Barnett Shale, Texas using coupled poroelastic model and surface deformation data

NASA Astrophysics Data System (ADS)

Zhai, G.; Shirzaei, M.

2017-12-01

Across the Barnett Shale, Texas a noticeable increase in seismic activity was observed during 2007 and 2015, which was accompanied by high volume injection at several nearby disposal wells. Many studies focused on the positive correlation between injection rate at individual wells and the adjacent seismicity, suggesting that seismicity is triggered or induced due to increased pore fluid pressure associated with fluid injection in hydraulically connected geological units. However, investigating temporal evolution of total volume of injected fluid and concurrent earthquakes in a larger area indicates more complex patterns, requiring a more comprehensive analysis of the spatiotemporal evolution of coupled poroelastic stress and pore fluid pressure. In this study, we created a coupled poroelastic model to simulate large scale spatiotemporal evolution of pore pressure, poroelastic stresses, and Coulomb failure stress in the Barnett Shale using injection time series of 96 high-volume injection wells spanning from 2007 to 2015. We additionally account for a layered poroelastic medium, where its parameters are set up using geological maps and seismic tomographic data sets. Fault orientations and relevant frictional properties are also extracted from published literatures. We further integrate observation of surface deformation obtained from interferometric processing of 16 ALOS L-Band SAR images to optimize rock hydraulic diffusivity and constrain the extent to which fluid may migrate. The preliminary modeling result shows that poroelastic stress is only 10% of pore pressure. However, the superimposition of these two effects is spatially and temporally responsible for the occurrence of earthquakes in the Barnett Shale. Also, not all area with increased Coulomb failure stress experiences elevated seismicity, suggesting possible heterogeneous background tectonic stresses, lacking pre-existing faults, and/or heterogeneous fault orientations.

Characterization of double continuum formulations of transport through pore-scale information

NASA Astrophysics Data System (ADS)

Porta, G.; Ceriotti, G.; Bijeljic, B.

2016-12-01

Information on pore-scale characteristics is becoming increasingly available at unprecedented levels of detail from modern visualization/data-acquisition techniques. These advancements are not completely matched by corresponding developments of operational procedures according to which we can engineer theoretical findings aiming at improving our ability to reduce the uncertainty associated with the outputs of continuum-scale models to be employed at large scales. We present here a modeling approach which rests on pore-scale information to achieve a complete characterization of a double continuum model of transport and fluid-fluid reactive processes. Our model makes full use of pore-scale velocity distributions to identify mobile and immobile regions. We do so on the basis of a pointwise (in the pore space) evaluation of the relative strength of advection and diffusion time scales, as rendered by spatially variable values of local Péclet numbers. After mobile and immobile regions are demarcated, we build a simplified unit cell which is employed as a representative proxy of the real porous domain. This model geometry is then employed to simplify the computation of the effective parameters embedded in the double continuum transport model, while retaining relevant information from the pore-scale characterization of the geometry and velocity field. We document results which illustrate the applicability of the methodology to predict transport of a passive tracer within two- and three-dimensional media upon comparison with direct pore-scale numerical simulation of transport in the same geometrical settings. We also show preliminary results about the extension of this model to fluid-fluid reactive transport processes. In this context, we focus on results obtained in two-dimensional porous systems. We discuss the impact of critical quantities required as input to our modeling approach to obtain continuum-scale outputs. We identify the key limitations of the proposed methodology and discuss its capability also in comparison with alternative approaches grounded, e.g., on nonlocal and particle-based approximations.

The Effect of the Pore Entrance on Particle Motion in Slit Pores: Implications for Ultrathin Membranes

PubMed Central

Delavari, Armin; Baltus, Ruth

2017-01-01

Membrane rejection models generally neglect the effect of the pore entrance on intrapore particle transport. However, entrance effects are expected to be particularly important with ultrathin membranes, where membrane thickness is typically comparable to pore size. In this work, a 2D model was developed to simulate particle motion for spherical particles moving at small Re and infinite Pe from the reservoir outside the pore into a slit pore. Using a finite element method, particles were tracked as they accelerated across the pore entrance until they reached a steady velocity in the pore. The axial position in the pore where particle motion becomes steady is defined as the particle entrance length (PEL). PELs were found to be comparable to the fluid entrance length, larger than the pore size and larger than the thickness typical of many ultrathin membranes. Results also show that, in the absence of particle diffusion, hydrodynamic particle–membrane interactions at the pore mouth result in particle “funneling” in the pore, yielding cross-pore particle concentration profiles focused at the pore centerline. The implications of these phenomena on rejection from ultrathin membranes are examined. PMID:28796197

A binomial modeling approach for upscaling colloid transport under unfavorable conditions: Emergent prediction of extended tailing

NASA Astrophysics Data System (ADS)

Hilpert, Markus; Rasmuson, Anna; Johnson, William P.

2017-07-01

Colloid transport in saturated porous media is significantly influenced by colloidal interactions with grain surfaces. Near-surface fluid domain colloids experience relatively low fluid drag and relatively strong colloidal forces that slow their downgradient translation relative to colloids in bulk fluid. Near-surface fluid domain colloids may reenter into the bulk fluid via diffusion (nanoparticles) or expulsion at rear flow stagnation zones, they may immobilize (attach) via primary minimum interactions, or they may move along a grain-to-grain contact to the near-surface fluid domain of an adjacent grain. We introduce a simple model that accounts for all possible permutations of mass transfer within a dual pore and grain network. The primary phenomena thereby represented in the model are mass transfer of colloids between the bulk and near-surface fluid domains and immobilization. Colloid movement is described by a Markov chain, i.e., a sequence of trials in a 1-D network of unit cells, which contain a pore and a grain. Using combinatorial analysis, which utilizes the binomial coefficient, we derive the residence time distribution, i.e., an inventory of the discrete colloid travel times through the network and of their probabilities to occur. To parameterize the network model, we performed mechanistic pore-scale simulations in a single unit cell that determined the likelihoods and timescales associated with the above colloid mass transfer processes. We found that intergrain transport of colloids in the near-surface fluid domain can cause extended tailing, which has traditionally been attributed to hydrodynamic dispersion emanating from flow tortuosity of solute trajectories.

Sequential peritoneal equilibration test: a new method for assessment and modelling of peritoneal transport.

PubMed

Galach, Magda; Antosiewicz, Stefan; Baczynski, Daniel; Wankowicz, Zofia; Waniewski, Jacek

2013-02-01

In spite of many peritoneal tests proposed, there is still a need for a simple and reliable new approach for deriving detailed information about peritoneal membrane characteristics, especially those related to fluid transport. The sequential peritoneal equilibration test (sPET) that includes PET (glucose 2.27%, 4 h) followed by miniPET (glucose 3.86%, 1 h) was performed in 27 stable continuous ambulatory peritoneal dialysis patients. Ultrafiltration volumes, glucose absorption, ratio of concentration in dialysis fluid to concentration in plasma (D/P), sodium dip (Dip D/P Sodium), free water fraction (FWF60) and the ultrafiltration passing through small pores at 60 min (UFSP60), were calculated using clinical data. Peritoneal transport parameters were estimated using the three-pore model (3p model) and clinical data. Osmotic conductance for glucose was calculated from the parameters of the model. D/P creatinine correlated with diffusive mass transport parameters for all considered solutes, but not with fluid transport characteristics. Hydraulic permeability (L(p)S) correlated with net ultrafiltration from miniPET, UFSP60, FWF60 and sodium dip. The fraction of ultrasmall pores correlated with FWF60 and sodium dip. The sequential PET described and interpreted mechanisms of ultrafiltration and solute transport. Fluid transport parameters from the 3p model were independent of the PET D/P creatinine, but correlated with fluid transport characteristics from PET and miniPET.

Fluid pressure and fault strength: insights from load-controlled experiments on carbonate-bearing rocks

NASA Astrophysics Data System (ADS)

Spagnuolo, E.; Violay, M.; Nielsen, S. B.; Di Toro, G.

2013-12-01

Fluid pressure Pf has been indicated as a major factor controlling natural (e.g., L'Aquila, Italy, 2009 Mw 6.3) and induced seismicity (e.g., Wilzetta, Oklahoma, 2011 Mw 5.7). The Terzaghi's principle states that the effective normal stress σeff= σn (1- α Pf ), with α the Biot coefficient and σn the normal stress, is reduced in proportion to Pf. A value of α=1 is often used by default; however, within a complex fault core of inhomogeneous permeability, α may vary in a yet poorly understood way. To shed light on this problem, we conducted experiments on carbonate-bearing rock samples (Carrara marble) in room humidity conditions and in the presence of pore fluids (drained conditions), where a pre-cut fault is loaded by shear stress τ in a rotary apparatus (SHIVA) under constant σn=15 MPa. Two types of tests were performed with fluids: (1) the fluid pressure was kept constant at Pf=5 MPa (close to hydrostatic conditions at a depth of 0.5 km) and the fault was driven to failure instability by gradually increasing τ; (2) the fluid pressure was kept at Pf=5 MPa and τ was increased until close to instability (τ = 7 MPa): at this point Pf was raised of 0.5 MPa every 10 s up to Pf =10 MPa to induce a main (failure) instability. Assuming α=1 and an effective peak strength (τp)eff=μp σeff at failure, the experiments reveal that: 1) (τp)eff is sensitive to the shear loading rate: fast loading rates (0.5 MPa every 20 s) induce higher peak shear-stress values than slow loading rates (0.5 MPa every 40 s). Such effect is not observed (minor or inexistent) in the absence of pore fluids. 2) Under fast loading rates the (τp)eff may surpass that measured in the absence of pore fluids under identical effective normal stress σeff. 3) An increase of Pf does not necessarily induce the main instability (within the time intervals studied here, i.e. up to ~10 s) even if the effective strength threshold is largely surpassed (e.g., (τp)eff=1.3 μp σeff). We interpret these results in terms of limited permeability of the fault slip zone which reduces α. Indeed result (3) may indicate that a Pf increase did not rapidly penetrate the slip zone because a seal (thin layer of wet ultrafine calcite gouge) formed during the slip preceding the main instability. On the other hand, shearing of the slip zone probably induces dilation (not measured because below resolution) in the slip zone and results in a decrease in pore pressure. Again, due to limited permeability, the drop in pore pressure within the slip zone does not have time to re-equilibrate with the imposed Pf, provided that the hold time is short (20 s) with respect to the diffusion time, but it may re-equilibrate under longer hold times (40 s). As a consequence the Biot coefficient depends on the time interval of observation, with α~0 at short time periods and α~1 at long time periods. This yields an approximate hydraulic diffusivity κ~10-8 m2 s-1 using κ=l2/td with the half length of the contact surface l=5 mm and td=30 s. Such diffusivity is compatible, for example, with a low porosity shale.

Fast intraslab fluid-flow events linked to pulses of high pore fluid pressure at the subducted plate interface

NASA Astrophysics Data System (ADS)

Taetz, Stephan; John, Timm; Bröcker, Michael; Spandler, Carl; Stracke, Andreas

2018-01-01

A better understanding of the subduction zone fluid cycle and its chemical-mechanical feedback requires in-depth knowledge about how fluids flow within and out of descending slabs. Relicts of fluid-flow systems in exhumed rocks of fossil subduction zones allow for identification of the general relationships between dehydration reactions, fluid pathway formation, the dimensions and timescales of distinct fluid flow events; all of which are required for quantitative models for fluid-induced subduction zone processes. Two types of garnet-quartz-phengite veins can be distinguished in an eclogite-facies mélange block from the Pouébo Eclogite Mélange, New Caledonia. These veins record synmetamorphic internal fluid release by mineral breakdown reactions (type I veins), and infiltration of an external fluid (type II veins) with the associated formation of a reaction selvage. The dehydration and fluid migration documented by the type I veins likely occurred on a timescale of 105-106 years, based on average subduction rates and metamorphic conditions required for mineral dehydration and fluid flow. The timeframe of fluid-rock interaction between the external fluid and the wall-rock of the type II veins is quantified using a continuous bulk-rock Li-diffusion profile perpendicular to a vein and its metasomatic selvage. Differences in Li concentration between the internal and external fluid reservoirs resulted in a distinct diffusion profile (decreasing Li concentration and increasing δ7 Li) as the reaction front propagated into the host rock. Li-chronometric constraints indicate that the timescales of fluid-rock interaction associated with type II vein formation are on the order of 1 to 4 months (0.150-0.08+0.14 years). The short-lived, pulse-like character of this process is consistent with the notion that fluid flow caused by oceanic crust dehydration at the blueschist-to-eclogite transition contributes to or even dominates episodic pore fluid pressure increases at the plate interface, which in turn, may trigger slip events reported from many subduction zones.

Robust determination of surface relaxivity from nuclear magnetic resonance DT2 measurements

NASA Astrophysics Data System (ADS)

Luo, Zhi-Xiang; Paulsen, Jeffrey; Song, Yi-Qiao

2015-10-01

Nuclear magnetic resonance (NMR) is a powerful tool to probe into geological materials such as hydrocarbon reservoir rocks and groundwater aquifers. It is unique in its ability to obtain in situ the fluid type and the pore size distributions (PSD). The T1 and T2 relaxation times are closely related to the pore geometry through the parameter called surface relaxivity. This parameter is critical for converting the relaxation time distribution into the PSD and so is key to accurately predicting permeability. The conventional way to determine the surface relaxivity ρ2 had required independent laboratory measurements of the pore size. Recently Zielinski et al. proposed a restricted diffusion model to extract the surface relaxivity from the NMR diffusion-T2 relaxation (DT2) measurement. Although this method significantly improved the ability to directly extract surface relaxivity from a pure NMR measurement, there are inconsistencies with their model and it relies on a number of preset parameters. Here we propose an improved signal model to incorporate a scalable LT and extend their method to extract the surface relaxivity based on analyzing multiple DT2 maps with varied diffusion observation time. With multiple diffusion observation times, the apparent diffusion coefficient correctly describes the restricted diffusion behavior in samples with wide PSDs, and the new method does not require predetermined parameters, such as the bulk diffusion coefficient and tortuosity. Laboratory experiments on glass beads packs with the beads diameter ranging from 50 μm to 500 μm are used to validate the new method. The extracted diffusion parameters are consistent with their known values and the determined surface relaxivity ρ2 agrees with the expected value within ±7%. This method is further successfully applied on a Berea sandstone core and yields surface relaxivity ρ2 consistent with the literature.

Fault reactivation by fluid injection considering permeability evolution in fault-bordering damage zones

NASA Astrophysics Data System (ADS)

Yang, Z.; Yehya, A.; Rice, J. R.; Yin, J.

2017-12-01

Earthquakes can be induced by human activity involving fluid injection, e.g., as wastewater disposal from hydrocarbon production. The occurrence of such events is thought to be, mainly, due to the increase in pore pressure, which reduces the effective normal stress and hence the strength of a nearby fault. Change in subsurface stress around suitably oriented faults at near-critical stress states may also contribute. We focus on improving the modeling and prediction of the hydro-mechanical response due to fluid injection, considering the full poroelastic effects and not solely changes in pore pressure in a rigid host. Thus we address the changes in porosity and permeability of the medium due to the changes in the local volumetric strains. Our results also focus on including effects of the fault architecture (low permeability fault core and higher permeability bordering damage zones) on the pressure diffusion and the fault poroelastic response. Field studies of faults have provided a generally common description for the size of their bordering damage zones and how they evolve along their direction of propagation. Empirical laws, from a large number of such observations, describe their fracture density, width, permeability, etc. We use those laws and related data to construct our study cases. We show that the existence of high permeability damage zones facilitates pore-pressure diffusion and, in some cases, results in a sharp increase in pore-pressure at levels much deeper than the injection wells, because these regions act as conduits for fluid pressure changes. This eventually results in higher seismicity rates. By better understanding the mechanisms of nucleation of injection-induced seismicity, and better predicting the hydro-mechanical response of faults, we can assess methodologies and injection strategies to avoid risks of high magnitude seismic events. Microseismic events occurring after the start of injection are very important indications of when injection should be stopped and how to avoid major events. Our work contributes to the assessment or mitigation of seismic hazard and risk, and our long-term target question is: How to not make an earthquake?

An earthquake instability model based on faults containing high fluid-pressure compartments

USGS Publications Warehouse

Lockner, D.A.; Byerlee, J.D.

1995-01-01

It has been proposed that large strike-slip faults such as the San Andreas contain water in seal-bounded compartments. Arguments based on heat flow and stress orientation suggest that in most of the compartments, the water pressure is so high that the average shear strength of the fault is less than 20 MPa. We propose a variation of this basic model in which most of the shear stress on the fault is supported by a small number of compartments where the pore pressure is relatively low. As a result, the fault gouge in these compartments is compacted and lithified and has a high undisturbed strength. When one of these locked regions fails, the system made up of the neighboring high and low pressure compartments can become unstable. Material in the high fluid pressure compartments is initially underconsolidated since the low effective confining pressure has retarded compaction. As these compartments are deformed, fluid pressure remains nearly unchanged so that they offer little resistance to shear. The low pore pressure compartments, however, are overconsolidated and dilate as they are sheared. Decompression of the pore fluid in these compartments lowers fluid pressure, increasing effective normal stress and shear strength. While this effect tends to stabilize the fault, it can be shown that this dilatancy hardening can be more than offset by displacement weakening of the fault (i.e., the drop from peak to residual strength). If the surrounding rock mass is sufficiently compliant to produce an instability, slip will propagate along the fault until the shear fracture runs into a low-stress region. Frictional heating and the accompanying increase in fluid pressure that are suggested to occur during shearing of the fault zone will act as additional destabilizers. However, significant heating occurs only after a finite amount of slip and therefore is more likely to contribute to the energetics of rupture propagation than to the initiation of the instability. We present results of a one-dimensional dynamic Burridge-Knopoff-type model to demonstrate various aspects of the fluid-assisted fault instability described above. In the numerical model, the fault is represented by a series of blocks and springs, with fault rheology expressed by static and dynamic friction. In addition, the fault surface of each block has associated with it pore pressure, porosity and permeability. All of these variables are allowed to evolve with time, resulting in a wide range of phenomena related to fluid diffusion, dilatancy, compaction and heating. These phenomena include creep events, diffusion-controlled precursors, triggered earthquakes, foreshocks, aftershocks, and multiple earthquakes. While the simulations have limitations inherent to 1-D fault models, they demonstrate that the fluid compartment model can, in principle, provide the rich assortment of phenomena that have been associated with earthquakes. ?? 1995 Birkha??user Verlag.

Shock Loading of Granular Ni/Al Composites. Part 1. Mechanics of Loading

DOE PAGES

Cherukara, Mathew J.; Germann, Timothy C.; Kober, Edward M.; ...

2014-10-16

We present molecular dynamics simulations of the thermomechanical response under shock loading of a granular material consisting of laminated Ni/Al grains. We observe two regimes: At low piston velocities (up ≲ 1km/s), the shock wave is diffuse, and the width of the shock front decreases with increasing piston velocity. Beyond a critical shock strength, however, the width remains relatively constant at approximately the mean grain radius. This change in behavior follows from an evolution of the mechanism of compaction with increasing insult strength. Furthermore, the mechanism evolves from plastic deformation-mediated pore collapse for relatively weak shocks, to solid extrusion andmore » fluid ejecta filling pores ahead of the shock front at intermediate strengths, and finally to atomic jetting into the pore for very strong shocks (up ≳ 2 km/s). High-energy fluid ejecta into pores leads to the formation of flow vorticity and can result in a large fraction of the input energy localizing into translational kinetic energy components including the formation of hot spots. This has implications for the mechanical mixing of Ni and Al in these reactive composites.« less

An investigation into the effects of pore connectivity on T2 NMR relaxation

NASA Astrophysics Data System (ADS)

Ghomeshi, Shahin; Kryuchkov, Sergey; Kantzas, Apostolos

2018-04-01

Nuclear Magnetic Resonance (NMR) is a powerful technique used to characterize fluids and flow in porous media. The NMR relaxation curves are closely related to pore geometry, and the inversion of the NMR relaxometry data is known to give useful information with regards to pore size distribution (PSD) through the relative amplitudes of the fluids stored in the small and large pores. While this information is crucial, the main challenge for the successful use of the NMR measurements is the proper interpretation of the measured signals. Natural porous media patterns consist of complex pore structures with many interconnected or "coupled" regions, as well as isolated pores. This connectivity along the throats changes the relaxation distribution and in order to properly interpret this data, a thorough understanding of the effects of pore connectivity on the NMR relaxation distribution is warranted. In this paper we address two main points. The first pertains to the fact that there is a discrepancy between the relaxation distribution obtained from experiments, and the ones obtained from solving the mathematical models of diffusion process in the digitized images of the pore space. There are several reasons that may attribute to this such as the lack of a proper incorporation of surface roughness into the model. However, here we are more interested in the effects of pore connectivity and to understand why the typical NMR relaxation distribution obtained from experiments are wider, while the numerical simulations predict that a wider NMR relaxation distribution may indicate poor connectivity. Secondly, by not taking into account the pore coupling effects, from our experience in interpreting the data, we tend to underestimate the pore volume of small pores and overestimate the amplitudes in the large pores. The role of pore coupling becomes even more prominent in rocks with small pore sizes such as for example in shales, clay in sandstones, and in the microstructures of carbonates.

New Laboratory Observations of Thermal Pressurization Weakening

NASA Astrophysics Data System (ADS)

Badt, N.; Tullis, T. E.; Hirth, G.

2017-12-01

Dynamic frictional weakening due to pore fluid thermal pressurization has been studied under elevated confining pressure in the laboratory, using a rotary-shear apparatus having a sample with independent pore pressure and confining pressure systems. Thermal pressurization is directly controlled by the permeability of the rocks, not only for the initiation of high-speed frictional weakening but also for a subsequent sequence of high-speed sliding events. First, the permeability is evaluated at different effective pressures using a method where the pore pressure drop and the flow-through rate are compared using Darcy's Law as well as a pore fluid oscillation method, the latter method also permitting measurement of the storage capacity. Then, the samples undergo a series of high-speed frictional sliding segments at a velocity of 2.5 mm/s, under an applied confining pressure and normal stress of 45 MPa and 50 MPa, respectively, and an initial pore pressure of 25 MPa. Finally the rock permeability and storage capacity are measured again to assess the evolution of the rock's pore fluid properties. For samples with a permeability of 10-20 m2 thermal pressurization promotes a 40% decrease in strength. However, after a sequence of three high-speed sliding events, the magnitude of weakening diminishes progressively from 40% to 15%. The weakening events coincide with dilation of the sliding interface. Moreover, the decrease in the weakening degree with progressive fast-slip events suggest that the hydraulic diffusivity may increase locally near the sliding interface during thermal pressurization-enhanced slip. This could result from stress- or thermally-induced damage to the host rock, which would perhaps increase both permeability and storage capacity, and so possibly decrease the susceptibility of dynamic weakening due to thermal pressurization in subsequent high-speed sliding events.

Control of optical transport parameters of 'porous medium – supercritical fluid' systems

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zimnyakov, D A; Ushakova, O V; Yuvchenko, S A

2015-11-30

The possibility of controlling optical transport parameters (in particular, transport scattering coefficient) of porous systems based on polymer fibres, saturated with carbon dioxide in different phase states (gaseous, liquid and supercritical) has been experimentally studied. An increase in the pressure of the saturating medium leads to a rise of its refractive index and, correspondingly, the diffuse-transmission coefficient of the system due to the decrease in the transport scattering coefficient. It is shown that, in the case of subcritical saturating carbon dioxide, the small-angle diffuse transmission of probed porous layers at pressures close to the saturated vapour pressure is determined bymore » the effect of capillary condensation in pores. The immersion effect in 'porous medium – supercritical fluid' systems, where the fluid pressure is used as a control parameter, is considered. The results of reconstructing the values of transport scattering coefficient of probed layers for different refractive indices of a saturating fluid are presented. (radiation scattering)« less

Reservoir Condition Pore-scale Imaging of Multiple Fluid Phases Using X-ray Microtomography

PubMed Central

Andrew, Matthew; Bijeljic, Branko; Blunt, Martin

2015-01-01

X-ray microtomography was used to image, at a resolution of 6.6 µm, the pore-scale arrangement of residual carbon dioxide ganglia in the pore-space of a carbonate rock at pressures and temperatures representative of typical formations used for CO2 storage. Chemical equilibrium between the CO2, brine and rock phases was maintained using a high pressure high temperature reactor, replicating conditions far away from the injection site. Fluid flow was controlled using high pressure high temperature syringe pumps. To maintain representative in-situ conditions within the micro-CT scanner a carbon fiber high pressure micro-CT coreholder was used. Diffusive CO2 exchange across the confining sleeve from the pore-space of the rock to the confining fluid was prevented by surrounding the core with a triple wrap of aluminum foil. Reconstructed brine contrast was modeled using a polychromatic x-ray source, and brine composition was chosen to maximize the three phase contrast between the two fluids and the rock. Flexible flow lines were used to reduce forces on the sample during image acquisition, potentially causing unwanted sample motion, a major shortcoming in previous techniques. An internal thermocouple, placed directly adjacent to the rock core, coupled with an external flexible heating wrap and a PID controller was used to maintain a constant temperature within the flow cell. Substantial amounts of CO2 were trapped, with a residual saturation of 0.203 ± 0.013, and the sizes of larger volume ganglia obey power law distributions, consistent with percolation theory. PMID:25741751

Fluid-Rock Characterization and Interactions in NMR Well Logging

DOE Office of Scientific and Technical Information (OSTI.GOV)

George J. Hirasaki; Kishore K. Mohanty

2005-09-05

The objective of this report is to characterize the fluid properties and fluid-rock interactions that are needed for formation evaluation by NMR well logging. The advances made in the understanding of NMR fluid properties are summarized in a chapter written for an AAPG book on NMR well logging. This includes live oils, viscous oils, natural gas mixtures, and the relation between relaxation time and diffusivity. Oil based drilling fluids can have an adverse effect on NMR well logging if it alters the wettability of the formation. The effect of various surfactants on wettability and surface relaxivity are evaluated for silicamore » sand. The relation between the relaxation time and diffusivity distinguishes the response of brine, oil, and gas in a NMR well log. A new NMR pulse sequence in the presence of a field gradient and a new inversion technique enables the T{sub 2} and diffusivity distributions to be displayed as a two-dimensional map. The objectives of pore morphology and rock characterization are to identify vug connectivity by using X-ray CT scan, and to improve NMR permeability correlation. Improved estimation of permeability from NMR response is possible by using estimated tortuosity as a parameter to interpolate between two existing permeability models.« less

A Binomial Modeling Approach for Upscaling Colloid Transport Under Unfavorable Attachment Conditions: Emergent Prediction of Nonmonotonic Retention Profiles

NASA Astrophysics Data System (ADS)

Hilpert, Markus; Johnson, William P.

2018-01-01

We used a recently developed simple mathematical network model to upscale pore-scale colloid transport information determined under unfavorable attachment conditions. Classical log-linear and nonmonotonic retention profiles, both well-reported under favorable and unfavorable attachment conditions, respectively, emerged from our upscaling. The primary attribute of the network is colloid transfer between bulk pore fluid, the near-surface fluid domain (NSFD), and attachment (treated as irreversible). The network model accounts for colloid transfer to the NSFD of downgradient grains and for reentrainment to bulk pore fluid via diffusion or via expulsion at rear flow stagnation zones (RFSZs). The model describes colloid transport by a sequence of random trials in a one-dimensional (1-D) network of Happel cells, which contain a grain and a pore. Using combinatorial analysis that capitalizes on the binomial coefficient, we derived from the pore-scale information the theoretical residence time distribution of colloids in the network. The transition from log-linear to nonmonotonic retention profiles occurs when the conditions underlying classical filtration theory are not fulfilled, i.e., when an NSFD colloid population is maintained. Then, nonmonotonic retention profiles result potentially both for attached and NSFD colloids. The concentration maxima shift downgradient depending on specific parameter choice. The concentration maxima were also shown to shift downgradient temporally (with continued elution) under conditions where attachment is negligible, explaining experimentally observed downgradient transport of retained concentration maxima of adhesion-deficient bacteria. For the case of zero reentrainment, we develop closed-form, analytical expressions for the shape, and the maximum of the colloid retention profile.

Molecular diffusion in disordered interfacial media as probed by pulsed field gradients and nuclear magnetic relaxation dispersion

NASA Astrophysics Data System (ADS)

Levitz, P.; Korb, J.-P.; Bryant, R. G.

1999-10-01

We address the question of probing the fluid dynamics in disordered interfacial media by Pulsed field gradient (PFG) and Magnetic relaxation dispersion (MRD) techniques. We show that the PFG method is useful to separate the effects of morphology from the connectivity in disordered macroporous media. We propose simulations of molecular dynamics and spectral density functions, J(ω), in a reconstructed mesoporous medium for different limiting conditions at the pore surface. An algebraic form is found for J(ω) in presence of a surface diffusion and a local exploration of the pore network. A logarithmic form of J(ω) is found in presence of a pure surface diffusion. We present magnetic relaxation dispersion experiments (MRD) for water and acetone in calibrated mesoporous media to support the main results of our simulations and theories. Nous présentons les avantages respectifs des méthodes de gradients de champs pulsés (PFG) et de relaxation magnétique nucléaire en champs cyclés (MRD) pour sonder la dynamique moléculaire dans les milieux interfaciaux désordonnés. La méthode PFG est utile pour séparer la morphologie et la connectivité dans des milieux macroporeux. Des simulations de diffusion moléculaire et de densité spectrale J(ω) en milieux mésoporeux sont présentées dans différentes conditions limites aux interfaces des pores. Nous trouvons une forme de dispersion algébrique de J(ω) pour une diffusion de surface assistée d'une exploration locale du réseau de pores et une forme logarithmique dans le cas d'une simple diffusion de surface. Les résultats expérimentaux de la méthode MRD pour de l'eau et de l'acétone dans des milieux mésoporeux calibrés supportent les résultats principaux de nos simulations et théories.

Aseismic Slip of a Thin Slab Due to a Fluid Source

NASA Astrophysics Data System (ADS)

Aubin, P. W.; Viesca, R. C.

2017-12-01

We explore the effects of an increase of pore pressure on the frictional interface along the base of a thin slab. The thin slab approximation corresponds to a layer overriding a substrate in which variations along the layer's length occur over distances much greater than the layer thickness. We consider deformation that may be in-plane or anti-plane, but approximately uniform in depth, such that spatial variations of displacement (and hence, slip) occur only along one direction parallel to the interface. Such a thin-sheet model may well represent the deformation of landslides and glacial ice streams, and also serves as a first-pass for fault systems, which, while better represented by elastic half-spaces in frictional contact, nonetheless show qualitatively similar behavior. We consider that the friction coefficient at the layer's interface remains (approximately) constant, and that aseismic slip is initiated by a (line) source of fluid at constant pressure, with one-dimensional diffusion parallel to the interface. As posed, the problem yields a self-similar expansion of slip, whose extent grows proportionally to (α * t)^(1/2) (where α is the hydraulic diffusivity) and can either lag behind or outpace the fluid diffusion front. The problem is controlled by a single parameter, accounting for the friction coefficient and the initial (pre-injection) states of stress and pore pressure. The problem solution consists of the self-similar slip profile and the coefficient of proportionality for the crack-front motion. Within the problem parameter range, two end-member scenarios result: one in which the initial level of shear stress on the interface is close to the value of the pre-injection strength (critically stressed) or another in which fluid pressure is just enough to induce slip (marginally pressurized). For the critically stressed and marginally pressurized cases, the aseismic slip front lies far ahead or far behind, respectively, the fluid diffusion front. We find closed-form solutions for both end-members, and in the former case, via matched asymptotics. These solutions provide a basis to solve the general problem, which we also solve numerically for comparison. The solutions also provide a starting point for examining the progression of slip and locking following the shutoff of the fluid source.

Lattice Boltzmann heat transfer model for permeable voxels

NASA Astrophysics Data System (ADS)

Pereira, Gerald G.; Wu, Bisheng; Ahmed, Shakil

2017-12-01

We develop a gray-scale lattice Boltzmann (LB) model to study fluid flow combined with heat transfer for flow through porous media where voxels may be partially solid (or void). Heat transfer in rocks may lead to deformation, which in turn can modulate the fluid flow and so has significant contribution to rock permeability. The LB temperature field is compared to a finite difference solution of the continuum partial differential equations for fluid flow in a channel. Excellent quantitative agreement is found for both Poiseuille channel flow and Brinkman flow. The LB model is then applied to sample porous media such as packed beds and also more realistic sandstone rock sample, and both the convective and diffusive regimes are recovered when varying the thermal diffusivity. It is found that while the rock permeability can be comparatively small (order milli-Darcy), the temperature field can show significant variation depending on the thermal convection of the fluid. This LB method has significant advantages over other numerical methods such as finite and boundary element methods in dealing with coupled fluid flow and heat transfer in rocks which have irregular and nonsmooth pore spaces.

Dehydration reactions, mass transfer and rock deformation relationships during subduction of Alpine metabauxites: insights from LIBS compositional profiles between metamorphic veins

NASA Astrophysics Data System (ADS)

Verlaguet, Anne; Brunet, Fabrice; Goffé, Bruno; Menut, Denis; Findling, Nathaniel; Poinssot, Christophe

2013-04-01

In subduction zones, the significant amounts of aqueous fluid released in the course of the successive dehydration reactions occurring during prograde metamorphism are expected to strongly influence the rock rheology, as well as kinetics of metamorphic reactions and mass transfer efficiency. Mineralized veins, ubiquitous in metamorphic rocks, can be seen as preserved witnesses of fluid and mass redistribution that partly accommodate the rock deformation (lateral segregation). However, the driving forces and mechanisms of mass transfer towards fluid-filled open spaces remain somewhat unclear. The aim of this study is to investigate the vein-forming processes and the modalities of mass transfer during local fluid-rock interactions, and their links with fluid production and rock deformation, with new insights from Laser Induced Breakdown Spectroscopy (LIBS) profiles. This study focuses on karstic pockets (metre scale) of Triassic metabauxites embedded in thick carbonate units, that have been isolated from large-scale fluid flow during HP-LT Alpine metamorphism (W. Vanoise, French Alps). These rocks display several generations of metamorphic veins containing various Al-bearing minerals, which give particular insights into mass transfer processes. It is proposed that the internally-derived fluid (~13 vol% produced by successive dehydration reactions) has promoted the opening of fluid-filled open spaces (euhedral habits of vein minerals) and served as medium for diffusive mass transfer from rock to vein. Based on mineralogical and textural features, two vein types can be distinguished: (1) some veins are filled with newly formed products of either prograde (chloritoid) or retrograde (chlorite) metamorphic reactions; in this case, fluid-filled open spaces seem to offer energetically favourable nucleation/growth sites; (2) the second vein type is filled with cookeite (Li-Al-rich chlorite) or pyrophyllite, that were present in the host rock prior to the vein formation. In this closed chemical system, mass transfer from rock to vein was achieved through the fluid, in a dissolution-transport-precipitation process, possibly stress-assisted. To investigate the modalities of mass transfer towards this second vein type, LIBS profiles were performed in the rock matrix, taking Li concentration as a proxy for cookeite distribution. Cookeite is highly concentrated (40-70 vol%) in regularly spaced veins, and the LIBS profiles show that cookeite is evenly distributed in the rock matrix comprised between two veins. The absence of diffusion profiles suggests that the characteristic diffusion length for Li, Al and Si is greater than or equal to the distance separating two cookeite veins (3-6 cm). This is in agreement with characteristic diffusion lengths calculated from both grain boundary and pore fluid diffusion coefficients, for the estimated duration of the peak of metamorphism. Concerning mass transfer driving forces, phyllosilicates have very different morphologies in the rock matrix (fibers) compared to veins (euhedral crystals): fluid-mineral interfacial energy may be maximal in the small matrix pores, which can maintain higher cookeite solubility than in fluid-filled open spaces. Therefore, as soon as veins open, chemical potential gradients may develop and drive cookeite transfer from rock matrix to veins.

Sructure and dynamics of fluids in micropous and mesoporous earth and engineered materials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Cole, David R; Mamontov, Eugene; Rother, Gernot

2009-01-01

The behavior of liquids in confined geometries (pores, fractures) typically differs, due to the effects of large internal surfaces and geometri-cal confinement, from their bulk behavior in many ways. Phase transitions (i.e., freezing and capillary condensation), sorption and wetting, and dy-namical properties, including diffusion and relaxation, may be modified, with the strongest changes observed for pores ranging in size from <2 nm to 50 nm the micro- and mesoporous regimes. Important factors influ-encing the structure and dynamics of confined liquids include the average pore size and pore size distribution, the degree of pore interconnection, and the character of the liquid-surfacemore » interaction. While confinement of liq-uids in hydrophobic matrices, such as carbon nanotubes, or near the sur-faces of mixed character, such as many proteins, has also been an area of rapidly growing interest, the confining matrices of interest to earth and ma-terials sciences usually contain oxide structural units and thus are hydro-philic. The pore size distribution and the degree of porosity and inter-connection vary greatly amongst porous matrices. Vycor, xerogels, aerogels, and rocks possess irregular porous structures, whereas mesopor-ous silicas (e.g., SBA-15, MCM-41, MCM-48), zeolites, and layered sys-tems, for instance clays, have high degrees of internal order. The pore type and size may be tailored by means of adjusting the synthesis regimen. In clays, the interlayer distance may depend on the level of hydration. Al-though studied less frequently, matrices such as artificial opals and chry-sotile asbestos represent other interesting examples of ordered porous structures. The properties of neutrons make them an ideal probe for com-paring the properties of bulk fluids with those in confined geometries. In this chapter, we provide a brief review of research performed on liquids confined in materials of interest to the earth and material sciences (silicas, aluminas, zeolites, clays, rocks, etc.), emphasizing those neutron scattering techniques which assess both structural modification and dynamical behav-ior. Quantitative understanding of the complex solid-fluid interactions under different thermodynamic situations will impact both the design of bet-ter substrates for technological applications (e.g., chromatography, fluid capture, storage and release, and heterogeneous catalysis) as well as our fundamental understanding of processes encountered in the environment (i.e., fluid and waste mitigation, carbon sequestration, etc.).« less

The effect of fluids on the frictional behavior of calcite gouge

NASA Astrophysics Data System (ADS)

Rempe, M.; Di Toro, G.; Mitchell, T. M.; Hirose, T.; Smith, S. A. F.; Renner, J.

2016-12-01

The presence of fluids in fault zones affects the faults' strength and the nucleation and propagation of earthquakes due to mechanical or physico-chemical weakening effects. To better understand the effect of pore fluids on the frictional behavior of gouge-bearing faults, a series of intermediate- to high-velocity experiments was conducted using the Phv rotary-shear apparatus (Kochi Core Center, Japan) equipped with a servo-controlled pore-fluid pressure system. Calcite gouge was sheared up to several meters displacement at room-humidity (dry) and water-saturated conditions. The pore-fluid factor, λ=pf/σn, ranged from 0.15 to 0.7 and the effective normal stress, σn,eff=σn-pf, from 1 to 12 MPa. Sheared samples were analyzed using scanning electron microscopy and Raman spectroscopy. The steady-state shear stress is lower for saturated than for dry gouges sliding at V=1 mm/s, possibly due to higher intergranular lubrication and/or accelerated subcritical crack growth, as evidenced also by the observed higher degree of compaction. At V=1 m/s, dry gouges show a pronounced strengthening phase preceding the onset of dynamic weakening; saturated gouges weaken abruptly. The higher λ, the lower the peak and steady-state shear stress, but -counterintuitively- the less localized deformation. Degree of weakening and localization might be influenced by insufficient drainage at high λ. In undrained experiments, the shear stress is slightly decreased likely due to thermal pressurization of the pore fluid, but the onset of dynamic weakening is not accelerated, indicating that dynamic weakening is due to more efficient mechanisms. For example, amorphous carbon may lubricate the slip surfaces of dry and saturated calcite gouges and cause dynamic weakening, yet Raman spectra only show the presence of disordered carbon on the principal slip surface. Furthermore, the presence of small recrystallized grains suggests that strain accommodation during steady-state slip might occur by non-frictional processes, such as grain-boundary sliding aided by diffusion creep.

Surface uplift and time-dependent seismic hazard due to fluid injection in eastern Texas.

PubMed

Shirzaei, Manoochehr; Ellsworth, William L; Tiampo, Kristy F; González, Pablo J; Manga, Michael

2016-09-23

Observations that unequivocally link seismicity and wastewater injection are scarce. Here we show that wastewater injection in eastern Texas causes uplift, detectable in radar interferometric data up to >8 kilometers from the wells. Using measurements of uplift, reported injection data, and a poroelastic model, we computed the crustal strain and pore pressure. We infer that an increase of >1 megapascal in pore pressure in rocks with low compressibility triggers earthquakes, including the 4.8-moment magnitude event that occurred on 17 May 2012, the largest earthquake recorded in eastern Texas. Seismic activity increased even while injection rates declined, owing to diffusion of pore pressure from earlier periods with higher injection rates. Induced seismicity potential is suppressed where tight confining formations prevent pore pressure from propagating into crystalline basement rocks. Copyright © 2016, American Association for the Advancement of Science.

The Scales Of Chemical Transport During Dehydration Reactions: The Roles Of Diffusion And Fluid Expulsion

NASA Astrophysics Data System (ADS)

Bedford, J. D.; Wheeler, J.; Fusseis, F.; Leclere, H.; Faulkner, D.

2016-12-01

The growth of new minerals in response to disequilibrium is the most fundamental metamorphic process. However, the kinetics controls on metamorphic reactions, including the transport of chemical components to growing minerals, are poorly understood as direct observation has been impossible both in nature and experiments. Dehydration reactions are common during prograde metamorphism and require both the expulsion of H2O and the transport of solid chemical components to the product mineral(s) in order for reaction to proceed. The expulsion of H2O is particularly important in understanding when fluid overpressures might develop, which can in turn lead to sudden failure (earthquakes) if the fluid is unable to drain. Using time-resolved (4D) synchrotron X-ray microtomography we have imaged a complete dehydration reaction and show how chemical transport of both the solid and fluid phases evolves during reaction. The reaction analysed is the dehydration of gypsum to form bassanite and H2O - an analogue for silicate dehydration but with much faster and controllable reaction rate. Like most dehydration reactions, the breakdown of gypsum is associated with a solid volume reduction which generates pore space. This new porosity wraps around growing bassanite grains, producing fluid-filled moats, across which transport of dissolved solutes to the growing grains occurs via diffusion. As moats grow in width, diffusion and hence reaction rate slow down. Each new grain-moat pair evolves in relative isolation from each other indicating that the chemical transport of solid components occurs over relatively short distances. This is despite the overall pore network becoming connected during the early stages of reaction, thus allowing efficient expulsion of excess H2O. A dehydrating system can therefore act as an open system with respect to the fluid phase and in contrast be restricted to the grain-scale with respect to the dissolved chemicals. This is in accord with observations that many rocks appear to have preserved their chemistry during dehydration (apart from the loss of H2O).

Pore fluids and the LGM ocean salinity-Reconsidered

NASA Astrophysics Data System (ADS)

Wunsch, Carl

2016-03-01

Pore fluid chlorinity/salinity data from deep-sea cores related to the salinity maximum of the last glacial maximum (LGM) are analyzed using estimation methods deriving from linear control theory. With conventional diffusion coefficient values and no vertical advection, results show a very strong dependence upon initial conditions at -100 ky. Earlier inferences that the abyssal Southern Ocean was strongly salt-stratified in the LGM with a relatively fresh North Atlantic Ocean are found to be consistent within uncertainties of the salinity determination, which remain of order ±1 g/kg. However, an LGM Southern Ocean abyss with an important relative excess of salt is an assumption, one not required by existing core data. None of the present results show statistically significant abyssal salinity values above the global average, and results remain consistent, apart from a general increase owing to diminished sea level, with a more conventional salinity distribution having deep values lower than the global mean. The Southern Ocean core does show a higher salinity than the North Atlantic one on the Bermuda Rise at different water depths. Although much more sophisticated models of the pore-fluid salinity can be used, they will only increase the resulting uncertainties, unless considerably more data can be obtained. Results are consistent with complex regional variations in abyssal salinity during deglaciation, but none are statistically significant.

Molecular Dynamic Simulation of Water Vapor and Determination of Diffusion Characteristics in the Pore

NASA Astrophysics Data System (ADS)

Nikonov, Eduard G.; Pavluš, Miron; Popovičová, Mária

2018-02-01

One of the varieties of pores, often found in natural or artificial building materials, are the so-called blind pores of dead-end or saccate type. Three-dimensional model of such kind of pore has been developed in this work. This model has been used for simulation of water vapor interaction with individual pore by molecular dynamics in combination with the diffusion equation method. Special investigations have been done to find dependencies between thermostats implementations and conservation of thermodynamic and statistical values of water vapor - pore system. The two types of evolution of water - pore system have been investigated: drying and wetting of the pore. Full research of diffusion coefficient, diffusion velocity and other diffusion parameters has been made.

NMR-based diffusion pore imaging.

PubMed

Laun, Frederik Bernd; Kuder, Tristan Anselm; Wetscherek, Andreas; Stieltjes, Bram; Semmler, Wolfhard

2012-08-01

Nuclear magnetic resonance (NMR) diffusion experiments offer a unique opportunity to study boundaries restricting the diffusion process. In a recent Letter [Phys. Rev. Lett. 107, 048102 (2011)], we introduced the idea and concept that such diffusion experiments can be interpreted as NMR imaging experiments. Consequently, images of closed pores, in which the spins diffuse, can be acquired. In the work presented here, an in-depth description of the diffusion pore imaging technique is provided. Image artifacts due to gradient profiles of finite duration, field inhomogeneities, and surface relaxation are considered. Gradients of finite duration lead to image blurring and edge enhancement artifacts. Field inhomogeneities have benign effects on diffusion pore images, and surface relaxation can lead to a shrinkage and shift of the pore image. The relation between boundary structure and the imaginary part of the diffusion weighted signal is analyzed, and it is shown that information on pore coherence can be obtained without the need to measure the phase of the diffusion weighted signal. Moreover, it is shown that quite arbitrary gradient profiles can be used for diffusion pore imaging. The matrices required for numerical calculations are stated and provided as supplemental material.

Pore Fluid Pressure Development in Compacting Fault Gouge in Theory, Experiments, and Nature

NASA Astrophysics Data System (ADS)

Faulkner, D. R.; Sanchez-Roa, C.; Boulton, C.; den Hartog, S. A. M.

2018-01-01

The strength of fault zones is strongly dependent on pore fluid pressures within them. Moreover, transient changes in pore fluid pressure can lead to a variety of slip behavior from creep to unstable slip manifested as earthquakes or slow slip events. The frictional properties of low-permeability fault gouge in nature and experiment can be affected by pore fluid pressure development through compaction within the gouge layer, even when the boundaries are drained. Here the conditions under which significant pore fluid pressures develop are analyzed analytically, numerically, and experimentally. Friction experiments on low-permeability fault gouge at different sliding velocities show progressive weakening as slip rate is increased, indicating that faster experiments are incapable of draining the pore fluid pressure produced by compaction. Experiments are used to constrain the evolution of the permeability and pore volume needed for numerical modeling of pore fluid pressure build up. The numerical results are in good agreement with the experiments, indicating that the principal physical processes have been considered. The model is used to analyze the effect of pore fluid pressure transients on the determination of the frictional properties, illustrating that intrinsic velocity-strengthening behavior can appear velocity weakening if pore fluid pressure is not given sufficient time to equilibrate. The results illustrate that care must be taken when measuring experimentally the frictional characteristics of low-permeability fault gouge. The contribution of compaction-induced pore fluid pressurization leading to weakening of natural faults is considered. Cyclic pressurization of pore fluid within fault gouge during successive earthquakes on larger faults may reset porosity and hence the capacity for compaction weakening.

Nuclear magnetic resonance diffusion pore imaging: Experimental phase detection by double diffusion encoding

NASA Astrophysics Data System (ADS)

Demberg, Kerstin; Laun, Frederik Bernd; Windschuh, Johannes; Umathum, Reiner; Bachert, Peter; Kuder, Tristan Anselm

2017-02-01

Diffusion pore imaging is an extension of diffusion-weighted nuclear magnetic resonance imaging enabling the direct measurement of the shape of arbitrarily formed, closed pores by probing diffusion restrictions using the motion of spin-bearing particles. Examples of such pores comprise cells in biological tissue or oil containing cavities in porous rocks. All pores contained in the measurement volume contribute to one reconstructed image, which reduces the problem of vanishing signal at increasing resolution present in conventional magnetic resonance imaging. It has been previously experimentally demonstrated that pore imaging using a combination of a long and a narrow magnetic field gradient pulse is feasible. In this work, an experimental verification is presented showing that pores can be imaged using short gradient pulses only. Experiments were carried out using hyperpolarized xenon gas in well-defined pores. The phase required for pore image reconstruction was retrieved from double diffusion encoded (DDE) measurements, while the magnitude could either be obtained from DDE signals or classical diffusion measurements with single encoding. The occurring image artifacts caused by restrictions of the gradient system, insufficient diffusion time, and by the phase reconstruction approach were investigated. Employing short gradient pulses only is advantageous compared to the initial long-narrow approach due to a more flexible sequence design when omitting the long gradient and due to faster convergence to the diffusion long-time limit, which may enable application to larger pores.

The action of water films at Å-scales in the Earth: Implications for the Nankai subduction system

NASA Astrophysics Data System (ADS)

Brown, Kevin M.; Poeppe, Dean; Josh, Matthew; Sample, James; Even, Emilie; Saffer, Demian; Tobin, Harold; Hirose, Takehiro; Kulongoski, J. T.; Toczko, Sean; Maeda, Lena; IODP Expedition 348 Shipboard Party

2017-04-01

Water properties change with confinement within nanofilms trapped between natural charged clay particles. We investigated nanofilm characteristics through high-stress laboratory compression tests in combination with analyses of expelled pore fluids. We utilized sediments obtained from deep drilling of the Nankai subduction zone at Site C0002 of the Integrated Ocean Drilling Program (IODP). We show that below 1-2 km, there should be widespread ultrafiltration of migrating fluids. Experiments to > ∼ 100 MPa normal compression collapse pores below a few ion monofilm thicknesses. A reduction towards a single condensing/dehydrating ion monofilm occurs as stresses rise >100-200 MPa and clay separations are reduced to <10-20 Å. Thus, porosity in high mineral surface area systems only consists of double and single monofilms at depths below a few km leaving little room for either bulk water or the deep biosphere. The resulting semipermeable properties result in variable segregation of ions and charged isotopes and water during active flow. The ultrafiltration and ion dehydration processes are coupled in that both require the partial immobilization of ions between the charged clay surfaces. The general effect is to increase salinities in residual pore fluids at depth and freshen fluids expelled during consolidation. Cessation of nanofilm collapse to a near constant ∼17 Å below 2 km depth at Nankai supports the contention for the onset of substantial geopressuring on the deeper seismogenic fault. The properties of monofilm water, thus, have considerable implications for the deep water properties of subduction zones generating major tremor and Mw 8+ earthquakes. Indeed, the combined effects of advective flow, ultrafiltration, diffusion, and diagenesis could provide a unifying explanation for the origins of overpressuring and pore water geochemical signals observed in many natural systems.

Capillary trapping quantification in sandstones using NMR relaxometry

NASA Astrophysics Data System (ADS)

Connolly, Paul R. J.; Vogt, Sarah J.; Iglauer, Stefan; May, Eric F.; Johns, Michael L.

2017-09-01

Capillary trapping of a non-wetting phase arising from two-phase immiscible flow in sedimentary rocks is critical to many geoscience scenarios, including oil and gas recovery, aquifer recharge and, with increasing interest, carbon sequestration. Here we demonstrate the successful use of low field 1H Nuclear Magnetic Resonance [NMR] to quantify capillary trapping; specifically we use transverse relaxation time [T2] time measurements to measure both residual water [wetting phase] content and the surface-to-volume ratio distribution (which is proportional to pore size] of the void space occupied by this residual water. Critically we systematically confirm this relationship between T2 and pore size by quantifying inter-pore magnetic field gradients due to magnetic susceptibility contrast, and demonstrate that our measurements at all water saturations are unaffected. Diffusion in such field gradients can potentially severely distort the T2-pore size relationship, rendering it unusable. Measurements are performed for nitrogen injection into a range of water-saturated sandstone plugs at reservoir conditions. Consistent with a water-wet system, water was preferentially displaced from larger pores while relatively little change was observed in the water occupying smaller pore spaces. The impact of cyclic wetting/non-wetting fluid injection was explored and indicated that such a regime increased non-wetting trapping efficiency by the sequential occupation of the most available larger pores by nitrogen. Finally the replacement of nitrogen by CO2 was considered; this revealed that dissolution of paramagnetic minerals from the sandstone caused by its exposure to carbonic acid reduced the in situ bulk fluid T2 relaxation time on a timescale comparable to our core flooding experiments. The implications of this for the T2-pore size relationship are discussed.

Estimating rates of authigenic carbonate precipitation in modern marine sediments

NASA Astrophysics Data System (ADS)

Mitnick, E. H.; Lammers, L. N.; DePaolo, D. J.

2015-12-01

The formation of authigenic carbonate (AC) in marine sediments provides a plausible explanation for large, long-lasting marine δ13C excursions that does not require extreme swings in atmospheric O2 or CO2. AC precipitation during diagenesis is driven by alkalinity production during anaerobic organic matter oxidation and is coupled to sulfate reduction. To evaluate the extent to which this process contributes to global carbon cycling, we need to relate AC production to the geochemical and geomicrobiological processes and ocean chemical conditions that control it. We present a method to estimate modern rates of AC precipitation using an inversion approach based on the diffusion-advection-reaction equation and sediment pore fluid chemistry profiles as a function of depth. SEM images and semi-quantitative elemental map analyses provide further constraints. Our initial focus is on ODP sites 807 and 1082. We sum the diffusive, advective, and reactive terms that describe changes in pore fluid Ca and Mg concentrations due to precipitation of secondary carbonate. We calculate the advective and diffusive terms from the first and second derivatives of the Ca and Mg pore fluid concentrations using a spline fit to the data. Assuming steady-state behavior we derive net AC precipitation rates of up to 8 x 10-4 mmol m-2 y-1 for Site 807 and 0.6 mmol m-2 y-1 for Site 1082. Site 1082 sediments contain pyrite, which increases in amount down-section towards the estimated peak carbonate precipitation rate, consistent with sulfate-reduction-induced AC precipitation. However, the presence of gypsum and barite throughout the sediment column implies incomplete sulfate reduction and merits further investigation of the biogeochemical reactions controlling authigenesis. Further adjustments to our method could account for the small but non-negligible fraction of groundmass with a CaSO4 signature. Our estimates demonstrate that AC formation may represent a sizeable flux in the modern global carbon cycle, on order of 1013 g C y-1. Further, it is likely to have played an even more impactful role in the Paleozoic and Precambrian, when lower surface O2 concentrations created reducing conditions favoring increased carbon burial and alkalinity production during diagenesis.

MuSET, A High Precision Logging Sensor For Downhole Spontaneous Electrical Potential.

NASA Astrophysics Data System (ADS)

Pezard, P. A.; Gautier, S.; Le Borgne, T.; Deltombe, J.

2008-12-01

MuSET has been designed by ALT and CNRS in the context of the EC ALIANCE research project. It is based on an existing multi-parameter borehole fluid sensor (p, T, Cw, pH, Eh) built by ALT. The new downhole geophysical tool aims to measure subsurface spontaneous electrical potentials (SP) in situ with great precision (< µV). For this, the device includes an unpolazirable Pb/PbCl2 electrode referred to a similar one at surface. Initial field testing in Montpellier (Languedoc, France), Ploemeur (Brittany, France) and Campos (Mallorca, Spain) took advantage of the set of field sites developed as part of ALIANCE then as part of the environmental research observatory (ORE) network for hydrogeology "H+". While Cretaceous marly limestone at Lavalette (Montpellier) proved to be almost exclusively the source of membrane potential, the clay-starved Miocene reefal carbonates of Campos generate a signal dominated by electrokinetic potential. This signal is generated due to nearby agricultural pumping, and associated strong horizontal flow. At the top of the salt to fresh water transtion, a discrepancy between the SP signal and the absence of vertical flow measured with a heat-pulse flowmeter hints at a capacity to detect the "fluid-junction", diffusion potential. At Ploemeur, the altered granite found in the vicinity of faults and fractures is also the source of a SP signal, mostly surface related while most fractures appear to be closed. In all, the MuSET demonstrates a capacity to identify several subsurface sources of natural electrical potential such as diffusion ones (membrane potential in the presence of clays, fickean processes due to pore fluid salinity gradients), or else the electrokinetic potential with pore fluid pressure gradients. While spontaneous electrical currents often loop out of the borehole, MuSET might be used as a radial electrical flowmeter once the diffusion components taken into account.

Mapping of accumulated nitrogen in the sediment pore water of a eutrophic lake in Iowa, USA

USGS Publications Warehouse

Iqbal, M.Z.; Fields, C.L.

2009-01-01

A large pool of nitrogen in the sediment pore fluid of a eutrophic lake in Iowa, USA, was mapped in this study. Previously, the lake had supported fishing and boating, but today it no longer supports its designated uses as a recreational water body. In the top 5 cm of the lake bottom, the pore water nitrogen ranges between 3.1 and 1,250 ??g/cm3 of sediments, with an average of 160.3 ??g/cm3. Vertically, nitrate concentrations were measured as 153 ??g/cm3 at 0-10 cm, 162 ??g/cm3 at 10-20 cm, and 32 ??g/cm3 at 20-30 cm. Nitrate mass distribution was quantified as 3.67 ?? 103 kg (65%) in the bottom sediments, 172 kg (3%) in suspended particulates, and 1.83 ?? 103 kg (32%) in the dissolved phase. Soil runoff nutrients arrive at the lake from the heavily fertilized lands in the watershed. Upon sedimentation, a large mass of nitrogen desorbs from mineral particles to the relatively immobile pore fluid. Under favorable conditions, this nitrogen diffuses back into the water column, thereby dramatically limiting the lake's capability to process incoming nutrients from farmlands. Consequently, a condition of oxygen deficiency disrupts the post-season biological activities in the lake. ?? 2008 Springer-Verlag.

Mobility Effect on Poroelastic Seismic Signatures in Partially Saturated Rocks With Applications in Time-Lapse Monitoring of a Heavy Oil Reservoir

NASA Astrophysics Data System (ADS)

Zhao, Luanxiao; Yuan, Hemin; Yang, Jingkang; Han, De-hua; Geng, Jianhua; Zhou, Rui; Li, Hui; Yao, Qiuliang

2017-11-01

Conventional seismic analysis in partially saturated rocks normally lays emphasis on estimating pore fluid content and saturation, typically ignoring the effect of mobility, which decides the ability of fluids moving in the porous rocks. Deformation resulting from a seismic wave in heterogeneous partially saturated media can cause pore fluid pressure relaxation at mesoscopic scale, thereby making the fluid mobility inherently associated with poroelastic reflectivity. For two typical gas-brine reservoir models, with the given rock and fluid properties, the numerical analysis suggests that variations of patchy fluid saturation, fluid compressibility contrast, and acoustic stiffness of rock frame collectively affect the seismic reflection dependence on mobility. In particular, the realistic compressibility contrast of fluid patches in shallow and deep reservoir environments plays an important role in determining the reflection sensitivity to mobility. We also use a time-lapse seismic data set from a Steam-Assisted Gravity Drainage producing heavy oil reservoir to demonstrate that mobility change coupled with patchy saturation possibly leads to seismic spectral energy shifting from the baseline to monitor line. Our workflow starts from performing seismic spectral analysis on the targeted reflectivity interface. Then, on the basis of mesoscopic fluid pressure diffusion between patches of steam and heavy oil, poroelastic reflectivity modeling is conducted to understand the shift of the central frequency toward low frequencies after the steam injection. The presented results open the possibility of monitoring mobility change of a partially saturated geological formation from dissipation-related seismic attributes.

Fluid Transport in Porous Media probed by Relaxation-Exchange NMR

NASA Astrophysics Data System (ADS)

Olaru, A. M.; Kowalski, J.; Sethi, V.; Blümich, B.

2011-12-01

The characterization of fluid transport in porous media represents a matter of high interest in fields like the construction industry, oil exploitation, and soil science. Moisture migration or flow at low rates, such as those occurring in soil during rain are difficult to characterize by classical high-field NMR velocimetry due to the dedicated hardware and elaborate techniques required for adequate signal encoding. The necessity of field studies raises additional technical problems, which can be solved only by the use of portable low-field NMR instruments. In this work we extend the use of low-field relaxation exchange experiments from the study of diffusive transport to that of advection. Relaxation exchange experiments were performed using a home-built Halbach magnet on model porous systems with controlled pore-size distributions and on natural porous systems (quartz sand with a broad pore-size distribution) exposed to unidirectional flow. Different flow rates leave distinctive marks on the exchange maps obtained by inverse Laplace transformation of the time domain results, due to the superposition of exchange, diffusion and inflow/outflow in multiple relaxation sites of the liquids in the porous media. In the case of slow velocities there is no loss of signal due to outflow, and the relaxation-exchange effects prevail, leading to a tilt of the diagonal distribution around a pivot point with increasing mixing time. The tilt suggests an asymmetry in the exchange between relaxation sites of large and small decay rates. Another observed phenomenon is the presence of a bigger number of exchange cross-peaks compared to the exchange maps obtained for the same systems in zero-flow conditions. We assume that this is due to enhanced exchange caused by the superposition of flow. For high velocities the outflow effects dominate and the relaxation-time distribution collapses towards lower values of the average relaxation times. In both cases the pore-size distribution has a strong effect on the results, the asymmetries being more obvious in the natural porous systems than in the glass bead packs used as models, while the enhanced exchange phenomenon appears predominantly in the maps obtained on the model systems. This is probably due to diffusion occurring in the presence of different internal field gradients. Shifts and tilts in the exchange maps can be simulated by solving the relaxation site-averaged Bloch-Torrey system forward in time and assuming an asymmetric closure for the transport, which might be realistic for preferential flow phenomena or for pore-size distributions with two or more clearly distinct pore size classes. When comparing the simulations results with the experimental data we observed a correspondence of signal collapse and translation towards lower relaxation times. The asymmetries could be qualitatively reproduced by making further assumptions on the pore structure, but further work is required to characterize and model the physical phenomenon behind. The results obtained reveal the possibility of characterizing advective fluid transport in porous systems by simple correlation experiments performed with inexpensive and mobile hardware.

The Pawnee Sequence: Poroelastic Effects from Injection in Osage County, Oklahoma

NASA Astrophysics Data System (ADS)

Barbour, A. J.; Rubinstein, J. L.

2016-12-01

Aggregate multi-year records of wastewater injection in Oklahoma show that the strongest change in injection within 20 km of the 2016 M5.8 Pawnee strike-slip earthquake was in Osage County, where injection rates increased rapidly in late-2012 by nearly a factor of three above previous levels. After this increase, rates there declined steadily over two years to an average rate characteristic of all other injection wells in Pawnee and Noble Counties, remaining relatively constant until the beginning of the earthquake sequence. Here we test if poroelastic effects associated with this injection-rate transient can help explain the relative timing between peak injection rates and the beginning of the Pawnee sequence. Although the alternative hypothesis that regional-scale faults and fractures in critically stressed rock serve as fast-pathways for fluid diffusion cannot be ruled out, it appears to be difficult to reconcile based solely on injection data and space-time patterns for this seismic sequence. We simulate the cylindrically symmetric, transient strain and pore pressure fields for an injection-source time function emulating the injection history in a layered half-space in accordance with linear poroelasticity. In the simulation domain, injection occurs at depths of 1300 - 1900 m, into a homogeneous basal sedimentary reservoir representing the Arbuckle Group, overlying a semi-infinite layer representing granitic basement; we determined the hydraulic, elastic, and poroelastic properties of these layers from published literature. At the mainshock hypocenter, this numerical model predicts a delay between peak injection rates and pore pressure increase that is strongly dependent on hydraulic diffusivity; however, the duration is also controlled by the bulk elastic properties and the undrained Skempton's coefficient of the rock. Furthermore, because of fluid-strain coupling, pore pressures in the basement rock decrease during this delay period, which would tend to stabilize temporarily a critically stressed fault. Even though pore pressure diffusion is the dominant mechanism at play, poroelastic effects do affect the relative timing assuming a reasonable set of material parameters, even though strain rates in the basement are relatively low compared to rates in the Arbuckle layer (and above).

Fluid‐driven seismicity response of the Rinconada fault near Paso Robles, California, to the 2003 M 6.5 San Simeon earthquake

USGS Publications Warehouse

Hardebeck, Jeanne L.

2012-01-01

The 2003 M 6.5 San Simeon, California, earthquake caused significant damage in the city of Paso Robles and a persistent cluster of aftershocks close to Paso Robles near the Rinconada fault. Given the importance of secondary aftershock triggering in sequences of large events, a concern is whether this cluster of events could trigger another damaging earthquake near Paso Robles. An epidemic‐type aftershock sequence (ETAS) model is fit to the Rinconada seismicity, and multiple realizations indicate a 0.36% probability of at least one M≥6.0 earthquake during the next 30 years. However, this probability estimate is only as good as the projection into the future of the ETAS model. There is evidence that the seismicity may be influenced by fluid pressure changes, which cannot be forecasted using ETAS. The strongest evidence for fluids is the delay between the San Simeon mainshock and a high rate of seismicity in mid to late 2004. This delay can be explained as having been caused by a pore pressure decrease due to an undrained response to the coseismic dilatation, followed by increased pore pressure during the return to equilibrium. Seismicity migration along the fault also suggests fluid involvement, although the migration is too slow to be consistent with pore pressure diffusion. All other evidence, including focal mechanisms and b‐value, is consistent with tectonic earthquakes. This suggests a model where the role of fluid pressure changes is limited to the first seven months, while the fluid pressure equilibrates. The ETAS modeling adequately fits the events after July 2004 when the pore pressure stabilizes. The ETAS models imply that while the probability of a damaging earthquake on the Rinconada fault has approximately doubled due to the San Simeon earthquake, the absolute probability remains low.

Lennard-Jones fluids in two-dimensional nano-pores. Multi-phase coexistence and fluid structure

NASA Astrophysics Data System (ADS)

Yatsyshin, Petr; Savva, Nikos; Kalliadasis, Serafim

2014-03-01

We present a number of fundamental findings on the wetting behaviour of nano-pores. A popular model for fluid confinement is a one-dimensional (1D) slit pore formed by two parallel planar walls and it exhibits capillary condensation (CC): a first-order phase transition from vapour to capillary-liquid (Kelvin shift). Capping such a pore at one end by a third orthogonal wall forms a prototypical two-dimensional (2D) pore. We show that 2D pores possess a wetting temperature such that below this temperature CC remains of first order, above it becomes a continuous phase transition manifested by a slab of capillary-liquid filling the pore from the capping wall. Continuous CC exhibits hysteresis and can be preceded by a first-order capillary prewetting transition. Additionally, liquid drops can form in the corners of the 2D pore (remnant of 2D wedge prewetting). The three fluid phases, vapour, capillary-liquid slab and corner drops, can coexist at the pore triple point. Our model is based on the statistical mechanics of fluids in the density functional formulation. The fluid-fluid and fluid-substrate interactions are dispersive. We analyze in detail the microscopic fluid structure, isotherms and full phase diagrams. Our findings also suggest novel ways to control wetting of nano-pores. We are grateful to the European Research Council via Advanced Grant No. 247031 for support.

Hydromechanics of Reservoir Induced Seismicity

NASA Astrophysics Data System (ADS)

Dura-Gomez, Inmaculada

Data from five reservoirs were analyzed to investigate the various factors and possible pore pressure thresholds associated with Reservoir Induced Seismicity (RIS). Data was obtained from the following reservoirs: Koyna and Warna Reservoirs in India, Itoiz Reservoir in the western Pyrenees, Spain, and Jocassee and Monticello Reservoirs in South Carolina, U.S.A. Koyna Reservoir is one out of four reservoirs in the world where M≥6.0 induced earthquakes have occurred, whereas Warna Reservoir accounts for one out of ten cases with 5.0≤M≤5.9 induced earthquakes. Induced seismicity in the Koyna-Warna region is associated with annual filling cycles in the two reservoirs, large water level changes (30 to 45 m) and the presence of regional scale fractures. The Koyna-Warna case includes 19 M≥5.0 earthquakes at non-repeating hypocenters. The calculation of excess pore pressures associated with these earthquakes suggests values >300 kPa or >600 kPa, before or after 1993 respectively. The need for larger pore pressures from 1993 suggests that M≥5 earthquakes were induced on stronger faults in the region. The exceedance of the previous water level maxima (stress memory) is the most important, although not determining factor in inducing these M≥5.0 earthquakes. Itoiz Reservoir is one of twenty nine reservoirs with 4.0≤M≤4.9 induced earthquakes. The analysis of the RIS associated with the Itoiz Reservoir impoundment, between January 2004 and the end of 2008, shows that that pore pressures diffuse away from Itoiz Reservoir through the carbonate megabreccia systems of the Early to Middle Eocene Hecho Group, and a series of near-vertical thrust faults above the gently dipping Gavarnie thrust. Excess diffused pore pressures destabilize saturated critically stressed seismogenic fractures where RIS takes place. In particular, M≥3.0 earthquakes in the region are associated with excess pore pressures of the order of 100 to 200 kPa. Jocassee and Monticello Reservoirs in South Carolina are among the best studied reservoir induced seismicity cases in the world, and have been associated with M<4 and M<3 earthquakes respectively. The analysis of some of these earthquakes emphasizes the contribution of the diffused pore pressures to the observed seismicity. In the case of Jocassee Reservoir, 2.0≤M≤3.0 earthquakes occurred in a homogenous, non-fractured crystalline rock and are associated with excess diffused pore pressures of the order of 600 to 700 kPa. Earthquakes associated with the impoundment of Monticello Reservoir occurred in a region of very complex geology, with many pre-existing local scale fractures. The calculation of the excess diffused pore pressures associated with 2.0≤M≤3.0 earthquakes yielded values of the order of 100 to 300 kPa. Synthesis of these data show that RIS occurs when excess pore pressures (which occur primarily by diffusion) reach threshold values needed to induce RIS. The occurrence of RIS and its magnitude are controlled by the filling history, availability of fluid filled saturated fractures and their hydrogeological properties.

A new methodology for determination of macroscopic transport parameters in drying porous media

NASA Astrophysics Data System (ADS)

Attari Moghaddam, A.; Kharaghani, A.; Tsotsas, E.; Prat, M.

2015-12-01

Two main approaches have been used to model the drying process: The first approach considers the partially saturated porous medium as a continuum and partial differential equations are used to describe the mass, momentum and energy balances of the fluid phases. The continuum-scale models (CM) obtained by this approach involve constitutive laws which require effective material properties, such as the diffusivity, permeability, and thermal conductivity which are often determined by experiments. The second approach considers the material at the pore scale, where the void space is represented by a network of pores (PN). Micro- or nanofluidics models used in each pore give rise to a large system of ordinary differential equations with degrees of freedom at each node of the pore network. In this work, the moisture transport coefficient (D), the pseudo desorption isotherm inside the network and at the evaporative surface are estimated from the post-processing of the three-dimensional pore network drying simulations for fifteen realizations of the pore space geometry from a given probability distribution. A slice sampling method is used in order to extract these parameters from PN simulations. The moisture transport coefficient obtained in this way is shown in Fig. 1a. The minimum of average D values demonstrates the transition between liquid dominated moisture transport region and vapor dominated moisture transport region; a similar behavior has been observed in previous experimental findings. A function is fitted to the average D values and then is fed into the non-linear moisture diffusion equation. The saturation profiles obtained from PN and CM simulations are shown in Fig. 1b. Figure 1: (a) extracted moisture transport coefficient during drying for fifteen realizations of the pore network, (b) average moisture profiles during drying obtained from PN and CM simulations.

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

DOE PAGES

Ackerman, David M.; Evans, James W.

2017-01-19

Here, we perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient D tr(x), at various positions x within themore » pore. D tr(x) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.« less

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

NASA Astrophysics Data System (ADS)

Ackerman, David M.; Evans, James W.

2017-01-01

We perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient Dtr(x ) , at various positions x within the pore. Dtr(x ) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.

Importance of Ion Packing on the Dynamics of Ionic Liquids during Micropore Charging.

PubMed

He, Yadong; Qiao, Rui; Vatamanu, Jenel; Borodin, Oleg; Bedrov, Dmitry; Huang, Jingsong; Sumpter, Bobby G

2016-01-07

Molecular simulations of the diffusion of EMIM(+) and TFSI(-) ions in slit-shaped micropores under conditions similar to those during charging show that in pores that accommodate only a single layer of ions, ions diffuse increasingly faster as the pore becomes charged (with diffusion coefficients even reaching ∼5 × 10(-9) m(2)/s), unless the pore becomes very highly charged. In pores wide enough to fit more than one layer of ions, ion diffusion is slower than in the bulk and changes modestly as the pore becomes charged. Analysis of these results revealed that the fast (or slow) diffusion of ions inside a micropore during charging is correlated most strongly with the dense (or loose) ion packing inside the pore. The molecular details of the ions and the precise width of the pores modify these trends weakly, except when the pore is so narrow that the ion conformation relaxation is strongly constrained by the pore walls.

Importance of Ion Packing on the Dynamics of Ionic Liquids during Micropore Charging

DOE PAGES

He, Yadong; Qiao, Rui; Vatamanu, Jenel; ...

2015-12-07

In molecular simulations of the diffusion of EMIM+ and TESI- ions in slit-shaped micropores under conditions similar to those during charging show that in pores that accommodate only a single layer of ions, ions diffuse increasingly faster as the pore becomes charged (with diffusion coefficients even reaching similar to 5 x 10 -9 m 2/s), unless the pore becomes very highly charged. In pores wide enough to fit more than one layer of ions, ion diffusion is slower than in the bulk and changes modestly as the pore becomes charged. Moreover, analysis of these results revealed that the fast (ormore » slow) diffusion of ions inside a micropore during charging is correlated most strongly with the dense (or loose) ion packing inside the pore. Finally, the molecular details of the ions and the precise width of the pores modify these trends weakly, except when the pore is so narrow that the ion conformation relaxation is strongly constrained by the pore walls.« less

The effect of initial resident fluid saturation on the interaction between resident and infiltrating fluids in porous media

NASA Astrophysics Data System (ADS)

Hsu, S. Y.; Chen, H.; Huang, Q. Z.; Lee, T. Y.; Chiu, Y.; Chang, L. C.; Lamorski, K.; Sławiński, C.; Tsao, C. W.

2017-12-01

The interplay between resident ("old") fluid already in the vadose zone and infiltrating ("new") fluid was examined with micromodel experiments. The geometric patterns of the micromodels are based on a pore doublet and a 2D pore geometry of a sand-packing soil scanned by Micro X-Ray CT. We studied the old and new fluid interaction during imbibition process subject to different evaporation times (different the initial old fluid saturations). The results found that, in the pore-doublet micromodel experiment, the old fluid was mixed and displaced by the new fluid, and an increase in the initial old fluid saturation led to a decrease in the amount of old fluid displaced by the new fluid. On the other hand, the most of the old fluid in the micromodel of 2D sand-packing pore geometry was displaced by and mixed with the new fluid. However, a small amount of the initial old fluid that occupied pore throats remained untouched by the new fluid due to the air blockage. The amount of untouched old fluid increased as the initial old fluid saturation decreased. Our finding reveals the effect of pore geometry and inital old fluid distribution on the interaction between resident and infiltrating fluids.

Hydromechanical coupling in fractured rock masses: mechanisms and processes of selected case studies

NASA Astrophysics Data System (ADS)

Zangerl, Christian

2015-04-01

Hydromechanical (HM) coupling in fractured rock play an important role when events including dam failures, landslides, surface subsidences due to water withdrawal or drainage, injection-induced earthquakes and others are analysed. Generally, hydromechanical coupling occurs when a rock mass contain interconnected pores and fractures which are filled with water and pore/fracture pressures evolves. In the on hand changes in the fluid pressure can lead to stress changes, deformations and failures of the rock mass. In the other hand rock mass stress changes and deformations can alter the hydraulic properties and fluid pressures of the rock mass. Herein well documented case studies focussing on surface subsidence due to water withdrawal, reversible deformations of large-scale valley flanks and failure as well as deformation processes of deep-seated rock slides in fractured rock masses are presented. Due to pore pressure variations HM coupling can lead to predominantly reversible rock mass deformations. Such processes can be considered by the theory of poroelasticity. Surface subsidence reaching magnitudes of few centimetres and are caused by water drainage into deep tunnels are phenomenas which can be assigned to processes of poroelasticity. Recently, particular focus was given on large tunnelling projects to monitor and predict surface subsidence in fractured rock mass in oder to avoid damage of surface structures such as dams of large reservoirs. It was found that surface subsidence due to tunnel drainage can adversely effect infrastructure when pore pressure drawdown is sufficiently large and spatially extended and differential displacements which can be amplified due to topographical effects e.g. valley closure are occurring. Reversible surface deformations were also ascertained on large mountain slopes and summits with the help of precise deformation measurements i.e. permanent GPS or episodic levelling/tacheometric methods. These reversible deformations are often in the range of millimetres to a very few centimetres and can be linked to annual groundwater fluctuations. Due to pore pressure variations HM coupling can influence seepage forces and effective stresses in the rock mass. Effective stress changes can adversely affect the stability and deformation behaviour of deep-seated rock slides by influencing the shear strength or the time dependent (viscous) material behaviour of the basal shear zone. The shear strength of active shear zones is often reasonably well described by Coulomb's law. In Coulomb's law the operative normal stresses to the shear surface/zone are effective stresses and hence pore pressures which should be taken into account reduces the shear strength. According to the time dependent material behaviour a few effective stress based viscous models exists which are able to consider pore pressures. For slowly moving rock slides HM coupling could be highly relevant when low-permeability clayey-silty shear zones (fault gouges) are existing. An important parameters therefore is the hydraulic diffusivity, which is controlled by the permeability and fluid-pore compressibility of the shear zone, and by fluid viscosity. Thus time dependent pore pressure diffusion in the shear zone can either control the stability condition or the viscous behaviour (creep) of the rock slide. Numerous cases studies show that HM coupling can effect deformability, shear strength and time dependent behaviour of fractured rock masses. A process-based consideration can be important to avoid unexpected impacts on infrastructures and to understand complex rock mass as well rock slide behaviour.

Experimental and theoretical studies of the effect of temperature on supercritical CO2 adsorption on illite

NASA Astrophysics Data System (ADS)

Joewondo, N.; Zhang, Y.; Prasad, M.

2016-12-01

Sequestration of carbon dioxide in shale has been a subject of interest as the result of the technological advancement in gas shale production. The process involves injection of CO2 to enhance methane recovery and storing CO2 in depleted shale reservoir at elevated pressures. To better understand both shale production and carbon storage one must study the physical phenomena acting at different scales that control the in situ fluid flow. Shale rocks are complex systems with heterogeneous structures and compositions. Pore structures of these systems are in nanometer scales and have significant gas storage capacity and surface area. Adsorption is prominent in nanometer sized pores due to the high attraction between gas molecules and the surface of the pores. Recent studies attempt to find correlation between storage capacity and the rock composition, particularly the clay content. This study, however, focuses on the study of supercritical adsorption of CO2 on pure clay sample. We have built an in-house manometric experimental setup that can be used to study both the equilibrium and kinetics of adsorption. The experiment is conducted at isothermal condition. The study of equilibrium of adsorption gives insight on the storage capacity of these systems, and the study of the kinetics of adsorption is essential in understanding the resistance to fluid transport. The diffusion coefficient, which can be estimated from the dynamic experimental results, is a parameter which quantify diffusion mobility, and is affected by many factors including pressure and temperature. The first part of this paper briefly discusses the study of both equilibrium and kinetics of the CO2 adsorption on illite. Both static and dynamic measurements on the system are compared to theoretical models available in the literature to estimate the storage capacity and the diffusion time constants. The main part of the paper discusses the effect of varying temperature on the static and dynamic experimental results.

Effect of pore geometry on the compressibility of a confined simple fluid

NASA Astrophysics Data System (ADS)

Dobrzanski, Christopher D.; Maximov, Max A.; Gor, Gennady Y.

2018-02-01

Fluids confined in nanopores exhibit properties different from the properties of the same fluids in bulk; among these properties is the isothermal compressibility or elastic modulus. The modulus of a fluid in nanopores can be extracted from ultrasonic experiments or calculated from molecular simulations. Using Monte Carlo simulations in the grand canonical ensemble, we calculated the modulus for liquid argon at its normal boiling point (87.3 K) adsorbed in model silica pores of two different morphologies and various sizes. For spherical pores, for all the pore sizes (diameters) exceeding 2 nm, we obtained a logarithmic dependence of fluid modulus on the vapor pressure. Calculation of the modulus at saturation showed that the modulus of the fluid in spherical pores is a linear function of the reciprocal pore size. The calculation of the modulus of the fluid in cylindrical pores appeared too scattered to make quantitative conclusions. We performed additional simulations at higher temperature (119.6 K), at which Monte Carlo insertions and removals become more efficient. The results of the simulations at higher temperature confirmed both regularities for cylindrical pores and showed quantitative difference between the fluid moduli in pores of different geometries. Both of the observed regularities for the modulus stem from the Tait-Murnaghan equation applied to the confined fluid. Our results, along with the development of the effective medium theories for nanoporous media, set the groundwork for analysis of the experimentally measured elastic properties of fluid-saturated nanoporous materials.

Irradiation behavior of the interaction product of U-Mo fuel particle dispersion in an Al matrix

NASA Astrophysics Data System (ADS)

Kim, Yeon Soo; Hofman, G. L.

2012-06-01

Irradiation performance of U-Mo fuel particles dispersed in Al matrix is stable in terms of fuel swelling and is suitable for the conversion of research and test reactors from highly enriched uranium (HEU) to low enriched uranium (LEU). However, tests of the fuel at high temperatures and high burnups revealed obstacles caused by the interaction layers forming between the fuel particle and matrix. In some cases, fission gas filled pores grow and interconnect in the interdiffusion layer resulting in fuel plate failure. Postirradiation observations are made to examine the behavior of the interdiffusion layers. The interdiffusion layers show a fluid-like behavior characteristic of amorphous materials. In the amorphous interdiffusion layers, fission gas diffusivity is high and the material viscosity is low so that the fission gas pores readily form and grow. Based on the observations, a pore formation mechanism is proposed and potential remedies to suppress the pore growth are also introduced.

Rock anelasticity due to patchy saturation and fabric heterogeneity: A double double-porosity model of wave propagation

NASA Astrophysics Data System (ADS)

Ba, Jing; Xu, Wenhao; Fu, Li-Yun; Carcione, José M.; Zhang, Lin

2017-03-01

Heterogeneity of rock's fabric can induce heterogeneous distribution of immiscible fluids in natural reservoirs, since the lithological variations (mainly permeability) may affect fluid migration in geological time scales, resulting in patchy saturation of fluids. Therefore, fabric and saturation inhomogeneities both affect wave propagation. To model the wave effects (attenuation and velocity dispersion), we introduce a double double-porosity model, where pores saturated with two different fluids overlap with pores having dissimilar compressibilities. The governing equations are derived by using Hamilton's principle based on the potential energy, kinetic energy, and dissipation functions, and the stiffness coefficients are determined by gedanken experiments, yielding one fast P wave and four slow Biot waves. Three examples are given, namely, muddy siltstones, clean dolomites, and tight sandstones, where fabric heterogeneities at three different spatial scales are analyzed in comparison with experimental data. In muddy siltstones, where intrapore clay and intergranular pores constitute a submicroscopic double-porosity structure, wave anelasticity mainly occurs in the frequency range (104-107 Hz), while in pure dolomites with microscopic heterogeneity of grain contacts and tight sandstones with mesoscopic heterogeneity of less consolidated sands, it occurs at 103-107 Hz and 101-103 Hz (seismic band), respectively. The predicted maximum quality factor of the fast compressional wave for the sandstone is the lowest (approximately 8), and that of the dolomite is the highest. The results of the diffusive slow waves are affected by the strong friction effects between solids and fluids. The model describes wave propagation in patchy-saturated rocks with fabric heterogeneity at different scales, and the relevant theoretical predictions agree well with the experimental data in fully and partially saturated rocks.

Nanoscale simulation of shale transport properties using the lattice Boltzmann method: Permeability and diffusivity

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Li; Zhang, Lei; Kang, Qinjun

Here, porous structures of shales are reconstructed using the markov chain monte carlo (MCMC) method based on scanning electron microscopy (SEM) images of shale samples from Sichuan Basin, China. Characterization analysis of the reconstructed shales is performed, including porosity, pore size distribution, specific surface area and pore connectivity. The lattice Boltzmann method (LBM) is adopted to simulate fluid flow and Knudsen diffusion within the reconstructed shales. Simulation results reveal that the tortuosity of the shales is much higher than that commonly employed in the Bruggeman equation, and such high tortuosity leads to extremely low intrinsic permeability. Correction of the intrinsicmore » permeability is performed based on the dusty gas model (DGM) by considering the contribution of Knudsen diffusion to the total flow flux, resulting in apparent permeability. The correction factor over a range of Knudsen number and pressure is estimated and compared with empirical correlations in the literature. We find that for the wide pressure range investigated, the correction factor is always greater than 1, indicating Knudsen diffusion always plays a role on shale gas transport mechanisms in the reconstructed shales. Specifically, we found that most of the values of correction factor fall in the slip and transition regime, with no Darcy flow regime observed.« less

Nanoscale simulation of shale transport properties using the lattice Boltzmann method: Permeability and diffusivity

DOE PAGES

Chen, Li; Zhang, Lei; Kang, Qinjun; ...

2015-01-28

Here, porous structures of shales are reconstructed using the markov chain monte carlo (MCMC) method based on scanning electron microscopy (SEM) images of shale samples from Sichuan Basin, China. Characterization analysis of the reconstructed shales is performed, including porosity, pore size distribution, specific surface area and pore connectivity. The lattice Boltzmann method (LBM) is adopted to simulate fluid flow and Knudsen diffusion within the reconstructed shales. Simulation results reveal that the tortuosity of the shales is much higher than that commonly employed in the Bruggeman equation, and such high tortuosity leads to extremely low intrinsic permeability. Correction of the intrinsicmore » permeability is performed based on the dusty gas model (DGM) by considering the contribution of Knudsen diffusion to the total flow flux, resulting in apparent permeability. The correction factor over a range of Knudsen number and pressure is estimated and compared with empirical correlations in the literature. We find that for the wide pressure range investigated, the correction factor is always greater than 1, indicating Knudsen diffusion always plays a role on shale gas transport mechanisms in the reconstructed shales. Specifically, we found that most of the values of correction factor fall in the slip and transition regime, with no Darcy flow regime observed.« less

Nanoscale simulation of shale transport properties using the lattice Boltzmann method: permeability and diffusivity

PubMed Central

Chen, Li; Zhang, Lei; Kang, Qinjun; Viswanathan, Hari S.; Yao, Jun; Tao, Wenquan

2015-01-01

Porous structures of shales are reconstructed using the markov chain monte carlo (MCMC) method based on scanning electron microscopy (SEM) images of shale samples from Sichuan Basin, China. Characterization analysis of the reconstructed shales is performed, including porosity, pore size distribution, specific surface area and pore connectivity. The lattice Boltzmann method (LBM) is adopted to simulate fluid flow and Knudsen diffusion within the reconstructed shales. Simulation results reveal that the tortuosity of the shales is much higher than that commonly employed in the Bruggeman equation, and such high tortuosity leads to extremely low intrinsic permeability. Correction of the intrinsic permeability is performed based on the dusty gas model (DGM) by considering the contribution of Knudsen diffusion to the total flow flux, resulting in apparent permeability. The correction factor over a range of Knudsen number and pressure is estimated and compared with empirical correlations in the literature. For the wide pressure range investigated, the correction factor is always greater than 1, indicating Knudsen diffusion always plays a role on shale gas transport mechanisms in the reconstructed shales. Specifically, we found that most of the values of correction factor fall in the slip and transition regime, with no Darcy flow regime observed. PMID:25627247

Localized fluid discharge in subduction zones: Insights from tension veins around an ancient megasplay fault (Nobeoka Thrust, SW Japan)

NASA Astrophysics Data System (ADS)

Otsubo, M.; Hardebeck, J.; Miyakawa, A.; Yamaguchi, A.; Kimura, G.

2017-12-01

Fluid-rock interactions along seismogenic faults are of great importance to understand fault mechanics. The fluid loss by the formation of mode I cracks (tension cracks) increases the fault strength and creates drainage asperities along the plate interface (Sibson, 2013, Tectonophysics). The Nobeoka Thrust, in southwestern Japan, is an on-land example of an ancient megasplay fault and provides an excellent record of deformation and fluid flow at seismogenic depths of a subduction zone (Kondo et al., 2005, Tectonics). We focus on (1) Pore fluid pressure loss, (2) Amount of fault strength recovery, and (3) Fluid circulation by the formation of mode I cracks in the post-seismic period around the fault zone of the Nobeoka Thrust. Many quartz veins that filled mode I crack at the coastal outcrops suggest a normal faulting stress regime after faulting of the Nobeoka Thrust (Otsubo et al., 2016, Island Arc). We estimated the decrease of the pore fluid pressure by the formation of the mode I cracks around the Nobeoka Thrust in the post-seismic period. When the pore fluid pressure exceeds σ3, veins filling mode I cracks are constructed (Jolly and Sanderson, 1997, Jour. Struct. Geol.). We call the pore fluid pressure that exceeds σ3 "pore fluid over pressure". The differential stress in the post-seismic period and the driving pore fluid pressure ratio P* (P* = (Pf - σ3) / (σ1 - σ3), Pf: pore fluid pressure) are parameters to estimate the pore fluid over pressure. In the case of the Nobeoka Thrust (P* = 0.4, Otsubo et al., 2016, Island Arc), the pore fluid over pressure is up to 20 MPa (assuming tensile strength = 10 MPa). 20 MPa is equivalent to <10% of the total pore fluid pressure around the Nobeoka Thrust (depth = 10 km, density = 2.7 kg/m3). When the pore fluid pressure decreases by 4%, the normalized pore pressure ratio λ* (λ* = (Pf - Ph) / (Pl - Ph), Pl: lithostatic pressure; Ph: hydrostatic pressure) changes from 0.95 to 0.86. In the case of the Nobeoka Thrust, the fault strength can increase by up to 10 MPa (assuming frictional coefficient = 0.6). 10 MPa is almost equivalent to the stress drop values in large trench type earthquakes. Hence, we suggest that the fluid loss caused by the formation of mode I cracks in the post-seismic period may play an important role by increasing frictional strength along the megasplay fault.

Role of Confinement on Adsorption and Dynamics of Ethane and an Ethane–CO 2 Mixture in Mesoporous CPG Silica

DOE PAGES

Patankar, Sumant; Gautam, Siddharth; Rother, Gernot; ...

2016-02-10

It was found that ethane is confined to mineral and organic pores in certain shale formations. Effects of confinement on structural and dynamic properties of ethane in mesoporous controlled pore glass (CPG) were studied by gravimetric adsorption and quasi-elastic neutron scattering (QENS) measurements. The obtained isotherms and scattering data complement each other by quantifying the relative strength of the solid–fluid interactions and the transport properties of the fluid under confinement, respectively. We used a magnetic suspension balance to measure the adsorption isotherms at two temperatures and over a range of pressures corresponding to a bulk density range of 0.01–0.35 g/cmmore » 3. Key confinement effects were highlighted through differences between isotherms for the two pore sizes. A comparison was made with previously published isotherms for CO 2 on the same CPG materials. Behavior of ethane in the smaller pore size was probed further using quasi-elastic neutron scattering. By extracting the self-diffusivity and residence time, we were able to study the effect of pressure and transition from gaseous to supercritical densities on the dynamics of confined ethane. Moreover, a temperature variation QENS study was also completed with pure ethane and a CO 2–ethane mixture. Activation energies extracted from the Arrhenius plots show the effects of CO 2 addition on ethane mobility.« less

Validation of model predictions of pore-scale fluid distributions during two-phase flow

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Lin, Qingyang; Gao, Ying; Raeini, Ali Q.; AlRatrout, Ahmed; Bijeljic, Branko; Blunt, Martin J.

2018-05-01

Pore-scale two-phase flow modeling is an important technology to study a rock's relative permeability behavior. To investigate if these models are predictive, the calculated pore-scale fluid distributions which determine the relative permeability need to be validated. In this work, we introduce a methodology to quantitatively compare models to experimental fluid distributions in flow experiments visualized with microcomputed tomography. First, we analyzed five repeated drainage-imbibition experiments on a single sample. In these experiments, the exact fluid distributions were not fully repeatable on a pore-by-pore basis, while the global properties of the fluid distribution were. Then two fractional flow experiments were used to validate a quasistatic pore network model. The model correctly predicted the fluid present in more than 75% of pores and throats in drainage and imbibition. To quantify what this means for the relevant global properties of the fluid distribution, we compare the main flow paths and the connectivity across the different pore sizes in the modeled and experimental fluid distributions. These essential topology characteristics matched well for drainage simulations, but not for imbibition. This suggests that the pore-filling rules in the network model we used need to be improved to make reliable predictions of imbibition. The presented analysis illustrates the potential of our methodology to systematically and robustly test two-phase flow models to aid in model development and calibration.

Oscillatory fluid flow in deformable tubes: Implications for pore-scale hydromechanics from comparing experimental observations with theoretical predictions.

PubMed

Kurzeja, Patrick; Steeb, Holger; Strutz, Marc A; Renner, Jörg

2016-12-01

Oscillatory flow of four fluids (air, water, two aqueous sodium-tungstate solutions) was excited at frequencies up to 250 Hz in tubes of two materials (steel, silicone) covering a wide range in length, diameter, and thickness. The hydrodynamical response was characterized by phase shift and amplitude ratio between pressures in an upstream (pressure excitation) and a downstream reservoir connected by the tubes. The resulting standing flow waves reflect viscosity-controlled diffusive behavior and inertia-controlled wave behavior for oscillation frequencies relatively low and high compared to Biot's critical frequency, respectively. Rigid-tube theories correspond well with the experimental results for steel tubes filled with air or water. The wave modes observed for silicone tubes filled with the rather incompressible liquids or air, however, require accounting for the solid's shear and bulk modulus to correctly predict speed of pressure propagation and deformation mode. The shear mode may be responsible for significant macroscopic attenuation in porous materials with effective frame-shear moduli lower than the bulk modulus of the pore fluid. Despite notable effects of the ratio of densities and of acoustic and shear velocity of fluid and solid, Biot's frequency remains an approximate indicator of the transition from the viscosity to the inertia controlled regime.

2D Simulations of Earthquake Cycles at a Subduction Zone Based on a Rate and State Friction Law -Effects of Pore Fluid Pressure Changes-

NASA Astrophysics Data System (ADS)

Mitsui, Y.; Hirahara, K.

2006-12-01

There have been a lot of studies that simulate large earthquakes occurring quasi-periodically at a subduction zone, based on the laboratory-derived rate-and-state friction law [eg. Kato and Hirasawa (1997), Hirose and Hirahara (2002)]. All of them assume that pore fluid pressure in the fault zone is constant. However, in the fault zone, pore fluid pressure changes suddenly, due to coseismic pore dilatation [Marone (1990)] and thermal pressurization [Mase and Smith (1987)]. If pore fluid pressure drops and effective normal stress rises, fault slip is decelerated. Inversely, if pore fluid pressure rises and effective normal stress drops, fault slip is accelerated. The effect of pore fluid may cause slow slip events and low-frequency tremor [Kodaira et al. (2004), Shelly et al. (2006)]. For a simple spring model, how pore dilatation affects slip instability was investigated [Segall and Rice (1995), Sleep (1995)]. When the rate of the slip becomes high, pore dilatation occurs and pore pressure drops, and the rate of the slip is restrained. Then the inflow of pore fluid recovers the pore pressure. We execute 2D earthquake cycle simulations at a subduction zone, taking into account such changes of pore fluid pressure following Segall and Rice (1995), in addition to the numerical scheme in Kato and Hirasawa (1997). We do not adopt hydrostatic pore pressure but excess pore pressure for initial condition, because upflow of dehydrated water seems to exist at a subduction zone. In our model, pore fluid is confined to the fault damage zone and flows along the plate interface. The smaller the flow rate is, the later pore pressure recovers. Since effective normal stress keeps larger, the fault slip is decelerated and stress drop becomes smaller. Therefore the smaller flow rate along the fault zone leads to the shorter earthquake recurrence time. Thus, not only the frictional parameters and the subduction rate but also the fault zone permeability affects the recurrence time of earthquake cycle. Further, the existence of heterogeneity in the permeability along the plate interface can bring about other slip behaviors, such as slow slip events. Our simulations indicate that, in addition to the frictional parameters, the permeability within the fault damage zone is one of essential parameters, which controls the whole earthquake cycle.

FINE PORE DIFFUSER FOULING: THE LOS ANGELES STUDIES

EPA Science Inventory

This report describes five fine pore diffuser evaluations conducted at three different wastewater treatment plants located in the greater Los Angeles area. The overall goal of the study was to evaluate the performance of fine pore diffusers using selected cleaning methods for ex...

On N. Park's Analytical solution for steady state density- and mixing regime—dependent solute transport in a vertical soil column

NASA Astrophysics Data System (ADS)

Thiele, Michael

1998-04-01

Recently, Park [1996] presented an analytical solution for stationary one-dimensional solute transport in a variable-density fluid flow through a vertical soil column. He used the widespread Bear-Scheidegger dispersion model describing solute mixing as a sum of molecular diffusion and velocity-proportional mechanical dispersion effects. His closed-form implicit concentration and pressure distributions thus allow for a discussion of the combined impact of molecular diffusion and mechanical dispersion in a variable-density environment. Whereas Park only considered the example of vanishing molecular diffusion in detail, both phenomena are taken into account simultaneously in the present study in order to elucidate their different influences on concentration distribution characteristics. The boundary value problem dealt with herein is based on an upward inflow of high-density fluid of constant solute concentration and corresponding outflow of a lower constant concentration fluid at the upper end of the column when dispersivity does not change along the flow path. The thickness of the transition zone between the two fluids appeared to strongly depend on the prevailing share of the molecular diffusion and mechanical dispersion mechanisms. The latter can be characterized by a molecular Peclet number Pe, which here is defined as the ratio of the column outflow velocity multiplied by a characteristic pore size and the molecular diffusion coefficient. For very small values of Pe, when molecular diffusion represents the exclusive mixing process, density differences have no impact on transition zone thicknesses. A relative density-;dependent thickness increases with flow velocities (increasing Pe values) very rapidly compared to the density-independent case, and after having passed a maximum decreases asymptotically to a constant value for the large Peclet number limit when mechanical dispersion is the only mixing mechanism. Hence the special transport problem analyzed gives further evidence for the importance of simultaneously considering molecular diffusion and mechanical dispersion in gravity-affected solute transport in porous media.

Calculation of effective transport properties of partially saturated gas diffusion layers

NASA Astrophysics Data System (ADS)

Bednarek, Tomasz; Tsotridis, Georgios

2017-02-01

A large number of currently available Computational Fluid Dynamics numerical models of Polymer Electrolyte Membrane Fuel Cells (PEMFC) are based on the assumption that porous structures are mainly considered as thin and homogenous layers, hence the mass transport equations in structures such as Gas Diffusion Layers (GDL) are usually modelled according to the Darcy assumptions. Application of homogenous models implies that the effects of porous structures are taken into consideration via the effective transport properties of porosity, tortuosity, permeability (or flow resistance), diffusivity, electric and thermal conductivity. Therefore, reliable values of those effective properties of GDL play a significant role for PEMFC modelling when employing Computational Fluid Dynamics, since these parameters are required as input values for performing the numerical calculations. The objective of the current study is to calculate the effective transport properties of GDL, namely gas permeability, diffusivity and thermal conductivity, as a function of liquid water saturation by using the Lattice-Boltzmann approach. The study proposes a method of uniform water impregnation of the GDL based on the "Fine-Mist" assumption by taking into account the surface tension of water droplets and the actual shape of GDL pores.

Acoustics of multiscale sorptive porous materials

NASA Astrophysics Data System (ADS)

Venegas, R.; Boutin, C.; Umnova, O.

2017-08-01

This paper investigates sound propagation in multiscale rigid-frame porous materials that support mass transfer processes, such as sorption and different types of diffusion, in addition to the usual visco-thermo-inertial interactions. The two-scale asymptotic method of homogenization for periodic media is successively used to derive the macroscopic equations describing sound propagation through the material. This allowed us to conclude that the macroscopic mass balance is significantly modified by sorption, inter-scale (micro- to/from nanopore scales) mass diffusion, and inter-scale (pore to/from micro- and nanopore scales) pressure diffusion. This modification is accounted for by the dynamic compressibility of the effective saturating fluid that presents atypical properties that lead to slower speed of sound and higher sound attenuation, particularly at low frequencies. In contrast, it is shown that the physical processes occurring at the micro-nano-scale do not affect the macroscopic fluid flow through the material. The developed theory is exemplified by introducing an analytical model for multiscale sorptive granular materials, which is experimentally validated by comparing its predictions with acoustic measurements on granular activated carbons. Furthermore, we provide empirical evidence supporting an alternative method for measuring sorption and mass diffusion properties of multiscale sorptive materials using sound waves.

Finite Element Analysis of Poroelastic Composites Undergoing Thermal and Gas Diffusion

NASA Technical Reports Server (NTRS)

Salamon, N. J. (Principal Investigator); Sullivan, Roy M.; Lee, Sunpyo

1995-01-01

A theory for time-dependent thermal and gas diffusion in mechanically time-rate-independent anisotropic poroelastic composites has been developed. This theory advances previous work by the latter two authors by providing for critical transverse shear through a three-dimensional axisymmetric formulation and using it in a new hypothesis for determining the Biot fluid pressure-solid stress coupling factor. The derived governing equations couple material deformation with temperature and internal pore pressure and more strongly couple gas diffusion and heat transfer than the previous theory. Hence the theory accounts for the interactions between conductive heat transfer in the porous body and convective heat carried by the mass flux through the pores. The Bubnov Galerkin finite element method is applied to the governing equations to transform them into a semidiscrete finite element system. A numerical procedure is developed to solve the coupled equations in the space and time domains. The method is used to simulate two high temperature tests involving thermal-chemical decomposition of carbon-phenolic composites. In comparison with measured data, the results are accurate. Moreover unlike previous work, for a single set of poroelastic parameters, they are consistent with two measurements in a restrained thermal growth test.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Garcia, Andres

Transport and reaction in zeolites and other porous materials, such as mesoporous silica particles, has been a focus of interest in recent years. This is in part due to the possibility of anomalous transport effects (e.g. single-file diffusion) and its impact in the reaction yield in catalytic processes. Computational simulations are often used to study these complex nonequilibrium systems. Computer simulations using Molecular Dynamics (MD) techniques are prohibitive, so instead coarse grained one-dimensional models with the aid of Kinetic Monte Carlo (KMC) simulations are used. Both techniques can be computationally expensive, both time and resource wise. These coarse-grained systems canmore » be exactly described by a set of coupled stochastic master equations, that describe the reaction-diffusion kinetics of the system. The equations can be written exactly, however, coupling between the equations and terms within the equations make it impossible to solve them exactly; approximations must be made. One of the most common methods to obtain approximate solutions is to use Mean Field (MF) theory. MF treatments yield reasonable results at high ratios of reaction rate k to hop rate h of the particles, but fail completely at low k=h due to the over-estimation of fluxes of particles within the pore. We develop a method to estimate fluxes and intrapore diffusivity in simple one- dimensional reaction-diffusion models at high and low k=h, where the pores are coupled to an equilibrated three-dimensional fluid. We thus successfully describe analytically these simple reaction-diffusion one-dimensional systems. Extensions to models considering behavior with long range steric interactions and wider pores require determination of multiple boundary conditions. We give a prescription to estimate the required parameters for these simulations. For one dimensional systems, if single-file diffusion is relaxed, additional parameters to describe particle exchange have to be introduced. We use Langevin Molecular Dynamics (MD) simulations to assess these parameters.« less

Flow of variably fluidized granular masses across three-dimensional terrain I. Coulomb mixture theory

USGS Publications Warehouse

Iverson, R.M.; Denlinger, R.P.

2001-01-01

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces, govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, threedimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Flow of variably fluidized granular masses across three-dimensional terrain: 1. Coulomb mixture theory

NASA Astrophysics Data System (ADS)

Iverson, Richard M.; Denlinger, Roger P.

2001-01-01

Rock avalanches, debris flows, and related phenomena consist of grain-fluid mixtures that move across three-dimensional terrain. In all these phenomena the same basic forces govern motion, but differing mixture compositions, initial conditions, and boundary conditions yield varied dynamics and deposits. To predict motion of diverse grain-fluid masses from initiation to deposition, we develop a depth-averaged, three-dimensional mathematical model that accounts explicitly for solid- and fluid-phase forces and interactions. Model input consists of initial conditions, path topography, basal and internal friction angles of solid grains, viscosity of pore fluid, mixture density, and a mixture diffusivity that controls pore pressure dissipation. Because these properties are constrained by independent measurements, the model requires little or no calibration and yields readily testable predictions. In the limit of vanishing Coulomb friction due to persistent high fluid pressure the model equations describe motion of viscous floods, and in the limit of vanishing fluid stress they describe one-phase granular avalanches. Analysis of intermediate phenomena such as debris flows and pyroclastic flows requires use of the full mixture equations, which can simulate interaction of high-friction surge fronts with more-fluid debris that follows. Special numerical methods (described in the companion paper) are necessary to solve the full equations, but exact analytical solutions of simplified equations provide critical insight. An analytical solution for translational motion of a Coulomb mixture accelerating from rest and descending a uniform slope demonstrates that steady flow can occur only asymptotically. A solution for the asymptotic limit of steady flow in a rectangular channel explains why shear may be concentrated in narrow marginal bands that border a plug of translating debris. Solutions for static equilibrium of source areas describe conditions of incipient slope instability, and other static solutions show that nonuniform distributions of pore fluid pressure produce bluntly tapered vertical profiles at the margins of deposits. Simplified equations and solutions may apply in additional situations identified by a scaling analysis. Assessment of dimensionless scaling parameters also reveals that miniature laboratory experiments poorly simulate the dynamics of full-scale flows in which fluid effects are significant. Therefore large geophysical flows can exhibit dynamics not evident at laboratory scales.

Biogeochemical cycles at the sulfate-methane transition zone (SMTZ) and geochemical characteristics of the pore fluids offshore southwestern Taiwan

NASA Astrophysics Data System (ADS)

Hu, Ching-Yi; Frank Yang, Tsanyao; Burr, George S.; Chuang, Pei-Chuan; Chen, Hsuan-Wen; Walia, Monika; Chen, Nai-Chen; Huang, Yu-Chun; Lin, Saulwood; Wang, Yunshuen; Chung, San-Hsiung; Huang, Chin-Da; Chen, Cheng-Hong

2017-11-01

In this study, we used pore water dissolved inorganic carbon (DIC), SO42-, Ca2+ and Mg2+ gradients at the sulfate-methane transition zone (SMTZ) to estimate biogeochemical fluxes for cored sediments collected offshore SW Taiwan. Net DIC flux changes (ΔDIC-Prod) were applied to determine the proportion of sulfate consumption by organic matter oxidation (heterotrophic sulfate reduction) and anaerobic oxidation of methane (AOM), and to determine reliable CH4 fluxes at the SMTZ. Our results show that SO42- profiles are mainly controlled by AOM rather than heterotrophic sulfate reduction. Refinement of CH4 flux estimates enhance our understanding of methane abundance from deep carbon reservoirs to the SMTZ. Concentrations of chloride (Cl-), bromide (Br-) and iodide (I-) dissolved in pore water were used to identify potential sources that control fluid compositions and the behavior of dissolved ions. Constant Cl- concentrations throughout ∼30 m sediment suggest no influence of gas hydrates for the compositions within the core. Bromide (Br-) and Iodine (I-) concentrations increase with sediment depth. The I-/Br- ratio appears to reflect organic matter degradation. SO42- concentrations decrease with sediment depth at a constant rate, and sediment depth profiles of Br- and I- concentrations suggests diffusion as the main transport mechanism. Therefore diffusive flux calculations are reasonable. Coring sites with high CH4 fluxes are more common in the accretionary wedge, amongst thrust faults and fractures, than in the passive continental margin offshore southwestern Taiwan. AOM reactions are a major sink for CH4 passing upward through the SMTZ and prevent high methane fluxes in the water column and to the atmosphere.

Relation Between Pore Size and the Compressibility of a Confined Fluid

PubMed Central

Gor, Gennady Y.; Siderius, Daniel W.; Rasmussen, Christopher J.; Krekelberg, William P.; Shen, Vincent K.; Bernstein, Noam

2015-01-01

When a fluid is confined to a nanopore, its thermodynamic properties differ from the properties of a bulk fluid, so measuring such properties of the confined fluid can provide information about the pore sizes. Here we report a simple relation between the pore size and isothermal compressibility of argon confined in these pores. Compressibility is calculated from the fluctuations of the number of particles in the grand canonical ensemble using two different simulation techniques: conventional grand-canonical Monte Carlo and grand-canonical ensemble transition-matrix Monte Carlo. Our results provide a theoretical framework for extracting the information on the pore sizes of fluid-saturated samples by measuring the compressibility from ultrasonic experiments. PMID:26590541

Glasslike Membrane Protein Diffusion in a Crowded Membrane.

PubMed

Munguira, Ignacio; Casuso, Ignacio; Takahashi, Hirohide; Rico, Felix; Miyagi, Atsushi; Chami, Mohamed; Scheuring, Simon

2016-02-23

Many functions of the plasma membrane depend critically on its structure and dynamics. Observation of anomalous diffusion in vivo and in vitro using fluorescence microscopy and single particle tracking has advanced our concept of the membrane from a homogeneous fluid bilayer with freely diffusing proteins to a highly organized crowded and clustered mosaic of lipids and proteins. Unfortunately, anomalous diffusion could not be related to local molecular details given the lack of direct and unlabeled molecular observation capabilities. Here, we use high-speed atomic force microscopy and a novel analysis methodology to analyze the pore forming protein lysenin in a highly crowded environment and document coexistence of several diffusion regimes within one membrane. We show the formation of local glassy phases, where proteins are trapped in neighbor-formed cages for time scales up to 10 s, which had not been previously experimentally reported for biological membranes. Furthermore, around solid-like patches and immobile molecules a slower glass phase is detected leading to protein trapping and creating a perimeter of decreased membrane diffusion.

Pore fluid pressure and the seismic cycle

NASA Astrophysics Data System (ADS)

French, M. E.; Zhu, W.; Hirth, G.; Belzer, B.

2017-12-01

In the brittle crust, the critical shear stress required for fault slip decreases with increasing pore fluid pressures according to the effective stress criterion. As a result, higher pore fluid pressures are thought to promote fault slip and seismogenesis, consistent with observations that increasing fluid pressure as a result of wastewater injection is correlated with increased seismicity. On the other hand, elevated pore fluid pressure is also proposed to promote slow stable failure rather than seismicity along some fault zones, including during slow slip in subduction zones. Here we review recent experimental evidence for the roles that pore fluid pressure and the effective stress play in controlling fault slip behavior. Using two sets of experiments on serpentine fault gouge, we show that increasing fluid pressure does decrease the shear stress for reactivation under brittle conditions. However, under semi-brittle conditions as expected near the base of the seismogenic zone, high pore fluid pressures are much less effective at reducing the shear stress of reactivation even though deformation is localized and frictional. We use an additional study on serpentinite to show that cohesive fault rocks, potentially the product of healing and cementation, experience an increase in fracture energy during faulting as fluid pressures approach lithostatic, which can lead to more stable failure. Structural observations show that the increased fracture energy is associated with a greater intensity of transgranular fracturing and delocalization of deformation. Experiments on several lithologies indicate that the stabilizing effect of fluid pressure occurs independent of rock composition and hydraulic properties. Thus, high pore fluid pressures have the potential to either enhance seismicity or promote stable faulting depending on pressure, temperature, and fluid pressure conditions. Together, the results of these studies indicate that pore fluid pressure promotes seismogenesis in the brittle shallow crust where fluid pressures are elevated but sub-lithostatic and promote slow, stable failure near seismic to aseismic transitions and under near-lithostatic fluid pressures.

Subcritical growth of natural hydraulic fractures

NASA Astrophysics Data System (ADS)

Garagash, D.

2014-12-01

Joints are the most common example of brittle tensile failure in the crust. Their genesis at depth is linked to the natural hydraulic fracturing, which requires pore fluid pressure in excess of the minimum in situ stress [Pollard and Aidyn, JSG1988]. Depending on the geological setting, high pore pressure can result form burial compaction of interbedded strata, diagenesis, or tectonics. Common to these loading scenarios is slow build-up of pore pressure over a geological timescale, until conditions for initiation of crack growth are met on favorably oriented/sized flaws. The flaws can vary in size from grain-size cracks in igneous rocks to a fossil-size flaws in clastic rock, and once activated, are inferred to propagate mostly subcritically [Segall JGR 1984; Olson JGR 1993]. Despite many observational studies of natural hydraulic fractures, the modeling attempts appear to be few [Renshaw and Harvey JGR 1994]. Here, we use boundary integral formulation for the pore fluid inflow from the permeable rock into a propagating joint [Berchenko et al. IJRMMS 1997] coupled with the criteria for subcritical propagation assisted by the environmental effects of pore fluid at the crack tip to solve for the evolution of a penny-shape joint, which, in interbedded rock, may eventually evolve to short-blade geometry (propagation confined to a bed). Initial growth is exceedingly slow, paced by the stress corrosion reaction kinetics at the crack tip. During this stage the crack is fully-drained (i.e. the fluid pressure in the crack is equilibrated with the ambient pore pressure). This "slow" stage is followed by a rapid acceleration, driven by the increase of the mechanical stress intensity factor with the crack length, towards the terminal joint velocity. We provide an analytical expression for the latter as a function of the rock diffusivity, net pressure loading at the initiation (or flaw lengthscale), and parameters describing resistance to fracture growth. Due to a much slower rate of the crack volume expansion of short-blade joints compared to that of penny-shape joints, the former would propagate much faster than the latter under otherwise identical conditions. Finally, we speculate about possible relation of the predicted patterns of joint development with morphology of joint fracture surfaces observed in sedimentary rock.

Dynamic Stability of the Rate, State, Temperature, and Pore Pressure Friction Model at a Rock Interface

NASA Astrophysics Data System (ADS)

Sinha, Nitish; Singh, Arun K.; Singh, Trilok N.

2018-05-01

In this article, we study numerically the dynamic stability of the rate, state, temperature, and pore pressure friction (RSTPF) model at a rock interface using standard spring-mass sliding system. This particular friction model is a basically modified form of the previously studied friction model namely the rate, state, and temperature friction (RSTF). The RSTPF takes into account the role of thermal pressurization including dilatancy and permeability of the pore fluid due to shear heating at the slip interface. The linear stability analysis shows that the critical stiffness, at which the sliding becomes stable to unstable or vice versa, increases with the coefficient of thermal pressurization. Critical stiffness, on the other hand, remains constant for small values of either dilatancy factor or hydraulic diffusivity, but the same decreases as their values are increased further from dilatancy factor (˜ 10^{ - 4} ) and hydraulic diffusivity (˜ 10^{ - 9} {m}2 {s}^{ - 1} ) . Moreover, steady-state friction is independent of the coefficient of thermal pressurization, hydraulic diffusivity, and dilatancy factor. The proposed model is also used for predicting time of failure of a creeping interface of a rock slope under the constant gravitational force. It is observed that time of failure decreases with increase in coefficient of thermal pressurization and hydraulic diffusivity, but the dilatancy factor delays the failure of the rock fault under the condition of heat accumulation at the creeping interface. Moreover, stiffness of the rock-mass also stabilizes the failure process of the interface as the strain energy due to the gravitational force accumulates in the rock-mass before it transfers to the sliding interface. Practical implications of the present study are also discussed.

Measuring temporal variability in pore-fluid chemistry to assess gas hydrate stability: development of a continuous pore-fluid array.

PubMed

Lapham, Laura L; Chanton, Jeffrey P; Martens, Christopher S; Higley, Paul D; Jannasch, Hans W; Woolsey, J Robert

2008-10-01

A specialized pore-fluid array (PFA) sampler was designed to collect and store pore fluids to monitor temporal changes of ions and gases in gas hydrate bearing sediments. We tested the hypothesis that pore-fluid chemistry records hydrate formation or decomposition events and reflects local seismic activity. The PFA is a seafloor probe that consists of an interchangeable instrument package that houses OsmoSamplers, long-term pore-fluid samplers, a specialized low-dead volume fluid coupler, and eight sample ports along a 10 m sediment probe shaft. The PFA was deployed at Mississippi Canyon 118, a Gulf of Mexico hydrate site. A 170 day record was acquired from the overlying water and 1.3 m below seafloor (mbsf). Fluids were measured for dissolved chloride, sulfate, and methane concentrations and dissolved inorganic carbon and methane stable carbon and deuterium isotope ratios. Chloride and sulfate did not change significantly over time, suggesting the absence of gas hydrate formation or decomposition events. Over the temporal record, methane concentrations averaged 4 mM at 1.3 mbsf, and methane was thermogenic in origin (delta13C-CH4 = -32.4 +/- 3.4 per thousand). The timing of an anomalous 14 mM methane spike coincided with a nearby earthquake (Mw = 5.8), consistent with the hypothesis that pore-fluid chemistry reflects seismic events.

Molecular simulation of simple fluids and polymers in nanoconfinement

NASA Astrophysics Data System (ADS)

Rasmussen, Christopher John

Prediction of phase behavior and transport properties of simple fluids and polymers confined to nanoscale pores is important to a wide range of chemical and biochemical engineering processes. A practical approach to investigate nanoscale systems is molecular simulation, specifically Monte Carlo (MC) methods. One of the most challenging problems is the need to calculate chemical potentials in simulated phases. Through the seminal work of Widom, practitioners have a powerful method for calculating chemical potentials. Yet, this method fails for dense and inhomogeneous systems, as well as for complex molecules such as polymers. In this dissertation, the gauge cell MC method, which had previously been successfully applied to confined simple fluids, was employed and extended to investigate nanoscale fluids in several key areas. Firstly, the process of cavitation (the formation and growth of bubbles) during desorption of fluids from nanopores was investigated. The dependence of cavitation pressure on pore size was determined with gauge cell MC calculations of the nucleation barriers correlated with experimental data. Additional computational studies elucidated the role of surface defects and pore connectivity in the formation of cavitation bubbles. Secondly, the gauge cell method was extended to polymers. The method was verified against the literature results and found significantly more efficient. It was used to examine adsorption of polymers in nanopores. These results were applied to model the dynamics of translocation, the act of a polymer threading through a small opening, which is implicated in drug packaging and delivery, and DNA sequencing. Translocation dynamics was studied as diffusion along the free energy landscape. Thirdly, we show how computer simulation of polymer adsorption could shed light on the specifics of polymer chromatography, which is a key tool for the analysis and purification of polymers. The quality of separation depends on the physico-chemical mechanisms of polymer/pore interaction. We considered liquid chromatography at critical conditions, and calculated the dependence of the partition coefficient on chain length. Finally, solvent-gradient chromatography was modeled using a statistical model of polymer adsorption. A model for predicting separation of complex polymers (with functional groups or copolymers) was developed for practical use in chromatographic separations.

In-pore exchange and diffusion of carbonate solvent mixtures in nanoporous carbon

DOE PAGES

Alam, Todd M.; Osborn Popp, Thomas M.

2016-06-04

High resolution magic angle spinning (HRMAS) 1H NMR spectroscopy has been used to resolve different surface and in-pore solvent environments of ethylene carbonate (EC) and dimethyl carbonate (DMC) mixtures absorbed within nanoporous carbon (NPC). Two dimensional (2D) 1H HRMAS NMR exchange measurements revealed that the inhomogeneous broadened in-pore resonances have pore-to-pore exchange rates on the millisecond timescale. Pulsed-field gradient (PFG) NMR diffusometry revealed the in-pore self-diffusion constants for both EC and DMC were reduced by up to a factor of five with respect to the diffusion in the non-absorbed solvent mixtures.

Cuticular gas exchange by Antarctic sea spiders.

PubMed

Lane, Steven J; Moran, Amy L; Shishido, Caitlin M; Tobalske, Bret W; Woods, H Arthur

2018-04-25

Many marine organisms and life stages lack specialized respiratory structures, like gills, and rely instead on cutaneous respiration, which they facilitate by having thin integuments. This respiratory mode may limit body size, especially if the integument also functions in support or locomotion. Pycnogonids, or sea spiders, are marine arthropods that lack gills and rely on cutaneous respiration but still grow to large sizes. Their cuticle contains pores, which may play a role in gas exchange. Here, we examined alternative paths of gas exchange in sea spiders: (1) oxygen diffuses across pores in the cuticle, a common mechanism in terrestrial eggshells, (2) oxygen diffuses directly across the cuticle, a common mechanism in small aquatic insects, or (3) oxygen diffuses across both pores and cuticle. We examined these possibilities by modeling diffusive oxygen fluxes across all pores in the body of sea spiders and asking whether those fluxes differed from measured metabolic rates. We estimated fluxes across pores using Fick's law parameterized with measurements of pore morphology and oxygen gradients. Modeled oxygen fluxes through pores closely matched oxygen consumption across a range of body sizes, which means the pores facilitate oxygen diffusion. Furthermore, pore volume scaled hypermetrically with body size, which helps larger species facilitate greater diffusive oxygen fluxes across their cuticle. This likely presents a functional trade-off between gas exchange and structural support, in which the cuticle must be thick enough to prevent buckling due to external forces but porous enough to allow sufficient gas exchange. © 2018. Published by The Company of Biologists Ltd.

Generalized Hydrodynamic Treatment of the Interplay between Restricted Transport and Catalytic Reactions in Nanoporous Materials

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ackerman, David M.; Wang, Jing; Evans, James W.

2012-05-30

Behavior of catalytic reactions in narrow pores is controlled by a delicate interplay between fluctuations in adsorption-desorption at pore openings, restricted diffusion, and reaction. This behavior is captured by a generalized hydrodynamic formulation of appropriate reaction-diffusion equations (RDE). These RDE incorporate an unconventional description of chemical diffusion in mixed-component quasi-single-file systems based on a refined picture of tracer diffusion for finite-length pores. The RDE elucidate the nonexponential decay of the steady-state reactant concentration into the pore and the non-mean-field scaling of the reactant penetration depth.

Generalized hydrodynamic treatment of the interplay between restricted transport and catalytic reactions in nanoporous materials.

PubMed

Ackerman, David M; Wang, Jing; Evans, James W

2012-06-01

Behavior of catalytic reactions in narrow pores is controlled by a delicate interplay between fluctuations in adsorption-desorption at pore openings, restricted diffusion, and reaction. This behavior is captured by a generalized hydrodynamic formulation of appropriate reaction-diffusion equations (RDE). These RDE incorporate an unconventional description of chemical diffusion in mixed-component quasi-single-file systems based on a refined picture of tracer diffusion for finite-length pores. The RDE elucidate the nonexponential decay of the steady-state reactant concentration into the pore and the non-mean-field scaling of the reactant penetration depth.

Pathways of fluid transport and reabsorption across the peritoneal membrane.

PubMed

Asghar, R B; Davies, S J

2008-05-01

The three-pore model of peritoneal fluid transport predicts that once the osmotic gradient has dissipated, fluid reabsorption will be due to a combination of small-pore reabsorption driven by the intravascular oncotic pressure, and an underlying disappearance of fluid from the cavity by lymphatic drainage. Our study measured fluid transport by these pathways in the presence and absence of an osmotic gradient. Paired hypertonic and standard glucose-dwell studies were performed using radio-iodinated serum albumin as an intraperitoneal volume marker and changes in intraperitoneal sodium mass to determine small-pore versus transcellular fluid transport. Disappearance of iodinated albumin was considered to indicate lymphatic drainage. Variability in transcellular ultrafiltration was largely explained by the rate of small-solute transport across the membrane. In the absence of an osmotic gradient, fluid reabsorption occurred via the small-pore pathway, the rate being proportional to the small-solute transport characteristics of the membrane. In most cases, fluid removal from the peritoneal cavity by this pathway was faster than by lymphatic drainage. Our study shows that the three-pore model describes the pathways of peritoneal fluid transport well. In the presence of high solute transport, poor transcellular ultrafiltration was due to loss of the osmotic gradient and an enhanced small-pore reabsorption rate after this gradient dissipated.

Effect of Polydispersity on Diffusion in Random Obstacle Matrices

NASA Astrophysics Data System (ADS)

Cho, Hyun Woo; Kwon, Gyemin; Sung, Bong June; Yethiraj, Arun

2012-10-01

The dynamics of tracers in disordered matrices is of interest in a number of diverse areas of physics such as the biophysics of crowding in cells and cell membranes, and the diffusion of fluids in porous media. To a good approximation the matrices can be modeled as a collection of spatially frozen particles. In this Letter, we consider the effect of polydispersity (in size) of the matrix particles on the dynamics of tracers. We study a two dimensional system of hard disks diffusing in a sea of hard disk obstacles, for different values of the polydispersity of the matrix. We find that for a given average size and area fraction, the diffusion of tracers is very sensitive to the polydispersity. We calculate the pore percolation threshold using Apollonius diagrams. The diffusion constant, D, follows a scaling relation D˜(ϕc-ϕm)μ-β for all values of the polydispersity, where ϕm is the area fraction and ϕc is the value of ϕm at the percolation threshold.

Effect of polydispersity on diffusion in random obstacle matrices.

PubMed

Cho, Hyun Woo; Kwon, Gyemin; Sung, Bong June; Yethiraj, Arun

2012-10-12

The dynamics of tracers in disordered matrices is of interest in a number of diverse areas of physics such as the biophysics of crowding in cells and cell membranes, and the diffusion of fluids in porous media. To a good approximation the matrices can be modeled as a collection of spatially frozen particles. In this Letter, we consider the effect of polydispersity (in size) of the matrix particles on the dynamics of tracers. We study a two dimensional system of hard disks diffusing in a sea of hard disk obstacles, for different values of the polydispersity of the matrix. We find that for a given average size and area fraction, the diffusion of tracers is very sensitive to the polydispersity. We calculate the pore percolation threshold using Apollonius diagrams. The diffusion constant, D, follows a scaling relation D~(φ(c)-φ(m))(μ-β) for all values of the polydispersity, where φ(m) is the area fraction and φ(c) is the value of φ(m) at the percolation threshold.

Evidence for Enhanced Matrix Diffusion in Geological Environment

NASA Astrophysics Data System (ADS)

Sato, Kiminori; Fujimoto, Koichiro; Nakata, Masataka; Shikazono, Naotatsu

2013-01-01

Molecular diffusion in rock matrix, called as matrix diffusion, has been appreciated as a static process for elemental migration in geological environment that has been acknowledged in the context of geological disposal of radioactive waste. However, incomprehensible enhancement of matrix diffusion has been reported at a number of field test sites. Here, the matrix diffusion of saline water at Horonobe, Hokkaido, Japan is highlighted directly probing angstrom-scale pores on a field scale up to 1 km by positron--positronium annihilation spectroscopy. The first application of positron--positronium annihilation spectroscopy to field-scale geophysical research reveals the slight variation of angstrom-scale pores influenced by saline water diffusion with complete accuracy. We found widely interconnected 3 Å pores, which offer the pathway of saline water diffusion with the highly enhanced effective matrix diffusion coefficient of 4× 10-6 cm2 s-1. The present findings provide unambiguous evidence that the angstrom-scale pores enhance effective matrix diffusion on a field scale in geological environment.

Decline of a Hydrothermal Vent Field - Escanaba Trough 12 Years Later

NASA Astrophysics Data System (ADS)

Zierenberg, R. A.; Clague, D. A.; Davis, A. S.; Lilley, M. D.; McClain, J. S.; Olson, E. S.; Ross, S. L.; Von Damm, K. L.

2001-12-01

Hydrothermal venting was discovered in Escanaba Trough, the southern sediment-covered portion of the Gorda Ridge, in 1988. Large pyrrhotite-rich massive sulfide mounds are abundant at each of the volcanic/intrusive centers that have been investigated in Escanaba Trough, but the only area of known hydrothermal venting is the NESCA site along the ridge axis at 41\\deg N. Hydrothermal fluids venting at 217\\deg C and 108\\deg C were sampled in 1988 on two sulfide mounds separated by about 275 m. The end-member fluid compositions were indistinguishable within analytical errors. Several sulfide mounds were observed in 1988 which had diffusely venting low temperature (< 20\\deg C) fluids that supported extensive vent communities dominated by fields of Ridgia. Nine holes were drilled in the NESCA area in 1996 on ODP Leg 169, including Hole 1036I that penetrated to basaltic basement at 405 m below sea floor (mbsf). Surveys of the area using the drill string camera located only one area of active venting at the same mound where 217\\deg C vent fluids were sampled from two active vents in 1988. Drill hole 1036A was spudded between the two active vents on this sulfide mound (approximately 4 and 8 m away) and penetrated to 115 mbsf. The NESCA site was revisited in 2000 using MBARI's R/V Western Flyer and ROV Tiburon. The hydrothermal vents appeared essentially identical to observations made from the drill string camera in 1996 despite the presence of a drill hole within meters of the two vents. The maximum vent temperature measured in 2000 was 212\\deg C. Fluid samples have major element and isotopic compositions very similar to those collected in 1988. The vent fluids have higher methane ( ~19 mmol/kg) than those from the geologically similar Middle Valley vent field, but lower values than those at Guaymas Basin. Drill hole 1036A was weakly venting, but the diffuse hydrothermal fluids could not be sampled with the equipment available. The walls of the drill hole were colonized by palm worms, limpets, and snails. Four other drill holes showed no hydrothermal flow nor visible evidence of down hole recharge. Mapping with Tiburon confirmed that the extent of hydrothermal venting at NESCA decreased dramatically since 1988. Formerly extensive colonies of Ridgia had vanished leaving no trace of their presence. Although hydrothermal venting has collapsed to a single mound, the temperature and composition of the fluids remained nearly unchanged. This is curious given that sediment pore fluids analyzed on Leg 169 included both high salinity and low salinity components of phase separated hydrothermal fluids in the shallow subsurface indicating that the hydrothermal field must have had a relatively recent (relative to the rate of pore fluid diffusion) high temperature history. Hydrothermal fluids presently venting at this site must be derived from an essentially homogeneous, approximately 215\\degC fluid reservoir that has declined in its fluid output on a decadal scale, but has not undergone significant changes in temperature and composition. Venting at the seafloor does not seem to have been affected by drilling in the hydrothermal field.

Enhanced transport of materials into enamel nanopores via electrokinetic flow.

PubMed

Gan, H Y; Sousa, F B; Carlo, H L; Maciel, P P; Macena, M S; Han, J

2015-04-01

The ability to infiltrate various molecules and resins into dental enamel is highly desirable in dentistry, yet transporting materials into dental enamel is limited by the nanometric scale of their pores. Materials that cannot be infiltrated into enamel by diffusion/capillarity are often considered molecules with sizes above a critical threshold, which are often considered to be larger than the pores of enamel. We challenge this notion by reporting the use of electrokinetic flow to transport solutions with molecules with sizes above a critical threshold-namely, an aqueous solution with a high refractive index (Thoulet's solution) and a curable fluid resin infiltrant (without acid etching)-deep into the normal enamel layer. Volume infiltration by Thoulet's solution is increased by 5- to 6-fold, and resin infiltration depths as large as 600 to 2,000 µm were achieved, in contrast to ~10 µm resulting from diffusion/capillarity. Incubation with demineralization solution for 192 h resulted in significant demineralization at noninfiltrated histologic points but not at resin infiltrated. These results open new avenues for the transport of materials in dental enamel. © International & American Associations for Dental Research 2015.

The impact of flow focusing on gas hydrate accumulations in overpressured marine sediments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nole, Michael; Daigle, Hugh; Cook, Ann

This study demonstrates the potential for flow focusing due to overpressuring in marine sedimentary environments to act as a significant methane transport mechanism from which methane hydrate can precipitate in large quantities in dipping sandstone bodies. Traditionally, gas hydrate accumulations in nature are discussed as resulting from either short-range diffusive methane migration or from long-range advective fluid transport sourced from depth. However, 3D simulations performed in this study demonstrate that a third migration mechanism, short-range advective transport, can provide a significant methane source that is unencumbered by limitations of the other two end-member mechanisms. Short-range advective sourcing is advantageous overmore » diffusion because it can convey greater amounts of methane to sands over shorter timespans, yet it is not necessarily limited by down-dip pore blocking in sands as is typical of updip advection from a deep source. These results are novel because they integrate pore size impacts on spatial solubility gradients, grid block properties that evolve through time, and methane sourcing through microbial methanogenesis into a holistic characterization of environments exposed to multiple methane hydrate sourcing mechanisms. We show that flow focusing toward sand bodies transports large quantities of methane, the magnitude of which are determined by the sand-clay solubility contrast, and generates larger quantities of hydrate in sands than a solely diffusive system; after depositing methane as hydrate, fluid exiting a sand body is depleted in methane and leaves a hydrate free region in its wake above the sand. Additionally, we demonstrate that in overpressured environments, hydrate growth is initially diffusively dominated before transitioning to an advection-dominated regime. The timescale and depth at which this transition takes place depends primarily on the rate of microbial metabolism and the sedimentation rate but only depends loosely on the degree of overpressuring.« less

INVESTIGATIONS INTO BIOFOULING PHENOMENA IN FINE PORE AERATION DEVICES

EPA Science Inventory

Microbiologically-based procedures were used to describe biofouling phenomena on fine pore aeration devices and to determine whether biofilm characteristics could be related to diffuser process performance parameters. Fine pore diffusers were obtained from five municipal wastewa...

Studying methane migration mechanisms at Walker Ridge, Gulf of Mexico, via 3D methane hydrate reservoir modeling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Nole, Michael; Daigle, Hugh; Mohanty, Kishore

We have developed a 3D methane hydrate reservoir simulator to model marine methane hydrate systems. Our simulator couples highly nonlinear heat and mass transport equations and includes heterogeneous sedimentation, in-situ microbial methanogenesis, the influence of pore size contrast on solubility gradients, and the impact of salt exclusion from the hydrate phase on dissolved methane equilibrium in pore water. Using environmental parameters from Walker Ridge in the Gulf of Mexico, we first simulate hydrate formation in and around a thin, dipping, planar sand stratum surrounded by clay lithology as it is buried to 295mbsf. We find that with sufficient methane beingmore » supplied by organic methanogenesis in the clays, a 200x pore size contrast between clays and sands allows for a strong enough concentration gradient to significantly drop the concentration of methane hydrate in clays immediately surrounding a thin sand layer, a phenomenon that is observed in well log data. Building upon previous work, our simulations account for the increase in sand-clay solubility contrast with depth from about 1.6% near the top of the sediment column to 8.6% at depth, which leads to a progressive strengthening of the diffusive flux of methane with time. By including an exponentially decaying organic methanogenesis input to the clay lithology with depth, we see a decrease in the aqueous methane supplied to the clays surrounding the sand layer with time, which works to further enhance the contrast in hydrate saturation between the sand and surrounding clays. Significant diffusive methane transport is observed in a clay interval of about 11m above the sand layer and about 4m below it, which matches well log observations. The clay-sand pore size contrast alone is not enough to completely eliminate hydrate (as observed in logs), because the diffusive flux of aqueous methane due to a contrast in pore size occurs slower than the rate at which methane is supplied via organic methanogenesis. Therefore, it is likely that additional mechanisms are at play, notably bound water activity reduction in clays. Three-dimensionality allows for inclusion of lithologic heterogeneities, which focus fluid flow and subsequently allow for heterogeneity in the methane migration mechanisms that dominate in marine sediments at a local scale. Incorporating recently acquired 3D seismic data from Walker Ridge to inform the lithologic structure of our modeled reservoir, we show that even with deep adjective sourcing of methane along highly permeable pathways, local hydrate accumulations can be sourced either by diffusive or advective methane flux; advectively-sourced hydrates accumulate evenly in highly permeable strata, while diffusively-sourced hydrates are characterized by thin strata-bound intervals with high clay-sand pore size contrasts.« less

Prediction of the filtrate particle size distribution from the pore size distribution in membrane filtration: Numerical correlations from computer simulations

NASA Astrophysics Data System (ADS)

Marrufo-Hernández, Norma Alejandra; Hernández-Guerrero, Maribel; Nápoles-Duarte, José Manuel; Palomares-Báez, Juan Pedro; Chávez-Rojo, Marco Antonio

2018-03-01

We present a computational model that describes the diffusion of a hard spheres colloidal fluid through a membrane. The membrane matrix is modeled as a series of flat parallel planes with circular pores of different sizes and random spatial distribution. This model was employed to determine how the size distribution of the colloidal filtrate depends on the size distributions of both, the particles in the feed and the pores of the membrane, as well as to describe the filtration kinetics. A Brownian dynamics simulation study considering normal distributions was developed in order to determine empirical correlations between the parameters that characterize these distributions. The model can also be extended to other distributions such as log-normal. This study could, therefore, facilitate the selection of membranes for industrial or scientific filtration processes once the size distribution of the feed is known and the expected characteristics in the filtrate have been defined.

Pore pressure development beneath the décollement at the Nankai subduction zone: Implications for plate boundary fault strength and sediment dewatering

NASA Astrophysics Data System (ADS)

Skarbek, Robert M.; Saffer, Demian M.

2009-07-01

Despite its importance for plate boundary fault processes, quantitative constraints on pore pressure are rare, especially within fault zones. Here, we combine laboratory permeability measurements from core samples with a model of loading and pore pressure diffusion to investigate pore fluid pressure evolution within underthrust sediment at the Nankai subduction zone. Independent estimates of pore pressure to ˜20 km from the trench, combined with permeability measurements conducted over a wide range of effective stresses and porosities, allow us to reliably simulate pore pressure development to greater depths than in previous studies and to directly quantify pore pressure within the plate boundary fault zone itself, which acts as the upper boundary of the underthrusting section. Our results suggest that the time-averaged excess pore pressure (P*) along the décollement ranges from 1.7-2.1 MPa at the trench to 30.2-35.9 MPa by 40 km landward, corresponding to pore pressure ratios of λb = 0.68-0.77. For friction coefficients of 0.30-0.40, the resulting shear strength along the décollement remains <12 MPa over this region. When noncohesive critical taper theory is applied using these values, the required pore pressure ratios within the wedge are near hydrostatic (λw = 0.41-0.59), implying either that pore pressure throughout the wedge is low or that the fault slips only during transient pulses of elevated pore pressure. In addition, simulated downward migration of minima in effective stress during drainage provides a quantitative explanation for down stepping of the décollement that is consistent with observations at Nankai.

Investigating axial diffusion in cylindrical pores using confocal single-particle fluorescence correlation spectroscopy.

PubMed

Chen, Fang; Neupane, Bhanu; Li, Peiyuan; Su, Wei; Wang, Gufeng

2016-08-01

We explored the feasibility of using confocal fluorescence correlation spectroscopy to study small nanoparticle diffusion in hundred-nanometer-sized cylindrical pores. By modeling single particle diffusion in tube-like confined three-dimensional space aligned parallel to the confocal optical axis, we showed that two diffusion dynamics can be observed in both original intensity traces and the autocorrelation functions (ACFs): the confined two-dimensional lateral diffusion and the unconfined one-dimensional (1D) axial diffusion. The separation of the axial and confined lateral diffusion dynamics provides an opportunity to study diffusions in different dimensions separately. We further experimentally studied 45 nm carboxylated polystyrene particles diffusing in 300 nm alumina pores. The experimental data showed consistency with the simulation. To extract the accurate axial diffusion coefficient, we found that a 1D diffusion model with a Lorentzian axial collection profile needs to be used to analyze the experimental ACFs. The diffusion of the 45 nm nanoparticles in polyethyleneglycol-passivated 300 nm pores slowed down by a factor of ∼2, which can be satisfactorily explained by hydrodynamic frictions. © 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Spatially resolved D-T(2) correlation NMR of porous media.

PubMed

Zhang, Yan; Blümich, Bernhard

2014-05-01

Within the past decade, 2D Laplace nuclear magnetic resonance (NMR) has been developed to analyze pore geometry and diffusion of fluids in porous media on the micrometer scale. Many objects like rocks and concrete are heterogeneous on the macroscopic scale, and an integral analysis of microscopic properties provides volume-averaged information. Magnetic resonance imaging (MRI) resolves this spatial average on the contrast scale set by the particular MRI technique. Desirable contrast parameters for studies of fluid transport in porous media derive from the pore-size distribution and the pore connectivity. These microscopic parameters are accessed by 1D and 2D Laplace NMR techniques. It is therefore desirable to combine MRI and 2D Laplace NMR to image functional information on fluid transport in porous media. Because 2D Laplace resolved MRI demands excessive measuring time, this study investigates the possibility to restrict the 2D Laplace analysis to the sum signals from low-resolution pixels, which correspond to pixels of similar amplitude in high-resolution images. In this exploratory study spatially resolved D-T2 correlation maps from glass beads and mortar are analyzed. Regions of similar contrast are first identified in high-resolution images to locate corresponding pixels in low-resolution images generated with D-T2 resolved MRI for subsequent pixel summation to improve the signal-to-noise ratio of contrast-specific D-T2 maps. This method is expected to contribute valuable information on correlated sample heterogeneity from the macroscopic and the microscopic scales in various types of porous materials including building materials and rock. Copyright © 2014 Elsevier Inc. All rights reserved.

Size Dependent Pore Formation in Germanium Nanowires Undergoing Reversible Delithiation Observed by In Situ TEM

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lu, Xiaotang; He, Yang; Mao, Scott X.

Germanium (Ge) nanowires coated with an amorphous silicon (Si) shell undergoing lithiation and delithiation were studied using in situ transmission electron microscopy (TEM). Delithiation creates pores in nanowires with diameters larger than ~25 nm, but not in smaller diameter nanowires. The formation of pores in Ge nanowires undergoing delithiation has been observed before in in situ TEM experiments, but there has been no indication that a critical diameter exists below which pores do not form. Pore formation occurs as a result of fast lithium diffusion compared to vacancy migration. We propose that a short diffusion path for vacancies to themore » nanowire surface plays a role in limiting pore formation even when lithium diffusion is fast.« less

A 2D Micromodel Study of Fines Migration and Clogging Behavior in Porous Media: Implications of Fines on Methane Extraction from Hydrate-Bearing Sediments

NASA Astrophysics Data System (ADS)

Cao, S. C.; Jang, J.; Waite, W. F.; Jafari, M.; Jung, J.

2017-12-01

Fine-grained sediment, or "fines," exist nearly ubiquitously in natural sediment, even in the predominantly coarse-grained sediments that host gas hydrates. Fines within these sandy sediments can play a crucial role during gas hydrate production activities. During methane extraction, several processes can alter the mobility and clogging potential of fines: 1) fluid flow as the formation is depressurized to release methane from hydrate; 2) pore-fluid chemistry shifts as pore-fluid brine freshens due to pure water released from dissociating hydrate; 3) the presence of a moving gas/water interface as gas evolves from dissociating hydrate and moves through the reservoir toward the production well. To evaluate fines migration and clogging behavior changes resulting from methane gas production and pore-water freshening during hydrate dissociation, 2D micromodel experiments have been conducted on a selection of pure fines, pore-fluids, and micromodel pore-throat sizes. Additionally, tests have been run with and without an invading gas phase (CO2) to test the significance of a moving meniscus on fines mobility and clogging. The endmember fine particles chosen for this research include silica silt, mica, calcium carbonate, diatoms, kaolinite, illite, and bentonite (primarily made of montmorillonite). The pore fluids include deionized water, sodium chloride brine (2M concentration), and kerosene. The microfluidic pore models, used as porous media analogs, were fabricated with pore-throat widths of 40, 60, and 100 µm. Results from this research show that in addition to the expected dependence of clogging on the ratio of particle-to-pore-throat size, pore-fluid chemistry is also a significant factor because the interaction between a particular type of fine and pore fluid influences that fine's capacity to cluster, clump together and effectively increase its particle "size" relative to the pore-throat width. The presence of a moving gas/fluid meniscus increases the clogging potential regardless of fine type as the advancing meniscus tends to gather and concentrate the fines. Results show the need to identify both the type and concentration of fines prior to evaluating whether a system's clogging potential will increase or decrease as pore waters freshen during methane extraction from hydrate.

Capillary filling rules and displacement mechanisms for spontaneous imbibition of CO2 for carbon storage and EOR using micro-model experiments and pore scale simulation

NASA Astrophysics Data System (ADS)

Chapman, E.; Yang, J.; Crawshaw, J.; Boek, E. S.

2012-04-01

In the 1980s, Lenormand et al. carried out their pioneering work on displacement mechanisms of fluids in etched networks [1]. Here we further examine displacement mechanisms in relation to capillary filling rules for spontaneous imbibition. Understanding the role of spontaneous imbibition in fluid displacement is essential for refining pore network models. Generally, pore network models use simple capillary filling rules and here we examine the validity of these rules for spontaneous imbibition. Improvement of pore network models is vital for the process of 'up-scaling' to the field scale for both enhanced oil recovery (EOR) and carbon sequestration. In this work, we present our experimental microfluidic research into the displacement of both supercritical CO2/deionised water (DI) systems and analogous n-decane/air - where supercritical CO2 and n-decane are the respective wetting fluids - controlled by imbibition at the pore scale. We conducted our experiments in etched PMMA and silicon/glass micro-fluidic hydrophobic chips. We first investigate displacement in single etched pore junctions, followed by displacement in complex network designs representing actual rock thin sections, i.e. Berea sandstone and Sucrosic dolomite. The n-decane/air experiments were conducted under ambient conditions, whereas the supercritical CO2/DI water experiments were conducted under high temperature and pressure in order to replicate reservoir conditions. Fluid displacement in all experiments was captured via a high speed video microscope. The direction and type of displacement the imbibing fluid takes when it enters a junction is dependent on the number of possible channels in which the wetting fluid can imbibe, i.e. I1, I2 and I3 [1]. Depending on the experiment conducted, the micro-models were initially filled with either DI water or air before the wetting fluid was injected. We found that the imbibition of the wetting fluid through a single pore is primarily controlled by the geometry of the pore body rather than the downstream pore throat sizes, contrary to the established capillary filling rules as used in current pore network models. Our experimental observations are confirmed by detailed lattice-Boltzmann pore scale computer simulations of fluid displacement in the same geometries. This suggests that capillary filling rules for imbibition as used in pore network models may need to be revised. [1] G. Lenormand, C. Zarcone and A. Sarr, J. Fluid Mech. 135 , 337-353 (1983).

Pore Pressure Diffusion as a possible mechanism for the Ag. Ioanis 2001 earthquake swarm activity (Gulf of Corinth, Central Greece).

NASA Astrophysics Data System (ADS)

Vallianatos, F.; Michas, G.; Papadakis, G.; Sammonds, P.

2012-04-01

The Gulf of Corinth rift (Central Greece) is one of the most seismotectonically active areas in Europe (Ambraseys and Jackson, 1990; 1997), with an important continental N-S extension of about 13 mm/yr and 6 mm/yr at the west and east part respectively (Clarke et al., 1997a). The seismicity of the area includes 5 main earthquakes of magnitude greater than 5.8 since 1960. In the western part of the rift, where the extension reaches its maximum value, earthquake swarms are often being observed (Bourouis and Cornet, 2009). Such an earthquake crisis has been occurred on 2001 at the southern margin of the west part of the rift. The crisis lasted about 100 days with a major event the Ag. Ioanis earthquake (4.3 Mw) on 8th of April 2001 (Pacchiani and Lyon-Caen, 2010). The possible relation between fluids flow and the observed earthquake swarms at the west part of the Gulf of Corinth rift has been discussed in the works of Bourouis and Cornet (2009) and Pacchiani and Lyon-Caen (2010). In the present work we examine the spatiotemporal properties of the Ag. Ioanis 2001 earthquake swarm, using data from the CRL network (http://crlab.eu/). We connect these properties to a mechanism due to pore pressure diffusion (Shapiro et al., 1997) and we estimate the hydraulic diffusivity and the permeability of the surrounding rocks. A back front of the seismicity (Parotidis et al., 2004) is also been observed, related to the migration of seismicity and the development of a quiescence region near the area of the initial pore pressure perturbation. Moreover, anisotropy of the hydraulic diffusivity has been observed, revealing the heterogeneity of the surrounding rocks and the fracture systems. This anisotropy is consistent in direction with the fault zone responsible for the Ag. Ioanis earthquake (Pacchiani and Lyon-Caen, 2010). Our results indicate that fluids flow and pore pressure perturbations are possible mechanisms for the initiation and the evolution of the Ag. Ioanis 2001 earthquake swarm activity and reveal the possible connection of the complex fracture network to the spatial evolution of seismicity on an active tectonic region as is the Gulf of Corinth rift. Acknowledgments. This work was supported in part by the THALES Program of the Ministry of Education of Greece and the European Union in the framework of the project entitled "Integrated understanding of Seismicity, using innovative Methodologies of Fracture mechanics along with Earthquake and non-extensive statistical physics - Application to the geodynamic system of the Hellenic Arc. SEISMO FEAR HELLARC". GM and GP wish to acknowledge the partial support of the Greek State Scholarships Foundation (IKY).

Adsorption behaviors of supercritical Lennard-Jones fluid in slit-like pores.

PubMed

Li, Yingfeng; Cui, Mengqi; Peng, Bo; Qin, Mingde

2018-05-18

Understanding the adsorption behaviors of supercritical fluid in confined space is pivotal for coupling the supercritical technology and the membrane separation technology. Based on grand canonical Monte Carlo simulations, the adsorption behaviors of a Lennard-Jones (LJ) fluid in slit-like pores at reduced temperatures over the critical temperature, T c *  = 1.312, are investigated; and impacts of the wall-fluid interactions, the pore width, and the temperature are taken into account. It is found that even if under supercritical conditions, the LJ fluid can undergo a "vapor-liquid phase transition" in confined space, i.e., the adsorption density undergoes a sudden increase with the bulk density. A greater wall-fluid attractive potential, a smaller pore width, and a lower temperature will bring about a stronger confinement effect. Besides, the adsorption pressure reaches a local minimum when the bulk density equals to a certain value, independent of the wall-fluid potential or pore width. The insights in this work have both practical and theoretical significances. Copyright © 2018 Elsevier Inc. All rights reserved.

Transport Mechanism of Guest Methane in Water-Filled Nanopores

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bui, Tai; Phan, Anh; Cole, David R.

We computed the transport of methane through 1 nm wide slit-shaped pores carved out of selected solid substrates using classical molecular dynamics simulations. The transport mechanism was elucidated via the implementation of the well-tempered metadynamics algorithm, which allowed for the quantification and visualization of the free energy landscape sampled by the guest molecule. Models for silica, magnesium oxide, alumina, muscovite, and calcite were used as solid substrates. Slit-shaped pores of width 1 nm were carved out of these materials and filled with liquid water. Methane was then inserted at low concentration. The results show that the diffusion of methane throughmore » the hydrated pores is strongly dependent on the solid substrate. While methane molecules diffuse isotropically along the directions parallel to the pore surfaces in most of the pores considered, anisotropic diffusion was observed in the hydrated calcite pore. The differences observed in the various pores are due to local molecular properties of confined water, including molecular structure and solvation free energy. The transport mechanism and the diffusion coefficients are dependent on the free energy barriers encountered by one methane molecule as it migrates from one preferential adsorption site to a neighboring one. It was found that the heterogeneous water distribution in different hydration layers and the low free energy pathways in the plane parallel to the pore surfaces yield the anisotropic diffusion of methane molecules in the hydrated calcite pore. Our observations contribute to an ongoing debate on the relation between local free energy profiles and diffusion coefficients and could have important practical consequences in various applications, ranging from the design of selective membranes for gas separations to the sustainable deployment of shale gas.« less

Transport Mechanism of Guest Methane in Water-Filled Nanopores

DOE PAGES

Bui, Tai; Phan, Anh; Cole, David R.; ...

2017-05-11

We computed the transport of methane through 1 nm wide slit-shaped pores carved out of selected solid substrates using classical molecular dynamics simulations. The transport mechanism was elucidated via the implementation of the well-tempered metadynamics algorithm, which allowed for the quantification and visualization of the free energy landscape sampled by the guest molecule. Models for silica, magnesium oxide, alumina, muscovite, and calcite were used as solid substrates. Slit-shaped pores of width 1 nm were carved out of these materials and filled with liquid water. Methane was then inserted at low concentration. The results show that the diffusion of methane throughmore » the hydrated pores is strongly dependent on the solid substrate. While methane molecules diffuse isotropically along the directions parallel to the pore surfaces in most of the pores considered, anisotropic diffusion was observed in the hydrated calcite pore. The differences observed in the various pores are due to local molecular properties of confined water, including molecular structure and solvation free energy. The transport mechanism and the diffusion coefficients are dependent on the free energy barriers encountered by one methane molecule as it migrates from one preferential adsorption site to a neighboring one. It was found that the heterogeneous water distribution in different hydration layers and the low free energy pathways in the plane parallel to the pore surfaces yield the anisotropic diffusion of methane molecules in the hydrated calcite pore. Our observations contribute to an ongoing debate on the relation between local free energy profiles and diffusion coefficients and could have important practical consequences in various applications, ranging from the design of selective membranes for gas separations to the sustainable deployment of shale gas.« less

Mouthpart conduit sizes of fluid-feeding insects determine the ability to feed from pores

DOE PAGES

Lehnert, Matthew S.; Bennett, Andrew; Reiter, Kristen E.; ...

2017-01-04

Fluid-feeding insects, such as butterflies, moths, and flies (20% of all animal species), are faced with the common selection pressure of having to remove and feed on trace amounts of fluids from porous surfaces. Insects able to acquire fluids that are confined to pores during drought conditions would have an adaptive advantage and increased fitness over other individuals. Here we performed feeding trials using solutions with magnetic nanoparticles to show that butterflies and flies have mouthparts adapted to pull liquids from porous surfaces using capillary action as the governing principle. In addition, the ability to feed on the liquids collectedmore » from pores depends on a relationship between the diameter of the mouthpart conduits and substrate pore size diameter; insects with mouthpart conduit diameters larger than the pores cannot successfully feed, thus there is a limiting substrate pore size from which each species can acquire liquids for fluid uptake. In conclusion, given that natural selection independently favored mouthpart architectures that support these methods of fluid uptake (Diptera and Lepidoptera share a common ancestor 280 mya that had chewing mouthparts), we suggest that the convergence of this mechanism advocates this as an optimal strategy for pulling trace amounts of fluids from porous surfaces.« less

Mouthpart conduit sizes of fluid-feeding insects determine the ability to feed from pores

DOE Office of Scientific and Technical Information (OSTI.GOV)

Lehnert, Matthew S.; Bennett, Andrew; Reiter, Kristen E.

Fluid-feeding insects, such as butterflies, moths, and flies (20% of all animal species), are faced with the common selection pressure of having to remove and feed on trace amounts of fluids from porous surfaces. Insects able to acquire fluids that are confined to pores during drought conditions would have an adaptive advantage and increased fitness over other individuals. Here we performed feeding trials using solutions with magnetic nanoparticles to show that butterflies and flies have mouthparts adapted to pull liquids from porous surfaces using capillary action as the governing principle. In addition, the ability to feed on the liquids collectedmore » from pores depends on a relationship between the diameter of the mouthpart conduits and substrate pore size diameter; insects with mouthpart conduit diameters larger than the pores cannot successfully feed, thus there is a limiting substrate pore size from which each species can acquire liquids for fluid uptake. In conclusion, given that natural selection independently favored mouthpart architectures that support these methods of fluid uptake (Diptera and Lepidoptera share a common ancestor 280 mya that had chewing mouthparts), we suggest that the convergence of this mechanism advocates this as an optimal strategy for pulling trace amounts of fluids from porous surfaces.« less

Surface sealing using self-assembled monolayers and its effect on metal diffusion in porous low-k dielectrics studied using monoenergetic positron beams

NASA Astrophysics Data System (ADS)

Uedono, Akira; Armini, Silvia; Zhang, Yu; Kakizaki, Takeaki; Krause-Rehberg, Reinhard; Anwand, Wolfgang; Wagner, Andreas

2016-04-01

Surface sealing effects on the diffusion of metal atoms in porous organosilicate glass (OSG) films were studied by monoenergetic positron beams. For a Cu(5 nm)/MnN(3 nm)/OSG(130 nm) sample fabricated with pore stuffing, C4F8 plasma etch, unstuffing, and a self-assembled monolayer (SAM) sealing process, it was found that pores with cubic pore side lengths of 1.1 and 3.1 nm coexisted in the OSG film. For the sample without the SAM sealing process, metal (Cu and Mn) atoms diffused from the top Cu/MnN layer into the OSG film and were trapped by the pores. As a result, almost all pore interiors were covered with those metals. For the sample damaged by an Ar/C4F8 plasma etch treatment before the SAM sealing process, SAMs diffused into the OSG film, and they were preferentially trapped by larger pores. The cubic pore side length in these pores containing self-assembled molecules was estimated to be 0.7 nm. Through this work, we have demonstrated that monoenergetic positron beams are a powerful tool for characterizing capped porous films and the trapping of atoms and molecules by pores.

Fractal Viscous Fingering in Fracture Networks

NASA Astrophysics Data System (ADS)

Boyle, E.; Sams, W.; Ferer, M.; Smith, D. H.

2007-12-01

We have used two very different physical models and computer codes to study miscible injection of a low- viscosity fluid into a simple fracture network, where it displaces a much-more viscous "defending" fluid through "rock" that is otherwise impermeable. The one code (NETfLow) is a standard pore level model, originally intended to treat laboratory-scale experiments; it assumes negligible mixing of the two fluids. The other code (NFFLOW) was written to treat reservoir-scale engineering problems; It explicitly treats the flow through the fractures and allows for significant mixing of the fluids at the interface. Both codes treat the fractures as parallel plates, of different effective apertures. Results are presented for the composition profiles from both codes. Independent of the degree of fluid-mixing, the profiles from both models have a functional form identical to that for fractal viscous fingering (i.e., diffusion limited aggregation, DLA). The two codes that solve the equations for different models gave similar results; together they suggest that the injection of a low-viscosity fluid into large- scale fracture networks may be much more significantly affected by fractal fingering than previously illustrated.

Transient poroelastic stress coupling between the 2015 M7.8 Gorkha, Nepal earthquake and its M7.3 aftershock

NASA Astrophysics Data System (ADS)

Tung, S.; Masterlark, T.; Dovovan, T.

2018-05-01

The large M7.3 aftershock occurred 17 days after the 2015 M7.8 Gorkha earthquake. We investigate if this sequence is mechanically favored by the mainshock via time-dependent fluid migration and pore pressure recovery. This study uses finite element models of fully-coupled poroelastic coseismic and postseismic behavior to simulate the evolving stress and pore-pressure fields. Using simulations of a reasonable permeability, the hypocenter was destabilized by an additional 0.15 MPa of Coulomb failure stress change (ΔCFS) and 0.17 MPa of pore pressure (Δp), the latter of which induced lateral and upward diffusive fluid flow (up to 2.76 mm/day) in the aftershock region. The M7.3 location is predicted next to a local maximum of Δp and a zone of positive ΔCFS northeast of Kathmandu. About 60% of the aftershocks occurred within zones having either Δp > 0 or ΔCFS > 0. Particularly in the eastern flank of the epicentral area, 83% of the aftershocks experienced postseismic fluid pressurization and 88% of them broke out with positive pore pressure, which are discernibly more than those with positive ΔCFS (71%). The transient scalar field of fluid pressurization provides a good proxy to predict aftershock-prone areas in space and time, because it does not require extraction of an assumed vector field from transient stress tensor fields as is the case for ΔCFS calculations. A bulk permeability of 8.32 × 10-18 m2 is resolved to match the transient response and the timing of the M7.3 rupture which occurred at the peak of the ΔCFS time-series. This estimate is consistent with the existing power-law permeability-versus-depth models, suggesting an intermediately-fractured upper crust coherent with the local geology of the central Himalayas. The contribution of poroelastic triggering is verified against different poroelastic moduli and surface flow-pressure boundaries, suggesting that a poroelastic component is essential to account for the time interval separating the mainshock and the M7.3 aftershock.

Microstructures of Randall's plaques and their interfaces with calcium oxalate monohydrate kidney stones reflect underlying mineral precipitation mechanisms.

PubMed

Sethmann, Ingo; Wendt-Nordahl, Gunnar; Knoll, Thomas; Enzmann, Frieder; Simon, Ludwig; Kleebe, Hans-Joachim

2017-06-01

Randall's plaques (RP) are preferred sites for the formation of calcium oxalate monohydrate (COM) kidney stones. However, although processes of interstitial calcium phosphate (CaP) plaque formation are not well understood, the potential of plaque microstructures as indicators of CaP precipitation conditions received only limited attention. We investigated RP-associated COM stones for structural details of the calcified tissues and microstructural features of plaque-stone interfaces as indicators of the initial processes of stone formation. Significantly increased CaP supersaturation can be expected for interstitial fluid, if reabsorbed ions from the tubular system continuously diffuse into the collagenous connective tissue. Densely packed, fine-grained CaP particles were found in dense textures of basement membranes while larger, laminated particles were scattered in coarse-meshed interstitial tissue, which we propose to be due to differential spatial confinements and restrictions of ion diffusion. Particle morphologies suggest an initial precipitation as metastable amorphous calcium phosphate (ACP). Morphologies and arrangements of first COM crystals at the RP-stone interface ranged from stacked euhedral platelets to skeletal morphologies and even porous, dendritic structures, indicating, in this order, increasing levels of COM supersaturation. Furthermore, these first COM crystals were often coated with CaP. On this basis, we propose that ions from CaP-supersaturated interstitial fluid may diffuse through porous RP into the urine, where a resulting local increase in COM supersaturation could trigger crystal nucleation and, hence, initiate stone formation. Ion-depleted fluid in persistent pores of initial COM layers may get replenished from interstitial fluid, leading to CaP precipitation in porous COM.

Pore Pressure and Stress Distributions Around a Hydraulic Fracture in Heterogeneous Rock

NASA Astrophysics Data System (ADS)

Gao, Qian; Ghassemi, Ahmad

2017-12-01

One of the most significant characteristics of unconventional petroleum bearing formations is their heterogeneity, which affects the stress distribution, hydraulic fracture propagation and also fluid flow. This study focuses on the stress and pore pressure redistributions during hydraulic stimulation in a heterogeneous poroelastic rock. Lognormal random distributions of Young's modulus and permeability are generated to simulate the heterogeneous distributions of material properties. A 3D fully coupled poroelastic model based on the finite element method is presented utilizing a displacement-pressure formulation. In order to verify the model, numerical results are compared with analytical solutions showing excellent agreements. The effects of heterogeneities on stress and pore pressure distributions around a penny-shaped fracture in poroelastic rock are then analyzed. Results indicate that the stress and pore pressure distributions are more complex in a heterogeneous reservoir than in a homogeneous one. The spatial extent of stress reorientation during hydraulic stimulations is a function of time and is continuously changing due to the diffusion of pore pressure in the heterogeneous system. In contrast to the stress distributions in homogeneous media, irregular distributions of stresses and pore pressure are observed. Due to the change of material properties, shear stresses and nonuniform deformations are generated. The induced shear stresses in heterogeneous rock cause the initial horizontal principal stresses to rotate out of horizontal planes.

Gint2D-T2 correlation NMR of porous media

NASA Astrophysics Data System (ADS)

Zhang, Yan; Blümich, Bernhard

2015-03-01

The internal magnetic field gradient induced in porous media by magnetic susceptibility differences at material interfaces impacts diffusion measurements in particular at high magnetic field and can be used to probe the pore structure. Insight about the relationship between pore space and internal gradient Gint can be obtained from 2D Laplace NMR experiments. When measuring distributions of transverse relaxation times T2 in fluid filled porous media, relaxation and diffusion in internal gradients arise simultaneously and data are often interpreted with the assumption that one or the other parameter be constant throughout the sample. To examine this assumption we measure correlations of the distributions of Gint2D and T2 by 2D Laplace NMR for three different kinds of samples, glass beads with different bead diameters saturated with water, glass beads filled with oil and water, and a wet mortar sample. For the first two samples the cases where either the internal gradient or diffusion dominates were examined separately in order to better understand the relationship between Gint and D. These results are useful for assessing the impact of internal gradients and diffusion in unknown samples, such as the mortar sample. The experiments were performed at different magnetic field strengths corresponding to 300 MHz and 700 MHz 1H Larmor frequency to identify the impact of the magnetic field on the internal gradient. Subsequently, spatially resolved Gint2D-T2 maps were obtained to study the sample heterogeneity.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.