Mechanics of dual-mode dilative failure in subaqueous sediment deposits

NASA Astrophysics Data System (ADS)

You, Yao; Flemings, Peter; Mohrig, David

2014-07-01

We introduce dual-mode dilative failure with flume experiments. Dual-mode dilative failure combines slow and steady release of sediments by breaching with periodic sliding, which rapidly releases an internally coherent wedge of sediments. It occurs in dilative sandy deposits. This periodic slope failure results from cyclic evolution of the excess pore pressure in the deposit. Sliding generates large, transient, negative excess pore pressure that strengthens the deposit and allows breaching to occur. During breaching, negative excess pore pressure dissipates, the deposit weakens, and ultimately sliding occurs once again. We show that the sliding frequency is proportional to the coefficient of consolidation. We find that thicker deposits are more susceptible to dual-mode dilative failure. Discovery of dual-mode dilative failure provides a new mechanism to consider when interpreting the sedimentary deposits linked to submarine slope failures.

Full paleostress tensor reconstruction using quartz veins of Panasqueira Mine, central Portugal; part I: Paleopressure determination

NASA Astrophysics Data System (ADS)

Jaques, Luís; Pascal, Christophe

2017-09-01

Paleostress tensor restoration methods are traditionally limited to reconstructing geometrical parameters and are unable to resolve stress magnitudes. Based on previous studies we further developed a methodology to restore full paleostress tensors. We concentrated on inversion of Mode I fractures and acquired data in Panasqueira Mine, Portugal, where optimal exposures of mineralized quartz veins can be found. To carry out full paleostress restoration we needed to determine (1) pore (paleo)pressure and (2) vein attitudes. The present contribution focuses specifically on the determination of pore pressure. To these aims we conducted an extensive fluid inclusion study to derive fluid isochores from the quartz of the studied veins. To constrain P-T conditions, we combined these isochores with crystallisation temperatures derived from geochemical analyses of coeval arsenopyrite. We also applied the sphalerite geobarometer and considered two other independent pressure indicators. Our results point to pore pressures of ∼300 MPa and formation depths of ∼10 km. Such formation depths are in good agreement with the regional geological evolution. The obtained pore pressure will be merged with vein inversion results, in order to achieve full paleostress tensor restoration, in a forthcoming companion paper.

The effect of stress on limestone permeability and its effective stress behavior

NASA Astrophysics Data System (ADS)

Meng, F.; Baud, P.; Ge, H.; Wong, T. F.

2017-12-01

The evolution of permeability and its effective stress behavior is related to inelastic deformation and failure mode. This was investigated in Indiana and Purbeck limestones with porosities of 18% and 13%, respectively. Hydrostatic and triaxial compression tests were conducted at room temperature on water-saturated samples at pore pressure of 5 MPa and confining pressures up to 90 MPa. Permeability was measured using steady flow at different stages of deformation. For Indiana limestone, under hydrostatic loading pore collapse initiated at critical pressure P* 55 MPa with an accelerated reduction of permeability by 1/2. At a confinement of 35 MPa and above, shear-enhanced compaction initiated at critical stress C*, beyond which permeability reduction up to a factor of 3 was observed. At a confinement of 15 MPa and below, dilatancy initiated at critical stress C', beyond which permeability continued to decrease, with a negative correlation between porosity and permeability changes. Purbeck limestone showed similar evolution of permeability. Microstructural and mercury porosimetry data showed that pore size distribution in both Indiana and Purbeck limestones is bimodal, with significant proportions of macropores and micropores. The effective stress behaviour of a limestone with dual porosity is different from the prediction for a microscopically homogeneous assemblage, in that its effective stress coefficients for permeability and porosity change may attain values significantly >1. Indeed this was confirmed by our measurements (at confining pressures of 7-15 MPa and pore pressures of 1-3 MPa) in samples that had not been deformed inelastically. We also investigated the behavior in samples hydrostatically and triaxially compacted to beyond the critical stresses P* and C*, respectively. Experimental data for these samples consistently showed effective stress coefficients for both permeability and porosity change with values <1. Thus the effective stress behavior in an inelastically compacted sample is fundamentally different, with attributes akin to that of a microscopically homogeneous assemblage. This is likely related to compaction from pervasive collapse of macropores, which would effectively homogenize the initially bimodal pore size distribution.

Structure, Mechanics and Flow Properties of Fractured Shale: Core-Scale Experimentation and In-situ Imaging

NASA Astrophysics Data System (ADS)

Abdelmalek, B. F.; Karpyn, Z.; Liu, S.

2014-12-01

Over the last several years, hydrocarbon exploitation and development in North America has been heavily centered on shale gas plays. However, the physical attributes of shales and their manifestation on transport properties and storage capacity remain poorly understood. Therefore, more experimentally based data are needed to fill the gaps in understanding both transport and storage of fluids in shale. The proposed work includes installation and testing of an experimental system which is capable of monitoring the dynamic evolution of shale core permeability under variable loading conditions and in coordination with X-ray microCT imaging. The goal of this study is to better understand and quantify fluid flow patterns and associated transport dynamics of fractured shale samples. The independent variables considered in this study are: mechanical loading and pore pressure. The mechanical response of shale core is captured for different loading paths. To best replicate the in-situ production scenario, the pore pressure is progressively depleted to mimic pressure decline. During the course of experimentation, permeability is estimated using the pulse-decay method under tri-axial stress boundary conditions. Simultaneously, X-ray microCT imaging is used with a tracer gas that is allowed to flow through the sample as an illuminating agent. In the presence of an illuminating agent, either Xenon or Krypton, the X-ray CT scanner can image fractures, global pathways and diffusional fronts in the matrix, as well as sorption sites that reflect heterogeneities in the sample and localized deformation. Anticipated results from these experiments will help quantify permeability evolution as a function of different loading conditions and pore pressure depletion. Also, the X-ray images will help visualize the change of flow patterns and the intensity of sorption as a function of mechanical loading and pore pressure.

The effect of the pore-fluid factor on strength and failure mechanism of Wilkeson sandstone

NASA Astrophysics Data System (ADS)

Kätker, A. K.; Rempe, M.; Renner, J.

2016-12-01

The effective stress law, σn,eff = σn - αpf, is a central tool in analysing phenomena related to hydromechanical coupling, such as fluid-induced seismicity or aftershock activity. The effective-stress coefficient α assumes different values for specific physical properties and may deviate from 1. The limited number of studies suggest that brittle compressive strength obeys an effective-stress law when effective drainage is achieved. Yet, open questions remain regarding, e.g., the role of the loading path. We performed suites of triaxial compression tests on samples of Wilkeson sandstone at a range of pore-fluid pressures but identical effective confining pressure (60, 100, and 120 MPa) maintaining the pore-fluid factor λ = pf / pc constant (0.05, 0.2, 0.4, 0.55) during the isostatic loading stage to ensure uniform loading paths. Samples were shortened with a strain rate of 4×10-7 s-1 yielding drained conditions. All tests were terminated at a total axial strain of 4.5% for comparability of microstructures. The tests also included continuous permeability determination and ultrasonic p-wave-velocity measurements to monitor microstructural evolution. Results from experiments conducted at peff = 100 MPa show that dry samples exhibit a higher peak strength and brittle failure while water-saturated samples tend to deform at lower stress by cataclastic flow indicating physico-chemical weakening. Regardless of pore-fluid factor, the saturated experiments exhibit similar peak and residual strength. Differences in failure mechanism (degree of macroscopic localization) and volumetric strain evolution are however noticed, albeit without systematic relation to pore-fluid factor. Microstructure analyses by optical and scanning electron microscopy revealed an evolution from localized shear zones in dry experiments and experiments with a low pore-fluid factor to rather distributed cataclastic flow for experiments with high pore fluid factors. Yet, mechanical and structural differences observed so far may result from sample-to-sample variability and the proximity of the experimental conditions to the brittle-ductile transition.

In situ pore-pressure evolution during dynamic CPT measurements in soft sediments of the western Baltic Sea

NASA Astrophysics Data System (ADS)

Seifert, Annedore; Stegmann, Sylvia; Mörz, Tobias; Lange, Matthias; Wever, Thomas; Kopf, Achim

2008-08-01

We present in situ strength and pore-pressure measurements from 57 dynamic cone penetration tests in sediments of Mecklenburg ( n = 51), Eckernförde ( n = 2) and Gelting ( n = 4) bays, western Baltic Sea, characterised by thick mud layers and partially free microbial gas resulting from the degradation of organic material. In Mecklenburg and Eckernförde bays, sediment sampling by nine gravity cores served sedimentological characterisation, analyses of geotechnical properties, and laboratory shear tests. At selected localities, high-resolution echo-sounder profiles were acquired. Our aim was to deploy a dynamic cone penetrometer (CPT) to infer sediment shear strength and cohesion of the sea bottom as a function of fluid saturation. The results show very variable changes in pore pressure and sediment strength during the CPT deployments. The majority of the CPT measurements ( n = 54) show initially negative pore-pressure values during penetration, and a delayed response towards positive pressures thereafter. This so-called type B pore-pressure signal was recorded in all three bays, and is typically found in soft muds with high water contents and undrained shear strengths of 1.6-6.4 kPa. The type B signal is further affected by displacement of sediment and fluid upon penetration of the lance, skin effects during dynamic profiling, enhanced consolidation and strength of individual horizons, the presence of free gas, and a dilatory response of the sediment. In Mecklenburg Bay, the remaining small number of CPT measurements ( n = 3) show a well-defined peak in both pore pressure and cone resistance during penetration, i.e. an initial marked increase which is followed by exponential pore-pressure decay during dissipation. This so-called type A pore-pressure signal is associated with normally consolidated mud, with indurated clay layers showing significantly higher undrained shear strength (up to 19 kPa). In Eckernförde and Gelting bays pore-pressure response type B is exclusively found, while in Mecklenburg Bay types A and B were detected. Despite the striking similarities in incremental density increase and shear strength behaviour with depth, gas occurrence and subtle variations in the coarse-grained fraction cause distinct pore-pressure curves. Gaseous muds interbedded with silty and sandy layers are most common in the three bays, and the potential effect of free gas (i.e. undersaturated pore space) on in situ strength has to be explored further.

Semipermeability Evolution of Wakkanai Mudstones During Isotropic Compression

NASA Astrophysics Data System (ADS)

Takeda, M.; Manaka, M.

2015-12-01

Precise identification of major processes that influence groundwater flow system is of fundamental importance for the performance assessment of waste disposal in subsurface. In the characterization of groundwater flow system, gravity- and pressure-driven flows have been conventionally assumed as dominant processes. However, recent studies have suggested that argillites can act as semipermeable membranes and they can cause chemically driven flow, i.e., chemical osmosis, under salinity gradients, which may generate erratic pore pressures in argillaceous formations. In order to identify the possibility that chemical osmosis is involved in erratic pore pressure generations in argillaceous formations, it is essential to measure the semipermeability of formation media; however, in the measurements of semipermeability, little consideration has been given to the stresses that the formation media would have experienced in past geologic processes. This study investigates the influence of stress history on the semipermeability of an argillite by an experimental approach. A series of chemical osmosis experiments were performed on Wakkanai mudstones to measure the evolution of semipermeability during loading and unloading confining pressure cycles. The osmotic efficiency, which represents the semipermeability, was estimated at each confining pressure. The results show that the osmotic efficiency increases almost linearly with increasing confining pressure; however, the increased osmotic efficiency does not recover during unloading unless the confining pressure is almost relieved. The observed unrecoverable change in osmotic efficiency may have an important implication on the evaluation of chemical osmosis in argillaceous formations that have been exposed to large stresses in past geologic processes. If the osmotic efficiency increased by the past stress can remain unchanged to date, the osmotic efficiency should be measured at the past highest stress rather than the current in-situ stress. Otherwise, the effect of chemical osmosis on the pore pressure generation would be underestimated.

Toward multiscale modelings of grain-fluid systems

NASA Astrophysics Data System (ADS)

Chareyre, Bruno; Yuan, Chao; Montella, Eduard P.; Salager, Simon

2017-06-01

Computationally efficient methods have been developed for simulating partially saturated granular materials in the pendular regime. In contrast, one hardly avoid expensive direct resolutions of 2-phase fluid dynamics problem for mixed pendular-funicular situations or even saturated regimes. Following previous developments for single-phase flow, a pore-network approach of the coupling problems is described. The geometry and movements of phases and interfaces are described on the basis of a tetrahedrization of the pore space, introducing elementary objects such as bridge, meniscus, pore body and pore throat, together with local rules of evolution. As firmly established local rules are still missing on some aspects (entry capillary pressure and pore-scale pressure-saturation relations, forces on the grains, or kinetics of transfers in mixed situations) a multi-scale numerical framework is introduced, enhancing the pore-network approach with the help of direct simulations. Small subsets of a granular system are extracted, in which multiphase scenario are solved using the Lattice-Boltzman method (LBM). In turns, a global problem is assembled and solved at the network scale, as illustrated by a simulated primary drainage.

Poromechanics of stick-slip frictional sliding and strength recovery on tectonic faults

DOE PAGES

Scuderi, Marco M.; Carpenter, Brett M.; Johnson, Paul A.; ...

2015-10-22

Pore fluids influence many aspects of tectonic faulting including frictional strength aseismic creep and effective stress during the seismic cycle. But, the role of pore fluid pressure during earthquake nucleation and dynamic rupture remains poorly understood. Here we report on the evolution of pore fluid pressure and porosity during laboratory stick-slip events as an analog for the seismic cycle. We sheared layers of simulated fault gouge consisting of glass beads in a double-direct shear configuration under true triaxial stresses using drained and undrained fluid conditions and effective normal stress of 5–10 MPa. Shear stress was applied via a constant displacementmore » rate, which we varied in velocity step tests from 0.1 to 30 µm/s. Here, we observe net pore pressure increases, or compaction, during dynamic failure and pore pressure decreases, or dilation, during the interseismic period, depending on fluid boundary conditions. In some cases, a brief period of dilation is attendant with the onset of dynamic stick slip. Our data show that time-dependent strengthening and dynamic stress drop increase with effective normal stress and vary with fluid conditions. For undrained conditions, dilation and preseismic slip are directly related to pore fluid depressurization; they increase with effective normal stress and recurrence time. Microstructural observations confirm the role of water-activated contact growth and shear-driven elastoplastic processes at grain junctions. These results indicate that physicochemical processes acting at grain junctions together with fluid pressure changes dictate stick-slip stress drop and interseismic creep rates and thus play a key role in earthquake nucleation and rupture propagation.« less

Deformation Failure Characteristics of Coal Body and Mining Induced Stress Evolution Law

PubMed Central

Wen, Zhijie; Wen, Jinhao; Shi, Yongkui; Jia, Chuanyang

2014-01-01

The results of the interaction between coal failure and mining pressure field evolution during mining are presented. Not only the mechanical model of stope and its relative structure division, but also the failure and behavior characteristic of coal body under different mining stages are built and demonstrated. Namely, the breaking arch and stress arch which influence the mining area are quantified calculated. A systematic method of stress field distribution is worked out. All this indicates that the pore distribution of coal body with different compressed volume has fractal character; it appears to be the linear relationship between propagation range of internal stress field and compressed volume of coal body and nonlinear relationship between the range of outburst coal mass and the number of pores which is influenced by mining pressure. The results provide theory reference for the research on the range of mining-induced stress and broken coal wall. PMID:24967438

Monitoring of Cyclic Steam Stimulation by Inversion of Surface Tilt Measurements

NASA Astrophysics Data System (ADS)

Maharramov, M.; Zoback, M. D.

2014-12-01

Temperature and pressure changes associated with the cyclic steam simulation (CSS) used in heavy oil production from sands are accompanied by significant deformation. Inversion of geomechanical data may provide a potentially powerful reservoir monitoring tool where geomechanical effects are significant. Induced pore pressure changes can be inverted from measurable surface deformations by solving an inverse problem of poroelasticity. In this work, we apply this approach to estimating pore pressure changes from surface tilt measurements at a heavy oil reservoir undergoing cyclic steam simulation. Steam was injected from November 2007 through January 2008. Surface tilt measurements were collected from 25 surface tilt stations during this period. The injection ran in two overlapping phases: Phase 1 ran from the beginning of the injection though mid-December, and Phase 2 overlapped with Phase 1 and ran through the beginning of January. During Phase 1 steam was injected in the western part of the reservoir, followed by injection in the eastern part in Phase 2. The pore pressure evolution was inverted from daily tilt measurements using regularized constrained least squares fitting, the results are shown on the plot. Estimated induced pore pressure change (color scale), observed daily incremental tilts (green arrows) and modeled daily incremental tilts (red arrows) are shown in three panels corresponding to two and five weeks of injection, and the end of injection period. DGPS measurements available for a single location were used as an additional inversion constraint. The results indicate that the pore pressure increase in the reservoir follows the same pattern as the steam injection, from west to east. This qualitative behaviour is independent of the amount of regularization, indirectly validating our inversion approach. Patches of lower pressure appear to be stable with regard to regularization and may provide valuable insight into the efficiency of steam injection. Inversion of pore pressure (and surface deformation) from tilts in this case is non-unique, and the DGPS measurement provided an important additional constraint. The method can be applied to inverting pore pressure changes from InSAR observations, and the latter can be expected to reduce limitations due to noise in tilt measurements.

One-dimensional pore pressure diffusion of different grain-fluid mixtures

NASA Astrophysics Data System (ADS)

von der Thannen, Magdalena; Kaitna, Roland

2015-04-01

During the release and the flow of fully saturated debris, non-hydrostatic fluid pressure can build up and probably dissipate during the event. This excess fluid pressure has a strong influence on the flow and deposition behaviour of debris flows. Therefore, we investigate the influence of mixture composition on the dissipation of non-hydrostatic fluid pressures. For this we use a cylindrical pipe of acrylic glass with installed pore water pressure sensors in different heights and measure the evolution of the pore water pressure over time. Several mixtures with variable content of fine sediment (silt and clay) and variable content of coarse sediment (with fixed relative fractions of grains between 2 and 32 mm) are tested. For the fines two types of clay (smectite and kaolinite) and loam (Stoober Lehm) are used. The analysis is based on the one-dimensional consolidation theory which uses a diffusion coefficient D to model the decay of excess fluid pressure over time. Starting from artificially induced super-hydrostatic fluid pressures, we find dissipation coefficients ranging from 10-5 m²/s for liquid mixtures to 10-8 m²/s for viscous mixtures. The results for kaolinite and smectite are quite similar. For our limited number of mixtures the effect of fines content is more pronounced than the effect of different amounts of coarse particles.

Shear enhanced compaction in a porous basalt from San Miguel Island, Azores

NASA Astrophysics Data System (ADS)

Loaiza, S.; Fortin, J.; Schubnel, A. J.; Vinciguerra, S.; Moreira, M.; Gueguen, Y.

2011-12-01

Basaltic rocks are the main component of the oceanic upper crust. This is of potential interest for water and geothermal resources, or for storage of CO2. The aim of our work is to investigate experimentally the mechanical behavior and the failure modes of porous basalt as well its permeability evolution during deformation. Cylindrical basalt samples, from the Azores, of 30 mm in diameter and 60 mm in length were deformed the triaxial cell at room temperature and at a constant axial strain rate of 10-5 s-1. The initial porosity of the sample was 18%. In our study, a set of experiments were performed at confining pressure in the range of 25-290 MPa. The samples were deformed under saturated conditions at a constant pore pressure of 5MPa. Two volumetric pumps kept the pore pressure constant, and the pore volume variations were recorded. The evolution of the porosity was calculated from the total volume variation inside the volumetric pumps. Permeability measurements were performed using the steady-state technique. Our result shows that two modes of deformation can be highlighted in this basalt. At low confining pressure (Pc < 50 MPa), the differential stress attains a peak before the sample undergoes strain softening; failure occurs by shear localization. The experiments performed at confining pressure higher than 50 MPa, show a totally different mode of deformation. In this second mode of deformation, an appreciable inelastic porosity reduction is observed. Comparing to the hydrostatic loading, the rock sample started to compact beyond a critical stress state; and from then, strain hardening, with stress drops are observed. Such a behavior is characteristic of the formation of compaction localization, due to grain crushing and pore collapse. In addition, this inelastic compaction is accompanied by a decrease of permeability, indicating that these compaction bands or zones act as barrier for fluid flow, in agreement with observations done in sandstone. Further studies, including Acoustic Emission locations and microstructural observations will be carried out in order to map the compaction bands or zones and confirm or infirm the formation of compaction localization, and the micromechanisms (pore collapse and grain crushing) taking place in this second mode of deformation.

Capillary Pressure of a Liquid Between Uniform Spheres Arranged in a Square-Packed Layer

NASA Technical Reports Server (NTRS)

Alexader, J. Iwan D.; Slobozhanin, Lev A.; Collicott, Steven H.

2004-01-01

The capillary pressure in the pores defined by equidimensional close-packed spheres is analyzed numerically. In the absence of gravity the menisci shapes are constructed using Surface Evolver code. This permits calculation the free surface mean curvature and hence the capillary pressure. The dependences of capillary pressure on the liquid volume constructed here for a set of contact angles allow one to determine the evolution of basic capillary characteristics under quasi-static infiltration and drainage. The maximum pressure difference between liquid and gas required for a meniscus passing through a pore is calculated and compared with that for hexagonal packing and with approximate solution given by Mason and Morrow [l]. The lower and upper critical liquid volumes that determine the stability limits for the equilibrium capillary liquid in contact with square packed array of spheres are tabulated for a set of contact angles.

Episodic Tremor and Slip Explained by Fluid-Enhanced Microfracturing and Sealing

NASA Astrophysics Data System (ADS)

Bernaudin, M.; Gueydan, F.

2018-04-01

Episodic tremor and slow-slip events at the deep extension of plate boundary faults illuminate seismic to aseismic processes around the brittle-ductile transition. These events occur in volumes characterized by overpressurized fluids and by near failure shear stress conditions. We present a new modeling approach based on a ductile grain size-sensitive rheology with microfracturing and sealing, which provides a mechanical and field-based explanation of such phenomena. We also model pore fluid pressure variation as a function of changes in porosity/permeability and strain rate-dependent fluid pumping. The fluid-enhanced dynamic evolution of microstructures defines cycles of ductile strain localization and implies increase in pore fluid pressure. We propose that slow-slip events are ductile processes related to transient strain localization, while nonvolcanic tremor corresponds to fracturing of the whole rock at the peak of pore fluid pressure. Our model shows that the availability of fluids and the efficiency of fluid pumping control the occurrence and the P-T conditions of episodic tremor and slip.

The 4D evolution of porosity during ongoing pressure-solution processes in NaCl using x-ray microtomography

NASA Astrophysics Data System (ADS)

Macente, Alice; Fusseis, Florian; Butler, Ian; Tudisco, Erika; Hall, Stephen; Andò, Edward

2016-04-01

Pressure-solution creep is a common deformation mechanism in the upper crust. It represents a mass transfer via dissolution-reprecipitation that critically affects the hydraulic properties of rocks. Successful management of safe radioactive storage sites in rock-salt deposits critically depends on an accurate knowledge of the hydro-mechanical behaviour of salt deposits. Despite numerous lab experiments that have been conducted, many aspects of pressure-solution are still poorly understood. There is little knowledge about the spatio-temporal evolution of porosity and permeability during pressure-solution creep. While rates of pressure-solution creep in silicates and carbonates are slow, which makes laboratory investigations of these materials impractical, compaction experiments have demonstrated that NaCl samples deform sufficiently fast to study pressure-solution creep in a lab environment at room temperature and modest loads. We present results from novel experiments that quantify the 4-dimensional (three spatial dimensions plus time) evolution of pressure-solution processes using in-situ x-ray microtomography. Our experiments are performed in custom made x-ray transparent presses. 5 mm diameter NaCl powder samples with a grain size of 250-300 μm are loaded dry into the press and pre-compacted to produce a starting aggregated material. The sample is then flooded with saturated NaCl solution and loaded uniaxially by means of a pneumatic actuator to a constant uniaxial stress. Different sample mixtures were tested, as well as different uniaxial loads. The resulting deformation of the samples is documented in 3-dimensional microtomographic datasets, acquired at regular time intervals. Image analysis allowed characterization of the microstructural evolution of the NaCl grains and the spatio-temporal distribution of porosity during ongoing mechanical and chemical compaction. The microtomography data have also been analysed with 3D Digital Image Correlation (3D-DIC or DVC) to quantify the fields of displacements in each direction, as well as volumetric and maximum shear strain fields. Following the approach described above, we have been able to quantify and characterize in 4D the evolution of pressure-solution creep and porosity distribution in relation to different sample materials and increasing uniaxial load. The presence of phyllosilicates (biotite) and more competent materials (glass beads) allowed pressure-solution to develop in a much shorter time compared to pure halite sample. The same trend is observed in samples experiencing bigger uniaxial loads (6.6 MPa v 1.6 MPa). We also found that, in the presence of phyllosilicates (biotite), pore size distribution clearly reflects the localisation of pressure-solution processes, as for natural stylolites. In the presence of glass beads, pressure-solution has a greater effect on the pore orientations rather than pore sizes. Our results extend the current understanding of the effect of pressure-solution creep on the mechanical and hydraulic properties of rocks, with implications for natural rock-salt, salt-based repository systems (nuclear and chemical waste storage) and salt mining.

Evolution of fracture permeability of ultramafic rocks undergoing serpentinization at hydrothermal conditions: An experimental study

NASA Astrophysics Data System (ADS)

Farough, A.; Moore, D. E.; Lockner, D. A.; Lowell, R. P.

2016-01-01

We performed flow-through laboratory experiments on five cylindrically cored samples of ultramafic rocks, in which we generated a well-mated through-going tensile fracture, to investigate evolution of fracture permeability during serpentinization. The samples were tested in a triaxial loading machine at a confining pressure of 50 MPa, pore pressure of 20 MPa, and temperature of 260°C, simulating a depth of 2 km under hydrostatic conditions. A pore pressure difference of up to 2 MPa was imposed across the ends of the sample. Fracture permeability decreased by 1-2 orders of magnitude during the 200-330 h experiments. Electron microprobe and SEM data indicated the formation of needle-shaped crystals of serpentine composition along the walls of the fracture, and chemical analyses of sampled pore fluids were consistent with dissolution of ferro-magnesian minerals. By comparing the difference between fracture permeability and matrix permeability measured on intact samples of the same rock types, we concluded that the contribution of the low matrix permeability to flow is negligible and essentially all of the flow is focused in the tensile fracture. The experimental results suggest that the fracture network in long-lived hydrothermal circulation systems can be sealed rapidly as a result of mineral precipitation, and generation of new permeability resulting from a combination of tectonic and crystallization-induced stresses is required to maintain fluid circulation.

Evolution of fracture permeability of ultramafic rocks undergoing serpentinization at hydrothermal conditions: An experimental study

USGS Publications Warehouse

Farough, Aida; Moore, Diane E.; Lockner, David A.; Lowell, R.P.

2016-01-01

We performed flow-through laboratory experiments on five cylindrically cored samples of ultramafic rocks, in which we generated a well-mated through-going tensile fracture, to investigate evolution of fracture permeability during serpentinization. The samples were tested in a triaxial loading machine at a confining pressure of 50 MPa, pore pressure of 20 MPa, and temperature of 260°C, simulating a depth of 2 km under hydrostatic conditions. A pore pressure difference of up to 2 MPa was imposed across the ends of the sample. Fracture permeability decreased by 1–2 orders of magnitude during the 200–330 h experiments. Electron microprobe and SEM data indicated the formation of needle-shaped crystals of serpentine composition along the walls of the fracture, and chemical analyses of sampled pore fluids were consistent with dissolution of ferro-magnesian minerals. By comparing the difference between fracture permeability and matrix permeability measured on intact samples of the same rock types, we concluded that the contribution of the low matrix permeability to flow is negligible and essentially all of the flow is focused in the tensile fracture. The experimental results suggest that the fracture network in long-lived hydrothermal circulation systems can be sealed rapidly as a result of mineral precipitation, and generation of new permeability resulting from a combination of tectonic and crystallization-induced stresses is required to maintain fluid circulation.

Significance of Dynamic Pore Pressure Variations - Comparison of Observations on Mud Volcanoes on the Costa Rica Margin and in the Gulf of Cadiz

NASA Astrophysics Data System (ADS)

Brueckmann, W.; Linke, P.; Pieper, M.; Hensen, C.; Tuerk, M.

2006-12-01

Research in the cooperative research center (SFB) 574 "Volatiles and Fluids in Subduction Zones" at the University Kiel focuses on volatile and fluid exchange processes at subduction zones. These have a significant impact on the long-term geochemical evolution of the hydrosphere and atmosphere. In the SFB 574 working area off Central America more than 120 mud volcanoes, mud diapirs and cold seeps have been identified and sampled. To better understand the internal dynamics of these structures and the temporal variability of fluid expulsion an in-situ tool for monitoring shallow pore pressure variations was devised. The tool (PWPL) monitors pore pressure variations along a 2m profile in the shallow subsurface using a stinger with 4 pressure ports. Positioned with a video-guided lander the stinger is gently pushed into the seafloor where it remains for several weeks or months in autonomous mode before being retrieved. While particular emphasis was placed on the convergent margin of Central America, mud volcanoes in other tectonic settings suitable for long-term observations of fluid flux are used for comparison. Here we will present data and interpretations from two mud volcanoes off Costa Rica and in the Gulf of Cadiz where we have conducted successful tests. Pore pressure data from short-term tests on Mound 11 on the continental slope off Costa Rica are compared with new results from a long-term (3-month) campaign on the Captain Arutjunov deep water mud volcano in the Gulf of Cadiz. Rates of fluid flow at both structures have been thoroughly characterized and quantified with geochemical methods providing a frame of reference for judging the significance of dynamic pore pressure variations.

Pore Pressure Evolution in Shallow Subduction Earthquake Sequences and Effects on Aseismic Slip Transients -- Numerical Modeling With Rate and State Friction

NASA Astrophysics Data System (ADS)

Liu, Y.; Rice, J. R.

2005-12-01

In 3D modeling of long tectonic loading and earthquake sequences on a shallow subduction fault [Liu and Rice, 2005], with depth-variable rate and state friction properties, we found that aseismic transient slip episodes emerge spontaneously with only a simplified representation of effects of metamorphic fluid release. That involved assumption of a constant in time but uniformly low effective normal stress in the downdip region. As suggested by observations in several major subduction zones [Obara, 2002; Rogers and Dragert, 2003; Kodaira et al, 2004], the presence of fluids, possibly released from dehydration reactions beneath the seismogenic zone, and their pressurization within the fault zone may play an important role in causing aseismic transients and associated non-volcanic tremors. To investigate the effects of fluids in the subduction zone, particularly on the generation of aseismic transients and their various features, we develop a more complete physical description of the pore pressure evolution (specifically, pore pressure increase due to supply from dehydration reactions and shear heating, decrease due to transport and dilatancy during slip), and incorporate that into the rate and state based 3D modeling. We first incorporated two important factors, dilatancy and shear heating, following Segall and Rice [1995, 2004] and Taylor [1998]. In the 2D simulations (slip varies with depth only), a dilatancy-stabilizing effect is seen which slows down the seismic rupture front and can prevent rapid slip from extending all the way to the trench, similarly to Taylor [1998]. Shear heating increases the pore pressure, and results in faster coseismic rupture propagation and larger final slips. In the 3D simulations, dilatancy also stabilizes the along-strike rupture propagation of both seismic and aseismic slips. That is, aseismic slip transients migrate along the strike faster with a shorter Tp (the characteristic time for pore pressure in the fault core to re-equilibrate with that of its surroundings). This is consistent with our previous simulations, which show that the aseismic transients migrate along the strike at a higher speed under a lower, constant in time, effective normal stress. As a combination of the two factors, we show the pore pressure evolution with drops (due to dilatancy during slip) and then rises (due to shear heating) on the fault over multiple time scales. We next plan to formulate, and merge with the slip-rupture analysis, fuller fluid release models based on phase equilibria and models of transport in which the average fault-parallel permeability is a decreasing function of the effective normal stress. The thrust fault zone, at seismogenic depths and slightly downdip, is represented in a conceptually similar manner to the well-studied major continental faults, assuming the fault core materials have a lower permeability than the neighboring damaged zone. Heat diffusion in the fault core and damaged zone will also be considered in the modeling. The simulation results may help to improve our understanding of the processes of the aseismic transients observed within a transform plate boundary along the SAF near Cholame, California [Nadeau and Dolenc, 2005].

Mechanical behavior and localized failure modes in a porous basalt from the Azores

NASA Astrophysics Data System (ADS)

Loaiza, S.; Fortin, J.; Schubnel, A.; Guéguen, Y.; Moreira, M.; Vinciguerra, S.

2012-04-01

Basaltic rocks are the main component of the oceanic upper crust. This is of potential interest for water and geothermal resources, or for storage of CO2. The aim of our work is to investigate experimentally the mechanical behavior and the failure modes of porous basalt as well as the permeability evolution during deformation. Cylindrical basalt samples, from the Azores, of 30 mm in diameter and 60 mm in length were deformed the triaxial cell of the Laboratoire de Geologie at the Ecole Normale Supérieure (Paris) at room temperature and at a constant axial strain rate of 10-5 s-1. The initial porosity of the sample was 18%. The Geodesign triaxial cell can reach 300MPa confining pressure; axial load is performed through a piston and can reach 900 MPa (for a 30mm diameter sample); maximum pore pressure is 100MPa (applied using two precision volumetric pumps). In our study, a set of experiments were performed at confining pressure in the range of 25-290 MPa. The samples were deformed under saturated conditions at a constant pore pressure of 5MPa. Two volumetric pumps kept the pore pressure constant, and the pore volume variations were recorded. The evolution of the porosity was calculated from the total volume variation inside the volumetric pumps. Permeability measurements were performed using the steady-state technique. Our result shows that two modes of deformation can be highlighted in this basalt. At low confining pressure (Pc < 50 MPa), the differential stress attains a peak before the sample undergoes strain softening; the failure of sample occurs by shear localization. Yet, the brittle regime is commonly observed in this low Pc range, the experiments performed at confining pressure higher than 50 MPa, show a totally different mode of deformation. In this second mode of deformation, an appreciable inelastic porosity reduction is observed. Comparing to the hydrostatic loading, the rock sample started to compact beyond a critical stress state; and from then, strain hardening, with stress drops are observed. Such a behavior is characteristic of the formation of compaction localization, due to grain crushing and pore collapse. In addition, this inelastic compaction is accompanied by a decrease of permeability, indicating that these compaction bands or zones act as barrier for fluid flow, in agreement with observations done in sandstone (Fortin et al., 2005). Further studies, including microstructural observations carried out by mapping the compaction bands or zones throughout a mosaic of SEM images at high resolution and acoustic emission recording will be carried in order to confirm the formation of compaction localization, and the micromechanisms (pore collapse and grain crushing) taking place in this second mode of deformation.

Poro-elasto-capillary wicking of cellulose sponges

PubMed Central

Kim, Do-Nyun

2018-01-01

We mundanely observe cellulose (kitchen) sponges swell while absorbing water. Fluid flows in deformable porous media, such as soils and hydrogels, are classically described on the basis of the theories of Darcy and poroelasticity, where the expansion of media arises due to increased pore pressure. However, the situation is qualitatively different in cellulosic porous materials like sponges because the pore expansion is driven by wetting of the surrounding cellulose walls rather than by increase of the internal pore pressure. We address a seemingly so simple but hitherto unanswered question of how fast water wicks into the swelling sponge. Our experiments uncover a power law of the wicking height versus time distinct from that for nonswelling materials. The observation using environmental scanning electron microscopy reveals the coalescence of microscale wall pores with wetting, which allows us to build a mathematical model for pore size evolution and the consequent wicking dynamics. Our study sheds light on the physics of water absorption in hygroscopically responsive multiscale porous materials, which have far more implications than everyday activities (for example, cleaning, writing, and painting) carried out with cellulosic materials (paper and sponge), including absorbent hygiene products, biomedical cell cultures, building safety, and cooking. PMID:29682606

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.

Permeability-porosity relationship for compaction of a low-permeability creeping material : Experimental evaluation using a single transient test

NASA Astrophysics Data System (ADS)

Ghabezloo, S.; Sulem, J.; Saint-Marc, J.

2009-04-01

It is well-known that there is no unique permeability-porosity relationship that can be applied to all porous materials. For a given evolution process that changes both permeability and porosity of a porous material, for example elastic or plastic compaction, microcracking or chemical alteration, it is usually assumed that there is an empirical relationship in the form of a power-law or exponential relationship between these parameters. The coefficients of these empirical relationships depend strongly on the properties of the material and of the evolution process. For the case of the power-law permeability-porosity relationship, a review of the literature shows that the exponent of this relation may be integer or non-integer, constant or variable, and the reported values of exponent vary between 1.1 and 25.4 for different materials and evolution processes, but no clear correlation between the exponenet and the petrophysical properties could be found. This wide variability of the permeability-porosity relationship highlights the necessity of experimental evaluation of this relationship for each material and evolution process. An experimental method is presented for the evaluation of a permeability-porosity relationship in a low-permeability porous material using the results of a single transient test. This method accounts for both elastic and non-elastic deformations of the sample during the test and is applied to a hardened class G oil well cement paste. An initial hydrostatic undrained loading is applied to the sample which generates an excess pore pressure, related to the applied hydrostatic stress by the Skempton coefficient of the material. The generated excess pore pressure is then released at one end of the sample while monitoring the pore pressure at the other end and the radial strain in the middle of the sample during the dissipation of the pore pressure. These measurements are back analysed using a finite-difference numerical scheme to evaluate the permeability and its evolution with porosity change. The stress-dependent character of the poroelastic parameters of the hardened cement paste (Ghabezloo et al., 2008) and also the creep of the material during the test add some particular aspects to the back-analysis, which makes this problem different from the classical solutions of transient permeability evaluation tests. The effect of creep of the sample during the test on the measured pore pressure and volume change is taken into account in the analysis. This approach permits to calibrate a power law permeability-porosity relationship for the tested hardened cement paste and also two parameters of a viscoelastic model for the creep of the material. The porosity sensitivity exponent of the power-law is evaluated equal to 11 and is shown to be mostly independent of the stress level and of the creep strains. The proposed method can be applied to different low permeability porous materials and for the case of non-creeping materials, the same type of analysis can be used to calibrate either a permeability-porosity or a permeability-effective stress relationship for the compaction of the tested material using a single transient test. References: 1.Ghabezloo S., Sulem J., Saint-Marc, J. (2008) Evaluation of a permeability-porosity relationship in a low permeability creeping material using a single transient test. Int J Rock Mech Min Sci, in press, DOI 10.1016/j.ijrmms.2008.10.003. 2.Ghabezloo, S., Sulem, J., Guédon, S., Martineau, F., Saint-Marc, J. (2008) Poromechanical behaviour of hardened cement paste under isotropic loading. Cement and Concrete Research, 38(12), 1424-1437.

Rate dependent deformation of porous sandstone across the brittle-ductile transition

NASA Astrophysics Data System (ADS)

Jefferd, M.; Brantut, N.; Mitchell, T. M.; Meredith, P. G.

2017-12-01

Porous sandstones transition from dilatant, brittle deformation at low pressure, to compactant, ductile deformation at high pressure. Both deformation modes are driven by microcracking, and are expected to exhibit a time dependency due to chemical interactions between the pore fluid and the rock matrix. In the brittle regime, time-dependent failure and brittle creep are well documented. However, much less is understood in the ductile regime. We present results from a series of triaxial deformation experiments, performed in the brittle-ductile transition zone of fluid saturated Bleurswiller sandstone (initial porosity = 23%). Samples were deformed at 40 MPa effective pressure, to 4% axial strain, under either constant strain rate (10-5 s-1) or constant stress (creep) conditions. In addition to stress, axial strain and pore volume change, P wave velocities and acoustic emission were monitored throughout. During constant stress tests, the strain rate initially decreased with increasing strain, before reaching a minimum and accelerating to a constant level beyond 2% axial strain. When plotted against axial strain, the strain rate evolution under constant stress conditions, mirrors the stress evolution during the constant strain rate tests; where strain hardening occurs prior to peak stress, which is followed by strain softening and an eventual plateau. In all our tests, the minimum strain rate during creep occurs at the same inelastic strain as the peak stress during constant strain tests, and strongly decreases with decreasing applied stress. The microstructural state of the rock, as interpreted from similar volumetric strain curves, as well as the P-wave velocity evolution and AE production rate, appears to be solely a function of the total inelastic strain, and is independent of the length of time required to reach said strain. We tested the sensitivity of fluid chemistry on the time dependency, through a series of experiments performed under similar stress conditions, but with chemically inert decane instead of water as the pore fluid. Under the same applied stress, decane saturated samples reached a minimum strain rate 2 orders of magnitude lower than the water saturated samples. This is consistent with a mechanism of subcritical crack growth driven by chemical interactions between the pore fluid and the rock.

Pore-scale simulations of drainage in granular materials: Finite size effects and the representative elementary volume

NASA Astrophysics Data System (ADS)

Yuan, Chao; Chareyre, Bruno; Darve, Félix

2016-09-01

A pore-scale model is introduced for two-phase flow in dense packings of polydisperse spheres. The model is developed as a component of a more general hydromechanical coupling framework based on the discrete element method, which will be elaborated in future papers and will apply to various processes of interest in soil science, in geomechanics and in oil and gas production. Here the emphasis is on the generation of a network of pores mapping the void space between spherical grains, and the definition of local criteria governing the primary drainage process. The pore space is decomposed by Regular Triangulation, from which a set of pores connected by throats are identified. A local entry capillary pressure is evaluated for each throat, based on the balance of capillary pressure and surface tension at equilibrium. The model reflects the possible entrapment of disconnected patches of the receding wetting phase. It is validated by a comparison with drainage experiments. In the last part of the paper, a series of simulations are reported to illustrate size and boundary effects, key questions when studying small samples made of spherical particles be it in simulations or experiments. Repeated tests on samples of different sizes give evolution of water content which are not only scattered but also strongly biased for small sample sizes. More than 20,000 spheres are needed to reduce the bias on saturation below 0.02. Additional statistics are generated by subsampling a large sample of 64,000 spheres. They suggest that the minimal sampling volume for evaluating saturation is one hundred times greater that the sampling volume needed for measuring porosity with the same accuracy. This requirement in terms of sample size induces a need for efficient computer codes. The method described herein has a low algorithmic complexity in order to satisfy this requirement. It will be well suited to further developments toward coupled flow-deformation problems in which evolution of the microstructure require frequent updates of the pore network.

Effects of stress paths on physical properties of sediments at the Nankai Trough subduction zone

NASA Astrophysics Data System (ADS)

Kitajima, H.; Saffer, D. M.

2011-12-01

Stress states are one of the most important factors governing deformation modes and fault strength. In subduction systems where tectonic stress is large, sediments are subjected to complicated stress conditions in time and space. Because direct measurements of stress are very limited, stress conditions at depths have been estimated by combining seismic reflection data with empirical relations between compressional-wave, porosity, and effective stress [Tsuji et al., 2008; Tobin and Saffer, 2009]. However, most of the empirical relations are derived from experiments conducted under isotropic conditions, and do not account for the more complicated stress states expected in active subduction-accretion complexes. In this study, we aim to derive relations between physical properties and stress states from triaxial deformation experiments on sediments. During the Integrated Ocean Drilling Program (IODP) Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) Expeditions 314, 315, 319, 322, and 333, core samples were recovered from shallow boreholes into the accretionary prism and two sites seaward of the deformation front (reference sites). We used core samples from reference sites (Sites C0011 and C0012) for this study because they represent input material for the subduction system, and have not been subjected to tectonic compression in the accretionary wedge. In our deformation tests, samples are loaded under a range of different stress paths including isotropic loading, triaxial compression, and triaxial extension by controlling axial stress (up to 100 MPa), confining pressure (up to 100 MPa), and pore pressure (0.5-28 MPa). During tests, all pressures, axial displacement, and pore volume change were monitored. Permeability, and ultrasonic velocity were also measured during the tests. Two experiments have been conducted on samples taken from the core 322-C0011B-19R-5 (Lower Shikoku Basin hemipelagic mudstone, initial porosity of 43 %). The first test was conducted under istotropic loading and unloading by (1) increase and decrease in confining pressure, and (2) decrease and increase in pore pressure. The evolution of physical properties depends on effective pressure regardless of whether confining pressure or pore pressure is controlled. As effective pressure increases from 0.2 to 30 MPa, porosity decreases from 43 to 18 %, permeability decreases from 1.1×10-18 to 4.1×10-20 m2, and compressional-wave velocity increases from 1.76 to 2.5 km/s, respectively. The same physical properties do not fully recover during unloading, which corresponds to overconsolidated or overpressured condition. The second test included various loading paths including triaxial compression and extension, and drained and undrained condition of pore pressure. The results indicate that the evolution of physical properties be dependent on both effective mean stress and differential stress. The experimental results suggest that it is important to consider consolidation state and loading paths. We will present more experimental results and derive relations between physical properties and stress states.

A mathematical multiscale model of bone remodeling, accounting for pore space-specific mechanosensation.

PubMed

Pastrama, Maria-Ioana; Scheiner, Stefan; Pivonka, Peter; Hellmich, Christian

2018-02-01

While bone tissue is a hierarchically organized material, mathematical formulations of bone remodeling are often defined on the level of a millimeter-sized representative volume element (RVE), "smeared" over all types of bone microstructures seen at lower observation scales. Thus, there is no explicit consideration of the fact that the biological cells and biochemical factors driving bone remodeling are actually located in differently sized pore spaces: active osteoblasts and osteoclasts can be found in the vascular pores, whereas the lacunar pores host osteocytes - bone cells originating from former osteoblasts which were then "buried" in newly deposited extracellular bone matrix. We here propose a mathematical description which considers size and shape of the pore spaces where the biological and biochemical events take place. In particular, a previously published systems biology formulation, accounting for biochemical regulatory mechanisms such as the rank-rankl-opg pathway, is cast into a multiscale framework coupled to a poromicromechanical model. The latter gives access to the vascular and lacunar pore pressures arising from macroscopic loading. Extensive experimental data on the biological consequences of this loading strongly suggest that the aforementioned pore pressures, together with the loading frequency, are essential drivers of bone remodeling. The novel approach presented here allows for satisfactory simulation of the evolution of bone tissue under various loading conditions, and for different species; including scenarios such as mechanical dis- and overuse of murine and human bone, or in osteocyte-free bone. Copyright © 2017 Elsevier Inc. All rights reserved.

Cavitation and pore blocking in nanoporous glasses.

PubMed

Reichenbach, C; Kalies, G; Enke, D; Klank, D

2011-09-06

In gas adsorption studies, porous glasses are frequently referred to as model materials for highly disordered mesopore systems. Numerous works suggest that an accurate interpretation of physisorption isotherms requires a complete understanding of network effects upon adsorption and desorption, respectively. The present article deals with nitrogen and argon adsorption at different temperatures (77 and 87 K) performed on a series of novel nanoporous glasses (NPG) with different mean pore widths. NPG samples contain smaller mesopores and significantly higher microporosity than porous Vycor glass or controlled pore glass. Since the mean pore width of NPG can be tuned sensitively, the evolution of adsorption characteristics with respect to a broadening pore network can be investigated starting from the narrowest nanopore width. With an increasing mean pore width, a H2-type hysteresis develops gradually which finally transforms into a H1-type. In this connection, a transition from a cavitation-induced desorption toward desorption controlled by pore blocking can be observed. Furthermore, we find concrete hints for a pore size dependence of the relative pressure of cavitation in highly disordered pore systems. By comparing nitrogen and argon adsorption, a comprehensive insight into adsorption mechanisms in novel disordered materials is provided. © 2011 American Chemical Society

Nano-Pore Size Analysis by SAXS Method of Cementitious Mortars Undergoing Delayed Ettringite Formation

NASA Astrophysics Data System (ADS)

Shekar, Yamini

This research investigates the nano-scale pore structure of cementitious mortars undergoing delayed ettringite formation (DEF) using small angle x-ray scattering (SAXS). DEF has been known to cause expansion and cracking during later ages (around 4000 days) in concrete that has been heat cured at temperatures of 70°C or above. Though DEF normally occurs in heat cured concrete, mass cured concrete can also experience DEF. Large crystalline pressures result in smaller pore sizes. The objectives of this research are: (1) to investigate why some samples expand early than later expansion, (2) to evaluate the effects of curing conditions and pore size distributions at high temperatures, and (3) to assess the evolution of the pore size distributions over time. The most important outcome of the research is the pore sizes obtained from SAXS were used in the development of a 3-stage model. From the data obtained, the pore sizes increase in stage 1 due to initial ettringite formation and in turn filling up the smallest pores. Once the critical pore size threshold is reached (around 20nm) stage 2 is formed due to cracking which tends to decrease in the pore sizes. Finally, in stage 3, the cracking continues, therefore increasing in the pore size.

Investigation on the Permeability Evolution of Gypsum Interlayer Under High Temperature and Triaxial Pressure

NASA Astrophysics Data System (ADS)

Tao, Meng; Yechao, You; Jie, Chen; Yaoqing, Hu

2017-08-01

The permeability of the surrounding rock is a critical parameter for the designing and assessment of radioactive waste disposal repositories in the rock salt. Generally, in the locations that are chosen for radioactive waste storage, the bedded rock salt is a sedimentary rock that contains NaCl and Na2SO4. Most likely, there are also layers of gypsum ( {CaSO}_{ 4} \\cdot 2 {H}_{ 2} {O)} present in the salt deposit. Radioactive wastes emit a large amount of heat and hydrogen during the process of disposal, which may result in thermal damage of the surrounding rocks and cause a great change in their permeability and tightness. Therefore, it is necessary to investigate the permeability evolution of the gypsum interlayer under high temperature and high pressure in order to evaluate the tightness and security of the nuclear waste repositories in bedded rock salt. In this study, a self-designed rock triaxial testing system by which high temperature and pressure can be applied is used; the μCT225kVFCB micro-CT system is also employed to investigate the permeability and microstructure of gypsum specimens under a constant hydrostatic pressure of 25 MPa, an increasing temperature (ranging from 20 to 650 °C), and a variable inlet gas pressure (1, 2, 4, 6 MPa). The experimental results show: (a) the maximum permeability measured during the whole experiment is less than 10-17 m2, which indicates that the gypsum interlayer has low permeability under high temperature and pressure that meet the requirements for radioactive waste repository. (b) Under the same temperature, the permeability of the gypsum specimen decreases at the beginning and then increases as the pore pressure elevates. When the inlet gas pressure is between 0 and 2 MPa, the Klinkenberg effect is very pronounced. Then, as the pore pressure increases, the movement behavior of gas molecules gradually changes from free motion to forced directional motion. So the role of free movement of gas molecules gradually reduced, which eventually leads to a decrease in permeability. When the inlet gas pressure is between 2 and 6 MPa, the Klinkenberg effect dribbles away, and the gas flow gradually obeys to the Darcy's law. Hence, the permeability increased with the increase in inlet gas pressure. (c) The curve of permeability versus temperature is divided into five stages based on its gradient. In the temperature range of 20-100 °C, the permeability of gypsum decreased slowly when the temperature decreased. From 100 to 200 °C, the permeability of gypsum increased dramatically when the temperature increased. However, a dramatic increase in permeability was observed from 200 to 450 °C. Subsequently, in the temperature range of 450-550 °C, due to closure of pores and fractures, the permeability of the specimens slowly lessened when the temperature increased. From 550 to 650 °C, the permeability of gypsum slightly increased when the temperature increased; (d) the micro-cracks and porosity obtained from the CT images show a high degree of consistency to the permeability evolution; (e) when compared to the permeability evolutions of sandstone, granite, and lignite, gypsum exhibits a stable evolution trend of permeability and has a much greater threshold temperature when its permeability increases sharply. The results of the paper may provide essential and valuable references for the design and construction of high-level radioactive wastes repository in bedded salt rock containing gypsum interlayers.

Constraints on Pore Pressure in Subduction Zones From Geotechnical Tests and Physical Properties Data

NASA Astrophysics Data System (ADS)

Saffer, D. M.; McKiernan, A. W.

2005-12-01

At subduction zones, as incoming sediments are either offscraped or underthrust at the trench, elevated pore pressures result from the combination of rapid loading and low permeability. Pore pressure within underthrust sediment is especially important for the mechanical strength of the plate boundary fault system, because the main décollement localizes immediately above this sediment, and at many subduction zones steps downward into it. Because the underthrust sediment undergoes progressive uniaxial (vertical) strain, quantitative estimates of in situ pore pressure can be obtained by several methods, including: (1) maximum past burial stress ( Pv'}) from laboratory consolidation tests on core samples, and (2) observed compaction trends in boreholes. These methods allow a detailed view of pore pressure and its variability down-section, providing insight into dewatering processes and the evolution of shear strength relevant to early development of the décollement. Geotechnical tests also provide independent measurement of the coefficient of consolidation ( Cv), compressibility ( mv), and permeability (k) of sediment samples, which can be used to parameterize forward models of pressure generation. Here, I discuss pore pressure estimates derived from (1) consolidation tests on core samples, and (2) observed porosity profiles, along transects where ODP drilling has sampled sediment at the Nankai, N. Barbados, and Costa Rican subduction zones. At all three margins, the two independent methods yield consistent results, and indicate development of significant overpressures that increase systematically with distance from the trench. The values are in good agreement with direct measurements in 2 instrumented boreholes at Barbados, maximum and minimum bounds from the known loading rate, and results of 2-D numerical models of fluid flow. Inferred pressures document nearly undrained conditions at the base of the section (excess pressures equal to the load emplaced by subduction burial), and partially drained conditions at the top (excess pressures of ~40% of the undrained response at Costa Rica, ~50-60% at Nankai, and ~90-100% at Barbados). The spatial pattern of excess pore pressure is most consistent with upward drainage to a highly permeable décollement, to distances of at least 5-10 km landward of the trench. When directly measured values of mv and k from laboratory geotechnical experiments are incorporated into simple 1-D models of vertical dewatering, simulated pore pressures are consistent with those inferred from consolidation tests and porosity data. Model results suggest that severe underconsolidation should persist for tens of km from the trench; notably, simulated underconsolidation is diminished by 20-30 km landward of the trench at Nankai, broadly coincident with the locations of both diminished seismic reflection amplitude observed at the décollement and the updip extent of coseismic slip. The consistent results achieved at these three margins indicate that: (1) geotechnical tests can provide viable estimates of in situ pore pressure, at least at shallow depths, and (2) laboratory-derived values of permeability and sediment compressibility may be representative of in situ properties, despite collection at small spatial scale and over short times. However, significant uncertainty exists in projecting models to greater depth using geotechnical parameters from shallow samples; more detailed laboratory investigations are clearly needed to better understand the roles of temperature, rate, and diagenetic effects.

An analysis of transient flow in upland watersheds: interactions between structure and process

Treesearch

David Lawrence Brown

1995-01-01

The physical structure and hydrological processes of upland watersheds interact in response to forcing functions such as rainfall, leading to storm runoff generation and pore pressure evolution. Transient fluid flow through distinct flow paths such as the soil matrix, macropores, saprolite, and bedrock may be viewed as a consequence of such interactions. Field...

Pore pressure migration during hydraulic stimulation due to permeability enhancement by low-pressure subcritical fracture slip

NASA Astrophysics Data System (ADS)

Mukuhira, Yusuke; Moriya, Hirokazu; Ito, Takatoshi; Asanuma, Hiroshi; Häring, Markus

2017-04-01

Understanding the details of pressure migration during hydraulic stimulation is important for the design of an energy extraction system and reservoir management, as well as for the mitigation of hazardous-induced seismicity. Based on microseismic and regional stress information, we estimated the pore pressure increase required to generate shear slip on an existing fracture during stimulation. Spatiotemporal analysis of pore pressure migration revealed that lower pore pressure migrates farther and faster and that higher pore pressure migrates more slowly. These phenomena can be explained by the relationship between fracture permeability and stress state criticality. Subcritical fractures experience shear slip following smaller increases of pore pressure and promote migration of pore pressure because of their enhanced permeability. The difference in migration rates between lower and higher pore pressures suggests that the optimum wellhead pressure is the one that can stimulate relatively permeable fractures, selectively. Its selection optimizes economic benefits and minimizes seismic risk.

From rock to magma and back again: The evolution of temperature and deformation mechanism in conduit margin zones

NASA Astrophysics Data System (ADS)

Heap, Michael J.; Violay, Marie; Wadsworth, Fabian B.; Vasseur, Jérémie

2017-04-01

Explosive silicic volcanism is driven by gas overpressure in systems that are inefficient at outgassing. The zone at the margin of a volcanic conduit-thought to play an important role in the outgassing of magma and therefore pore pressure changes and explosivity-is the boundary through which heat is exchanged from the hot magma to the colder country rock. Using a simple heat transfer model, we first show that the isotherm for the glass transition temperature (whereat the glass within the groundmass transitions from a glass to an undercooled liquid) moves into the country rock when the magma within the conduit can stay hot, or into the conduit when the magma is quasi-stagnant and cools (on the centimetric scale over days to months). We then explore the influence of a migrating viscous boundary on compactive deformation micromechanisms in the conduit margin zone using high-pressure (effective pressure of 40 MPa), high-temperature (up to 800 °C) triaxial deformation experiments on porous andesite. Our experiments show that the micromechanism facilitating compaction in andesite is localised cataclastic pore collapse at all temperatures below the glass transition of the amorphous groundmass glass Tg (i.e., rock). In this regime, porosity is only reduced within the bands of crushed pores; the porosity outside the bands remains unchanged. Further, the strength of andesite is a positive function of temperature below the threshold Tg due to thermal expansion driven microcrack closure. The micromechanism driving compaction above Tg (i.e., magma) is the distributed viscous flow of the melt phase. In this regime, porosity loss is distributed and is accommodated by the widespread flattening and closure of pores. We find that viscous flow is much more efficient at reducing porosity than cataclastic pore collapse, and that it requires stresses much lower than those required to form bands of crushed pores. Our study therefore highlights that temperature excursions can result in a change in deformation micromechanism that drastically alters the mechanical and hydraulic properties of the material within the conduit margin zone, with possible implications for pore pressure augmentation and explosive behaviour.

Pore Pressure Pulse Drove the 2012 Emilia (Italy) Series of Earthquakes

NASA Astrophysics Data System (ADS)

Pezzo, Giuseppe; De Gori, Pasquale; Lucente, Francesco Pio; Chiarabba, Claudio

2018-01-01

The 2012 Emilia earthquakes sequence is the first debated case in Italy of destructive event possibly induced by anthropic activity. During this sequence, two main earthquakes occurred separated by 9 days on contiguous thrust faults. Scientific commissions engaged by the Italian government reported complementary scenarios on the potential trigger mechanism ascribable to exploitation of a nearby oil field. In this study, we combine a refined geodetic source model constrained by precise aftershock locations and an improved tomographic model of the area to define the geometrical relation between the activated faults and investigate possible triggering mechanisms. An aftershock decay rate that deviates from the classical Omori-like pattern and Vp/Vs changes along the fault system suggests that natural pore pressure pulse drove the space-time evolution of seismicity and the activation of the second main shock.

Experimental proof of the existence of a Widom line based on peculiarities of the behavior of hydrogen in nanoporous silicate at -45°C and atmospheric pressure

NASA Astrophysics Data System (ADS)

Bordonskii, G. S.; Gurulev, A. A.

2017-04-01

We have experimentally studied the thermal and microwave properties of a nanoporous medium (silica gel) with hydrogen-filled pores. On cooling down to about -45°C at atmospheric pressure, the system exhibited chemical transformations, a first-order phase transition with heat evolution, and a sharp change in the power of microwave radiation at 34 GHz transmitted through a sample. It is concluded that this point on the phase diagram corresponds to a point on the Widom line featuring sharp increase in fluctuations of the entropy and density of supercooled water formed during hydrogen interaction with the surface of pores in silica gel. These results confirm the existence of a second critical point of water, from which the Widom line originates.

Transient Deformation Patterns in Response to Quaternary Glacial Advance-Retreat Across the Offshore St. Elias Mountains, southern Alaska

NASA Astrophysics Data System (ADS)

Worthington, L. L.; Clary, W. A.; Daigle, H.; Koons, P. O.; Gulick, S. P. S.; Jaeger, J. M.

2016-12-01

The southern Alaska margin, home to the St. Elias Mountains, the highest coastal mountain range on Earth experiencing the highest erosion rates on Earth, provides a superb setting for evaluating competing influences of rheological and climate control on orogen development. Previous studies have recognized this potential, but conclusions were limited due to the absence of information on the time-dependent behavior of climate and rheological processes. These limitations can now be surpassed due to 1) the recent availability of high-precision age constraints on the structural and stratigraphic evolution of offshore sediments and structures and 2) geotechnical information on the extent of dewatering and related spatial changes in the material properties of these sediments. We correlate emerging results from Integrated Ocean Drilling Program (IODP) Expedition 341 Sites U1420 and U1421 with regional seismic data across the continental shelf and slope to determine the spatial and temporal evolution of thrusting in response to Yakutat-North American convergence. Our mapping shows that the pattern of faulting changed from distributed across the shelf to highly localized away from the primary glacial depocenter over the course of one glacial cycle. Core samples suggest that the glacially derived sediment is overpressured, with pore pressures possibly reaching >90% of lithostatic stress. Elevated pore pressures develop rapidly in response to focused glaciomarine sedimentation, in addition to direct ice loading, and may induce a transient state of wedge reorganization manifested as a change in localization of deformation. This relationship suggests that the additive response of pore pressure variations over glacial cycles throughout the Pleistocene and Holocene result in constant reorganization of deformation style and location.

Overpressure, Flow Focusing, Compaction and Slope Stability on the continental slope: Insights from IODP Expedition 308

NASA Astrophysics Data System (ADS)

Flemings, P. B.

2010-12-01

Integrated Ocean Drilling Program Expepedition 308 used direct measurements of pore pressure, analysis of hydromechanical properties, and geological analysis to illuminate how sedimentation, flow focusing, overpressure, and slope stability couple beneath the seafloor on the deepwater continental slope in the Gulf of Mexico. We used pore pressure penetrometers to measure severe overpressures (60% of the difference between lithostatic stress and hydrostatic pressure) that extend from the seafloor for 100’s of meters. We ran uniaxial consolidation experiments on whole core and found that although permeability is relatively high near the seafloor, the sediments are highly compressible. As a result, the coefficient of consolidation (the hydraulic diffusivity) is remarkably constant over a large range of effective stresses. This behavior accounts for the high overpressure that begins near the seafloor and extends to depth. Forward modeling suggests that flow is driven laterally along a permeable unit called the Blue Unit. Calculations suggest that soon after deposition, lateral flow lowered the effective stress and triggered the submarine landslides that we observe. Later in the evolution of this system, overpressure may have pre-conditioned the slope to failure by earthquakes. Results from IODP Expedition 308 illustrate how pore pressure and sedimentation control the large-scale form of continental margins, how submarine landslides form, and provide strategies for designing stable drilling programs.

An experimental study of the effects of adsorbing and non-adsorbing gases on friction and permeability evolution in clay-rich fault gouge

NASA Astrophysics Data System (ADS)

Lisabeth, H. P.; Zoback, M. D.

2017-12-01

Understanding the flow of fluids through fractures in clay-rich rocks is fundamental to a number of geoengineering enterprises, including development of unconventional hydrocarbon resources, nuclear waste storage and geological carbon sequestration. High clay content tends to make rocks plastic, low-porosity and anisotropic. In addition, some gasses adsorb to clay mineral surfaces, resulting in swelling and concomitant changes in physical properties. These complexities can lead to coupled behaviors that render prediction of fluid behavior in the subsurface difficult. We present the results of a suite of triaxial experiments on binary mixtures of quartz and illite grains to separate and quantify the effects of hydrostatic pressure, differential stress, clay content and gas chemistry on the evolution of mechanical and hydraulic characteristics of the gouge material during deformation. Tests are run on saw-cut samples prepared with gouge at 20 MPa confining pressure, 10 MPa pore pressure and at room temperature. Argon or carbon dioxide is used as pore fluid. Sample permeability, stress and strain are monitored continuously during hydrostatic and axial deformation. We find that pressure and shearing both lead to reductions in permeability. Adsorbing gas leads to swelling and promotes permeability reduction, but appears to have no effect on frictional properties. These results indicate that the seal integrity of clay-rich caprocks may not be compromised by shear deformation, and that depletion and shear deformation of unconventional reservoirs is expected to result in production declines.

The physics of debris flows

NASA Astrophysics Data System (ADS)

Iverson, Richard M.

1997-08-01

Recent advances in theory and experimentation motivate a thorough reassessment of the physics of debris flows. Analyses of flows of dry, granular solids and solid-fluid mixtures provide a foundation for a comprehensive debris flow theory, and experiments provide data that reveal the strengths and limitations of theoretical models. Both debris flow materials and dry granular materials can sustain shear stresses while remaining static; both can deform in a slow, tranquil mode characterized by enduring, frictional grain contacts; and both can flow in a more rapid, agitated mode characterized by brief, inelastic grain collisions. In debris flows, however, pore fluid that is highly viscous and nearly incompressible, composed of water with suspended silt and clay, can strongly mediate intergranular friction and collisions. Grain friction, grain collisions, and viscous fluid flow may transfer significant momentum simultaneously. Both the vibrational kinetic energy of solid grains (measured by a quantity termed the granular temperature) and the pressure of the intervening pore fluid facilitate motion of grains past one another, thereby enhancing debris flow mobility. Granular temperature arises from conversion of flow translational energy to grain vibrational energy, a process that depends on shear rates, grain properties, boundary conditions, and the ambient fluid viscosity and pressure. Pore fluid pressures that exceed static equilibrium pressures result from local or global debris contraction. Like larger, natural debris flows, experimental debris flows of ˜10 m³ of poorly sorted, water-saturated sediment invariably move as an unsteady surge or series of surges. Measurements at the base of experimental flows show that coarse-grained surge fronts have little or no pore fluid pressure. In contrast, finer-grained, thoroughly saturated debris behind surge fronts is nearly liquefied by high pore pressure, which persists owing to the great compressibility and moderate permeability of the debris. Realistic models of debris flows therefore require equations that simulate inertial motion of surges in which high-resistance fronts dominated by solid forces impede the motion of low-resistance tails more strongly influenced by fluid forces. Furthermore, because debris flows characteristically originate as nearly rigid sediment masses, transform at least partly to liquefied flows, and then transform again to nearly rigid deposits, acceptable models must simulate an evolution of material behavior without invoking preternatural changes in material properties. A simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure. These equations can describe a spectrum of debris flow behaviors intermediate between those of wet rock avalanches and sediment-laden water floods. With appropriate pore pressure distributions the equations yield numerical solutions that successfully predict unsteady, nonuniform motion of experimental debris flows.

The physics of debris flows

USGS Publications Warehouse

Iverson, R.M.

1997-01-01

Recent advances in theory and experimentation motivate a thorough reassessment of the physics of debris flows. Analyses of flows of dry, granular solids and solid-fluid mixtures provide a foundation for a comprehensive debris flow theory, and experiments provide data that reveal the strengths and limitations of theoretical models. Both debris flow materials and dry granular materials can sustain shear stresses while remaining static; both can deform in a slow, tranquil mode characterized by enduring, frictional grain contacts; and both can flow in a more rapid, agitated mode characterized by brief, inelastic grain collisions. In debris flows, however, pore fluid that is highly viscous and nearly incompressible, composed of water with suspended silt and clay, can strongly mediate intergranular friction and collisions. Grain friction, grain collisions, and viscous fluid flow may transfer significant momentum simultaneously. Both the vibrational kinetic energy of solid grains (measured by a quantity termed the granular temperature) and the pressure of the intervening pore fluid facilitate motion of grains past one another, thereby enhancing debris flow mobility. Granular temperature arises from conversion of flow translational energy to grain vibrational energy, a process that depends on shear rates, grain properties, boundary conditions, and the ambient fluid viscosity and pressure. Pore fluid pressures that exceed static equilibrium pressures result from local or global debris contraction. Like larger, natural debris flows, experimental debris flows of ???10 m3 of poorly sorted, water-saturated sediment invariably move as an unsteady surge or series of surges. Measurements at the base of experimental flows show that coarse-grained surge fronts have little or no pore fluid pressure. In contrast, finer-grained, thoroughly saturated debris behind surge fronts is nearly liquefied by high pore pressure, which persists owing to the great compressibility and moderate permeability of the debris. Realistic models of debris flows therefore require equations that simulate inertial motion of surges in which high-resistance fronts dominated by solid forces impede the motion of low-resistance tails more strongly influenced by fluid forces. Furthermore, because debris flows characteristically originate as nearly rigid sediment masses, transform at least partly to liquefied flows, and then transform again to nearly rigid deposits, acceptable models must simulate an evolution of material behavior without invoking preternatural changes in material properties. A simple model that satisfies most of these criteria uses depth-averaged equations of motion patterned after those of the Savage-Hutter theory for gravity-driven flow of dry granular masses but generalized to include the effects of viscous pore fluid with varying pressure. These equations can describe a spectrum of debris flow behaviors intermediate between those of wet rock avalanches and sediment-laden water floods. With appropriate pore pressure distributions the equations yield numerical solutions that successfully predict unsteady, nonuniform motion of experimental debris flows.

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.

Pore water evolution in sandstones of the Groundhog Coalfield, northern Bowser Basin, British Columbia

NASA Astrophysics Data System (ADS)

Cookenboo, H. O.; Bustin, R. M.

1999-01-01

The succession of sandstone cements in chert and volcanic lithic arenites and wackes from the northern Bowser Basin of British Columbia comprises a record of diagenesis in shallow marine, deltaic, and coastal plain siliciclastic sediments that pass through the oil window and reach temperatures near the onset of metamorphism. The succession of cements is consistent with seawater in the sandstones mixing with acid waters derived from dewatering of interbedded organic rich muds. Sandstone cement paragenesis includes seven discrete cement stages. From earliest to latest the cement stages are: (1) pore-lining chlorite; (2) pore-lining to pore-filling illite; (3) pore-filling kaolinite; (4) oil migration through some of the remaining connected pores; (5) chlorite dissolution; (6) quartz cement; and (7) calcite cement. These seven cement stages are interpreted as a record of the evolution of pore waters circulating through the sandstones after burial. The earliest cement stages, as well as the depositional environments, are compatible with seawater as the initial pore fluid. Seawater composition changed during transport through the sandstones, first by loss of Mg 2+ and Fe 2+ during chlorite precipitation (stage 1). Dewatering of interbedded organic-rich mudstones probably added Mg 2+ and Fe 2+ to partially buffer the loss of these cations to chlorite. Acids produced during breakdown of organic matter are presumed to have mixed into sandstone pore fluids due to further compaction of the muds, leading to reduction of initial alkalinity. Reduction in alkalinity, in turn, favours change from chlorite to illite precipitation (stage 2), and finally to kaolinite (stage 3). Pore waters likely reached their peak acidity at the time of oil migration (stage 4). Chlorite dissolution (stage 5) and quartz precipitation (stage 6) occurred when pores were filled by these hydrocarbon-bearing and presumably acidic fluids. Fluid inclusions in fracture-filling quartz cements contain petroleum, high-pressure methane, and methane-rich aqueous solutions. Homogenization temperatures from primary two-phase inclusions are consistent with quartz cementation during progressive heating between approximately 100 and 200°C. Following quartz precipitation, alkaline pore waters were re-established, as evidenced by late-stage calcite cement (stage 7).

Elastic wave velocity and acoustic emission monitoring during Gypsum dehydration under triaxial stress conditions

NASA Astrophysics Data System (ADS)

Brantut, N.; David, E. C.; Héripré, E.; Schubnel, A. J.; Zimmerman, R. W.; Gueguen, Y.

2010-12-01

Dehydration experiments were performed on natural Gypsum polycrystal samples coming from Volterra, Italy in order to study contemporaneously the evolution of P and S elastic wave velocities and acoustic emission (AE) triggering. During these experiments, temperature was slowly raised at 0.15 degrees C per minute under constant stress conditions. Two experiments were realized under quasi-hydrostatic stress (15 and 55 MPa respectively). The third experiment was realized under constant triaxial stress (σ3=45MPa, σ1=75MPa). All three were drained (10MPa constant pore pressure). In each experiments, both P and S wave velocities reduced drastically (as much as approx. 50% in the low confining pressure case) at the onset of dehydration. Importantly, the Vp/Vs ratio also decreased. Shortly after the onset of decrease in P and S wave velocities, the dehydration reaction was also accompanied by bursts of AEs. Time serie locations of the AEs show that they initiated from the pore pressure port, ie from where the pore fluid could easily be drained, and then slowly migrated within the sample. In each experiments, the AE rate could be positively correlated to the reaction rate, inferred from pore volumetry. In such a way, the AE rate reached a peak when the reaction was the fastest. Focal mechanism analysis of the largest AEs showed they had a large volumetric component in compaction, confirming that AEs were indeed related to pore closure and/or collapse. In addition, the AE rate also increased with confinement, ie when a larger amount of compaction was observed. Interestingly, when under differential stress conditions, AE focal mechanisms were mainly in shear. Additional dehydration experiments performed within an environmental scanning electron microscope under low vacuum highlight that, in drained conditions at least, the reaction seems to take place in two phases. First, cracks are being opened along cleavage planes within a single gypsum crystal, which allows for the fluid to escape. Second, the solid volume shrinks and pore collapse can occur. Such a scenario is also consistent with our in-situ analysis under pressure. Finally, a differential effective medium theory approach is used to invert crack density and crack average aspect ratio from elastic wave velocity measurements. Coupling this to Biot-Gassman equation, we can correct for some of the dispersion effects (mainly squirt flow) between the ultrasonic (MHz) and the seismic frequency (Hz) ranges. When doing so, we observe, that, under low confining pressures and in drained conditions at least, the evolution of elastic wave velocities is dominated by the effect due to nucleation of low aspect ratio crack during dehydration. Our results thus seem to point out that, because dehydration reaction are accompanied with crack nucleation, the signature of these reactions in nature, should, in fact, possibly be that of a low Vp/Vs ratio, contrarily to what has been instinctively assumed until now.

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.

Evolution of Micro-Pores in a Single-Crystal Nickel-Based Superalloy During Solution Heat Treatment

NASA Astrophysics Data System (ADS)

Li, Xiangwei; Wang, Li; Dong, Jiasheng; Lou, Langhong; Zhang, Jian

2017-06-01

Evolution of micro-pores in a third-generation single-crystal nickel-based superalloy during solution heat treatment at 1603 K (1330 °C) was investigated by X-ray computed tomography. 3D information including morphology, size, number, and volume fraction of micro-pores formed during solidification (S-pores) and solution (H-pores) was analyzed. The growth behaviors of both S-pores and H-pores can be related to the vacancy formation and diffusion during heat treatment.

The contribution of air-fluidization to the mobility of rapid flowslides involving fine particles

NASA Astrophysics Data System (ADS)

Stilmant, Frédéric; Dewals, Benjamin; Archambeau, Pierre; Erpicum, Sébastien; Pirotton, Michel

2016-04-01

Air-fluidization can be the origin of the long runout of gravitational flows involving fine particles such as ash. An excessive air pore pressure dramatically reduces the friction angle of the material as long as this pressure has not been dissipated, which occurs during the flow. This phenomenon can be modelled thanks to the 2D depth-averaged equations of mass and momentum conservation and an additional transport equation for basal pore pressure evolution (Iverson and Denlinger, 2001). In this contribution, we discuss the application of this model in relation to recent experimental results on air-fluidized flows by Roche et al. (2008) and Roche (2012). The experimental results were used to set a priori the value of the diffusion coefficient in the model, taking into account the difference of scale between the experiments and real-world applications. We also compare the model predictions against detailed observations of a well-documented historical event, the collapse of a fly-ash heap in Belgium (Stilmant et al., 2015). In particular, we analyse the influence of the different components of the model on the results (pore pressure dissipation vs. pore pressure generation). The diffusion coefficient which characterizes the dissipation of air pore pressure is found sufficiently low for maintaining a fluidized flow over hundreds of meters. The study concludes that an air-fluidization theory is consistent with the field observations. These findings are particularly interesting as they seem not in line with the mainstream acceptation in landslide modelling that air generally plays a secondary role (e.g., Legros, 2002). References Iverson, R.M., Denlinger, R.P., 2001. Flow of variably fluidized granular masses across three-dimensional terrain - 1. Coulomb mixture theory. J. Geophys. Res. 106, 537 552. Legros, F., 2002. The mobility of long-runout landslides. Eng. Geol. 63, 301-331. Roche, O., 2012. Depositional processes and gas pore pressure in pyroclastic flows: an experimental perspective. Bull. Volcanol. 74, 1807-1820. Roche, O., Montserrat, S., Niño, Y., Tamburrino, A., 2008. Experimental observations of water-like behavior of initially fluidized, dam break granular flows and their relevance for the propagation of ash-rich pyroclastic flows. J. Geophys. Res. 113, B12203. Stilmant, F., Pirotton, M., Archambeau, P., Erpicum, S., & Dewals, B. (2015). Can the collapse of a fly ash heap develop into an air-fluidized flow? - Reanalysis of the Jupille accident (1961). Geomorphology, 228, 746-755.

Compaction and Permeability Reduction of Castlegate Sandstone under Pore Pressure Cycling

NASA Astrophysics Data System (ADS)

Bauer, S. J.

2014-12-01

We investigate time-dependent compaction and permeability changes by cycling pore pressure with application to compressed air energy storage (CAES) in a reservoir. Preliminary experiments capture the impacts of hydrostatic stress, pore water pressure, pore pressure cycling, chemical, and time-dependent considerations near a borehole in a CAES reservoir analog. CAES involves creating an air bubble in a reservoir. The high pressure bubble serves as a mechanical battery to store potential energy. When there is excess grid energy, bubble pressure is increased by air compression, and when there is energy needed on the grid, stored air pressure is released through turbines to generate electricity. The analog conditions considered are depth ~1 km, overburden stress ~20 MPa and a pore pressure ~10MPa. Pore pressure is cycled daily or more frequently between ~10 MPa and 6 MPa, consistent with operations of a CAES facility at this depth and may continue for operational lifetime (25 years). The rock can vary from initially fully-to-partially saturated. Pore pressure cycling changes the effective stress.Jacketed, room temperature tap water-saturated samples of Castlegate Sandstone are hydrostatically confined (20 MPa) and subjected to a pore pressure resulting in an effective pressure of ~10 MPa. Pore pressure is cycled between 6 to 10 MPa. Sample displacement measurements yielded determinations of volumetric strain and from water flow measurements permeability was determined. Experiments ran for two to four weeks, with 2 to 3 pore pressure cycles per day. The Castlegate is a fluvial high porosity (>20%) primarily quartz sandstone, loosely calcite cemented, containing a small amount of clay.Pore pressure cycling induces compaction (~.1%) and permeability decreases (~20%). The results imply that time-dependent compactive processes are operative. The load path, of increasing and decreasing pore pressure, may facilitate local loosening and grain readjustments that results in the compaction and permeability decreases observed. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Dept. of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.SAND2014-16586A

Earthquakes, fluid pressures and rapid subduction zone metamorphism

NASA Astrophysics Data System (ADS)

Viete, D. R.

2013-12-01

High-pressure/low-temperature (HP/LT) metamorphism is commonly incomplete, meaning that large tracts of rock can remain metastable at blueschist- and eclogite-facies conditions for timescales up to millions of years [1]. When HP/LT metamorphism does take place, it can occur over extremely short durations (<<1 Myr) [1-2]. HP/LT metamorphism must be associated with processes that allow large volumes of rock to remain unaffected over long periods of time, but then suddenly undergo localized metamorphism. Existing models for HP/LT metamorphism have focussed on the role of fluids in providing heat for metamorphism [2] or catalyzing metamorphic reactions [1]. Earthquakes in subduction zone settings can occur to depths of 100s of km. Metamorphic dehydration and the associated development of elevated pore pressures in HP/LT metamorphic rocks has been identified as a cause of earthquake activity at such great depths [3-4]. The process of fracturing/faulting significantly increases rock permeability, causing channelized fluid flow and dissipation of pore pressures [3-4]. Thus, deep subduction zone earthquakes are thought to reflect an evolution in fluid pressure, involving: (1) an initial increase in pore pressure by heating-related dehydration of subduction zone rocks, and (2) rapid relief of pore pressures by faulting and channelized flow. Models for earthquakes at depth in subduction zones have focussed on the in situ effects of dehydration and then sudden escape of fluids from the rock mass following fracturing [3-4]. On the other hand, existing models for rapid and incomplete metamorphism in subduction zones have focussed only on the effects of heating and/or hydration with the arrival of external fluids [1-2]. Significant changes in pressure over very short timescales should result in rapid mineral growth and/or disequilibrium texture development in response to overstepping of mineral reaction boundaries. The repeated process of dehydration-pore pressure development-earthquake-pore pressure relief could conceivably produce a record of episodic HP/LT metamorphism driven by rapid pressure pulses. A new hypothesis is presented for the origins of HP/LT metamorphism: that HP/LT metamorphism is driven by effective pressure pulses caused by localized, earthquake-related modifications to fluid pressures in the subducted slab. In other words, HP/LT metamorphism marks abrupt changes in stress state within the subducted slab, driven by earthquake rupture and fluid flow, and involving a rapid return toward lithostatic pressure from effective pressures well below lithostatic. References: 1. Bjørnerud, MG, Austrheim, H & Lund, MG, 2002. Processes leading to eclogitization (densification) of subducted and tectonically buried crust. Journal of Geophysical Research 107, 2252. 2. Camacho, A, Lee, JKW, Hensen, BJ & Braun, J, 2005. Short-lived orogenic cycles and the eclogitization of cold crust by spasmodic hot fluids. Nature 435, 1191-1196. 3. Green, HW & Houston, H, 1995. The mechanics of deep earthquakes. Annual Reviews of Earth and Planetary Sciences 23, 169-213. 4. Hacker, BR, Peacock, SM, Abers, GA & Holloway, SD, 2003. Subduction factory 2. Are intermediate-depth earthquakes in subducting slabs linked to metamorphic dehydration reactions?. Journal of Geophysical Research 108, 2030.

Numerical Modeling of Initial Slip and Poroelastic Effects of the 2012 Costa Rica Earthquake Using GPS Data

NASA Astrophysics Data System (ADS)

McCormack, K. A.; Hesse, M. A.; Stadler, G.

2015-12-01

Remote sensing and geodetic measurements are providing a new wealth of spatially distributed, time-series data that have the ability to improve our understanding of co-seismic rupture and post-seismic processes in subduction zones. We formulate a Bayesian inverse problem to infer the slip distribution on the plate interface using an elastic finite element model and GPS surface deformation measurements. We present an application to the co-seismic displacement during the 2012 earthquake on the Nicoya Peninsula in Costa Rica, which is uniquely positioned close to the Middle America Trench and directly over the seismogenic zone of the plate interface. The results of our inversion are then used as an initial condition in a coupled poroelastic forward model to investigate the role of poroelastic effects on post-seismic deformation and stress transfer. From this study we identify a horseshoe-shaped rupture area with a maximum slip of approximately 2.5 meters surrounding a locked patch that is likely to release stress in the future. We model the co-seismic pore pressure change as well as the pressure evolution and resulting deformation in the months after the earthquake. The results of the forward model indicate that earthquake-induced pore pressure changes dissipate quickly near the surface, resulting in relaxation of the surface in the seven to ten days following the earthquake. Near the subducting slab interface, pore pressure changes are an order of magnitude larger and may persist for many months after the earthquake.

Quantification of subsurface pore pressure through IODP drilling

NASA Astrophysics Data System (ADS)

Saffer, D. M.; Flemings, P. B.

2010-12-01

It is critical to understand the magnitude and distribution of subsurface pore fluid pressure: it controls effective stress and thus mechanical strength, slope stability, and sediment compaction. Elevated pore pressures also drive fluid flows that serve as agents of mass, solute, and heat fluxes. The Ocean Drilling Program (ODP) and Integrated Ocean Drilling Program (IODP) have provided important avenues to quantify pore pressure in a range of geologic and tectonic settings. These approaches include 1) analysis of continuous downhole logs and shipboard physical properties data to infer compaction state and in situ pressure and stress, 2) laboratory consolidation testing of core samples collected by drilling, 3) direct downhole measurements using pore pressure probes, 3) pore pressure and stress measurements using downhole tools that can be deployed in wide diameter pipe recently acquired for riser drilling, and 4) long-term monitoring of formation pore pressure in sealed boreholes within hydraulically isolated intervals. Here, we summarize key advances in quantification of subsurface pore pressure rooted in scientific drilling, highlighting with examples from subduction zones, the Gulf of Mexico, and the New Jersey continental shelf. At the Nankai, Costa Rican, and Barbados subduction zones, consolidation testing of cores samples, combined with analysis of physical properties data, indicates that even within a few km landward of the trench, pore pressures in and below plate boundary décollement zones reach a significant fraction of the lithostatic load (λ*=0.25-0.91). These results document a viable and quantifiable mechanism to explain the mechanical weakness of subduction décollements, and are corroborated by a small number of direct measurements in sealed boreholes and by inferences from seismic reflection data. Recent downhole measurements conducted during riser drilling using the modular formation dynamics tester wireline tool (MDT) in a forearc basin ~50 km from the trench document hydrostatic pore pressures in the basin fill down to ~1500 mbsf, and illustrate a promising technique for obtaining pore pressure and stress magnitude. In the Gulf of Mexico, we used pore pressure penetrometers to measure severe overpressures (λ*=0.7); a comprehensive program of consolidation testing on recovered core samples confirms elevated pore pressures due to rapid sedimentation, reflecting disequilibrium compaction. Similarly, along the New Jersey continental shelf, analysis of porosity data from downhole logs and augmented by geotechnical testing of cores demonstrates elevated pore pressures in the shallow subsurface. In both offshore New Jersey and the Gulf of Mexico, integration of direct measurements, geotechnical testing, and hydrodynamic modeling illustrate how flow is focused along permeable layers to reduce effective stress and drive submarine landslides. In sum, pore pressure observations made through the ODP and IODP provide insight into how pore pressure controls the large-scale form of passive and active continental margins, how submarine landslides form, and provide strategies for engineering deep boreholes.

Seismic attributes and advanced computer algorithm to predict formation pore pressure: Qalibah formation of Northwest Saudi Arabia

NASA Astrophysics Data System (ADS)

Nour, Abdoulshakour M.

Oil and gas exploration professionals have long recognized the importance of predicting pore pressure before drilling wells. Pre-drill pore pressure estimation not only helps with drilling wells safely but also aids in the determination of formation fluids migration and seal integrity. With respect to the hydrocarbon reservoirs, the appropriate drilling mud weight is directly related to the estimated pore pressure in the formation. If the mud weight is lower than the formation pressure, a blowout may occur, and conversely, if it is higher than the formation pressure, the formation may suffer irreparable damage due to the invasion of drilling fluids into the formation. A simple definition of pore pressure is the pressure of the pore fluids in excess of the hydrostatic pressure. In this thesis, I investigated the utility of advance computer algorithm called Support Vector Machine (SVM) to learn the pattern of high pore pressure regime, using seismic attributes such as Instantaneous phase, t*Attenuation, Cosine of Phase, Vp/Vs ratio, P-Impedance, Reflection Acoustic Impedance, Dominant frequency and one well attribute (Mud-Weigh) as the learning dataset. I applied this technique to the over pressured Qalibah formation of Northwest Saudi Arabia. The results of my research revealed that in the Qalibah formation of Northwest Saudi Arabia, the pore pressure trend can be predicted using SVM with seismic and well attributes as the learning dataset. I was able to show the pore pressure trend at any given point within the geographical extent of the 3D seismic data from which the seismic attributes were derived. In addition, my results surprisingly showed the subtle variation of pressure within the thick succession of shale units of the Qalibah formation.

Heat of adsorption, adsorption stress, and optimal storage of methane in slit and cylindrical carbon pores predicted by classical density functional theory.

PubMed

Hlushak, Stepan

2018-01-03

Temperature, pressure and pore-size dependences of the heat of adsorption, adsorption stress, and adsorption capacity of methane in simple models of slit and cylindrical carbon pores are studied using classical density functional theory (CDFT) and grand-canonical Monte-Carlo (MC) simulation. Studied properties depend nontrivially on the bulk pressure and the size of the pores. Heat of adsorption increases with loading, but only for sufficiently narrow pores. While the increase is advantageous for gas storage applications, it is less significant for cylindrical pores than for slits. Adsorption stress and the average adsorbed fluid density show oscillatory dependence on the pore size and increase with bulk pressure. Slit pores exhibit larger amplitude of oscillations of the normal adsorption stress with pore size increase than cylindrical pores. However, the increase of the magnitude of the adsorption stress with bulk pressure increase is more significant for cylindrical than for slit pores. Adsorption stress appears to be negative for a wide range of pore sizes and external conditions. The pore size dependence of the average delivered density of the gas is analyzed and the optimal pore sizes for storage applications are estimated. The optimal width of slit pore appears to be almost independent of storage pressure at room temperature and pressures above 10 bar. Similarly to the case of slit pores, the optimal radius of cylindrical pores does not exhibit much dependence on the storage pressure above 15 bar. Both optimal width and optimal radii of slit and cylindrical pores increase as the temperature decreases. A comparison of the results of CDFT theory and MC simulations reveals subtle but important differences in the underlying fluid models employed by the approaches. The differences in the high-pressure behaviour between the hard-sphere 2-Yukawa and Lennard-Jones models of methane, employed by the CDFT and MC approaches, respectively, result in an overestimation of the heat of adsorption by the CDFT theory at higher loadings. However, both adsorption stress and adsorption capacity appear to be much less sensitive to the differences between the models and demonstrate excellent agreement between the theory and the computer experiment.

Evolution of a natural debris flow: In situ measurements of flow dynamics, video imagery, and terrestrial laser scanning

USGS Publications Warehouse

McCoy, S.W.; Kean, J.W.; Coe, J.A.; Staley, D.M.; Wasklewicz, T.A.; Tucker, G.E.

2010-01-01

Many theoretical and laboratory studies have been undertaken to understand debris-flow processes and their associated hazards. However, complete and quantitative data sets from natural debris flows needed for confirmation of these results are limited. We used a novel combination of in situ measurements of debris-flow dynamics, video imagery, and pre- and postflow 2-cm-resolution digital terrain models to study a natural debris-flow event. Our field data constrain the initial and final reach morphology and key flow dynamics. The observed event consisted of multiple surges, each with clear variation of flow properties along the length of the surge. Steep, highly resistant, surge fronts of coarse-grained material without measurable pore-fluid pressure were pushed along by relatively fine-grained and water-rich tails that had a wide range of pore-fluid pressures (some two times greater than hydrostatic). Surges with larger nonequilibrium pore-fluid pressures had longer travel distances. A wide range of travel distances from different surges of similar size indicates that dynamic flow properties are of equal or greater importance than channel properties in determining where a particular surge will stop. Progressive vertical accretion of multiple surges generated the total thickness of mapped debris-flow deposits; nevertheless, deposits had massive, vertically unstratified sedimentological textures. ?? 2010 Geological Society of America.

Anisotropy of permeability of reservoir rocks over Miaoli area, NW Taiwan.

NASA Astrophysics Data System (ADS)

Bo-Siang, Xiong; Loung-Yie, Tsai

2013-04-01

The amount of the CO2 has risen since the Industrial Evolution. In order to reduce the amount of CO2 in atmosphere, CO2 sequestration is considered to be the most effective way. In recent years, research about subsurface storage of CO2 into geological formations has increased rapidly. Assessment of storage capability is needed before selecting a site for sequestration. Porosity and permeability are important assessment factors for CO2 sequestration in reservoir rocks. In order to improve the assessment, reservoir rock properties are important and need to be evaluated in advance. Porosity of sandstone is controlled by texture and degree of cementation, whereas permeability is controlled by pore-throat size, pore types and connectivity of pore throat. Sandstones of Miocene to Pleistocene in Miaoli area, NW Taiwan, were collected in this study. YOKO2 porosity/permeability detector is used to measure their permeability perpendicular and parallel to bedding planes under 3 to 60MPa confining pressure with Helium as media. Optical microscope and scanning electron microscope (SEM) were then used to observe the mineral composition, lithology, texture and pore type of sandstones, so as to explore the influence of rock properties on porosity and anisotropy of permeability, as well as the storage potential for CO2 sequestration in the future. The experimental results show that most of the horizontal permeability exceeds the vertical permeability and the anisotropy increases with increasing confining pressure. Mineral composition of sandstones studied were mainly quartz and lithic with little feldspar content. The pore types were mainly primary pores and micropores in this study. The correlation between quantity of macropores and permeability were higher than total porosity and permeability, mainly due to total porosity contains micropores which contribute little to permeability.

A new confined high pressure rotary shear apparatus: preliminary results

NASA Astrophysics Data System (ADS)

Faulkner, D.; Coughlan, G.; Bedford, J. D.

2017-12-01

The frictional properties of fault zone materials, and their evolution during slip, are of paramount importance for determining the earthquake mechanics of large tectonic faults. Friction is a parameter that is difficult to determine from seismological methods so much of our understanding comes from experiment. Rotary shear apparatuses have been widely used in experimental studies to elucidate the frictional properties of faults under realistic earthquake slip velocities (0.1-10 m/s) and displacements (>20 m). However one technical limitation of rotary shear experiments at seismic slip rates has been the lack of confinement. This has led to a limit on the normal stress (due to the strength of the forcing blocks) and also a lack of control of measurements of the pore fluid pressure. Here we present the first preliminary results from a rotary shear apparatus that has been developed to attempt to address this issue. The new fully confined ring shear apparatus has a fast-acting servo-hydraulic confining pressure system of up to 200 MPa and a servo-controlled upstream and downstream pore pressure system of up to 200 MPa. Displacement rates of 0.01μ/s to 2 m/s can be achieved. Fault gouge samples can therefore be sheared at earthquake speed whilst being subject to pressures typically associated with the depth of earthquake nucleation.

Examination of Anisotropy Using Amplitude Variation with Angle and Azimuth (AVAZ) in the Woodford Shale, Anadarko Basin, Oklahoma

NASA Astrophysics Data System (ADS)

Bailey, Austin

Amplitude Variation with Angle and Azimuth (AVAZ) is a method that examines the azimuthal change in seismic amplitude to calculate the anisotropy of a horizontally transverse isotropic (HTI) formation. Anisotropy is generally indicative of heterogeneity in the rock fabric, be it fractures, crack-like pores, or local stress changes. The aim of this study as a whole is to examine the relationship between AVAZ anisotropy magnitude from seismic data and pore pressure gradient from wells. Pore pressure is an important reservoir metric that is often used to understand the production variations within a hydrocarbon reservoir. Predicting pore pressure from seismic data can be extremely useful in not only estimating production, but also in predicting the completion and development strategies that may be most effective. However, seismic-based pore pressure prediction methods have not evolved much in the past decade, with the industry standard to rely on the Bowers (1995) or Eaton (1987) method of converting seismic velocities to pore pressure volumes. These methods may fall short as a predictive tool in many cases, due to their lack of spatial resolution and dependency on a stable velocity model, which may not always be available. Therefore, this study was begun in order to examine if an alternative method of detecting pore pressure variations could be found using AVAZ. The AVAZ methodology was applied to a merged 3D seismic dataset in the Anadarko Basin, Oklahoma provided by Cimarex Energy, in order to examine the Woodford Shale. The Woodford has been a key player in hydrocarbon production from the Anadarko Basin for decades, mainly serving as a source rock until the mid-2000's during the "unconventional revolution''. Throughout its extent, the Woodford Formation shows significant heterogeneity due to both the structure and faults of the basin, as well as changes in the rock fabric. This study aims to use the AVAZ methodology to examine heterogeneity in the Woodford and to relate its anisotropy to pore pressure. Before examining the AVAZ effect in the seismic data, forward modeling from well logs was completed to conceptualize a relationship between pore pressure and anisotropy. Theoretically, at higher pore pressures the reservoir fluid may be effectively propping the fractures open, thus having a greater effect on any pressure wave traveling through the fluid. At lower pore pressure, the overburden pressure dominates the fluid-filled fractures and closes them down. Therefore, at higher pore pressure the AVAZ anisotropy would be greater than at lower pore pressure. The forward modeling from dipole sonic well logs confirms this conceptual model by showing a positive relationship between pore pressure and AVAZ anisotropy. Before the results of the AVAZ workflow were obtained, a variety of pre-processing steps and quality controls were done on the merged 3D seismic dataset. Although the pore pressure - anisotropy relationship appears robust in modeling, the AVAZ results from the seismic data do not appear to correlate with pore pressure. It is likely that acquisition-related artifacts in the seismic data, as well as small magnitude of change in pore pressure, contribute to this lack of correlation. However, further interpretation of the AVAZ volumes shows local stress variations near faults as well as a potential secondary stress trend striking to the north-east. Such information has implications for completion and overall development of the Woodford as an unconventional resource play.

Multi-scale characterization of pore evolution in a combustion metamorphic complex, Hatrurim basin, Israel: Combining (ultra) small-angle neutron scattering and image analysis

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

Wang, Hsiu-Wen; Anovitz, Lawrence; Burg, Avihu

Backscattered scanning electron micrograph and ultra small- and small-angle neutron scattering data have been combined to provide statistically meaningful data on the pore/grain structure and pore evolution of combustion metamorphic complexes from the Hatrurim basin, Israel. Three processes, anti-sintering roughening, alteration of protolith (dehydration, decarbonation, and oxidation) and crystallization of high-temperature minerals, occurred simultaneously, leading to significant changes in observed pore/grain structures. Pore structures in the protoliths, and in lowand high-grade metamorphic rocks show surface (Ds) and mass (Dm) pore fractal geometries with gradual increases in both Ds and Dm values as a function of metamorphic grade. This suggests thatmore » increases in pore volume and formation of less branching pore networks are accompanied by a roughening of pore/grain interfaces. Additionally, pore evolution during combustion metamorphism is also characterized by reduced contributions from small-scale pores to the cumulative porosity in the high-grade rocks. At high temperatures, small-scale pores may be preferentially closed by the formation of high-temperature minerals, producing a rougher morphology with increasing temperature. Alternatively, large-scale pores may develop at the expense of small-scale pores. These observations (pore fractal geometry and cumulative porosity) indicate that the evolution of pore/grain structures is correlated with the growth of high-temperature phases and is a consequence of the energy balance between pore/grain surface energy and energy arising from heterogeneous phase contacts. The apparent pore volume density further suggests that the localized time/temperature development of the high-grade Hatrurim rocks is not simply an extension of that of the low-grade rocks. The former likely represents the "hot spots (burning foci)" in the overall metamorphic terrain while the latter may represent contact aureoles.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

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

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of GPFs including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtration on a singlemore » channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. The microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE PAGES

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla; ...

2018-01-03

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

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

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Characterization of excess pore pressures at the toe of the Nankai accretionary complex, Ocean Drilling Program sites 1173, 1174, and 808: Results of one-dimensional modeling

NASA Astrophysics Data System (ADS)

Gamage, K.; Screaton, E.

2006-04-01

Elevated fluid pore pressures play a critical role in the development of accretionary complexes, including the development of the décollement zone. In this study, we used measured permeabilities of core samples from Ocean Drilling Program (ODP) Leg 190 to develop a permeability-porosity relationship for hemipelagic sediments at the toe of the Nankai accretionary complex. This permeability-porosity relationship was used in a one-dimensional loading and fluid flow model to simulate excess pore pressures and porosities. Simulated excess pore pressure ratios (as a fraction of lithostatic pressure-hydrostatic pressure) using the best fit permeability-porosity relationship were lower than predicted from previous studies. We then tested sensitivity of excess pore pressure ratios in the underthrust sediments to bulk permeability, lateral stress in the prism, and a hypothetical low-permeability barrier at the décollement. Our results demonstrated significant increase in pore pressures below the décollement with lower bulk permeability, such as obtained by using the lower boundary of permeability-porosity data, or when a low-permeability barrier is added at the décollement. In contrast, pore pressures in the underthrust sediments demonstrated less sensitivity to added lateral stresses in the prism, although the profile of the excess pore pressure ratio is affected. Both simulations with lateral stress and a low-permeability barrier at the décollement resulted in sharp increases in porosity at the décollement, similar to that observed in measured porosities. Furthermore, in both scenarios, maximum excess pore pressure ratios were found at the décollement, suggesting that either of these factors would contribute to stable sliding along the décollement.

Proton Gradients as a Key Physical Factor in the Evolution of the Forced Transport Mechanism Across the Lipid Membrane.

PubMed

Strbak, Oliver; Kanuchova, Zuzana; Krafcik, Andrej

2016-11-01

A critical phase in the transition from prebiotic chemistry to biological evolution was apparently an asymmetric ion flow across the lipid membrane. Due to imbalance in the ion flow, the early lipid vesicles could selectively take the necessary molecules from the environment, and release the side-products from the vesicle. Natural proton gradients played a definitively crucial role in this process, since they remain the basis of energy transfer in the present-day cells. On the basis of this supposition, and the premise of the early vesicle membrane's impermeability to protons, we have shown that the emergence of the proton gradient in the lipid vesicle could be a key physical factor in the evolution of the forced transport mechanism (pore formation and active transport) across the lipid bilayer. This driven flow of protons across the membrane is the result of the electrochemical proton gradient and osmotic pressures on the integrity of the lipid vesicle. At a critical number of new lipid molecules incorporated into the vesicle, the energies associated with the creation of the proton gradient exceed the bending stiffness of the lipid membrane, and overlap the free energy of the lipid bilayer pore formation.

Coupled hydrogeological and geomechanical modelling for the analysis of large slope instabilities.

NASA Astrophysics Data System (ADS)

Laloui, Lyesse; Ferrari, Alessio; Bonnard, Christophe

2010-05-01

Slowly-moving landslides (average velocity between 2 and 10 cm/year) are quite frequent in mountainous or hilly areas and they may display occasional crises, generally due to exceptional climatic conditions. The hazard related to these events cannot be analysed in terms of probability analysis, as the number of recorded past events is generally very small and climate changes could significantly modify the environmental setting. Quantitative relationships relating climatic condition fluctuations and sliding area velocity must then be pursued by taking into account the most relevant physical processes involved in the landslide behaviours. Conventional stability analyses are unable to deal with such questions because they do not allow the velocity fields to be determined. With regard to the behaviour of large slope instabilities, a methodology is presented which aims to describe the behaviour of slow-moving landslides by means of a coupled hydrogeological and geomechanical modelling framework. As it is well known, the evolution of the pore water pressure within the landslide body is often recognized as the main cause for the occurrence of displacement accelerations. In this sense the interaction among the hydrological and the mechanical responses must be considered to analyse the landslide behaviour, with the aim of quantitatively relating pore water pressure variations and movements. For a given case study, pore water pressure evolutions in space and time are obtained from a duly calibrated finite element hydrogeological model, which can take into account the role of several key factors such as infiltration, preferential flows and vegetation. Computed groundwater pressures resulting from the hydrogeological simulations are introduced as nodal forces in a finite element geomechanical model in order to calculate stress evolutions and displacements. The use of advanced constitutive models based on the generalised effective stress concept allows taking into account specific behavioural features such as the effects of the changes in the degree of saturation, associated to the fluctuation of the groundwater level. The geomechanical model is calibrated comparing computed and measured displacements in relevant points of the slope. When appropriate, the outcomes from the geomechanical model can be used in an iterative way to update the hydrogeological model settings. In this way it is possible to simulate the evolution of critical factors (such as permeability or retention properties of the involved materials) associated to the cumulated displacements. Once calibrated, the coupled models can be used to assess the landslide behaviour under different scenarios, including modified climatic conditions and the implementation of mitigation measures. Applications to relevant case studies are presented in order to demonstrate the adequacy and the usefulness of the proposed modelling framework.

Poroelastic Response to the 2012 Costa Rica Earthquake and the Effects on Geodetic Surface Deformation and Groundwater Fluxes

NASA Astrophysics Data System (ADS)

McCormack, K. A.; Hesse, M.

2016-12-01

Remote sensing and geodetic measurements are providing a new wealth of spatially distributed, time-series data that have the ability to improve our understanding of co-seismic rupture and post-seismic processes in subduction zones. Following a large earthquake, large-scale deformation is influenced by a myriad of post-seismic processes occurring on different spatial and temporal scales. These include continued slip on the fault plane (after-slip), a poroelastic response due to the movement of over-pressurized groundwater and viscoelastic relaxation of the underlying mantle. Often, the only means of observing these phenomena are through surface deformation measurements - either GPS or InSAR. Such tools measure the combined result of all these processes, which makes studying the effects of any single process difficult. For the 2012 Mw 7.6 Costa Rica Earthquake, we formulate a Bayesian inverse problem to infer the slip distribution on the plate interface using an elastic finite element model and GPS surface deformation measurements. From this study we identify a horseshoe-shaped rupture area surrounding a locked patch that is likely to release stress in the future. The results of our inversion are then used as an initial condition in a coupled poroelastic forward model to investigate the role of poroelastic effects on post-seismic deformation and stress transfer. We model the co-seismic pore pressure change as well as the pressure evolution and resulting deformation in the months after the earthquake. The surface permeability field is constrained by pump-test data from 526 groundwater wells throughout the study area. The results of the forward model indicate that earthquake-induced pore pressure changes dissipate quickly in most areas near the surface, resulting in relaxation of the surface in the seven to twenty days following the earthquake. Near the subducting slab interface, pore pressure changes can be an order of magnitude larger and may persist for many months after the earthquake. Dissipation of earthquake-induced pore pressure in deeper, low permeability areas manifests as surface deformation over a much longer timescale - on the order of months - which may influence the interpretation of longer timescale post-seismic deformation as purely viscoelastic relaxation.

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?

Temperature and Pressure from Collapsing Pores in HMX

NASA Astrophysics Data System (ADS)

Hardin, D. Barrett

2017-06-01

The thermal and mechanical response of collapsing voids in HMX is analyzed. In this work, the focus is simulating the temperature and pressure fields arising from isolated, idealized pores as they collapse in the presence of a shock. HMX slabs are numerically generated which contain a single pore, isolated from the boundaries to remove all wave reflections. In order to understand the primary pore characteristics leading to temperature rise, a series of 2D, plane strain simulations are conducted on HMX slabs containing both cylindrical and elliptical pores of constant size equal to the area of a circular pore with a 1 micron diameter. Each of these pore types is then subjected to shock pressures ranging from a weak shock that is unable to fully collapse the pore to a strong shock which overwhelms the tendency for localization. Results indicate that as shock strength increases, pore collapse phenomenology for a cylindrical pore transitions from a mode dominated by localized melt cracking to an idealized hydrodynamic pore collapse. For the case of elliptical pores, the orientation causing maximum temperature and pressure rise is found. The relative heating in elliptical pores is then quantified as a function of pore orientation and aspect ratio for a pore of a given area. Distribution A: Distribution unlimited. (96TW 2017-0036).

Time-dependent permeability evolution in compacting volcanic fracture systems and implications for gas overpressure

NASA Astrophysics Data System (ADS)

Farquharson, Jamie I.; Wadsworth, Fabian B.; Heap, Michael J.; Baud, Patrick

2017-06-01

Volcanic eruptions are driven by the ascent of volatile-laden magma. The capacity of a volcano system to outgas these volatiles-its permeability-controls the explosive potential, and fractures at volcanic conduit margins play a crucial role in tempering eruption explosivity by acting as outgassing pathways. However, these fractures are often filled with hot volcanic debris that welds and compacts over time, meaning that these permeable pathways have a finite lifetime. While numerous studies emphasize that permeability evolution is important for regulating pressure in shallow volcanic systems, how and when this occurs remains an outstanding question in volcanology. In this contribution, we show that different pressure evolution regimes can be expected across a range of silicic systems as a function of the width and distribution of fractures in the system, the timescales over which they can outgas (a function of depth and temperature), and the permeability of the host material. We define outgassing, diffusive relaxation, and pressure increase regimes, which are distinguished by comparing the characteristic timescales over which they operate. Moreover, we define a critical permeability threshold, which determines (in concert with characteristic timescales of diffusive mass exchange between the pore and melt phases) whether systems fracture and outgas efficiently, or if a volcano will be prone to pressure increases, incomplete healing, and explosive failure.

Numerical simulation on hydromechanical coupling in porous media adopting three-dimensional pore-scale model.

PubMed

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view.

Numerical Simulation on Hydromechanical Coupling in Porous Media Adopting Three-Dimensional Pore-Scale Model

PubMed Central

Liu, Jianjun; Song, Rui; Cui, Mengmeng

2014-01-01

A novel approach of simulating hydromechanical coupling in pore-scale models of porous media is presented in this paper. Parameters of the sandstone samples, such as the stress-strain curve, Poisson's ratio, and permeability under different pore pressure and confining pressure, are tested in laboratory scale. The micro-CT scanner is employed to scan the samples for three-dimensional images, as input to construct the model. Accordingly, four physical models possessing the same pore and rock matrix characteristics as the natural sandstones are developed. Based on the micro-CT images, the three-dimensional finite element models of both rock matrix and pore space are established by MIMICS and ICEM software platform. Navier-Stokes equation and elastic constitutive equation are used as the mathematical model for simulation. A hydromechanical coupling analysis in pore-scale finite element model of porous media is simulated by ANSYS and CFX software. Hereby, permeability of sandstone samples under different pore pressure and confining pressure has been predicted. The simulation results agree well with the benchmark data. Through reproducing its stress state underground, the prediction accuracy of the porous rock permeability in pore-scale simulation is promoted. Consequently, the effects of pore pressure and confining pressure on permeability are revealed from the microscopic view. PMID:24955384

Detecting Pore Fluid Pressure Changes by Using the Vp/Vs Ratio

NASA Astrophysics Data System (ADS)

Vanorio, T.; Mavko, G.

2006-12-01

A central problem in studies aimed at predicting the dynamic behavior of faults is monitoring and quantifying fluid changes in areas prone to overpressure. Experimental and modeling studies show the Vp/Vs ratio to be a good determinant of the saturation state of a rock formation as well as of its inner pore pressure condition. Dectecting pore pressure changes depends, among other causes, on the reliability of laboratory data to calibrate the in-situ measured velocities. Ideally, laboratory experiments performed under controlled conditions would identify the fundamental mechanisms responsible for changes in the measured acoustic properties. However, technical limitations in the laboratory together with the assumptions driving the experimental and modeling approaches rise spouriuos mechanisms which hinder our present understanding of the actual role of high pore pressure on the elastic and poroelastic parameters. Critical issues unclude: a) the frequencies used in the laboratory are responsible for high-frequency fluid effects which induce velocity dispersion. As a result, both the effective stress parameter and velocities (and their pressure-dependence) estimated from high- frequency ultrasonic data are different from those applicable to crustal low frequency wave propagation; b) laboratory measurements made at dry, drained conditions are assumed to mimic those in gas pressured rocks. However, in dry, drained conditions, no pore pressure is exerted in the pore space, and the pore gas is infinitely compressible; c) when using room-dry, drained measurements as the baseline to model pressured rock formations, the unloading path (i.e. decreasing confining pressure) is supposed to mimic the inflationary path due to pore pressure increase. Doing so, it is assumed that the amount of crack opening due to pore pressure is equal to that of crack closure caused by the overburden stress and thus, the effective stress coefficient is implicitely assumed equal to 1. To minimize the assumptions and limitations described above, we designed a laboratory experiment which used gas as pore fluid medium. Experimental results show that in gas-pressured saturated rocks the Vp/Vs ratio, while remaining lower than values reported for liquid saturation conditions, increases with decreasing differential pressure, similarly to the trend observed in liquid saturated rocks.

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.

Protein osmotic pressure gradients and microvascular reflection coefficients.

PubMed

Drake, R E; Dhother, S; Teague, R A; Gabel, J C

1997-08-01

Microvascular membranes are heteroporous, so the mean osmotic reflection coefficient for a microvascular membrane (sigma d) is a function of the reflection coefficient for each pore. Investigators have derived equations for sigma d based on the assumption that the protein osmotic pressure gradient across the membrane (delta II) does not vary from pore to pore. However, for most microvascular membranes, delta II probably does vary from pore to pore. In this study, we derived a new equation for sigma d. According to our equation, pore-to-pore differences in delta II increase the effect of small pores and decrease the effect of large pores on the overall membrane osmotic reflection coefficient. Thus sigma d for a heteroporous membrane may be much higher than previously derived equations indicate. Furthermore, pore-to-pore delta II differences increase the effect of plasma protein osmotic pressure to oppose microvascular fluid filtration.

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.

Quantification of a maximum injection volume of CO2 to avert geomechanical perturbations using a compositional fluid flow reservoir simulator

NASA Astrophysics Data System (ADS)

Jung, Hojung; Singh, Gurpreet; Espinoza, D. Nicolas; Wheeler, Mary F.

2018-02-01

Subsurface CO2 injection and storage alters formation pressure. Changes of pore pressure may result in fault reactivation and hydraulic fracturing if the pressure exceeds the corresponding thresholds. Most simulation models predict such thresholds utilizing relatively homogeneous reservoir rock models and do not account for CO2 dissolution in the brine phase to calculate pore pressure evolution. This study presents an estimation of reservoir capacity in terms of allowable injection volume and rate utilizing the Frio CO2 injection site in the coast of the Gulf of Mexico as a case study. The work includes laboratory core testing, well-logging data analyses, and reservoir numerical simulation. We built a fine-scale reservoir model of the Frio pilot test in our in-house reservoir simulator IPARS (Integrated Parallel Accurate Reservoir Simulator). We first performed history matching of the pressure transient data of the Frio pilot test, and then used this history-matched reservoir model to investigate the effect of the CO2 dissolution into brine and predict the implications of larger CO2 injection volumes. Our simulation results -including CO2 dissolution- exhibited 33% lower pressure build-up relative to the simulation excluding dissolution. Capillary heterogeneity helps spread the CO2 plume and facilitate early breakthrough. Formation expansivity helps alleviate pore pressure build-up. Simulation results suggest that the injection schedule adopted during the actual pilot test very likely did not affect the mechanical integrity of the storage complex. Fault reactivation requires injection volumes of at least about sixty times larger than the actual injected volume at the same injection rate. Hydraulic fracturing necessitates much larger injection rates than the ones used in the Frio pilot test. Tested rock samples exhibit ductile deformation at in-situ effective stresses. Hence, we do not expect an increase of fault permeability in the Frio sand even in the presence of fault reactivation.

HPHT reservoir evolution: a case study from Jade and Judy fields, Central Graben, UK North Sea

NASA Astrophysics Data System (ADS)

di Primio, Rolando; Neumann, Volkmar

2008-09-01

3D basin modelling of a study area in Quadrant 30, UK North Sea was performed in order to elucidate the burial, thermal, pressure and hydrocarbon generation, migration and accumulation history in the Jurassic and Triassic high pressure high temperature sequences. Calibration data, including reservoir temperatures, pressures, petroleum compositional data, vitrinite reflectance profiles and published fluid inclusion data were used to constrain model predictions. The comparison of different pressure generating processes indicated that only when gas generation is taken into account as a pressure generating mechanism, both the predicted present day as well as palaeo-pressure evolution matches the available calibration data. Compositional modelling of hydrocarbon generation, migration and accumulation also reproduced present and palaeo bulk fluid properties such as the reservoir fluid gas to oil ratios. The reconstruction of the filling histories of both reservoirs indicates that both were first charged around 100 Ma ago and contained initially a two-phase system in which gas dominated volumetrically. Upon burial reservoir fluid composition evolved to higher GORs and became undersaturated as a function of increasing pore pressure up to the present day situation. Our results indicate that gas compositions must be taken into account when calculating the volumetric effect of gas generation on overpressure.

Evaluating the Poroelastic Effect on Anisotropic, Organic-Rich, Mudstone Systems

NASA Astrophysics Data System (ADS)

Suarez-Rivera, Roberto; Fjær, Erling

2013-05-01

Understanding the poroelastic effect on anisotropic organic-rich mudstones is of high interest and value for evaluating coupled effects of rock deformation and pore pressure, during drilling, completion and production operations in the oilfield. These applications include modeling and prevention of time-dependent wellbore failure, improved predictions of fracture initiation during hydraulic fracturing operations (Suarez-Rivera et al. Presented at the Canadian Unconventional Resources Conference held in Calgary, Alberta, Canada, 15-17 November 2011. CSUG/SPE 146998 2011), improved understanding of the evolution of pore pressure during basin development, including subsidence and uplift, and the equilibrated effective in situ stress (Charlez, Rock mechanics, vol 2 1997; Katahara and Corrigan, Pressure regimes in sedimentary basins and their prediction: AAPG Memoir, vol 76, pp 73-78 2002; Fjær et al. Petroleum related rock mechanics. 2nd edn 2008). In isotropic rocks, the coupled poro-elastic deformations of the solid framework and the pore fluids are controlled by the Biot and Skempton coefficients. These are the two fundamental properties that relate the rock framework and fluid compressibility and define the magnitude of the poroelastic effect. In transversely isotropic rocks, one desires to understand the variability of these coefficients along the directions parallel and longitudinal to the principal directions of material symmetry (usually the direction of bedding). These types of measurements are complex and uncommon in low-porosity rocks, and particularly problematic and scarce in tight shales. In this paper, we discuss a methodology for evaluating the Biot's coefficient, its variability along the directions parallel and perpendicular to bedding as a function of stress, and the homogenized Skempton coefficient, also as a function of stress. We also predict the pore pressure change that results during undrained compression. Most importantly, we provide values of transverse and longitudinal Biot's coefficients and the homogenized Skempton coefficient for two important North American, gas-producing, organic-rich mudstones. These results could be used for petroleum-related applications.

A new approximation for pore pressure accumulation in marine sediment due to water waves

NASA Astrophysics Data System (ADS)

Jeng, D.-S.; Seymour, B. R.; Li, J.

2007-01-01

The residual mechanism of wave-induced pore water pressure accumulation in marine sediments is re-examined. An analytical approximation is derived using a linear relation for pore pressure generation in cyclic loading, and mistakes in previous solutions (Int. J. Numer. Anal. Methods Geomech. 2001; 25:885-907; J. Offshore Mech. Arctic Eng. (ASME) 1989; 111(1):1-11) are corrected. A numerical scheme is then employed to solve the case with a non-linear relation for pore pressure generation. Both analytical and numerical solutions are verified with experimental data (Laboratory and field investigation of wave-sediment interaction. Joseph H. Defrees Hydraulics Laboratory, School of Civil and Environmental Engineering, Cornell University, Ithaca, NY, 1983), and provide a better prediction of pore pressure accumulation than the previous solution (J. Offshore Mech. Arctic Eng. (ASME) 1989; 111(1):1-11). The parametric study concludes that the pore pressure accumulation and use of full non-linear relation of pore pressure become more important under the following conditions: (1) large wave amplitude, (2) longer wave period, (3) shallow water, (4) shallow soil and (5) softer soils with a low consolidation coefficient. Copyright

Solute mixing regulates heterogeneity of mineral precipitation in porous media: Effect of Solute Mixing on Precipitation

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

Cil, Mehmet B.; Xie, Minwei; Packman, Aaron I.

Synchrotron X-ray microtomography was used to track the spatiotemporal evolution of mineral precipitation and the consequent alteration of the pore structure. Column experiments were conducted by injecting CaCl2 and NaHCO3 solutions into granular porous media either as a premixed supersaturated solution (external mixing) or as separate solutions that mixed within the specimen (internal mixing). The two mixing modes produced distinct mineral growth patterns. While internal mixing promoted transverse heterogeneity with precipitation at the mixing zone, external mixing favored relatively homogeneous precipitation along the flow direction. The impact of precipitation on pore water flow and permeability was assessed via 3-D flowmore » simulations, which indicated anisotropic permeability evolution for both mixing modes. Under both mixing modes, precipitation decreased the median pore size and increased the skewness of the pore size distribution. Such similar pore-scale evolution patterns suggest that the clogging of individual pores depends primarily on local supersaturation state and pore geometry.« less

Pore pressure control on faulting behavior in a block-gouge system

NASA Astrophysics Data System (ADS)

Yang, Z.; Juanes, R.

2016-12-01

Pore fluid pressure in a fault zone can be altered by natural processes (e.g., mineral dehydration and thermal pressurization) and industrial operations involving subsurface fluid injection/extraction for the development of energy and water resources. However, the effect of pore pressure change on the stability and slip motion of a preexisting geologic fault remain poorly understood; yet they are critical for the assessment of seismic risk. In this work, we develop a micromechanical model to investigate the effect of pore pressure on faulting behavior. The model couples pore network fluid flow and mechanics of the solid grains. We conceptualize the fault zone as a gouge layer sandwiched between two blocks; the block material is represented by a group of contact-bonded grains and the gouge is composed of unbonded grains. A pore network is extracted from the particulate pack of the block-gouge system with pore body volumes and pore throat conductivities calculated rigorously based on the geometry of the local pore space. Pore fluid exerts pressure force onto the grains, the motion of which is solved using the discrete element method (DEM). The model updates the pore network regularly in response to deformation of the solid matrix. We study the fault stability in the presence of a pressure inhomogeneity (gradient) across the gouge layer, and compare it with the case of homogeneous pore pressure. We consider both normal and thrust faulting scenarios with a focus on the onset of shear failure along the block-gouge interfaces. Numerical simulations show that the slip behavior is characterized by intermittent dynamics, which is evident in the number of slipping contacts at the block-gouge interfaces and the total kinetic energy of the gouge particles. Numerical results also show that, for the case of pressure inhomogeneity, the onset of slip occurs earlier for the side with higher pressure, and that this onset appears to be controlled by the maximum pressure of both sides of the fault. We conclude that the stability of the fault should be evaluated separately for both sides of the gouge layer, a result that sheds new light on the use of the effective stress principle and the Coulomb failure criterion in evaluating the stability of a complex fault zone.

Preliminary Investigation on the Behavior of Pore Air Pressure During Rainfall Infiltration

NASA Astrophysics Data System (ADS)

Ashraf Mohamad Ismail, Mohd; Min, Ng Soon; Hasliza Hamzah, Nur; Hazreek Zainal Abidin, Mohd; Madun, Aziman; Tajudin, Saiful Azhar Ahmad

2018-04-01

This paper focused on the preliminary investigation of pore air pressure behaviour during rainfall infiltration in order to substantiate the mechanism of rainfall induced slope failure. The actual behaviour or pore air pressure during infiltration is yet to be clearly understood as it is regularly assumed as atmospheric. Numerical modelling of one dimensional (1D) soil column was utilized in this study to provide a preliminary insight of this highlighted uncertainty. Parametric study was performed by using rainfall intensities of 1.85 x 10-3m/s and 1.16 x 10-4m/s applied on glass beads to simulate intense and modest rainfall conditions. Analysis results show that the high rainfall intensity causes more development of pore air pressure compared to low rainfall intensity. This is because at high rainfall intensity, the rainwater cannot replace the pore air smoothly thus confining the pore air. Therefore, the effect of pore air pressure has to be taken into consideration particularly during heavy rainfall.

Fluid pressure development beneath the décollement at the Nankai subduction zone: its implications for slow earthquakes

NASA Astrophysics Data System (ADS)

Hirose, T.; Kamiya, N.; Yamamoto, Y.; Heuer, V.; Inagaki, F.; Kubo, Y.

2017-12-01

Pore fluid pressure along a fault zone is very important for understanding earthquake generation processes in subduction zones. However, quantitative constraints on the pore pressure are quite limited. Here we report two estimates of the pore pressure developed within the underthrust sediments in the Nankai Trough off Cape Muroto, Japan, using the shipboard data obtained during IODP Expedition 370 (Heuer et al., 2017). First estimates are based on the depth trend of porosity data in the lower Shikoku Basin (LSB) facies, in which the décollement zone has propagated. Porosities in the LSB facies generally decrease with depth, but turn to increase by 5-7% below the décollement zone at 760 mbsf. Deeper than 830 mbsf, porosities resume a general compaction trend. By applying the method followed by Screaton et al. (2002) in which the downward porosity-increase is reflected by an excess pore pressure, we estimated the highest excess pore pressure of 4.2 MPa (λ* = 0.4: a ratio of excess pore pressure to effective overburden stress) at 1020 mbsf within the underthrust sediments. Another estimate is based on the analysis of upwelling drilling-mud flow from the borehole, which is a direct evidence the development of overpressure. We assumed that the borehole penetrated a disc-shaped high pore pressure zone with 10 m thickness and the steady-state flow. Then the pore pressure for a given radius of the disc-shaped zone, which is necessary for explaining the observed flow rate, was calculated using Darcy's law. The calculation yields that the pore pressure exceeded by 2-4 MPa above hydrostatic in case of the 10-13 m2 permeability and the 100-1000 m radius of the disc-shaped zone. Our analysis indicates a significant development of excess pore pressure beneath the décollement zone, most likely at the depth of 1020 mbsf where the highest overpressure was estimated from the downhole porosity trend and also an anomaly in relative hydrocarbon gas concentrations. Friction experiments by Sawai et al. (2016) show that a transition from stable to unstable slip behavior appears with increasing pore fluid pressure that is a prerequisite for the generation of slow earthquakes. Thus, slow earthquakes that occurred off Cape Muroto (Obara & Kato, 2016) can be attributed with the observed significant overpressure beneath the décollement.

Parameterization of the Porous-Material Model for Sand with Different Levels of Water Saturation

DTIC Science & Technology

2008-01-01

equation of state defines pressures dependence on mass density and internal-energy density (and in the case of anisotropic materials, on deviatoric ...strain). The strength and failure equations define the evolutions of the deviatoric stress in the elastic regime, elastic–plastic regime, and in the...via the so-called ‘‘pore pressure’’) [6]. Furthermore, the deformation of soil is controlled by the effective stress since the water and gas do not

A heating experiment in the argillites in the Meuse/Haute-Marne underground research laboratory

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

Wileveau, Yannick; Su, Kun; Ghoreychi, Mehdi

2007-07-01

A heating experiment named TER is being conducted with the objectives to identify the thermal properties, as well as to enhance the knowledge on THM processes in the Callovo-Oxfordian clay at the Meuse/Haute Marne Underground Research Laboratory (France). The in situ experiment has being switched on from early 2006. The heater, 3 m length, is designed to inject the power in the undisturbed zone at 6 m from the gallery wall. A heater packer is inflated in a metallic tubing. During the experiment, numerous sensors are emplaced in the surrounding rock and are experienced to monitor the evolution in temperature,more » pore-water pressure and deformation. The models and numerical codes applied should be validated by comparing the modeling results with the measurements. In parallel, some lab testing have been achieved in order to compare the results given with two different scales (cm up to meter scale). In this paper, we present a general description of the TER experiment with installation of the heater equipment and the surrounding instrumentation. Details of the in situ measurements of temperature, pore-pressure and strain evolutions are given for the several heating and cooling phases. The thermal conductivity and some predominant parameters in THM processes (as linear thermal expansion coefficient and permeability) will be discussed. (authors)« less

A depth-averaged debris-flow model that includes the effects of evolving dilatancy. I. physical basis

USGS Publications Warehouse

Iverson, Richard M.; George, David L.

2014-01-01

To simulate debris-flow behaviour from initiation to deposition, we derive a depth-averaged, two-phase model that combines concepts of critical-state soil mechanics, grain-flow mechanics and fluid mechanics. The model's balance equations describe coupled evolution of the solid volume fraction, m, basal pore-fluid pressure, flow thickness and two components of flow velocity. Basal friction is evaluated using a generalized Coulomb rule, and fluid motion is evaluated in a frame of reference that translates with the velocity of the granular phase, vs. Source terms in each of the depth-averaged balance equations account for the influence of the granular dilation rate, defined as the depth integral of ∇⋅vs. Calculation of the dilation rate involves the effects of an elastic compressibility and an inelastic dilatancy angle proportional to m−meq, where meq is the value of m in equilibrium with the ambient stress state and flow rate. Normalization of the model equations shows that predicted debris-flow behaviour depends principally on the initial value of m−meq and on the ratio of two fundamental timescales. One of these timescales governs downslope debris-flow motion, and the other governs pore-pressure relaxation that modifies Coulomb friction and regulates evolution of m. A companion paper presents a suite of model predictions and tests.

Two innovative pore pressure calculation methods for shallow deep-water formations

NASA Astrophysics Data System (ADS)

Deng, Song; Fan, Honghai; Liu, Yuhan; He, Yanfeng; Zhang, Shifeng; Yang, Jing; Fu, Lipei

2017-11-01

There are many geological hazards in shallow formations associated with oil and gas exploration and development in deep-water settings. Abnormal pore pressure can lead to water flow and gas and gas hydrate accumulations, which may affect drilling safety. Therefore, it is of great importance to accurately predict pore pressure in shallow deep-water formations. Experience over previous decades has shown, however, that there are not appropriate pressure calculation methods for these shallow formations. Pore pressure change is reflected closely in log data, particularly for mudstone formations. In this paper, pore pressure calculations for shallow formations are highlighted, and two concrete methods using log data are presented. The first method is modified from an E. Philips test in which a linear-exponential overburden pressure model is used. The second method is a new pore pressure method based on P-wave velocity that accounts for the effect of shallow gas and shallow water flow. Afterwards, the two methods are validated using case studies from two wells in the Yingqiong basin. Calculated results are compared with those obtained by the Eaton method, which demonstrates that the multi-regression method is more suitable for quick prediction of geological hazards in shallow layers.

A characterization of the coupled evolution of grain fabric and pore space using complex networks: Pore connectivity and optimized flows in the presence of shear bands

NASA Astrophysics Data System (ADS)

Russell, Scott; Walker, David M.; Tordesillas, Antoinette

2016-03-01

A framework for the multiscale characterization of the coupled evolution of the solid grain fabric and its associated pore space in dense granular media is developed. In this framework, a pseudo-dual graph transformation of the grain contact network produces a graph of pores which can be readily interpreted as a pore space network. Survivability, a new metric succinctly summarizing the connectivity of the solid grain and pore space networks, measures material robustness. The size distribution and the connectivity of pores can be characterized quantitatively through various network properties. Assortativity characterizes the pore space with respect to the parity of the number of particles enclosing the pore. Multiscale clusters of odd parity versus even parity contact cycles alternate spatially along the shear band: these represent, respectively, local jamming and unjamming regions that continually switch positions in time throughout the failure regime. Optimal paths, established using network shortest paths in favor of large pores, provide clues on preferential paths for interstitial matter transport. In systems with higher rolling resistance at contacts, less tortuous shortest paths thread through larger pores in shear bands. Notably the structural patterns uncovered in the pore space suggest that more robust models of interstitial pore flow through deforming granular systems require a proper consideration of the evolution of in situ shear band and fracture patterns - not just globally, but also inside these localized failure zones.

Thermally-driven Coupled THM Processes in Shales

NASA Astrophysics Data System (ADS)

Rutqvist, J.

2017-12-01

Temperature changes can trigger strongly coupled thermal-hydrological-mechanical (THM) processes in shales that are important to a number of subsurface energy applications, including geologic nuclear waste disposal and hydrocarbon extraction. These coupled processes include (1) direct pore-volume couplings, by thermal expansion of trapped pore-fluid that triggers instantaneous two-way couplings between pore fluid pressure and mechanical deformation, and (2) indirect couplings in terms of property changes, such as changes in mechanical stiffness, strength, and permeability. Direct pore-volume couplings have been studied in situ during borehole heating experiments in shale (or clay stone) formations at Mont Terri and Bure underground research laboratories in Switzerland and France. Typically, the temperature changes are accompanied with a rapid increase in pore pressure followed by a slower decrease towards initial (pre-heating) pore pressure. Coupled THM modeling of these heater tests shows that the pore pressure increases because the thermal expansion coefficient of the fluid is much higher than that of the porous clay stone. Such thermal pressurization induces fluid flow away from the pressurized area towards areas of lower pressure. The rate of pressure increase and magnitude of peak pressure depends on the rate of heating, pore-compressibility, and permeability of the shale. Modeling as well as laboratory experiments have shown that if the pore pressure increase is sufficiently large it could lead to fracturing of the shale or shear slip along pre-existing bedding planes. Another set of data and observations have been collected associated with studies related to concentrated heating and cooling of oil-shales and shale-gas formations. Heating may be used to enhance production from tight oil-shale, whereas thermal stimulation has been attempted for enhanced shale-gas extraction. Laboratory experiments on shale have shown that strength and elastic deformation modulus decreases with temperature while the rate creep deformations increase with temperature. Such temperature dependency also affects the well stability and zonal sealing across shale layers.

Abnormal formation velocities and applications to pore pressure prediction

NASA Astrophysics Data System (ADS)

Liu, Libin; Shen, Guoqiang; Wang, Zhentao; Yang, Hongwei; Han, Hongwei; Cheng, Yuanfeng

2018-06-01

The pore pressure is a vital concept to the petroleum industry and cannot be ignored by either reservoir engineers or geoscientists. Based on theoretical analyses of effective stresses and the grain packing model, a new equation is proposed for predicting pore pressures from formation velocity data. The predictions agree well with both measured pressures and estimations using Eaton's empirical equation, but the application of the new equation to seismic data is simple and convenient. One application example shows that the identification of sweet spots is much easier using pore pressure data than with inverted seismic velocity data. In another application example using field seismic data, a distribution of overpressured strata is revealed, which is a crucial clue for petroleum generation and accumulation. Still, the accuracy of pore pressure prediction is hardly always guaranteed, mainly owing to the complexity of the real geology and the suitability of specific assumptions about the underlying rock physics.

Dynamic three-dimensional pore-scale imaging of reaction in a carbonate at reservoir conditions.

PubMed

Menke, Hannah P; Bijeljic, Branko; Andrew, Matthew G; Blunt, Martin J

2015-04-07

Quantifying CO2 transport and average effective reaction rates in the subsurface is essential to assess the risks associated with underground carbon capture and storage. We use X-ray microtomography to investigate dynamic pore structure evolution in situ at temperatures and pressures representative of underground reservoirs and aquifers. A 4 mm diameter Ketton carbonate core is injected with CO2-saturated brine at 50 °C and 10 MPa while tomographic images are taken at 15 min intervals with a 3.8 μm spatial resolution over a period of 2(1/2) h. An approximate doubling of porosity with only a 3.6% increase in surface area to volume ratio is measured from the images. Pore-scale direct simulation and network modeling on the images quantify an order of magnitude increase in permeability and an appreciable alteration of the velocity field. We study the uniform reaction regime, with dissolution throughout the core. However, at the pore scale, we see variations in the degree of dissolution with an overall reaction rate which is approximately 14 times lower than estimated from batch measurements. This work implies that in heterogeneous rocks, pore-scale transport of reactants limits dissolution and can reduce the average effective reaction rate by an order of magnitude.

Characterization of reactive flow-induced evolution of carbonate rocks using digital core analysis - part 2: Calculation of the evolution of percolation and transport properties.

PubMed

Qajar, Jafar; Arns, Christoph H

2017-09-01

Percolation of reactive fluids in carbonate rocks affects the rock microstructure and hence changes the rock macroscopic properties. In Part 1 paper, we examined the voxel-wise evolution of microstructure of the rock in terms of mineral dissolution/detachment, mineral deposition, and unchanged regions. In the present work, we investigate the relationships between changes in two characteristic transport properties, i.e. permeability and electrical conductivity and two critical parameters of the pore phase, i.e. the fraction of the pore space connecting the inlet and outlet faces of the core sample and the critical pore-throat diameter. We calculate the aforementioned properties on the images of the sample, wherein a homogeneous modification of pore structure occurred in order to ensure the representativeness of the calculated transport properties at the core scale. From images, the evolution of pore connectivity and the potential role of micropores on the connectivity are quantified. It is found that the changing permeability and electrical conductivity distributions along the core length are generally in good agreement with the longitudinal evolution of macro-connected macroporosity and the critical pore-throat diameter. We incorporate microporosity into critical length and permeability calculations and show how microporosity locally plays a role in permeability. It is shown that the Katz-Thompson model reasonably predicts the post-alteration permeability in terms of pre-alteration simulated parameters. This suggests that the evolution of permeability and electrical conductivity of the studied complex carbonate core are controlled by the changes in the macro-connected macroporosity as well as the smallest pore-throats between the connected macropores. Copyright © 2017 Elsevier B.V. All rights reserved.

Elements of an improved model of debris‐flow motion

USGS Publications Warehouse

Iverson, Richard M.

2009-01-01

A new depth‐averaged model of debris‐flow motion describes simultaneous evolution of flow velocity and depth, solid and fluid volume fractions, and pore‐fluid pressure. Non‐hydrostatic pore‐fluid pressure is produced by dilatancy, a state‐dependent property that links the depth‐averaged shear rate and volumetric strain rate of the granular phase. Pore‐pressure changes caused by shearing allow the model to exhibit rate‐dependent flow resistance, despite the fact that the basal shear traction involves only rate‐independent Coulomb friction. An analytical solution of simplified model equations shows that the onset of downslope motion can be accelerated or retarded by pore‐pressure change, contingent on whether dilatancy is positive or negative. A different analytical solution shows that such effects will likely be muted if downslope motion continues long enough, because dilatancy then evolves toward zero, and volume fractions and pore pressure concurrently evolve toward steady states.

Thermal separation of soil particles from thermal conductivity measurement under various air pressures.

PubMed

Lu, Sen; Ren, Tusheng; Lu, Yili; Meng, Ping; Zhang, Jinsong

2017-01-05

The thermal conductivity of dry soils is related closely to air pressure and the contact areas between solid particles. In this study, the thermal conductivity of two-phase soil systems was determined under reduced and increased air pressures. The thermal separation of soil particles, i.e., the characteristic dimension of the pore space (d), was then estimated based on the relationship between soil thermal conductivity and air pressure. Results showed that under both reduced and increased air pressures, d estimations were significantly larger than the geometrical mean separation of solid particles (D), which suggested that conductive heat transfer through solid particles dominated heat transfer in dry soils. The increased air pressure approach gave d values lower than that of the reduced air pressure method. With increasing air pressure, more collisions between gas molecules and solid surface occurred in micro-pores and intra-aggregate pores due to the reduction of mean free path of air molecules. Compared to the reduced air pressure approach, the increased air pressure approach expressed more micro-pore structure attributes in heat transfer. We concluded that measuring thermal conductivity under increased air pressure procedures gave better-quality d values, and improved soil micro-pore structure estimation.

Relationship between compressional-wave velocity and porosity of sediments along subduction plate interface

NASA Astrophysics Data System (ADS)

Yamaguchi, M.; Hashimoto, Y.

2012-12-01

Evolution of physical properties of sediments along subduction interface has effects on wedge strength, wedge geometry, dewatering and dehydration processes, and seismic behavior. Sediments have initially more than 70% of porosity prior to subduction. Through underthrusting and accretion, porosity of sediments decreases by compaction and cementation to be lithified sediments. The purpose of this study is to understand evolution of physical properties from a state before subduction to a state within a wedge using a relationship between compressional-wave velocity and porosity. In this study, we obtained new data for sediments from a reference site in IODP NanTroSEIZE, Expedition 333. In addition to that, we have complied velocity-porosity relationships for the samples and also for previous studies from NanTroSEIZE (off Kumano) (Hashimoto et al., 2010, 2011), ODP Leg 190 (off Shikoku) (Hoffman and Tobin, 2004) and ODP Leg 170 (off Costa Rica) (Gettemy and Tobin, 2003). Velocity measurement procedure in this study to obtain new data is as following: Two pumps were used to control pore fluid pressure and confining pressure. The pore pressure of 1000kPa was kept under drained conditions. Confining (effective) pressure was increased stepwise in the measurements. Velocity measurements were conducted under isotropic pressure conditions. Confining pressure was pressurized in tens seconds and kept for more than 8 hours for next step to obtain equilibrium conditions between effective pressure and sediments strain. Lead zirconate titanate (PZT) shear wave transducers (500kHz) were used in a source-receiver pair to measure wave speed. Porosity and P-wave velocity ranges about 27 - 75% and 1.4 - 2.2 km/s in this study, respectively. In the comparison in Vp-porosity relationships between sedimetns from reference sites and others, sediments were classified into two, simply compacted sediments (reference site and slope sediments) and wedge sediments. Different trends in Vp-porosity relationships were observed for the classified sediments. For compacted sediments, Vp-porosity relationships are along the global empirical relationships (Erickson and Jarrard 1988) and almost within the area between normal and highly compaction curves. On the other hand, some of Vp-porosity relationships for wedge sediments represent trends with higher velocity at a porosity. Such trend was observed for wedge sediments from Site C0001 and C0004. Those higher Vp trend in Vp-porosity relationship for wedge sediments can be explained by shear strain of sediments and/or cementation. Even though the velocity measurements was conducted under hydrostatic condition, we examined the void ratio-porosity curve as a kind of compaction curve. On the basis of the curves, break points were observed at the pressure which corresponds to the effective pressure assuming the hydrostatic pore fluid pressure. The result suggests that the sediments were under condition of normal compaction. Some of void ratio-porosity curve represent a evidence of weak cement which can correspond with anomaly in porosity-depth curve in the shallow portion of the reference sites.

Capillary pressure-saturation relationships for porous granular materials: Pore morphology method vs. pore unit assembly method

NASA Astrophysics Data System (ADS)

Sweijen, Thomas; Aslannejad, Hamed; Hassanizadeh, S. Majid

2017-09-01

In studies of two-phase flow in complex porous media it is often desirable to have an estimation of the capillary pressure-saturation curve prior to measurements. Therefore, we compare in this research the capability of three pore-scale approaches in reproducing experimentally measured capillary pressure-saturation curves. To do so, we have generated 12 packings of spheres that are representative of four different glass-bead packings and eight different sand packings, for which we have found experimental data on the capillary pressure-saturation curve in the literature. In generating the packings, we matched the particle size distributions and porosity values of the granular materials. We have used three different pore-scale approaches for generating the capillary pressure-saturation curves of each packing: i) the Pore Unit Assembly (PUA) method in combination with the Mayer and Stowe-Princen (MS-P) approximation for estimating the entry pressures of pore throats, ii) the PUA method in combination with the hemisphere approximation, and iii) the Pore Morphology Method (PMM) in combination with the hemisphere approximation. The three approaches were also used to produce capillary pressure-saturation curves for the coating layer of paper, used in inkjet printing. Curves for such layers are extremely difficult to determine experimentally, due to their very small thickness and the presence of extremely small pores (less than one micrometer in size). Results indicate that the PMM and PUA-hemisphere method give similar capillary pressure-saturation curves, because both methods rely on a hemisphere to represent the air-water interface. The ability of the hemisphere approximation and the MS-P approximation to reproduce correct capillary pressure seems to depend on the type of particle size distribution, with the hemisphere approximation working well for narrowly distributed granular materials.

Extrusion of transmitter, water and ions generates forces to close fusion pore.

PubMed

Tajparast, M; Glavinović, M I

2009-05-01

During exocytosis the fusion pore opens rapidly, then dilates gradually, and may subsequently close completely, but what controls its dynamics is not well understood. In this study we focus our attention on forces acting on the pore wall, and which are generated solely by the passage of transmitter, ions and water through the open fusion pore. The transport through the charged cylindrical nano-size pore is simulated using a coupled system of Poisson-Nernst-Planck and Navier-Stokes equations and the forces that act radially on the wall of the fusion pore are then estimated. Four forces are considered: a) inertial force, b) pressure, c) viscotic force, and d) electrostatic force. The inertial and viscotic forces are small, but the electrostatic force and the pressure are typically significant. High vesicular pressure tends to open the fusion pore, but the pressure induced by the transport of charged particles (glutamate, ions), which is predominant when the pore wall charge density is high tends to close the pore. The electrostatic force, which also depends on the charge density on the pore wall, is weakly repulsive before the pore dilates, but becomes attractive and pronounced as the pore dilates. Given that the vesicular concentration of free transmitter can change rapidly due to the release, or owing to the dissociation from the gel matrix, we evaluated how much and how rapidly a change of the vesicular K(+)-glutamate(-) concentration affects the concentration of glutamate(-) and ions in the pore and how such changes alter the radial force on the wall of the fusion pore. A step-like rise of the vesicular K(+)-glutamate(-) concentration leads to a chain of events. Pore concentration (and efflux) of both K(+) and glutamate(-) rise reaching their new steady-state values in less than 100 ns. Interestingly within a similar time interval the pore concentration of Na(+) also rises, whereas that of Cl(-) diminishes, although their extra-cellular concentration does not change. Finally such changes affect also the water movement. Water efflux changes bi-phasically, first increasing before decreasing to a new, but lower steady-state value. Nevertheless, even under such conditions an overall approximate neutrality of the pore is maintained remarkably well, and the electrostatic, but also inertial, viscotic and pressure forces acting on the pore wall remain constant. In conclusion the extrusion of the vesicular content generates forces, primarily the force due to the electro-kinetically induced pressure and electrostatic force (both influenced by the pore radius and even more by the charge density on the pore wall), which tend to close the fusion pore.

Analysis on the overpressure characterization with respect to depositional environment facies: Case Study in Miri Area, Baram Delta

NASA Astrophysics Data System (ADS)

Mhd Hanapiah, N.; Yusoff, W. I. Wan; Zakariah, M. N. A.

2017-10-01

Overpressure studies in oil and gas exploration and production are carried out in order to mitigate any losses that could happen while drilling. These concerns can be addressed by enhancing the understanding of overpressure characterization in the fields. This research emphasizes in determining the pore pressure trend in Miri area to assist pore pressure prediction for future hydrocarbon exploration and production. Generally, pore pressure trends are related to mechanisms that contribute to the overpressure generation. In the region predominant overpressure are disequilibrium compaction within the prodelta shales meanwhile in outer shelf overpressure generation controlled by fluid expansion in deltaic sequence of inner shelf area. The objective of this research is to analyze the pore pressure profile of wells for determining vertical trends of pore pressure for various depositional environment facies of Miri area. Integration of rock physics and pore pressure analysis and relating the trends to environment depositional environment facies within shale underlying sand interval. Analysis done shows that overpressure top is characterize by depositional environment facies within shale underlying sand interval.

Condensation in Nanoporous Packed Beds.

PubMed

Ally, Javed; Molla, Shahnawaz; Mostowfi, Farshid

2016-05-10

In materials with tiny, nanometer-scale pores, liquid condensation is shifted from the bulk saturation pressure observed at larger scales. This effect is called capillary condensation and can block pores, which has major consequences in hydrocarbon production, as well as in fuel cells, catalysis, and powder adhesion. In this study, high pressure nanofluidic condensation studies are performed using propane and carbon dioxide in a colloidal crystal packed bed. Direct visualization allows the extent of condensation to be observed, as well as inference of the pore geometry from Bragg diffraction. We show experimentally that capillary condensation depends on pore geometry and wettability because these factors determine the shape of the menisci that coalesce when pore filling occurs, contrary to the typical assumption that all pore structures can be modeled as cylindrical and perfectly wetting. We also observe capillary condensation at higher pressures than has been done previously, which is important because many applications involving this phenomenon occur well above atmospheric pressure, and there is little, if any, experimental validation of capillary condensation at such pressures, particularly with direct visualization.

Regulation of landslide motion by dilatancy and pore pressure feedback

USGS Publications Warehouse

Iverson, R.M.

2005-01-01

A new mathematical model clarifies how diverse styles and rates of landslide motion can result from regulation of Coulomb friction by dilation or contraction of water-saturated basal shear zones. Normalization of the model equations shows that feedback due to coupling between landslide motion, shear zone volume change, and pore pressure change depends on a single dimensionless parameter ??, which, in turn, depends on the dilatancy angle ?? and the intrinsic timescales for pore pressure generation and dissipation. If shear zone soil contracts during slope failure, then ?? 0, and negative feedback permits slow, steady landslide motion to occur while positive pore pressure is supplied by rain infiltration. Steady state slip velocities v0 obey v0 = -(K/??) p*e, where K is the hydraulic conductivity and p*e is the normalized (dimensionless) negative pore pressure generated by dilation. If rain infiltration and attendant pore pressure growth continue unabated, however, their influence ultimately overwhelms the stabilizing influence of negative p*e. Then, unbounded landslide acceleration occurs, accentuated by an instability that develops if ?? diminishes as landslide motion proceeds. Nonetheless, numerical solutions of the model equations show that slow, nearly steady motion of a clay-rich landslide may persist for many months as a result of negative pore pressure feedback that regulates basal Coulomb friction. Similarly stabilized motion is less likely to occur in sand-rich landslides that are characterized by weaker negative feedback.

Characterization of pore structure in cement-based materials using pressurization-depressurization cycling mercury intrusion porosimetry (PDC-MIP)

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

Zhou Jian, E-mail: Jian.Zhou@tudelft.n; Ye Guang, E-mail: g.ye@tudelft.n; Magnel Laboratory for Concrete Research, Department of Structural Engineering, Ghent University, Technologiepark-Zwijnaarde 904 B-9052, Ghent

2010-07-15

Numerous mercury intrusion porosimetry (MIP) studies have been carried out to investigate the pore structure in cement-based materials. However, the standard MIP often results in an underestimation of large pores and an overestimation of small pores because of its intrinsic limitation. In this paper, an innovative MIP method is developed in order to provide a more accurate estimation of pore size distribution. The new MIP measurements are conducted following a unique mercury intrusion procedure, in which the applied pressure is increased from the minimum to the maximum by repeating pressurization-depressurization cycles instead of a continuous pressurization followed by a continuousmore » depressurization. Accordingly, this method is called pressurization-depressurization cycling MIP (PDC-MIP). By following the PDC-MIP testing sequence, the volumes of the throat pores and the corresponding ink-bottle pores can be determined at every pore size. These values are used to calculate pore size distribution by using the newly developed analysis method. This paper presents an application of PDC-MIP on the investigation of the pore size distribution in cement-based materials. The experimental results of PDC-MIP are compared with those measured by standard MIP. The PDC-MIP is further validated with the other experimental methods and numerical tool, including nitrogen sorption, backscanning electron (BSE) image analysis, Wood's metal intrusion porosimetry (WMIP) and the numerical simulation by the cement hydration model HYMOSTRUC3D.« less

Numerical analysis of pressure and porosity evolution in lava domes during periodic degassing conditions

NASA Astrophysics Data System (ADS)

Hyman, D.; Bursik, M. I.; Pitman, E. B.

2017-12-01

The collapse or explosive breakup of growing and degassing lava domes presents a significant hazard due to the generation of dense, mobile pyroclastic flows as well as the wide dispersal of dense ballistic blocks. Lava dome stability is in large part governed by the balance of transport and storage of gas within the pore space. Because pore pressurization reduces the effective stress within a dome, the transient distribution of elevated gas pressure is critically important to understanding dome break up. We combine mathematical and numerical analyses to gain a better understanding of the temporal variation in gas flow and storage within the dome system. In doing so, we develop and analyze new governing equations describing nonlinear gas pressure diffusion in a deforming dome with an evolving porosity field. By relating porosity, permeability, and pressure, we show that the flux of gas through a dome is highly sensitive to the porosity distribution and viscosity of the lava, as well as the timescale and magnitude of the gas supply. The numerical results suggest that the diffusion of pressure and porosity variations play an integral role in the cyclic growth and destruction of small domes.The nearly continuous cycles of lava dome growth, pressurization, and failure that have characterized the last two decades of eruptive history at Volcán Popocatépetl, Mexico provide excellent natural data with which to compare new models of transient dome pressurization. At Popocatépetl, periodic pressure increases brought on by changes in gas supply into the base of the dome may play a role in its cyclic growth and destruction behavior. We compare our model of cyclic pressurization with lava dome survival data from Popocatépetl. We show that transient changes in pore pressure explain how small lava domes evolve to a state of criticality before explosion or collapse. Additionally, numerical analyses presented here suggest that short-term oscillations cannot arise within the dome, and must be the result of an oscillating supply of gas into the dome. The oscillating gas supply may result from alternating gas-rich and gas-poor regions of rising magma, so-called "porosity waves" within the conduit. These internal pressure fluctuations lead to periodic reductions in the stress required to fracture the dome and induce explosion.

Earthquake-Induced Liquefaction of Confined Soil Zones: A Centrifuge Study

DTIC Science & Technology

1993-12-09

8217. 4.2 Pore pressure transducers (PPT) Pore pressures in the saturated soil were monitored by Druck PDCR 81 pore pressure transducers. This type of pore...PPT5406 -20 MTn-23.3 2 3 4 5 6 7 a 2 10 II 12 Mo, 34 .8 PPT6270 20~ Mlm=-I .2 1 2 3 4 5 6 8 a 1 0 II 12 4 M....7 0 PPT6263 A 3d Moa20.S PPT6260 7 -0 1 2 3

Numerical Modeling of the 2014 Oso, Washington, Landslide.

NASA Astrophysics Data System (ADS)

George, D. L.; Iverson, R. M.

2014-12-01

Numerical simulations of alternative scenarios that could have transpired during the Oso, Washington, landslide of 22 March 2014 provide insight into factors responsible for the landslide's devastating high-speed runout.We performed these simulations using D-Claw, a numerical model we recently developed to simulate landslide and debris-flow motion from initiation to deposition. D-Claw solves a hyperbolic system of five partial differential equations that describe simultaneous evolution of the thickness,solid volume fraction, basal pore-fluid pressure, and two components of momentum of the moving mass. D-Claw embodies the concept ofstate-dependent dilatancy, which causes the solid volume fraction m to evolve toward a value that is equilibrated to the ambient stress state andshear rate. As the value of m evolves, basal pore-fluid pressure coevolves,and thereby causes an evolution in frictional resistance to motion. Our Oso simulations considered alternative scenarios in which values of all model parameters except the initial solid volume fraction m0 were held constant.These values were: basal friction angle = 36°; static critical-state solidvolume fraction = 0.64; initial sediment permeability = 10-8 m2; pore-fluid density = 1100 kg/m3; sediment grain density = 2700 kg/m3; pore-fluid viscosity = 0.005 Pa-s; and dimensionless sediment compressibility coefficient = 0.03. Simulations performed using these values and m0 = 0.62 produced widespread landslide liquefaction, runaway acceleration, andlandslide runout distances, patterns and speeds similar to those observed or inferred for the devastating Oso event. Alternative simulations that usedm0 = 0.64 produced a much slower landslide that did not liquefy and that traveled only about 100 m before stopping. This relatively benign behavioris similar to that of several landslides at the Oso site prior to 2014. Our findings illustrate a behavioral bifurcation that is highly sensitive to the initial solid volume fraction. They suggest that the destructiveness of the2014 Oso event may have resulted in part from prior slope deformation that produced a dilated sediment state that made the sediment susceptible to contraction and liquefaction as it began to fail on March 22.

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

Cook, J.M.; Sheppard, M.C.; Houwen, O.H.

Previous work on shale mechanical properties has focused on the slow deformation rates appropriate to wellbore deformation. Deformation of shale under a drill bit occurs at a very high rate, and the failure properties of the rock under these conditions are crucial in determining bit performance and in extracting lithology and pore-pressure information from drilling parameters. Triaxial tests were performed on two nonswelling shales under a wide range of strain rates and confining and pore pressures. At low strain rates, when fluid is relatively free to move within the shale, shale deformation and failure are governed by effective stress ormore » pressure (i.e., total confining pressure minus pore pressure), as is the case for ordinary rock. If the pore pressure in the shale is high, increasing the strain rate beyond about 0.1%/sec causes large increases in the strength and ductility of the shale. Total pressure begins to influence the strength. At high stain rates, the influence of effective pressure decreases, except when it is very low (i.e., when pore pressure is very high); ductility then rises rapidly. This behavior is opposite that expected in ordinary rocks. This paper briefly discusses the reasons for these phenomena and their impact on wellbore and drilling problems.« less

Upscaling pore pressure-dependent gas permeability in shales

NASA Astrophysics Data System (ADS)

Ghanbarian, Behzad; Javadpour, Farzam

2017-04-01

Upscaling pore pressure dependence of shale gas permeability is of great importance and interest in the investigation of gas production in unconventional reservoirs. In this study, we apply the Effective Medium Approximation, an upscaling technique from statistical physics, and modify the Doyen model for unconventional rocks. We develop an upscaling model to estimate the pore pressure-dependent gas permeability from pore throat size distribution, pore connectivity, tortuosity, porosity, and gas characteristics. We compare our adapted model with six data sets: three experiments, one pore-network model, and two lattice-Boltzmann simulations. Results showed that the proposed model estimated the gas permeability within a factor of 3 of the measurements/simulations in all data sets except the Eagle Ford experiment for which we discuss plausible sources of discrepancies.

Impoundment of the Zipingpu reservoir and triggering of the 2008 Mw 7.9 Wenchuan earthquake, China

PubMed Central

Tao, Wei; Masterlark, Timothy; Ronchin, Erika

2015-01-01

Abstract Impoundment of the Zipingpu reservoir (ZR), China, began in September 2005 and was followed 2.7 years later by the 2008 Mw 7.9 Wenchuan earthquake (WE) rupturing the Longmen Shan Fault (LSF), with its epicenter ~12 km away from the ZR. Based on the poroelastic theory, we employ three‐dimensional finite element models to simulate the evolution of stress and pore pressure due to reservoir impoundment, and its effect on the Coulomb failure stress on the LSF. The results indicate that the reservoir impoundment formed a pore pressure front that slowly propagated through the crust with fluid diffusion. The reservoir loading induced either moderate or no increase of the Coulomb failure stress at the hypocenter prior to the WE. The Coulomb failure stress, however, grew ~9.3–69.1 kPa in the depth range of 1–8 km on the LSF, which may have advanced tectonic loading of the fault system by ~60–450 years. Due to uncertainties of fault geometry and hypocenter location of the WE, it is inconclusive whether impoundment of the ZR directly triggered the WE. However, a small event at the hypocenter could have triggered large rupture elsewhere on fault, where the asperities were weakened by the ZR. The microseismicity around the ZR also showed an expanding pattern from the ZR since its impoundment, likely associated with diffusion of a positive pore pressure pulse. These results suggest a poroelastic triggering effect (even if indirectly) of the WE due to the impoundment of the ZR. PMID:27812436

Diagenesis and porosity evolution of tight sand reservoirs in Carboniferous Benxi Formation, Southeast Ordos Basin

NASA Astrophysics Data System (ADS)

Hu, Peng; Yu, Xinghe; Shan, Xin; Su, Dongxu; Wang, Jiao; Li, Yalong; Shi, Xin; Xu, Liqiang

2016-04-01

The Ordos Basin, situated in west-central China, is one of the oldest and most important fossil-fuel energy base, which contains large reserves of coal, oil and natural gas. The Upper Palaeozoic strata are widely distributed with rich gas-bearing and large natural gas resources, whose potential is tremendous. Recent years have witnessed a great tight gas exploration improvement of the Upper Paleozoic in Southeastern Ordos basin. The Carboniferous Benxi Formation, mainly buried more than 2,500m, is the key target strata for hydrocarbon exploration, which was deposited in a barrier island and tidal flat environment. The sandy bars and flats are the favorable sedimentary microfacies. With an integrated approach of thin-section petrophysics, constant velocity mercury injection test, scanning electron microscopy and X-ray diffractometry, diagenesis and porosity evolution of tight sand reservoirs of Benxi Formation were analyzed in detail. The result shows that the main lithology of sandstone in this area is dominated by moderately to well sorted quartz sandstone. The average porosity and permeability is 4.72% and 1.22mD. The reservoirs of Benxi Formation holds a variety of pore types and the pore throats, with obvious heterogeneity and poor connection. Based on the capillary pressure curve morphological characteristics and parameters, combined with thin section and phycical property data, the reservoir pore structure of Benxi Formation can be divided into 4 types, including mid pore mid throat type(I), mid pore fine throat type(II), small pore fine throat type(III) and micro pro micro throat type(Ⅳ). The reservoirs primarily fall in B-subsate of middle diagenesis and late diagenesis, which mainly undergo compaction, cmentation, dissolution and fracturing process. Employing the empirical formula of different sorting for unconsolideated sandstone porosity, the initial sandstone porosity is 38.32% on average. Quantitative evaluation of the increase and decrease of porosity caused by different diagenesis reveals that mechanical compaction and chemical cementation are the main mechanisms for destroying primary pores, which contribute 19.61% and 8.75% to the loss of primary posoity, respectively. Dissolution of volcanic fragments and feldspar increased reservoir porosity by 4.14%. The pores were occluded by late minerals and carbonate cements, resulting in a reduction of 9.38%. Overall, the dual influence of compaction and cementation is the key of the key, controlling formation of tight gas sandstone reservoirs. Keywords: diagenesis, porosity evolution, tight sandstone, Benxi Formation, Southeast Ordos Basin Acknowledgements: We greatfully acknowledge Yanchang Petroleum for providing the samples and data access and for permission to publish this work. The first author, Peng Hu, would like to thank the support from Prof. Xinghe Yu.

The influence of temperature on brittle creep in sandstones

NASA Astrophysics Data System (ADS)

Heap, M. J.; Baud, P.; Meredith, P. G.; Vinciguerra, S.

2009-04-01

The characterization of time-dependent brittle rock deformation is fundamental to understanding the long-term evolution and dynamics of the Earth's upper crust. The presence of water promotes time-dependent deformation through environment-dependent stress corrosion cracking that allows rocks to deform at stresses far below their short-term failure stress. Here we report results from an experimental study of the influence of an elevated temperature on time-dependent brittle creep in water-saturated samples of Darley Dale (initial porosity of 13%), Bentheim (23%) and Crab Orchard (4%) sandstones. We present results from both conventional creep experiments (or ‘static fatigue' tests) and stress-stepping creep experiments performed under 20°C and 75°C and an effective confining pressure of 30 MPa (50 MPa confining pressure and a 20 MPa pore fluid pressure). The evolution of crack damage was monitored throughout each experiment by measuring the three proxies for damage (1) axial strain (2) pore volume change and (3) the output of AE energy. Conventional creep experiments have demonstrated that, for any given applied differential stress, the time-to-failure is dramatically reduced and the creep strain rate is significantly increased by application of an elevated temperature. Stress-stepping creep experiments have allowed us to investigate the influence of temperature in detail. Results from these experiments show that the creep strain rate for Darley Dale and Bentheim sandstones increases by approximately 3 orders of magnitude, and for Crab Orchard sandstone increases by approximately 2 orders of magnitude, as temperature is increased from 20°C to 75°C at a fixed effective differential stress. We discuss these results in the context of the different mineralogical and microstructural properties of the three rock types and the micro-mechanical and chemical processes operating on them.

Pore-Scale X-ray Micro-CT Imaging and Analysis of Oil Shales

NASA Astrophysics Data System (ADS)

Saif, T.

2015-12-01

The pore structure and the connectivity of the pore space during the pyrolysis of oil shales are important characteristics which determine hydrocarbon flow behaviour and ultimate recovery. We study the effect of temperature on the evolution of pore space and subsequent permeability on five oil shale samples: (1) Vernal Utah United States, (2) El Lajjun Al Karak Jordan, (3) Gladstone Queensland Australia (4) Fushun China and (5) Kimmerdige United Kingdom. Oil Shale cores of 5mm in diameter were pyrolized at 300, 400 and 500 °C. 3D imaging of 5mm diameter core samples was performed at 1μm voxel resolution using X-ray micro computed tomography (CT) and the evolution of the pore structures were characterized. The experimental results indicate that the thermal decomposition of kerogen at high temperatures is a major factor causing micro-scale changes in the internal structure of oil shales. At the early stage of pyrolysis, micron-scale heterogeneous pores were formed and with a further increase in temperature, the pores expanded and became interconnected by fractures. Permeability for each oil shale sample at each temperature was computed by simulation directly on the image voxels and by pore network extraction and simulation. Future work will investigate different samples and pursue insitu micro-CT imaging of oil shale pyrolysis to characterize the time evolution of the pore space.

Temporal and Spatial Pore Water Pressure Distribution Surrounding a Vertical Landfill Leachate Recirculation Well

PubMed Central

Kadambala, Ravi; Townsend, Timothy G.; Jain, Pradeep; Singh, Karamjit

2011-01-01

Addition of liquids into landfilled waste can result in an increase in pore water pressure, and this in turn may increase concerns with respect to geotechnical stability of the landfilled waste mass. While the impact of vertical well leachate recirculation on landfill pore water pressures has been mathematically modeled, measurements of these systems in operating landfills have not been reported. Pressure readings from vibrating wire piezometers placed in the waste surrounding a liquids addition well at a full-scale operating landfill in Florida were recorded over a 2-year period. Prior to the addition of liquids, measured pore pressures were found to increase with landfill depth, an indication of gas pressure increase and decreasing waste permeability with depth. When liquid addition commenced, piezometers located closer to either the leachate injection well or the landfill surface responded more rapidly to leachate addition relative to those far from the well and those at deeper locations. After liquid addition stopped, measured pore pressures did not immediately drop, but slowly decreased with time. Despite the large pressures present at the bottom of the liquid addition well, much smaller pressures were measured in the surrounding waste. The spatial variation of the pressures recorded in this study suggests that waste permeability is anisotropic and decreases with depth. PMID:21655145

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.

2017-01-01

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix. PMID:28772465

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N

2017-01-25

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix.

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).

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.; Gao, Shengyan

2015-01-01

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir. PMID:26310236

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N; Gao, Shengyan

2015-08-27

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir.

Entrainment of bed sediment by debris flows: results from large-scale experiments

USGS Publications Warehouse

Reid, Mark E.; Iverson, Richard M.; Logan, Matthew; LaHusen, Richard G.; Godt, Jonathan W.; Griswold, Julie P.

2011-01-01

When debris flows grow by entraining sediment, they can become especially hazardous owing to increased volume, speed, and runout. To investigate the entrainment process, we conducted eight largescale experiments in the USGS debris-flow flume. In each experiment, we released a 6 m3 water-saturated debris flow across a 47-m long, ~12-cm thick bed of partially saturated sediment lining the 31º flume. Prior to release, we used low-intensity overhead sprinkling and real-time monitoring to control the bed-sediment wetness. As each debris flow descended the flume, we measured the evolution of flow thickness, basal total normal stress, basal pore-fluid pressure, and sediment scour depth. When debris flows traveled over relatively dry sediment, net scour was minimal, but when debris flows traveled over wetter sediment (volumetric water content > 0.22), debris-flow volume grew rapidly and flow speed and runout were enhanced. Data from scour sensors showed that entrainment occurred by rapid (5-10 cm/s), progressive scour rather than by mass failure at depth. Overriding debris flows rapidly generated high basal pore-fluid pressures when they loaded and deformed bed sediment, and in wetter beds these pressures approached lithostatic levels. Reduction of intergranular friction within the bed sediment thereby enhanced scour efficiency, entrainment, and runout.

Study of pore pressure reaction on hydraulic fracturing

NASA Astrophysics Data System (ADS)

Trimonova, Mariia; Baryshnikov, Nikolay; Turuntaev, Sergey; Zenchenko, Evgeniy; Zenchenko, Petr

2017-04-01

We represent the results of the experimental study of the hydraulic fracture propagation influence on the fluid pore pressure. Initial pore pressure was induced by injection and production wells. The experiments were carried out according to scaling analysis based on the radial model of the fracture. All required geomechanical and hydrodynamical properties of a sample were derived from the scaling laws. So, gypsum was chosen as a sample material and vacuum oil as a fracturing fluid. The laboratory setup allows us to investigate the samples of cylindrical shape. It can be considered as an advantage in comparison with standard cubic samples, because we shouldn't consider the stress field inhomogeneity induced by the corners. Moreover, we can set 3D-loading by this setting. Also the sample diameter is big enough (43cm) for placing several wells: the fracturing well in the center and injection and production wells on two opposite sides of the central well. The experiment consisted of several stages: a) applying the horizontal pressure; b) applying the vertical pressure; c) water solution injection in the injection well with a constant pressure; d) the steady state obtaining; e) the oil injection in the central well with a constant rate. The pore pressure was recorded in the 15 points along bottom side of the sample during the whole experiment. We observe the pore pressure change during all the time of the experiment. First, the pore pressure changed due to water injection. Then we began to inject oil in the central well. We compared the obtained experimental data on the pore pressure changes with the solution of the 2D single-phase equation of pore-elasticity, and we found significant difference. The variation of the equation parameters couldn't help to resolve the discrepancy. After the experiment, we found that oil penetrated into the sample before and after the fracture initiation. This fact encouraged us to consider another physical process - the oil-water displacement. Have taken into account the phenomenon, we could find the parameter values for the best matching the experimental data with the analytical one. After such a comparison, we could estimate the permeability variation in the different directions due to changes in the pore pressure during fracturing. Thus it was found that for the correct solution of hydrodynamic problems in relation with hydraulic fracturing (for example, to estimate the production rate of the fractured well) one should take into account the change of the permeability in the vicinity of the fracture and solve nonlinear pore-elasticity problem.

Analogue modelling of caprock failure and sediment mobilisation due to pore fluid overpressure in shallow reservoirs

NASA Astrophysics Data System (ADS)

Warsitzka, Michael; Kukowski, Nina; May, Franz

2017-04-01

Injection of CO2 in geological formations may cause excess pore fluid pressure by enhancing the fluid volume in the reservoir rock and by buoyancy-driven flow. If sediments in the reservoir and the caprock are undercompacted, pore fluid overpressure can lead to hydro-fractures in the caprock and fluidisation of sediments. Eventually, these processes trigger the formation of pipe structures, gas chimneys, gas domes or sand injections. Generally, such structures serve as high permeable pathways for fluid migration through a low-permeable seal layer and have to be considered in risk assessment or modelling of caprock integrity of CO2 storage sites. We applied scaled analogue experiments to characterise and quantify mechanisms determining the onset and migration of hydro-fractures in a low-permeable, cohesive caprock and fluidisation of unconsolidated sediments of the reservoir layer. The caprock is simulated by different types of cohesive powder. The reservoir layer consists of granulates with small particle density. Air injected through the base of the experiment and additionally through a single needle valve reaching into the analogue material is applied to generate fluid pressure within the materials. With this procedure, regional fluid pressure increase or a point-like local fluid pressure increase (e.g. injection well), respectively, can be simulated. The deformation in the analogue materials is analysed with a particle tracking imaging velocimetry technique. Pressure sensors at the base of the experiment and in the needle valve record the air pressure during an experimental run. The structural evolution observed in the experiments reveal that the cohesive cap rock first forms a dome-like anticline. Extensional fractures occur at the hinges of the anticline. A further increase of fluid pressure causes a migration of this fractures towards the surface, which is followed by intrusion of reservoir material into the fractures and the collapse of the anticline. The breakthrough of the fractures at the surface is accompanied by a significant drop of air pressure at the base of the analogue materials. The width of the dome shaped uplift is narrower and the initiating fluid pressure in the needle valve is lower, if the fluid pressure at the base of the experiment is larger. The experimental outcomes help to evaluate if the injection of CO2 into a reservoir potentially provokes initiation or reactivation of fractures and sediment mobilisation structures.

Time evolution of pore system in lime - Pozzolana composites

NASA Astrophysics Data System (ADS)

Doleželová, Magdaléna; Čáchová, Monika; Scheinherrová, Lenka; Keppert, Martin

2017-11-01

The lime - pozzolana mortars and plasters are used in restoration works on building cultural heritage but these materials are also following the trend of energy - efficient solutions in civil engineering. Porosity and pore size distribution is one of crucial parameters influencing engineering properties of porous materials. The pore size distribution of lime based system is changing in time due to chemical processes occurring in the material. The present paper describes time evolution of pore system in lime - pozzolana composites; the obtained results are useful in prediction of performance of lime - pozzolana systems in building structures.

The change in orientation of subsidiary shears near faults containing pore fluid under high pressure

USGS Publications Warehouse

Byerlee, J.

1992-01-01

Byerlee, J., 1992. The change in orientation of subsidiary shears near faults containing pore fluid under high pressure. In: T. Mikumo, K. Aki, M. Ohnaka, L.J. Ruff and P.K.P. Spudich (Editors), Earthquake Source Physics and Earthquake Precursors. Tectonophysics, 211: 295-303. The mechanical effects of a fault containing near-lithostatic fluid pressure in which fluid pressure decreases monotonically from the core of the fault zone to the adjacent country rock is considered. This fluid pressure distribution has mechanical implications for the orientation of subsidiary shears around a fault. Analysis shows that the maximum principal stress is oriented at a high angle to the fault in the country rock where the pore pressure is hydrostatic, and rotates to 45?? to the fault within the fault zone where the pore pressure is much higher. This analysis suggests that on the San Andreas fault, where heat flow constraints require that the coefficient of friction for slip on the fault be less than 0.1, the pore fluid pressure on the main fault is 85% of the lithostatic pressure. The observed geometry of the subsidiary shears in the creeping section of the San Andreas are broadly consistent with this model, with differences that may be due to the heterogeneous nature of the fault. ?? 1992.

Electrohydraulic shock wave generation as a means to increase intrinsic permeability of mortar

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

Maurel, O.; Reess, T.; Matallah, M.

2010-12-15

This article discusses the influence of compressive shock waves on the permeability of cementitious materials. Shock waves are generated in water by Pulsed Arc Electrohydraulic Discharges (PAED). The practical aim is to increase the intrinsic permeability of the specimens. The maximum pressure amplitude of the shock wave is 250 MPa. It generates damage in the specimens and the evolution of damage is correlated with the intrinsic permeability of the mortar. A threshold of pressure is observed. From this threshold, the increase of permeability is linear in a semi-log plot. The influence of repeated shocks on permeability is also discussed. Qualitativemore » X Ray Tomography illustrates the evolution of the microstructure of the material leading to the increase of permeability. Comparative results from mercury intrusion porosimetry (MIP) show that the micro-structural damage process starts at the sub-micrometric level and that the characteristic size of pores of growing volume increases.« less

Modelling the influence of pore size on the response of materials to infrared lasers An application to human enamel

NASA Astrophysics Data System (ADS)

Vila Verde, A.; Ramos, Marta M. D.

2005-07-01

We present an analytical model for a ceramic material (hydroxyapatite, HA) containing nanometre-scale water pores, and use it to estimate the pressure at the pore as a function of temperature at the end of a single 0.35 μs laser pulse by Er:YAG (2.94 μm) and CO 2 (10.6 μm) lasers. Our results suggest that the pressure at the pore is directly related to pore temperature, and that very high pressures can be generated simply by the thermal expansion of liquid water. Since the temperature reached in the pores at the end of the laser pulse is a strong function of pore size for Er:YAG lasers, but is independent of pore size for CO 2 lasers, our present results provide a possible explanation for the fact that human dental enamel threshold ablation fluences vary more for Er:YAG lasers than for CO 2 lasers. This suggests that experimentalists should analyse their results accounting for factors, like age or type of tooth, that may change the pore size distribution in their samples.

Evolution of Abnormally Low Pressure at Bravo Dome and its Implications for Carbon Capture and Storage (CCS)

NASA Astrophysics Data System (ADS)

Akhbari, D.; Hesse, M. A.

2015-12-01

Carbon capture and storage allows reductions of the rapidly rising CO2 from fossil fuel-based power generation, if large storage rates and capacities can be achieved. The injection of large fluid volumes at high rates leads to a build-up of pore-pressure in the storage formation that may induce seismicity and compromise the storage security. Many natural CO2 fields in midcontinent US, in contrast, are under-pressured rather than over-pressured suggesting that natural processes reduce initial over-pressures and generate significant under-pressures. The question is therefore to understand the sequence of process(es) that allow the initial over-pressure to be eliminated and the under-pressure to be maintained over geological periods of time. We therefore look into pressure evolution in Bravo Dome, one of the largest natural CO2 accumulations in North America, which stores 1.3 Gt of CO2. Bravo Dome is only 580-900 m deep and is divided into several compartments with near gas-static pressure (see Figure). The pre-production gas pressures in the two main compartments that account for 70% of the mass of CO2 stored at Bravo Dome are more than 6 MPa below hydrostatic pressure. Here we show that the under-pressure in the Bravo Dome CO2 reservoir is maintained by hydrological compartmentalization over millennial timescales and generated by a combination of processes including cooling, erosional unloading, limited leakage into overlying formations, and CO2 dissolution into brine. Herein, we introduce CO2 dissolution into brine as a new process that reduce gas pressure in a compartmentalized reservoir and our results suggest that it may contribute significantly to reduce the initial pressure build-up due to injection. Bravo Dome is the first documented case of pressure drop due to CO2 dissolution. To have an accurate prediction of pressure evolution in Bravo Dome, our models must include geomechanics and thermodynamics for the reservoir while they account for the pressure changes due to the CO2 dissolution.

Effect of confinement in nano-porous materials on the solubility of a supercritical gas

NASA Astrophysics Data System (ADS)

Hu, Yaofeng; Huang, Liangliang; Zhao, Shuangliang; Liu, Honglai; Gubbins, Keith E.

2016-11-01

By combining Gibbs Ensemble Monte Carlo simulations and density functional theory, we investigate the influence of confinement in a slit-shaped carbon pore on the solubility of a supercritical solute gas in a liquid solvent. In the cases studied here, competing adsorption of the solvent and solute determines whether the solubility is enhanced or suppressed for larger pores. We find that the solubility in the confined system is strongly dependent on pore width, and that molecular packing effects are important for small pore widths. In addition, the solubility decreases on increase in the temperature, as for the bulk mixture, but the rate of decrease is greater in the pore due to a decrease in the partial molar enthalpy of the solute in the pore; this effect becomes greater as pore width is decreased. The solubility is increased on increasing the bulk pressure of the gas in equilibrium with the pore, and obeys Henry's law at lower pressures. However, the Henry constant differs significantly from that for the bulk mixture, and the range of pressure over which Henry's law applies is reduced relative to that for the bulk mixture. The latter observation indicates that solute-solute interactions become more important in the pore than for the bulk at a given bulk pressure. Finally, we note that different authors use different definitions of the solubility in pores, leading to some confusion over the reported phenomenon of 'oversolubility'. We recommend that solubility be defined as the overall mole fraction of solute in the pores, since it takes into account the increase in density of the solvent in the pores, and avoids ambiguity in the definition of the pore volume.

The effects of pressure, temperature, and pore water on velocities in Westerly granite. [for seismic wave propagation

NASA Technical Reports Server (NTRS)

Spencer, J. W., Jr.; Nur, A. M.

1976-01-01

A description is presented of an experimental assembly which has been developed to conduct concurrent measurements of compressional and shear wave velocities in rocks at high temperatures and confining pressures and with independent control of the pore pressure. The apparatus was used in studies of the joint effects of temperature, external confining pressure, and internal pore water on sonic velocities in Westerly granite. It was found that at a given temperature, confining pressure has a larger accelerating effect on compressional waves in dry rock, whereas at a given confining pressure, temperature has a larger retarding effect on shear waves.

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.

Measurements of unjacketed moduli of porous rock

NASA Astrophysics Data System (ADS)

Tarokh, A.; Makhnenko, R. Y.; Labuz, J.

2017-12-01

Coupling of stress and pore pressure appears in a number of applications dealing with subsurface (sedimentary) rock, including petroleum exploration and waste storage. Poroelastic analyses consider the compressibility of the solid constituents forming the rock, and often times solid bulk modulus Ks is assumed to be the same as the dominant mineral bulk modulus. In fact, there are two different parameters describing solid compressibility of a porous rock: the unjacketed bulk modulus Ks' and the unjacketed pore modulus Ks". Experimental techniques are developed to measure the two poroelastic parameters of fluid-saturated porous rock under the unjacketed condition. In an unjacketed experiment, the rock without a membrane is loaded by the fluid in a pressure vessel. The confining fluid permeates the connected pore space throughout the interior of the rock. Therefore, changes in mean stress P will produce equal changes in pore pressure p, i.e. ΔP = Δp. The test can also be performed with a jacketed rock specimen by applying equal increments of mean stress and pore pressure. The unjacketed bulk modulus, Ks', is obtained by measuring the bulk strain with resistive strain gages. The unjacketed pore modulus, Ks", the pore volume counterpart to Ks', is a measure of the change in pore pressure per unit pore volume strain under the unjacketed condition. Several indirect estimates of Ks" have been reported but limitations of these approaches do not provide an accurate value. We present direct measurements of Ks" with detailed calibration on the system volumetric response. The results indicate that for Dunnville sandstone Ks' and Ks" are equal while for Berea sandstone, a difference between the two moduli exists, which is explained by the presence of non-connected pores. The experiments also strongly suggest that both Ks' and Ks" are independent of effective stress.

Coupled Electro-Hydrodynamic Effects of Electro-Osmosis from Pore Scale to Darcy Scale

NASA Astrophysics Data System (ADS)

Schotting, R.; Joekar-Niasar, V.; Leijnse, A.

2011-12-01

Electro-osmosis is "movement of a fluid under the effect of an electric field in a porous medium". This phenomenon has many applications in civil engineering (slope stabilization, dewatering), environmental engineering (soil remediation, sludge dewatering), chemical engineering (micro- or nano- mixers), medical engineering (drug delivery), etc. The key factor in electro-osmosis is the competition between the electrochemical and hydrodynamic forces as well as the coupling between the solid surface and the electrolyte properties. The objective of this research is to understand the influence of pore-scale heterogeneities of surface properties on the Darcy-scale behavior. We develop novel analytical solutions for the flow and transport of electrolyte including electro-hydrodynamic forces in a single micro-channel. We propose the complete analytical solution for monovalent electrolyte at full range overlapping double layers, and nonlinear electric field, including the Donan effect in transport of ions. These pore-scale formulations are numerically upscaled to obtain the Darcy-scale behavior. Our results show the contribution of electro-osmotic, chemical-osmotic and hydrodynamic components of the flow equation on pressure field evolution and multi-directional flow field at Darcy scale.

Virial series expansion and Monte Carlo studies of equation of state for hard spheres in narrow cylindrical pores

NASA Astrophysics Data System (ADS)

Mon, K. K.

2018-05-01

In this paper, the virial series expansion and constant pressure Monte Carlo method are used to study the longitudinal pressure equation of state for hard spheres in narrow cylindrical pores. We invoke dimensional reduction and map the model into an effective one-dimensional fluid model with interacting internal degrees of freedom. The one-dimensional model is extensive. The Euler relation holds, and longitudinal pressure can be probed with the standard virial series expansion method. Virial coefficients B2 and B3 were obtained analytically, and numerical quadrature was used for B4. A range of narrow pore widths (2 Rp) , Rp<(√{3 }+2 ) /4 =0.9330 ... (in units of the hard sphere diameter) was used, corresponding to fluids in the important single-file formations. We have also computed the virial pressure series coefficients B2', B3', and B4' to compare a truncated virial pressure series equation of state with accurate constant pressure Monte Carlo data. We find very good agreement for a wide range of pressures for narrow pores. These results contribute toward increasing the rather limited understanding of virial coefficients and the equation of state of hard sphere fluids in narrow cylindrical pores.

An evaluation of factors influencing pore pressure in accretionary complexes: Implications for taper angle and wedge mechanics

USGS Publications Warehouse

Saffer, D.M.; Bekins, B.A.

2006-01-01

At many subduction zones, accretionary complexes form as sediment is off-scraped from the subducting plate. Mechanical models that treat accretionary complexes as critically tapered wedges of sediment demonstrate that pore pressure controls their taper angle by modifying basal and internal shear strength. Here, we combine a numerical model of groundwater flow with critical taper theory to quantify the effects of sediment and de??collement permeability, sediment thickness, sediment partitioning between accretion and underthrusting, and plate convergence rate on steady state pore pressure. Our results show that pore pressure in accretionary wedges can be viewed as a dynamically maintained response to factors which drive pore pressure (source terms) and those that limit flow (permeability and drainage path length). We find that sediment permeability and incoming sediment thickness are the most important factors, whereas fault permeability and the partitioning of sediment have a small effect. For our base case model scenario, as sediment permeability is increased, pore pressure decreases from near-lithostatic to hydrostatic values and allows stable taper angles to increase from ??? 2.5?? to 8??-12.5??. With increased sediment thickness in our models (from 100 to 8000 m), increased pore pressure drives a decrease in stable taper angle from 8.4??-12.5?? to 15?? to <4??) with increased sediment thickness (from <1 to 7 km). One key implication is that hydrologic properties may strongly influence the strength of the crust in a wide range of geologic settings. Copyright 2006 by the American Geophysical Union.

Preparation of high porosity xerogels by chemical surface modification.

DOEpatents

Deshpande, Ravindra; Smith, Douglas M.; Brinker, C. Jeffrey

1996-01-01

This invention provides an extremely porous xerogel dried at vacuum-to-below supercritical pressures but having the properties of aerogels which are typically dried at supercritical pressures. This is done by reacting the internal pore surface of the wet gel with organic substances in order to change the contact angle of the fluid meniscus in the pores during drying. Shrinkage of the gel (which is normally prevented by use of high autoclave pressures, such that the pore fluid is at temperature and pressure above its critical values) is avoided even at vacuum or ambient pressures.

Long Term Observations of Subsurface Pore Pressure in the Kumano Basin and Upper Accretionary Wedge along the NanTroSIEZE Transect, offshore Japan: Signals from the 2011 Tohoku Earthquake

NASA Astrophysics Data System (ADS)

Zhang, Y.; Saffer, D. M.

2013-12-01

Subsurface pore pressure as a sensitive measure of strain and formation properties has provided insights into the wide range of fault slip behaviors, contributing to the understanding of fault and earthquake mechanics. Pore pressures from off shore borehole observatory are especially important, as 1) they are the only detectable signals of small and slow events; 2) they provide our only access to the outer forearc, where the tsunami hazards are triggered by the fault slip. As part of the Nankai Trough Seismogenic Zone Experiment (NanTroSEIZE) a suite of borehole sensors were installed as part of a long-term borehole observatory at IODP Site C0002, during IODP Expedition # 332 in December of 2010. The observatory includes a broadband seismometer, short period geophones, a volumetric strainmeter, temperature sensors, an accelerometer, and formation pore pressure monitoring at two depths: one in the mudstones of the Kumano Basin in an interval spanning 757-780 meters below seafloor (mbsf), and a second in the uppermost accretionary wedge in an interval from 937 - 980 mbsf. Here, we report on pore pressure records acquired at a sampling frequency of 1/60 Hz, spanning the period from December 2010 to January 2013, which were recovered in early 2013. We observe a clear hydraulic signal from March 11, 2011 Tohoku earthquake and aftershocks, including both dynamic pore pressure changes during passage of surface waves and shifts in formation pressure following the event. Pressure exhibit an increase of ~3 kPa in the upper sediment screened interval following the earthquake, and decrease by ~5 kPa in the accretionary prism interval. Both of the offset changes persist through the end of the data recording. These pore pressure changes may reflect static stress changes from the earthquake, or local site effects related to shaking. We also observe a clear increase in formation pore pressures associated with drilling operations at nearby holes in November and December 2012. These inadvertent two-well tests provide information about formation hydraulic properties at the ~20-50 m scale.

Experimental study of dynamic effective stress coefficient for ultrasonic velocities of Bakken cores

NASA Astrophysics Data System (ADS)

Ma, X.; Zoback, M. D.

2016-12-01

We have performed a series of exhaustive experiments to measure the effective stress coefficient (α) of the tight cores from the Bakken shale oil play. Five distinct, bedding-normal cores from a vertical well were tested, covering the sequences of Lodgepole, Middle Bakken, and Three Forks. The scope of this laboratory study is two-fold: (1) to obtain the dynamic effective stress coefficient for ultrasonic velocities; (2) to characterize the poromechanical properties in relation to rock's mineral composition and microstructure. The experiments were carried out as follows: Argon-saturated specimen (1-inch length, 1-inch diameter) was subjected to hydrostatic confining pressure under drained conditions. Pore pressure was regulated as Argon was injected into both ends of the specimen. We drilled multiple non-through-going boreholes (1-mm diameter) in the specimen to facilitate pore pressure equilibrium, without compromising its integrity. The specimen was put through a loading path to experience confining pressure and pore pressure up to 70 and 60 MPa, respectively. P- and S- wave velocities were measured and used to calculate the rock's dynamic effective stress coefficient. Results of all five cores unanimously show that the dynamic a is a function of both confining and pore pressures, regardless of the wave type and loading path. When the simple effective stress is low, α is close to unity; however, α consistently increases as the simple effective stress rises and can reach as much as 3 when the latter reaches 60 MPa. This trend is rather surprising as it is diametrically the opposite of what was observed for the static α. A possible explanation is that high-frequency wave-induced pore pressure increment may have not remained equilibrated throughout the pore space, especially in very thin cracks, according to the squirt model. This phenomenon can be enhanced when the bulk modulus of pore fluid (gas typically considered to be `soft' and `non-viscous') increases with pore pressure and becomes comparable to the crack stiffness.

Modeling postshock evolution of large electropores

NASA Astrophysics Data System (ADS)

Neu, John C.; Krassowska, Wanda

2003-02-01

The Smoluchowski equation (SE), which describes the evolution of pores created by electric shocks, cannot be applied to modeling large and long-lived pores for two reasons: (1) it does not predict pores of radius above 20 nm without also predicting membrane rupture; (2) it does not predict postshock growth of pores. This study proposes a model in which pores are coupled by membrane tension, resulting in a nonlinear generalization of SE. The predictions of the model are explored using examples of homogeneous (all pore radii r are equal) and heterogeneous (0⩽r⩽rmax) distributions of pores. Pores in a homogeneous population either shrink to zero or assume a stable radius corresponding to the minimum of the bilayer energy. For a heterogeneous population, such a stable radius does not exist. All pores, except rmax, shrink to zero and rmax grows to infinity. However, the unbounded growth of rmax is not physical because the number of pores per cell decreases in time and the continuum model loses validity. When the continuum formulation is replaced by the discrete one, the model predicts the coarsening process: all pores, except rmax, shrink to zero and rmax assumes a stable radius. Thus, the model with tension-coupled pores does not predict membrane rupture and the predicted postshock growth of pores is consistent with experimental evidence.

Evolution of the pore structure during the early stages of the alkali-activation reaction: An in situ small-angle neutron scattering investigation

DOE PAGES

White, Claire E.; Olds, Daniel P.; Hartl, Monika; ...

2017-02-01

The long-term durability of cement-based materials is influenced by the pore structure and associated permeability at the sub-micrometre length scale. With the emergence of new types of sustainable cements in recent decades, there is a pressing need to be able to predict the durability of these new materials, and therefore nondestructive experimental techniques capable of characterizing the evolution of the pore structure are increasingly crucial for investigating cement durability. Here, small-angle neutron scattering is used to analyze the evolution of the pore structure in alkali-activated materials over the initial 24 h of reaction in order to assess the characteristic poremore » sizes that emerge during these short time scales. By using a unified fitting approach for data modeling, information on the pore size and surface roughness is obtained for a variety of precursor chemistries and morphologies (metakaolin- and slag-based pastes). Furthermore, the impact of activator chemistry is elucidatedviathe analysis of pastes synthesized using hydroxide- and silicate-based activators. It is found that the main aspect influencing the size of pores that are accessible using small-angle neutron scattering analysis (approximately 10–500 Å in diameter) is the availability of free silica in the activating solution, which leads to a more refined pore structure with smaller average pore size. Furthermore, as the reaction progresses the gel pores visible using this scattering technique are seen to increase in size.« less

Distributed optical fibre sensing for early detection of shallow landslides triggering.

PubMed

Schenato, Luca; Palmieri, Luca; Camporese, Matteo; Bersan, Silvia; Cola, Simonetta; Pasuto, Alessandro; Galtarossa, Andrea; Salandin, Paolo; Simonini, Paolo

2017-10-31

A distributed optical fibre sensing system is used to measure landslide-induced strains on an optical fibre buried in a large scale physical model of a slope. The fibre sensing cable is deployed at the predefined failure surface and interrogated by means of optical frequency domain reflectometry. The strain evolution is measured with centimetre spatial resolution until the occurrence of the slope failure. Standard legacy sensors measuring soil moisture and pore water pressure are installed at different depths and positions along the slope for comparison and validation. The evolution of the strain field is related to landslide dynamics with unprecedented resolution and insight. In fact, the results of the experiment clearly identify several phases within the evolution of the landslide and show that optical fibres can detect precursory signs of failure well before the collapse, paving the way for the development of more effective early warning systems.

Overpressures in the Uinta Basin, Utah: Analysis using a three-dimensional basin evolution model

NASA Astrophysics Data System (ADS)

McPherson, Brian J. O. L.; Bredehoeft, John D.

2001-04-01

High pore fluid pressures, approaching lithostatic, are observed in the deepest sections of the Uinta basin, Utah. Geologic observations and previous modeling studies suggest that the most likely cause of observed overpressures is hydrocarbon generation. We studied Uinta overpressures by developing and applying a three-dimensional, numerical model of the evolution of the basin. The model was developed from a public domain computer code, with addition of a new mesh generator that builds the basin through time, coupling the structural, thermal, and hydrodynamic evolution. Also included in the model are in situ hydrocarbon generation and multiphase migration. The modeling study affirmed oil generation as an overpressure mechanism, but also elucidated the relative roles of multiphase fluid interaction, oil density and viscosity, and sedimentary compaction. An important result is that overpressures by oil generation create conditions for rock fracturing, and associated fracture permeability may regulate or control the propensity to maintain overpressures.

Examining the effects of microstructure and loading on the shock initiation of HMX with mesoscale simulations

NASA Astrophysics Data System (ADS)

Springer, H. Keo; Tarver, Craig; Bastea, Sorin

2015-06-01

We perform reactive mesoscale simulations to study shock initiation in HMX over a range of pore morphologies and sizes, porosities, and loading conditions in order to improve our understanding of structure-performance relationships. These relationships are important because they guide the development of advanced macroscale models incorporating hot spot mechanisms and the optimization of novel energetic material microstructures. Mesoscale simulations are performed using the multiphysics hydrocode, ALE3D. Spherical, elliptical, polygonal, and crack-like pore geometries 0.1, 1, 10, and 100 microns in size and 2, 5, 10, and 14% porosity are explored. Loading conditions are realized with shock pressures of 6, 10, 20, 38, and 50 GPa. A Cheetah-based tabular model, including temperature-dependent heat capacity, is used for the unreacted and the product equation-of-state. Also, in-line Cheetah is used to probe chemical species evolution. The influence of microstructure and shock loading on shock-to-detonation-transition run distance, reaction rate and product gas species evolution are discussed. This work performed under the auspices of the U.S. DOE by LLNL under Contract DE-AC52-07NA27344. This work is funded by the Joint DoD-DOE Munitions Program.

Characteristics of nuclepore filters with large pore size—I. Physical properties

NASA Astrophysics Data System (ADS)

John, W.; Hering, S.; Reischl, G.; Sasaki, G.; Goren, S.

Measurements of pore diameter, pore density and filter thickness have been made on Nuclepore filters of 5, 8 and 12 μm pore size. The areal distribution of the pores is random, as verified by total hole counts and by counts of overlapping holes. Filter thicknesses decrease with increasing pore diameter. The Hagen-Poiseuille formula accounts for less than half of the measured pressure drop across 12 μm pore size filters. A new calculation, including a term for the pressure drop external to the filter, accounts quantitatively for the observations. There are sufficient variations among filter batches to require knowledge of the filter parameters for each batch to ensure accurate measurements using these filters.

Factors affecting plant growth in membrane nutrient delivery

NASA Technical Reports Server (NTRS)

Dreschel, T. W.; Wheeler, R. M.; Sager, J. C.; Knott, W. M.

1990-01-01

The development of the tubular membrane plant growth unit for the delivery of water and nutrients to roots in microgravity has recently focused on measuring the effects of changes in physical variables controlling solution availability to the plants. Significant effects of membrane pore size and the negative pressure used to contain the solution were demonstrated. Generally, wheat grew better in units with a larger pore size but equal negative pressure and in units with the same pore size but less negative pressure. Lettuce also exhibited better plant growth at less negative pressure.

Pore Characterization of Shale Rock and Shale Interaction with Fluids at Reservoir Pressure-Temperature Conditions Using Small-Angle Neutron Scattering

NASA Astrophysics Data System (ADS)

Ding, M.; Hjelm, R.; Watkins, E.; Xu, H.; Pawar, R.

2015-12-01

Oil/gas produced from unconventional reservoirs has become strategically important for the US domestic energy independence. In unconventional realm, hydrocarbons are generated and stored in nanopores media ranging from a few to hundreds of nanometers. Fundamental knowledge of coupled thermo-hydro-mechanical-chemical (THMC) processes that control fluid flow and propagation within nano-pore confinement is critical for maximizing unconventional oil/gas production. The size and confinement of the nanometer pores creates many complex rock-fluid interface interactions. It is imperative to promote innovative experimental studies to decipher physical and chemical processes at the nanopore scale that govern hydrocarbon generation and mass transport of hydrocarbon mixtures in tight shale and other low permeability formations at reservoir pressure-temperature conditions. We have carried out laboratory investigations exploring quantitative relationship between pore characteristics of the Wolfcamp shale from Western Texas and the shale interaction with fluids at reservoir P-T conditions using small-angle neutron scattering (SANS). We have performed SANS measurements of the shale rock in single fluid (e.g., H2O and D2O) and multifluid (CH4/(30% H2O+70% D2O)) systems at various pressures up to 20000 psi and temperature up to 150 oF. Figure 1 shows our SANS data at different pressures with H2O as the pressure medium. Our data analysis using IRENA software suggests that the principal changes of pore volume in the shale occurred on smaller than 50 nm pores and pressure at 5000 psi (Figure 2). Our results also suggest that with increasing P, more water flows into pores; with decreasing P, water is retained in the pores.

Pore Formation During Solidification of Aluminum: Reconciliation of Experimental Observations, Modeling Assumptions, and Classical Nucleation Theory

NASA Astrophysics Data System (ADS)

Yousefian, Pedram; Tiryakioğlu, Murat

2018-02-01

An in-depth discussion of pore formation is presented in this paper by first reinterpreting in situ observations reported in the literature as well as assumptions commonly made to model pore formation in aluminum castings. The physics of pore formation is reviewed through theoretical fracture pressure calculations based on classical nucleation theory for homogeneous and heterogeneous nucleation, with and without dissolved gas, i.e., hydrogen. Based on the fracture pressure for aluminum, critical pore size and the corresponding probability of vacancies clustering to form that size have been calculated using thermodynamic data reported in the literature. Calculations show that it is impossible for a pore to nucleate either homogeneously or heterogeneously in aluminum, even with dissolved hydrogen. The formation of pores in aluminum castings can only be explained by inflation of entrained surface oxide films (bifilms) under reduced pressure and/or with dissolved gas, which involves only growth, avoiding any nucleation problem. This mechanism is consistent with the reinterpretations of in situ observations as well as the assumptions made in the literature to model pore formation.

CO2 breakthrough pressure and permeability for unsaturated low-permeability sandstone of the Ordos Basin

NASA Astrophysics Data System (ADS)

Zhao, Yan; Yu, Qingchun

2017-07-01

With rising threats from greenhouse gases, capture and injection of CO2 into suitable underground formations is being considered as a method to reduce anthropogenic emissions of CO2 to the atmosphere. As the injected CO2 will remain in storage for hundreds of years, the safety of CO2 geologic sequestration is a major concern. The low-permeability sandstone of the Ordos Basin in China is regarded as both caprock and reservoir rock, so understanding the breakthrough pressure and permeability of the rock is necessary. Because part of the pore volume experiences a non-wetting phase during the CO2 injection and migration process, the rock may be in an unsaturated condition. And if accidental leakage occurs, CO2 will migrate up into the unsaturated zone. In this study, breakthrough experiments were performed at various degrees of water saturation with five core samples of low-permeability sandstone obtained from the Ordos Basin. The experiments were conducted at 40 °C and pressures of >8 MPa to simulate the geological conditions for CO2 sequestration. The results indicate that the degree of water saturation and the pore structure are the main factors affecting the rock breakthrough pressure and permeability, since the influence of calcite dissolution and clay mineral swelling during the saturation process is excluded. Increasing the average pore radius or most probable pore radius leads to a reduction in the breakthrough pressure and an increase by several orders of magnitude in scCO2 effective permeability. In addition, the breakthrough pressure rises and the scCO2 effective permeability decreases when the water saturation increases. However, when the average pore radius is greater than 0.151 μm, the degree of water saturation will has a little effect on the breakthrough pressure. On this foundation, if the most probable pore radius of the core sample reaches 1.760 μm, the breakthrough pressure will not be impacted by the increasing water saturation. We establish correlations between (1) the breakthrough pressure and average pore radius or most probable pore radius, (2) the breakthrough pressure and scCO2 effective permeability, (3) the breakthrough pressure and water saturation, and (4) the scCO2 effective permeability and water saturation. This study provides practical information for further studies of CO2 sequestration as well as the caprock evaluation.

Poromicromechanics reveals that physiological bone strains induce osteocyte-stimulating lacunar pressure.

PubMed

Scheiner, Stefan; Pivonka, Peter; Hellmich, Christian

2016-02-01

Mechanical loads which are macroscopically acting onto bony organs, are known to influence the activities of biological cells located in the pore spaces of bone, in particular so the signaling and production processes mediated by osteocytes. The exact mechanisms by which osteocytes are actually able to "feel" the mechanical loading and changes thereof, has been the subject of numerous studies, and, while several hypotheses have been brought forth over time, this topic has remained a matter of debate. Relaxation times reported in a recent experimental study of Gardinier et al. (Bone 46(4):1075-1081, 2010) strongly suggest that the lacunar pores are likely to experience, during typical physiological load cycles, not only fluid transport, but also undrained conditions. The latter entail the buildup of lacunar pore pressures, which we here quantify by means of a thorough multiscale modeling approach. In particular, the proposed model is based on classical poroelasticity theory, and able to account for multiple pore spaces. First, the model reveals distinct nonlinear dependencies of the resulting lacunar (and vascular) pore pressures on the underlying bone composition, highlighting the importance of a rigorous multiscale approach for appropriate computation of the aforementioned pore pressures. Then, the derived equations are evaluated for macroscopic (uniaxial as well as hydrostatic) mechanical loading of physiological magnitude. The resulting model-predicted pore pressures agree very well with the pressures that have been revealed, by means of in vitro studies, to be of adequate magnitude for modulating the responses of biological cells, including osteocytes. This underlines that osteocytes may respond to many types of loading stimuli at the same time, in particular so to fluid flow and hydrostatic pressure.

Microseismicity Induced by Fluid Pressure Drop (Laboratory Study)

NASA Astrophysics Data System (ADS)

Turuntaev, Sergey; Zenchenko, Evgeny; Melchaeva, Olga

2013-04-01

Pore pressure change in saturated porous rocks may result in its fracturing (Maury et Fourmaintraux, 1993) and corresponding microseismic event occurrences. Microseismicity due to fluid injection is considered in numerous papers (Maxwell, 2010, Shapiro et al., 2005). Another type of the porous medium fracturing is related with rapid pore pressure drop at some boundary. The mechanism of such fracturing was considered by (Khristianovich, 1985) as a model of sudden coal blowing and by (Alidibirov, Panov, 1998) as a model of volcano eruptions. If the porous saturated medium has a boundary where it directly contacted with fluid under the high pressure (in a hydraulic fracture or in a borehole), and the pressure at that boundary is dropped, the conditions for tensile cracks can be achieved at some distance from the boundary. In the paper, the results of experimental study of saturated porous sample fracturing due to pore pressure rapid drop are discussed. The samples (82 mm high, ∅60 mm) were made of quartz sand, which was cemented by "liquid glass" glue with mass fraction 1%. The sample (porosity 35%, uniaxial unconfined compression strength 2.5 MPa) was placed in a mould and saturated by oil. The upper end of the sample contacted with the mould upper lid, the lower end contacted with fluid. The fluid pressure was increased to 10 MPa and then discharged through the bottom nipple. The pressure increases/drops were repeated 30-50 times. Pore pressure and acoustic emission (AE) were registered by transducers mounted into upper and bottom lids of the mould. It was found, that AE sources (corresponded to microfracturing) were spreading from the open end to the closed end of the sample, and that maximal number of AE events was registered at some distance from the opened end. The number of AE pulses increased with every next pressure drop, meanwhile the number of pulses with high amplitudes diminished. It was found that AE maximal rate corresponded to the fluid pressure gradient maximal values. The model of AE relation with the pore pressure gradient was considered based on the following assumptions: AE event occurred when the pore pressure gradient reaches some critical value; the critical value varies and can be described by Weibull distribution. Permeability variation during the fluid pressure drop was estimated by means of fluid pressure data and pore-elastic equation solution for small time intervals (0.01 sec). The study showed possibility to solve both a direct problem of microseismicity variation relation with fluid pressure changes and an inverse problem of defining permeability by registering microseismic activity variation in particular volume of porous medium alongside with pore pressure measurements at some point.

Axisymmetric deformation of a poroelastic layer overlying an elastic half-space due to surface loading

NASA Astrophysics Data System (ADS)

Rani, Sunita; Rani, Sunita

2017-11-01

The axisymmetric deformation of a homogeneous, isotropic, poroelastic layer of uniform thickness overlying a homogeneous, isotropic, elastic half-space due to surface loads has been obtained. The fluid and the solid constituents of the porous layer are compressible and the permeability in vertical direction is different from its permeability in horizontal direction. The displacements and pore-pressure are taken as basic state variables. An analytical solution for the pore-pressure, displacements and stresses has been obtained using the Laplace-Hankel transform technique. The case of normal disc loading is discussed in detail. Diffusion of pore-pressure is obtained in the space-time domain. The Laplace inversion is evaluated using the fixed Talbot algorithm and the Hankel inversion using the extended Simpson's rule. Two different models of the Earth have been considered: continental crust model and oceanic crust model. For continental crust model, the layer is assumed to be of Westerly Granite and for the oceanic crust model of Hanford Basalt. The effect of the compressibilities of the fluid as well as solid constituents and anisotropy in permeability has been studied on the diffusion of pore-pressure. Contour maps have been plotted for the diffusion of pore-pressure for both models. It is observed that the pore-pressure changes to compression for the continental crust model with time, which is not true for the oceanic crust.

Variations of permeability and pore size distribution of porous media with pressure.

PubMed

Chen, Quan; Kinzelbach, Wolfgang; Ye, Chaohui; Yue, Yong

2002-01-01

Porosity and permeability of porous and fractured geological media decrease with the exploitation of formation fluids such as petroleum, natural gas, or ground water. This may result in ground subsidence and a decrease of recovery of petroleum, natural gas, or ground water. Therefore, an evaluation of the behavior of permeability and porosity under formation fluid pressure changes is important to petroleum and ground water industries. This study for the first time establishes a method, which allows for the measurement of permeability, porosity, and pore size distribution of cores simultaneously. From the observation of the pore size distribution by low-field nuclear magnetic resonance (NMR) relaxation time spectrometry the mechanisms of pressure-dependent porosity and permeability change can be derived. This information cannot be obtained by traditional methods. As the large-size pores or fractures contribute significantly to the permeability, their change consequently leads to a large permeability change. The contribution of fractures to permeability is even larger than that of pores. Thus, the permeability of the cores with fractures decreased more than that of cores without fractures during formation pressure decrease. Furthermore, it did not recover during formation pressure increase. It can be concluded that in fractures, mainly plastic deformation takes place, while matrix pores mainly show elastic deformation. Therefore, it is very important to keep an appropriate formation fluid pressure during the exploitation of ground water and petroleum in a fractured formation.

Pore Structure Model for Predicting Elastic Wavespeeds in Fluid-Saturated Sandstones

NASA Astrophysics Data System (ADS)

Zimmerman, R. W.; David, E. C.

2011-12-01

During hydrostatic compression, in the elastic regime, ultrasonic P and S wave velocities measured on rock cores generally increase with pressure, and reach asymptotic values at high pressures. The pressure dependence of seismic velocities is generally thought to be due to the closure of compliant cracks, in which case the high-pressure velocities must reflect only the influence of the non-closable, equant "pores". Assuming that pores can be represented by spheroids, we can relate the elastic properties to the pore structure using an effective medium theory. Moreover, the closure pressure of a thin crack-like pore is directly proportional to its aspect ratio. Hence, our first aim is to use the pressure dependence of seismic velocities to invert the aspect ratio distribution. We use a simple analytical algorithm developed by Zimmerman (Compressibility of Sandstones, 1991), which can be used for any effective medium theory. Previous works have used overly restrictive assumptions, such as assuming that the stiff pores are spherical, or that the interactions between pores can be neglected. Here, we assume that the rock contains an exponential distribution of crack aspect ratios, and one family of stiff pores having an aspect ratio lying somewhere between 0.01 and 1. We develop our model in two versions, using the Differential Scheme, and the Mori-Tanaka scheme. The inversion is done using data obtained in dry experiments, since pore fluids have a strong effect on velocities and tend to mask the effect of the pore geometry. This avoids complicated joint inversion of dry and wet data, such as done by Cheng and Toksoz (JGR, 1979). Our results show that for many sets of data on sandstones, we can fit very well the dry velocities. Our second aim is to predict the saturated velocities from our pore structure model, noting that at a given differential stress, the pore structure should be the same as for a dry test. Our results show that the Biot-Gassmann predictions always underpredict the rock stiffness and that, for ultrasonic measurements performed at high frequencies (~MHz), it is more accurate to use the results from effective medium theories, which implicitly assume that the fluid is trapped in the pores. Hence, we use the aspect ratio distribution inverted from dry data, but this time introducing fluid in the pores. For a good number of experimental data on sandstones, our predictions for the saturated velocities match well the experimental data. This validates the use of a spheroidal model for pores. The results are only very weakly dependent on the choice of the effective medium theory. We conclude that our method, which remain relatively simple, is a useful tool to extract the pore aspect ratio distribution, as well as predicting the saturated velocities for sandstones.

Modeling Thermal Pressurization Around Shallow Dikes Using Temperature-Dependent Hydraulic Properties: Implications for Deformation Around Intrusions

NASA Astrophysics Data System (ADS)

Townsend, Meredith R.

2018-01-01

Pressurization and flow of groundwater around igneous intrusions depend in part on the hydraulic diffusivity of the host rocks and processes that enhance diffusivity, such as fracturing, or decrease diffusivity, such as mineral precipitation during chemical alteration. Characterizing and quantifying the coupled effects of alteration, pore pressurization, and deformation have significant implications for deformation around intrusions, geothermal energy, contact metamorphism, and heat transfer at mid-ocean ridges. Fractures around dikes at Ship Rock, New Mexico, indicate that pore pressures in the host rocks exceeded hydrostatic conditions by at least 15 MPa following dike emplacement. Hydraulic measurements and petrographic analysis indicate that mineral precipitation clogged the pores of the host rock, reducing porosity from 0.25 to <0.10 and reducing permeability by 5 orders of magnitude. Field data from Ship Rock are used to motivate and constrain numerical models for thermal pore fluid pressurization adjacent to a meter-scale dike, using temperature-dependent hydraulic properties in the host rock as a proxy for porosity loss by mineral precipitation during chemical alteration. Reduction in permeability by chemical alteration has a negligible effect on pressurization. However, reduction in porosity by mineral precipitation increases fluid pressure by constricting pore volume and is identified as a potentially significant source of pressure. A scaling relationship is derived to determine when porosity loss becomes important; if permeability is low enough, pressurization by porosity loss outweighs pressurization by thermal expansion of fluids.

Theoretical Analysis of Pore Pressure Diffusion in Some Basic Rock Mechanics Experiments

NASA Astrophysics Data System (ADS)

Braun, Philipp; Ghabezloo, Siavash; Delage, Pierre; Sulem, Jean; Conil, Nathalie

2018-05-01

Non-homogeneity of the pore pressure field in a specimen is an issue for characterization of the thermo-poromechanical behaviour of low-permeability geomaterials, as in the case of the Callovo-Oxfordian claystone ( k < 10-20 m2), a possible host rock for deep radioactive waste disposal in France. In tests with drained boundary conditions, excess pore pressure can result in significant errors in the measurement of material parameters. Analytical solutions are presented for the change in time of the pore pressure field in a specimen submitted to various loading paths and different rates. The pore pressure field in mechanical and thermal undrained tests is simulated with a 1D finite difference model taking into account the dead volume of the drainage system of the triaxial cell connected to the specimen. These solutions provide a simple and efficient tool for the estimation of the conditions that must hold for reliable determination of material parameters and for optimization of various test conditions to minimize the experimental duration, while keeping the measurement errors at an acceptable level.

Numerical simulation of pore pressure changes in levee under flood conditions

NASA Astrophysics Data System (ADS)

Stanisz, Jacek; Borecka, Aleksandra; Pilecki, Zenon; Kaczmarczyk, Robert

2017-11-01

The article discusses the potential for using numerical simulation to assess the development of deformation and pore pressure changes in a levee as a result of the increase and decrease of the flood wave. The simulation made in FLAC 2D did not take into account the filter-erosion deformation associated with seepage in the levee. The simulations were carried out for a field experimental storage consisting of two combined levees, which was constructed with the help of homogeneous cohesive materials with different filtration coefficients. Calculated and measured pore pressure changes were analysed at 4 monitoring points. The water level was increased to 4 m in 96 hours and decreased in 120 hours. The characteristics of the calculated and measured pore pressure changes over time were similar. The maximum values of the calculated and measured pore pressure were almost identical. The only differences were the greater delay of the experimental levee response to changes in water level increase compared to the response of the numerical model. These differences were probably related to filtering-erosion effects during seepage in the levee.

Laser-assisted formation of micropores and nanobubbles in sclera promote stable normalization of intraocular pressure

NASA Astrophysics Data System (ADS)

Baum, Olga; Wachsmann-Hogiu, Sebastian; Milner, Thomas; Sobol, Emil

2017-06-01

Pores in sclera enhance uveoscleral water outflow and can normalize intraocular pressure in glaucomatous eyes. The aims of this study are to demonstrate laser-induced formation of pores with a dendritic structure and to answer the questions: How is a pore system stable and can laser treatment provide a long-lasting pressure stabilization effect? Effect of 1.56 µm laser radiation on porcine eye sclera was studied using atomic force microscopy and super resolution structured irradiation microscopy with fluorescent markers. Results suggest that the pores with a complex spatial configuration can arise as a result of laser irradiation and that laser-generated stable gas nanobubbles coated with calcium ions allow pore stabilization in the sclera. Our results support a laser based approach for treatment of glaucoma.

Constraints on the stress state of the San Andreas Fault with analysis based on core and cuttings from San Andreas Fault Observatory at Depth (SAFOD) drilling phases 1 and 2

USGS Publications Warehouse

Tembe, S.; Lockner, D.; Wong, T.-F.

2009-01-01

Analysis of field data has led different investigators to conclude that the San Andreas Fault (SAF) has either anomalously low frictional sliding strength (?? 0.6). Arguments for the apparent weakness of the SAF generally hinge on conceptual models involving intrinsically weak gouge or elevated pore pressure within the fault zone. Some models assert that weak gouge and/or high pore pressure exist under static conditions while others consider strength loss or fluid pressure increase due to rapid coseismic fault slip. The present paper is composed of three parts. First, we develop generalized equations, based on and consistent with the Rice (1992) fault zone model to relate stress orientation and magnitude to depth-dependent coefficient of friction and pore pressure. Second, we present temperature-and pressure-dependent friction measurements from wet illite-rich fault gouge extracted from San Andreas Fault Observatory at Depth (SAFOD) phase 1 core samples and from weak minerals associated with the San Andreas Fault. Third, we reevaluate the state of stress on the San Andreas Fault in light of new constraints imposed by SAFOD borehole data. Pure talc (?????0.1) had the lowest strength considered and was sufficiently weak to satisfy weak fault heat flow and stress orientation constraints with hydrostatic pore pressure. Other fault gouges showed a systematic increase in strength with increasing temperature and pressure. In this case, heat flow and stress orientation constraints would require elevated pore pressure and, in some cases, fault zone pore pressure in excess of vertical stress. Copyright 2009 by the American Geophysical Union.

Constraints on the stress state of the San Andreas fault with analysis based on core and cuttings from SAFOD drilling phases I and II

USGS Publications Warehouse

Lockner, David A.; Tembe, Cheryl; Wong, Teng-fong

2009-01-01

Analysis of field data has led different investigators to conclude that the San Andreas Fault (SAF) has either anomalously low frictional sliding strength (m < 0.2) or strength consistent with standard laboratory tests (m > 0.6). Arguments for the apparent weakness of the SAF generally hinge on conceptual models involving intrinsically weak gouge or elevated pore pressure within the fault zone. Some models assert that weak gouge and/or high pore pressure exist under static conditions while others consider strength loss or fluid pressure increase due to rapid coseismic fault slip. The present paper is composed of three parts. First, we develop generalized equations, based on and consistent with the Rice (1992) fault zone model to relate stress orientation and magnitude to depth-dependent coefficient of friction and pore pressure. Second, we present temperature- and pressure-dependent friction measurements from wet illite-rich fault gouge extracted from San Andreas Fault Observatory at Depth (SAFOD) phase 1 core samples and from weak minerals associated with the San Andreas Fault. Third, we reevaluate the state of stress on the San Andreas Fault in light of new constraints imposed by SAFOD borehole data. Pure talc (m0.1) had the lowest strength considered and was sufficiently weak to satisfy weak fault heat flow and stress orientation constraints with hydrostatic pore pressure. Other fault gouges showed a systematic increase in strength with increasing temperature and pressure. In this case, heat flow and stress orientation constraints would require elevated pore pressure and, in some cases, fault zone pore pressure in excess of vertical stress.

Dynamic Heterogeneous Multiscale Filtration Model: Probing Micro- and Macroscopic Filtration Characteristics of Gasoline Particulate Filters.

PubMed

Gong, Jian; Viswanathan, Sandeep; Rothamer, David A; Foster, David E; Rutland, Christopher J

2017-10-03

Motivated by high filtration efficiency (mass- and number-based) and low pressure drop requirements for gasoline particulate filters (GPFs), a previously developed heterogeneous multiscale filtration (HMF) model is extended to simulate dynamic filtration characteristics of GPFs. This dynamic HMF model is based on a probability density function (PDF) description of the pore size distribution and classical filtration theory. The microstructure of the porous substrate in a GPF is resolved and included in the model. Fundamental particulate filtration experiments were conducted using an exhaust filtration analysis (EFA) system for model validation. The particulate in the filtration experiments was sampled from a spark-ignition direct-injection (SIDI) gasoline engine. With the dynamic HMF model, evolution of the microscopic characteristics of the substrate (pore size distribution, porosity, permeability, and deposited particulate inside the porous substrate) during filtration can be probed. Also, predicted macroscopic filtration characteristics including particle number concentration and normalized pressure drop show good agreement with the experimental data. The resulting dynamic HMF model can be used to study the dynamic particulate filtration process in GPFs with distinct microstructures, serving as a powerful tool for GPF design and optimization.

Root induced changes of effective 1D hydraulic properties in a soil column.

PubMed

Scholl, P; Leitner, D; Kammerer, G; Loiskandl, W; Kaul, H-P; Bodner, G

Roots are essential drivers of soil structure and pore formation. This study aimed at quantifying root induced changes of the pore size distribution (PSD). The focus was on the extent of clogging vs. formation of pores during active root growth. Parameters of Kosugi's lognormal PSD model were determined by inverse estimation in a column experiment with two cover crops (mustard, rye) and an unplanted control. Pore dynamics were described using a convection-dispersion like pore evolution model. Rooted treatments showed a wider range of pore radii with increasing volumes of large macropores >500 μm and micropores <2.5 μm, while fine macropores, mesopores and larger micropores decreased. The non-rooted control showed narrowing of the PSD and reduced porosity over all radius classes. The pore evolution model accurately described root induced changes, while structure degradation in the non-rooted control was not captured properly. Our study demonstrated significant short term root effects with heterogenization of the pore system as dominant process of root induced structure formation. Pore clogging is suggested as a partial cause for reduced pore volume. The important change in micro- and large macropores however indicates that multiple mechanic and biochemical processes are involved in root-pore interactions.

Effect of thermal pressurization on dynamic rupture propagation under depth-dependent stress

NASA Astrophysics Data System (ADS)

Urata, Y.; Kuge, K.; Kase, Y.

2009-12-01

Fluid and pore pressure evolution can affect dynamic propagation of earthquake ruptures owing to thermal pressurization (e.g., Mase and Smith, 1985). We investigate dynamic rupture propagation with thermal pressurization on a fault subjected to depth-dependent stress, on the basis of 3-D numerical simulations for spontaneous dynamic ruptures. We put a vertical strike-slip rectangular fault in a semi-infinite, homogenous, and elastic medium. The length and width of the fault are 8 and 3 km, respectively. We assume a depth-dependent stress estimated by Yamashita et al. (2004). The numerical algorithm is based on the finite-difference method by Kase and Kuge (2001). A rupture is initiated by increasing shear stress in a small patch at the bottom of the fault, and then proceeds spontaneously, governed by a slip-weakening law with the Coulomb failure criteria. Coefficients of friction and Dc are homogeneous on the fault. On a fault with thermal pressurization, we allow effective normal stress to vary with pore pressure change due to frictional heating by the formulation of Bizzarri and Cocco (2006). When thermal pressurization does not work, tractions drop in the same way everywhere and rupture velocity is subshear except near the free surface. Due to thermal pressurization, dynamic friction on the fault decreases and is heterogeneous not only vertically but horizontally, slip increases, and rupture velocity along the strike direction becomes supershear. As a result, plural peaks of final slip appear, as observed in the case of undrained dip-slip fault by Urata et al. (2008). We found in this study that the early stage of rupture growth under the depth-dependent stress is affected by the location of an initial crack. When a rupture is initiated at the center of the fault without thermal pressurization, the rupture cannot propagate and terminates. Thermal pressurization can help such a powerless rupture to keep propagating.

Inertial effects during irreversible meniscus reconfiguration in angular pores

NASA Astrophysics Data System (ADS)

Ferrari, Andrea; Lunati, Ivan

2014-12-01

In porous media, the dynamics of the invading front between two immiscible fluids is often characterized by abrupt reconfigurations caused by local instabilities of the interface. As a prototype of these phenomena we consider the dynamics of a meniscus in a corner as it can be encountered in angular pores. We investigate this process in detail by means of direct numerical simulations that solve the Navier-Stokes equations in the pore space and employ the Volume of Fluid method (VOF) to track the evolution of the interface. We show that for a quasi-static displacement, the numerically calculated surface energy agrees well with the analytical solutions that we have derived for pores with circular and square cross sections. However, the spontaneous reconfigurations are irreversible and cannot be controlled by the injection rate: they are characterized by the amount of surface energy that is spontaneously released and transformed into kinetic energy. The resulting local velocities can be orders of magnitude larger than the injection velocity and they induce damped oscillations of the interface that possess their own time scales and depend only on fluid properties and pore geometry. In complex media (we consider a network of cubic pores) reconfigurations are so frequent and oscillations last long enough that increasing inertial effects leads to a different fluid distribution by influencing the selection of the next pore to be invaded. This calls into question simple pore-filling rules based only on capillary forces. Also, we demonstrate that inertial effects during irreversible reconfigurations can influence the work done by the external forces that is related to the pressure drop in Darcy's law. This suggests that these phenomena have to be considered when upscaling multiphase flow because local oscillations of the menisci affect macroscopic quantities and modify the constitutive relationships to be used in macro-scale models. These results can be extrapolated to other interface instabilities that are at the origin of fast pore-scale events, such as Haines jumps, snap-off and coalescence.

Research of CO2 and N2 Adsorption Behavior in K-Illite Slit Pores by GCMC Method

PubMed Central

Chen, Guohui; Lu, Shuangfang; Zhang, Junfang; Xue, Qingzhong; Han, Tongcheng; Xue, Haitao; Tian, Shansi; Li, Jinbu; Xu, Chenxi; Pervukhina, Marina; Clennell, Ben

2016-01-01

Understanding the adsorption mechanisms of CO2 and N2 in illite, one of the main components of clay in shale, is important to improve the precision of the shale gas exploration and development. We investigated the adsorption mechanisms of CO2 and N2 in K-illite with varying pore sizes at the temperature of 333, 363 and 393 K over a broad range of pressures up to 30 MPa using the grand canonical Monte Carlo (GCMC) simulation method. The simulation system is proved to be reasonable and suitable through the discussion of the impact of cation dynamics and pore wall thickness. The simulation results of the excess adsorption amount, expressed per unit surface area of illite, is in general consistency with published experimental results. It is found that the sorption potential overlaps in micropores, leading to a decreasing excess adsorption amount with the increase of pore size at low pressure, and a reverse trend at high pressure. The excess adsorption amount increases with increasing pressure to a maximum and then decreases with further increase in the pressure, and the decreasing amount is found to increase with the increasing pore size. For pores with size greater larger than 2 nm, the overlap effect disappears. PMID:27897232

Measurement of variation in soil solute tracer concentration across a range of effective pore sizes

USGS Publications Warehouse

Harvey, Judson W.

1993-01-01

Solute transport concepts in soil are based on speculation that solutes are distributed nonuniformly within large and small pores. Solute concentrations have not previously been measured across a range of pore sizes and examined in relation to soil hydrological properties. For this study, modified pressure cells were used to measure variation in concentration of a solute tracer across a range of pore sizes. Intact cores were removed from the site of a field tracer experiment, and soil water was eluted from 10 or more discrete classes of pore size. Simultaneous changes in water content and unsaturated hydraulic conductivity were determined on cores using standard pressure cell techniques. Bromide tracer concentration varied by as much as 100% across the range of pore sizes sampled. Immediately following application of the bromide tracer on field plots, bromide was most concentrated in the largest pores; concentrations were lower in pores of progressively smaller sizes. After 27 days, bromide was most dilute in the largest pores and concentrations were higher in the smaller pores. A sharp, threefold decrease in specific water capacity during elution indicated separation of two major pore size classes at a pressure of 47 cm H2O and a corresponding effective pore diameter of 70 μm. Variation in tracer concentration, on the other hand, was spread across the entire range of pore sizes investigated in this study. A two-porosity characterization of the transport domain, based on water retention criteria, only broadly characterized the pattern of variation in tracer concentration across pore size classes during transport through a macroporous soil.

Prediction of porosity of food materials during drying: Current challenges and directions.

PubMed

Joardder, Mohammad U H; Kumar, C; Karim, M A

2017-07-18

Pore formation in food samples is a common physical phenomenon observed during dehydration processes. The pore evolution during drying significantly affects the physical properties and quality of dried foods. Therefore, it should be taken into consideration when predicting transport processes in the drying sample. Characteristics of pore formation depend on the drying process parameters, product properties and processing time. Understanding the physics of pore formation and evolution during drying will assist in accurately predicting the drying kinetics and quality of food materials. Researchers have been trying to develop mathematical models to describe the pore formation and evolution during drying. In this study, existing porosity models are critically analysed and limitations are identified. Better insight into the factors affecting porosity is provided, and suggestions are proposed to overcome the limitations. These include considerations of process parameters such as glass transition temperature, sample temperature, and variable material properties in the porosity models. Several researchers have proposed models for porosity prediction of food materials during drying. However, these models are either very simplistic or empirical in nature and failed to consider relevant significant factors that influence porosity. In-depth understanding of characteristics of the pore is required for developing a generic model of porosity. A micro-level analysis of pore formation is presented for better understanding, which will help in developing an accurate and generic porosity model.

Ultrasonic Inspection and Fatigue Evaluation of Critical Pore Size in Welds.

DTIC Science & Technology

1981-09-01

Boiler and Pressure Vessel Code ) 20...Five porosity levels were produced that parallelled ASME boiler and pressure vessel code specification (Section VIII). Appendix IV of the pressure...Figure 2 shows porosity charts (ASME Boiler and Pressure Vessel Code ) which classify and designate the number and size of pores in any six inch length

ROLE OF PRESSURE IN SMECTITE DEHYDRATION - EFFECTS ON GEOPRESSURE AND SMECTITE-TO-ILLITE TRANSFORMATION.

USGS Publications Warehouse

Colten-Bradley, Virginia

1987-01-01

Evaluation of the effects of pressure on the temperature of interlayer water loss (dehydration) by smectites under diagenetic conditions indicates that smectites are stable as hydrated phases in the deep subsurface. Hydraulic and differential pressure conditions affect dehydration differently. The temperature of dehydration increase with pore fluid pressure and interlayer water density. The temperatures of dehydration increase with pore fluid pressure and interlayer water density. The temperatures of dehydration under differential-presssure conditions are inversely related to pressure and interlayer water density. The model presented assumes the effects of pore fluid composition and 2:1 layer reactivity to be negligible. Agreement between theoretical and experimental results validate this assumption. Additional aspects of the subject are discussed.

Effect of Pore Pressure on Slip Failure of an Impermeable Fault: A Coupled Micro Hydro-Geomechanical Model

NASA Astrophysics Data System (ADS)

Yang, Z.; Juanes, R.

2015-12-01

The geomechanical processes associated with subsurface fluid injection/extraction is of central importance for many industrial operations related to energy and water resources. However, the mechanisms controlling the stability and slip motion of a preexisting geologic fault remain poorly understood and are critical for the assessment of seismic risk. In this work, we develop a coupled hydro-geomechanical model to investigate the effect of fluid injection induced pressure perturbation on the slip behavior of a sealing fault. The model couples single-phase flow in the pores and mechanics of the solid phase. Granular packs (see example in Fig. 1a) are numerically generated where the grains can be either bonded or not, depending on the degree of cementation. A pore network is extracted for each granular pack with pore body volumes and pore throat conductivities calculated rigorously based on geometry of the local pore space. The pore fluid pressure is solved via an explicit scheme, taking into account the effect of deformation of the solid matrix. The mechanics part of the model is solved using the discrete element method (DEM). We first test the validity of the model with regard to the classical one-dimensional consolidation problem where an analytical solution exists. We then demonstrate the ability of the coupled model to reproduce rock deformation behavior measured in triaxial laboratory tests under the influence of pore pressure. We proceed to study the fault stability in presence of a pressure discontinuity across the impermeable fault which is implemented as a plane with its intersected pore throats being deactivated and thus obstructing fluid flow (Fig. 1b, c). We focus on the onset of shear failure along preexisting faults. We discuss the fault stability criterion in light of the numerical results obtained from the DEM simulations coupled with pore fluid flow. The implication on how should faults be treated in a large-scale continuum model is also presented.

Framework flexibility of ZIF-8 under liquid intrusion: discovering time-dependent mechanical response and structural relaxation.

PubMed

Sun, Yueting; Li, Yibing; Tan, Jin-Chong

2018-04-18

The structural flexibility of a topical zeolitic imidazolate framework with sodalite topology, termed ZIF-8, has been elucidated through liquid intrusion under moderate pressures (i.e. tens of MPa). By tracking the evolution of water intrusion pressure under cyclic conditions, we interrogate the role of the gate-opening mechanism controlling the size variation of the pore channels of ZIF-8. Interestingly, we demonstrate that its channel deformation is recoverable through structural relaxation over time, hence revealing the viscoelastic mechanical response in ZIF-8. We propose a simple approach employing a glycerol-water solution mixture, which can significantly enhance the sensitivity of intrusion pressure for the detection of structural deformation in ZIF-8. By leveraging the time-dependent gate-opening phenomenon in ZIF-8, we achieved a notable improvement (50%) in energy dissipation during multicycle mechanical deformation experiments.

Thermo-Hydro-Micro-Mechanical 3D Modeling of a Fault Gouge During Co-seismic Slip

NASA Astrophysics Data System (ADS)

Papachristos, E.; Stefanou, I.; Sulem, J.; Donze, F. V.

2017-12-01

A coupled Thermo-Hydro-Micro-Mechanical (THMM) model based on the Discrete Elements method (DEM) is presented for studying the evolving fault gouge properties during pre- and co-seismic slip. Modeling the behavior of the fault gouge at the microscale is expected to improve our understanding on the various mechanisms that lead to slip weakening and finally control the transition from aseismic to seismic slip.The gouge is considered as a granular material of spherical particles [1]. Upon loading, the interactions between particles follow a frictional behavior and explicit dynamics. Using regular triangulation, a pore network is defined by the physical pore space between the particles. The network is saturated by a compressible fluid, and flow takes place following Stoke's equations. Particles' movement leads to pore deformation and thus to local pore pressure increase. Forces exerted from the fluid onto the particles are calculated using mid-step velocities. The fluid forces are then added to the contact forces resulting from the mechanical interactions before the next step.The same semi-implicit, two way iterative coupling is used for the heat-exchange through conduction.Simple tests have been performed to verify the model against analytical solutions and experimental results. Furthermore, the model was used to study the effect of temperature on the evolution of effective stress in the system and to highlight the role of thermal pressurization during seismic slip [2, 3].The analyses are expected to give grounds for enhancing the current state-of-the-art constitutive models regarding fault friction and shed light on the evolution of fault zone propertiesduring seismic slip.[1] Omid Dorostkar, Robert A Guyer, Paul A Johnson, Chris Marone, and Jan Carmeliet. On the role of fluids in stick-slip dynamics of saturated granular fault gouge using a coupled computational fluid dynamics-discrete element approach. Journal of Geophysical Research: Solid Earth, 122(5):3689-3700, 2017.[2] James R Rice. Heating and weakening of faults during earthquake slip. Journal of Geophysical Research: Solid Earth, 111(B5), 2006.[3] Jean Sulem, Ioannis Stefanou, and Emmanuil Veveakis. Stability analysis of undrained adiabatic shearing of a rock layer with cosserat microstructure. Granular Matter, 13(3):261-268,2011.

Untangling Topographic and Climatic Forcing of Earthflow Motion

NASA Astrophysics Data System (ADS)

Finnegan, N. J.; Nereson, A. L.

2017-12-01

Earthflows commonly form in steep river canyons and are argued to initiate from rapid incision that destabilizes hill slope toes. At the same time, earthflows are known to exhibit a temporal pattern of movement that is correlated with seasonal precipitation and associated changes in effective stress. In this contribution, we use infinite slope analysis to illuminate the relative roles of topographic slope and climate (via its control on pore fluid pressure) in influencing earthflow motion at Oak Ridge earthflow, near San Jose, CA. To this end, we synthesize two years of shallow (2.7 m depth) pore fluid pressure data and continuous GPS-derived velocities with an 80-year record of historical deformation derived from tracking of trees and rocks on orthophotos along much of the 1.4 km length and 400 m relief of the earthflow. Multiple lines of evidence suggest that motion of Oak Ridge earthflow occurs as frictional sliding along a discrete failure surface, as argued for other earthflows. Spatial patterns of sliding velocity along the earthflow show the same sensitivity to topographic slope for five discrete periods of historical sliding, accelerating by roughly an order of magnitude along a 20 degree increase in earthflow gradient. In contrast, during the 2016-2017 winter, velocity increased much more rapidly for an equivalent increase in driving stress due to pore-fluid pressure rise at our GPS antenna. During this time period, Oak Ridge earthflow moved approximately 30 cm and we observed a relatively simple, non-linear relationship between GPS-derived sliding velocity and shallow pore fluid pressure. Rapid sliding in 2016-2017 (> 0.6 cm/day) occurred exclusively during the week following a large winter storm event that raised pore pressures to seasonal highs within only 1-2 days of the storm peak. These observations suggests that a mechanism, such as dilatant strengthening, acts to stabilize velocities for a given value of pore fluid pressure in the landslide mass. They also suggest that earthflow motion is more sensitive to pore-fluid pressure forcing than to topographic forcing and challenge the view that attenuation of pore fluid pressure with depth renders large landslides relatively insensitive to high frequency climate variability.

Poroelastic stress-triggering of the 2005 M8.7 Nias earthquake by the 2004 M9.2 Sumatra-Andaman earthquake

USGS Publications Warehouse

Hughes, K.L.H.; Masterlark, Timothy; Mooney, W.D.

2010-01-01

The M9.2 Sumatra-Andaman earthquake (SAE) occurred three months prior to the M8.7 Nias earthquake (NE). We propose that the NE was mechanically triggered by the SAE, and that poroelastic effects were a major component of this triggering. This study uses 3D finite element models (FEMs) of the Sumatra-Andaman subduction zone (SASZ) to predict the deformation, stress, and pore pressure fields of the SAE. The coseismic slip distribution for the SAE is calibrated to near-field GPS data using FEM-generated Green's Functions and linear inverse methods. The calibrated FEM is then used to predict the postseismic poroelastic contribution to stress-triggering along the rupture surface of the NE, which is adjacent to the southern margin of the SAE. The coseismic deformation of the SAE, combined with the rheologic configuration of the SASZ produces two transient fluid flow regimes having separate time constants. SAE coseismic pore pressures in the relatively shallow forearc and volcanic arc regions (within a few km depth) dissipate within one month after the SAE. However, pore pressures in the oceanic crust of the down-going slab persist several months after the SAE. Predictions suggest that the SAE initially induced MPa-scale negative pore pressure near the hypocenter of the NE. This pore pressure slowly recovered (increased) during the three-month interval separating the SAE and NE due to lateral migration of pore fluids, driven by coseismic pressure gradients, within the subducting oceanic crust. Because pore pressure is a fundamental component of Coulomb stress, the MPa-scale increase in pore pressure significantly decreased stability of the NE fault during the three-month interval after the SAE and prior to rupture of the NE. A complete analysis of stress-triggering due to the SAE must include a poroelastic component. Failure to include poroelastic mechanics will lead to an incomplete model that cannot account for the time interval between the SAE and NE. Our transient poroelastic model explains both the spatial and temporal characteristics of triggering of the NE by the SAE. ?? 2010 Elsevier B.V.

Methane gas hydrate effect on sediment acoustic and strength properties

USGS Publications Warehouse

Winters, W.J.; Waite, W.F.; Mason, D.H.; Gilbert, L.Y.; Pecher, I.A.

2007-01-01

To improve our understanding of the interaction of methane gas hydrate with host sediment, we studied: (1) the effects of gas hydrate and ice on acoustic velocity in different sediment types, (2) effect of different hydrate formation mechanisms on measured acoustic properties (3) dependence of shear strength on pore space contents, and (4) pore pressure effects during undrained shear.A wide range in acoustic p-wave velocities (Vp) were measured in coarse-grained sediment for different pore space occupants. Vp ranged from less than 1 km/s for gas-charged sediment to 1.77–1.94 km/s for water-saturated sediment, 2.91–4.00 km/s for sediment with varying degrees of hydrate saturation, and 3.88–4.33 km/s for frozen sediment. Vp measured in fine-grained sediment containing gas hydrate was substantially lower (1.97 km/s). Acoustic models based on measured Vp indicate that hydrate which formed in high gas flux environments can cement coarse-grained sediment, whereas hydrate formed from methane dissolved in the pore fluid may not.The presence of gas hydrate and other solid pore-filling material, such as ice, increased the sediment shear strength. The magnitude of that increase is related to the amount of hydrate in the pore space and cementation characteristics between the hydrate and sediment grains. We have found, that for consolidation stresses associated with the upper several hundred meters of sub-bottom depth, pore pressures decreased during shear in coarse-grained sediment containing gas hydrate, whereas pore pressure in fine-grained sediment typically increased during shear. The presence of free gas in pore spaces damped pore pressure response during shear and reduced the strengthening effect of gas hydrate in sands.

Fluid displacement fronts in porous media: pore scale interfacial jumps, pressure bursts and acoustic emissions

NASA Astrophysics Data System (ADS)

Moebius, Franziska; Or, Dani

2014-05-01

The macroscopically smooth and regular motion of fluid fronts in porous media is composed of numerous rapid pore-scale interfacial jumps and pressure bursts that involve intense interfacial energy release in the form of acoustic emissions. The characteristics of these pore scale events affect residual phase entrapment and transport properties behind the front. We present experimental studies using acoustic emission technique (AE), rapid imaging, and liquid pressure measurements to characterize these processes during drainage and imbibition in simple porous media. Imbibition and drainage produce different AE signatures (AE amplitudes obey a power law). For rapid drainage, AE signals persist long after cessation of front motion reflecting fluid redistribution and interfacial relaxation. Imaging revealed that the velocity of interfacial jumps often exceeds front velocity by more than 50 fold and is highly inertial component (Re>1000). Pore invasion volumes reduced deduced from pressure fluctuations waiting times (for constant withdrawal rates) show remarkable agreement with geometrically-deduced pore volumes. Discrepancies between invaded volumes and geometrical pores increase with increasing capillary numbers due to constraints on evacuation opportunity times and simultaneous invasion events. A mechanistic model for interfacial motions in a pore-throat network was developed to investigate interfacial dynamics focusing on the role of inertia. Results suggest that while pore scale dynamics were sensitive to variations in pore geometry and boundary conditions, inertia exerted only a minor effect on phase entrapment. The study on pore scale invasion events paints a complex picture of rapid and inertial motions and provides new insights on mechanisms at displacement fronts that are essential for improved macroscopic description of multiphase flows in porous media.

Pore-Lining Composition and Capillary Breakthrough Pressure of Mudstone Caprocks: Sealing Efficiency of Geologic CO2 Storage Sites

NASA Astrophysics Data System (ADS)

Heath, J. E.; Dewers, T. A.; McPherson, B. J.; Kotula, P. G.

2010-12-01

Subsurface containment of CO2 is predicated on effective caprock sealing. Many previous studies have relied on macroscopic measurements of capillary breakthrough pressure and other petrophysical properties without direct examination of solid phases that line pore networks and directly contact fluids. However, pore-lining phases strongly contribute to sealing behavior through interfacial interactions among CO2, brine, and the mineral or non-mineral phases. Our high resolution (i.e., sub-micron) examination of the composition of pore-lining phases of several continental and marine mudstones indicates that sealing efficiency (i.e., breakthrough pressure) is governed by pore shapes and pore-lining phases that are not identifiable except through direct characterization of pores. Bulk X-ray diffraction data does not indicate which phases line the pores and may be especially lacking for mudstones with organic material. Organics can line pores and may represent once-mobile phases that modify the wettability of an originally clay-lined pore network. For shallow formations (i.e., < ~800 m depth), interfacial tension and contact angles result in breakthrough pressures that may be as high as those needed to fracture the rock—thus, in the absence of fractures, capillary sealing efficiency is indicated. Deeper seals have poorer capillary sealing if mica-like wetting dominates the wettability. We thank the U.S. Department of Energy’s National Energy Technology Laboratory and the Office of Basic Energy Sciences, and the Southeast and Southwest Carbon Sequestration Partnerships for supporting this work. Sandia National Laboratories is a multi-program laboratory operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Company, for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-AC04-94AL85000.

Observations of wave-induced pore pressure gradients and bed level response on a surf zone sandbar

NASA Astrophysics Data System (ADS)

Anderson, Dylan; Cox, Dan; Mieras, Ryan; Puleo, Jack A.; Hsu, Tian-Jian

2017-06-01

Horizontal and vertical pressure gradients may be important physical mechanisms contributing to onshore sediment transport beneath steep, near-breaking waves in the surf zone. A barred beach was constructed in a large-scale laboratory wave flume with a fixed profile containing a mobile sediment layer on the crest of the sandbar. Horizontal and vertical pore pressure gradients were obtained by finite differences of measurements from an array of pressure transducers buried within the upper several centimeters of the bed. Colocated observations of erosion depth were made during asymmetric wave trials with wave heights between 0.10 and 0.98 m, consistently resulting in onshore sheet flow sediment transport. The pore pressure gradient vector within the bed exhibited temporal rotations during each wave cycle, directed predominantly upward under the trough and then rapidly rotating onshore and downward as the wavefront passed. The magnitude of the pore pressure gradient during each phase of rotation was correlated with local wave steepness and relative depth. Momentary bed failures as deep as 20 grain diameters were coincident with sharp increases in the onshore-directed pore pressure gradients, but occurred at horizontal pressure gradients less than theoretical critical values for initiation of the motion for compact beds. An expression combining the effects of both horizontal and vertical pore pressure gradients with bed shear stress and soil stability is used to determine that failure of the bed is initiated at nonnegligible values of both forces.