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

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.

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.

Gas hydrate property measurements in porous sediments with resonant ultrasound spectroscopy

NASA Astrophysics Data System (ADS)

McGrail, B. P.; Ahmed, S.; Schaef, H. T.; Owen, A. T.; Martin, P. F.; Zhu, T.

2007-05-01

Resonant ultrasound spectroscopy was used to characterize a natural geological core sample obtained from the Mallik 5L-38 gas hydrate research well at high pressure and subambient temperatures. Using deuterated methane gas to form gas hydrate in the core sample, it was discovered that resonance amplitudes are correlated with the fraction of the pore space occupied by the gas hydrate crystals. A pore water freezing model was developed that utilizes the known pore size distribution and pore water chemistry to predict gas hydrate saturation as a function of pressure and temperature. The model showed good agreement with the experimental measurements and demonstrated that pore water chemistry is the most important factor controlling equilibrium gas hydrate saturations in these sediments when gas hydrates are formed artificially in laboratory pressure vessels. With further development, the resonant ultrasound technique can provide a rapid, nondestructive, field portable means of measuring the equilibrium P-T properties and dissociation kinetics of gas hydrates in porous media, determining gas hydrate saturations, and may provide new insights into the nature of gas hydrate formation mechanisms in geologic materials.

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.

High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland

NASA Astrophysics Data System (ADS)

Terakawa, Toshiko; Miller, Stephen A.; Deichmann, Nicholas

2012-07-01

We analyzed 118 well-constrained focal mechanisms to estimate the pore fluid pressure field of the stimulated region during the fluid injection experiment in Basel, Switzerland. This technique, termed focal mechanism tomography (FMT), uses the orientations of slip planes within the prevailing regional stress field as an indicator of the fluid pressure along the plane at the time of slip. The maximum value and temporal change of excess pore fluid pressures are consistent with the known history of the wellhead pressure applied at the borehole. Elevated pore fluid pressures were concentrated within 500 m of the open hole section, which are consistent with the spatiotemporal evolution of the induced microseismicity. Our results demonstrate that FMT is a robust approach, being validated at the meso-scale of the Basel stimulation experiment. We found average earthquake triggering excess pore fluid pressures of about 10 MPa above hydrostatic. Overpressured fluids induced many small events (M < 3) along faults unfavorably oriented relative to the tectonic stress pattern, while the larger events tended to occur along optimally oriented faults. This suggests that small-scale hydraulic networks, developed from the high pressure stimulation, interact to load (hydraulically isolated) high strength bridges that produce the larger events. The triggering pore fluid pressures are substantially higher than that predicted from a linear pressure diffusion process from the source boundary, and shows that the system is highly permeable along flow paths that allow fast pressure diffusion to the boundaries of the stimulated region.

MONITORING OF PORE WATER PRESSURE AND WATER CONTENT AROUND A HORIZONTAL DRIFT THROUGH EXCAVATION - MEASUREMENT AT THE 140m GALLERY IN THE HORONOBE URL -

NASA Astrophysics Data System (ADS)

Yabuuchi, Satoshi; Kunimaru, Takanori; Kishi, Atsuyasu; Komatsu, Mitsuru

Japan Atomic Energy Agency has been conducting the Horonobe Underground Research Laboratory (URL) project in Horonobe, Hokkaido, as a part of the research and development program on geological disposal of high-level radioactive waste. Pore water pressure and water content around a horizontal drift in the URL have been monitored for over 18 months since before the drift excavation was started. During the drift excavation, both pore water pressure and water content were decreasing. Pore water pressure has been still positive though it continued to decrease with its gradient gradually smaller after excavation, while water content turned to increase about 6 months after the completion of the excavation. It turned to fall again about 5 months later. An unsaturated zone containing gases which were dissolved in groundwater may have been formed around the horizontal drift.

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

Elastic and viscoelastic model of the stress history of sedimentary rocks

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

Warpinski, N.R.

A model has been developed to calculate the elastic and viscoelastic stresses which develop in rocks at depth due to burial, uplift and diagenesis. This model includes the effect of the overburden load, tectonic or geometric strains, thermal strains, varying material properties, pore pressure variations, and viscoeleastic relaxation. Calculations for some simple examples are given to show the contributions of the individual stress components due to gravity, tectonics, thermal effects and pore pressure. A complete stress history for Mesaverde rocks in the Piceance basin is calculated based on available burial history, thermal history and expected pore pressure, material property andmore » tectonic strain variations through time. These calculations show the importance of including material property changes and viscoelastic effects. 15 refs., 48 figs.« less

Surge dynamics coupled to pore-pressure evolution in debris flows

USGS Publications Warehouse

Savage, S.B.; Iverson, R.M.; ,

2003-01-01

Temporally and spatially varying pore-fluid pressures exert strong controls on debris-flow motion by mediating internal and basal friction at grain contacts. We analyze these effects by deriving a one-dimensional model of pore-pressure diffusion explicitly coupled to changes in debris-flow thickness. The new pore-pressure equation is combined with Iverson's (1997) extension of the depth-averaged Savage-Hutter (1989, 1991) granular avalanche equations to predict motion of unsteady debris-flow surges with evolving pore-pressure distributions. Computational results illustrate the profound effects of pore-pressure diffusivities on debris-flow surge depths and velocities. ?? 2003 Millpress,.

Pore-scale Analysis of the effects of Contact Angle Hysteresis on Blob Mobilization in a Pore Doublet

NASA Astrophysics Data System (ADS)

Hsu, Shao-Yiu; Glantz, Roland; Hilpert, Markus

2011-11-01

The mobilization of residual oil blobs in porous media is of major interest to the petroleum industry. We studied the Jamin effect, which hampers the blob mobilization, experimentally in a pore doublet model and explain the Jamin effect through contact angle hysteresis. A liquid blob was trapped in one of the tubes of the pore doublet model and then subjected to different pressure gradients. We measured the contact angles (in 2D and 3D) as well as the mean curvatures of the blob. Due to gravity effects and hysteresis, the contact angles of the blob were initially (zero pressure gradient) non-uniform and exhibited a pronounced altitude dependence. As the pressure gradient was increased, the contact angles became more uniform and the altitude dependence of the contact angle decreased. At the same time, the mean curvature of the drainage interface increased, and the mean curvature of the imbibition interface decreased. The pressure drops across the pore model, which we inferred with our theory from the measured contact angles and mean curvatures, were in line with the directly measured pressure data. We not only show that a trapped blob can sustain a finite pressure gradient but also develop methods to measure the contact angles and mean curvatures in 3D.

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.

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.

Three mechanisms model of shale gas in real state transport through a single nanopore

NASA Astrophysics Data System (ADS)

Li, Dongdong; Zhang, Yanyu; Sun, Xiaofei; Li, Peng; Zhao, Fengkai

2018-02-01

At present, the apparent permeability models of shale gas consider only the viscous flow and Knudsen diffusion of free gas, but do not take into account the influence of surface diffusion. Moreover, it is assumed that shale gas is in ideal state. In this paper, shale gas is assumed in real state, a new apparent permeability model for shale gas transport through a single nanopore is developed that captures many important migration mechanisms, such as viscous flow and Knudsen diffusion of free gas, surface diffusion of adsorbed gas. According to experimental data, the accuracy of apparent permeability model was verified. What’s more, the effects of pressure and pore radius on apparent permeability, and the effects on the permeability fraction of viscous flow, Knudsen diffusion and surface diffusion were analysed, separately. Finally, the results indicate that the error of the developed model in this paper was 3.02%, which is less than the existing models. Pressure and pore radius seriously affect the apparent permeability of shale gas. When the pore radius is small or pressure is low, the surface diffusion cannot be ignored. When the pressure and the pore radius is big, the viscous flow occupies the main position.

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

Quaternary Sediment Accumulation in the Aleutian Trench: Implications for Dehydration Reaction Progress and Pore Pressure Development Offshore Alaska

NASA Astrophysics Data System (ADS)

Meridth, L. N.; Screaton, E.; Jaeger, J. M.; James, S. R.; Villaseñor, T. G.

2015-12-01

Sediment inputs to subduction zones impart a significant control on diagenetic reaction progress, fluid production and pore pressure development and thus affect hydrologic and tectonic behavior during subduction. Intensified glaciation following the mid-Pleistocene transition increased sediment flux to the Gulf of Alaska. Rapid sediment accumulation (>1 km/my) in the Aleutian Trench increases overburden and should accelerate dehydration of hydrous sedimentary components by elevating temperatures in the incoming sediment column. These processes have the potential to generate fluid overpressures in the mud-dominated, low permeability sediments deposited on the incoming plate, offshore SE Alaska. Mineralogical analyses on incoming sediments from Deep Sea Drilling Project Leg 18 and Integrated Ocean Drilling Program Expedition 341 show that both smectite and Opal-A are present as hydrous mineral phases. A 1-D numerical model was developed to track dehydration reaction progress and pore pressures in the incoming sediment column from the abyssal plain to the Aleutian Trench. Simulated temperatures in the incoming column increase due to the insulating effect of trench sediments. As a result, trench sedimentation causes smectite dehydration to begin and Opal-A dehydration to nearly reach completion at the deformation front. Simulated excess pore pressures in the proto-decollement zone increase from nearly hydrostatic to almost half of lithostatic due to the rapid deposition of trench sediments. The 1-D modeling results were incorporated into a 2-D model that follows the underthrust column at the deformation front into the subduction zone. Simulated results of the 2-D flow model illustrate the effects of lateral flow on pore pressure distribution following subduction.

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.

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.

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.

Effect of excess pore pressure on the long runout of debris flows over low gradient channels: A case study of the Dongyuege debris flow in Nu River, China

NASA Astrophysics Data System (ADS)

Zhou, Zhen-Hua; Ren, Zhe; Wang, Kun; Yang, Kui; Tang, Yong-Jun; Tian, Lin; Xu, Ze-Min

2018-05-01

Debris flows with long reaches are one of the major natural hazards to human life and property on alluvial fans, as shown by the debris flow that occurred in the Dongyuege (DYG) Gully in August 18, 2010, and caused 96 deaths. The travel distance and the runout distance of the DYG large-scale tragic debris flow were 11 km and 9 km, respectively. In particular, the runout distance over the low gradient channel (channel slope < 5°) upstream of the depositional fan apex reached up to 3.3 km. The build-up and maintenance of excess pore pressure in the debris-flow mass might have played a crucial role in the persistence and long runout of the bouldery viscous debris flow. Experiments to measure pore pressure and pore water escape have been carried out by reconstituting the debris flow bodies with the DYG debris flow deposit. The slurrying of the debris is governed by solid volumetric concentration (SVC), and the difference between the lower SVC limit and the upper SVC limit can be defined as debris flow index (Id). Peak value (Kp) and rate of dissipation (R) of relative excess pore pressure are dependent on SVC. Further, the SVC that gives the lowest rate of dissipation is regarded as the optimum SVC (Cvo). The dissipation response of excess pore pressure can be characterized by the R value under Cvo at a given moment (i.e., 0.5 h, 1 h or 2 h later after peak time). The results reveal that a relatively high level of excess pore pressure developed within the DYG debris-flow mass and had a strong persistence capability. Further research shows that the development, peak value and dissipation of excess pore pressure in a mixture of sediment and water are related to the maximum grain size (MGS), gradation and mineralogy of clay-size particles of the sediment. The layer-lattice silicates in clay particles can be the typical clay minerals, including kaolinite, montmorillonite and illite, and also the unrepresentative clay minerals such as muscovite and chlorite. Moreover, small woody debris can also contribute to the slurrying of sediments and maintenance of debris flows in well vegetated mountainous areas and the boulders suspended in debris flows can elevate excess pore pressure and extend debris-flow mobility. The parameters, including Id, Kp, R and etc., are affected by the intrinsic properties of debris. They, therefore, can reflect the slurrying susceptibility of sediments, and can also be applied to the research on the occurrence mechanisms and risk assessment of other debris flows.

Spectral element modelling of seismic wave propagation in visco-elastoplastic media including excess-pore pressure development

NASA Astrophysics Data System (ADS)

Oral, Elif; Gélis, Céline; Bonilla, Luis Fabián; Delavaud, Elise

2017-12-01

Numerical modelling of seismic wave propagation, considering soil nonlinearity, has become a major topic in seismic hazard studies when strong shaking is involved under particular soil conditions. Indeed, when strong ground motion propagates in saturated soils, pore pressure is another important parameter to take into account when successive phases of contractive and dilatant soil behaviour are expected. Here, we model 1-D seismic wave propagation in linear and nonlinear media using the spectral element numerical method. The study uses a three-component (3C) nonlinear rheology and includes pore-pressure excess. The 1-D-3C model is used to study the 1987 Superstition Hills earthquake (ML 6.6), which was recorded at the Wildlife Refuge Liquefaction Array, USA. The data of this event present strong soil nonlinearity involving pore-pressure effects. The ground motion is numerically modelled for different assumptions on soil rheology and input motion (1C versus 3C), using the recorded borehole signals as input motion. The computed acceleration-time histories show low-frequency amplification and strong high-frequency damping due to the development of pore pressure in one of the soil layers. Furthermore, the soil is found to be more nonlinear and more dilatant under triaxial loading compared to the classical 1C analysis, and significant differences in surface displacements are observed between the 1C and 3C approaches. This study contributes to identify and understand the dominant phenomena occurring in superficial layers, depending on local soil properties and input motions, conditions relevant for site-specific studies.

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.

Pressure prediction in non-uniaxial settings based on field data and geomechanical modeling: a well example

NASA Astrophysics Data System (ADS)

Lockhart, L. P.; Flemings, P. B.; Nikolinakou, M. A.; Heidari, M.

2016-12-01

We apply a new pressure prediction approach that couples sonic velocity data, geomechanical modeling, and a critical state soil model to estimate pore pressure from wellbore data adjacent to a salt body where the stress field is complex. Specifically, we study pressure and stress in front of the Mad Dog salt body, in the Gulf of Mexico. Because of the loading from the salt, stresses are not uniaxial; the horizontal stress is elevated, leading to higher mean and shear stresses. For the Mad Dog field, we develop a relationship between velocity and equivalent effective stress, in order to account for both the mean and shear stress effect on pore pressure. We obtain this equivalent effective stress using a geomechanical model of the Mad Dog field. We show that the new approach improves pressure prediction in areas near salt where mean and shear stress are different than the control well. Our methodology and results show that pore pressure is driven by a combination of mean stress and shear stress, and highlight the importance of shear-induced pore pressures. Furthermore, the impact of our study extends beyond salt bodies; the methodology and gained insights are applicable to geological environments around the world with a complex geologic history, where the stress state is not uniaxial (fault zones, anticlines, synclines, continental margins, etc.).

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.

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.

Extension of the thermal porosimetry method to high gas pressure for nanoporosimetry estimation

NASA Astrophysics Data System (ADS)

Jannot, Y.; Degiovanni, A.; Camus, M.

2018-04-01

Standard pore size determination methods like mercury porosimetry, nitrogen sorption, microscopy, or X-ray tomography are not suited to highly porous, low density, and thus very fragile materials. For this kind of materials, a method based on thermal characterization has been developed in a previous study. This method has been used with air pressure varying from 10-1 to 105 Pa for materials having a thermal conductivity less than 0.05 W m-1 K-1 at atmospheric pressure. It enables the estimation of pore size distribution between 100 nm and 1 mm. In this paper, we present a new experimental device enabling thermal conductivity measurement under gas pressure up to 106 Pa, enabling the estimation of the volume fraction of pores having a 10 nm diameter. It is also demonstrated that the main thermal conductivity models (parallel, series, Maxwell, Bruggeman, self-consistent) lead to the same estimation of the pore size distribution as the extended parallel model (EPM) presented in this paper and then used to process the experimental data. Three materials with thermal conductivities at atmospheric pressure ranging from 0.014 W m-1 K-1 to 0.04 W m-1 K-1 are studied. The thermal conductivity measurement results obtained with the three materials are presented, and the corresponding pore size distributions between 10 nm and 1 mm are presented and discussed.

Variability of in situ sediment strength and pore pressure behavior of tidal estuary surface sediments

NASA Astrophysics Data System (ADS)

Lucking, Greg; Stark, Nina; Lippmann, Thomas; Smyth, Stephen

2017-10-01

Tidal estuaries feature spatially and temporally varying sediment dynamics and characteristics. Particularly, the variability of geotechnical sediment parameters is still poorly understood, limiting the prediction of long-term sediment stability and dynamics. This paper presents results from an in situ investigation of surficial sediments (≤50 cm) in a tidal estuary in New Hampshire (USA), using a portable free fall penetrometer. The aim is to investigate variations in sediment strength and pore pressure behavior with regard to sediment type and seabed morphology. The study also provides a detailed analysis of high velocity impact pore pressure data to derive information about sediment type and permeability. The penetrometer was deployed 227 times, and the findings are correlated to 78 sediment samples. Differences in sediment strength and type were found when transitioning from tidal flats to the deeper channels. Finer-grained sediments located predominantly on the tidal flats appeared well consolidated with noticeable and spatially consistent sediment strength (reflected in an estimate of quasi-static bearing capacity qsbcmax 10 kPa). Sediments with higher sand content (>75%) showed more variations in strength relating to differences in gradation, and likely represent loose and poorly consolidated sands (qsbcmax 10-55 kPa). The rate at which the recorded excess pore pressures approached equilibrium after penetration was classified and related to sediment type. The data indicate that the development of excess pore pressures upon impact and during penetration may provide additional insight into the nature and layering of bed material, such as identifying a desiccated or over-consolidated dilative surficial layer. In summary, with varying sediment grain size distributions, bulk densities and morphology, sediment strength and pore pressure behavior can vary significantly within a tidal estuary.

An Experimental Investigation of the Risk of Triggering Geological Disasters by Injection under Shear Stress

PubMed Central

Liu, Yixin; Xu, Jiang; Peng, Shoujian

2016-01-01

Fluid injection has been applied in many fields, such as hazardous waste deep well injection, forced circulation in geothermal fields, hydraulic fracturing, and CO2 geological storage. However, current research mainly focuses on geological data statistics and the dominating effects of pore pressure. There are only a few laboratory-conditioned studies on the role of drilling boreholes and the effect of injection pressure on the borehole wall. Through experimental phenomenology, this study examines the risk of triggering geological disasters by fluid injection under shear stress. We developed a new direct shear test apparatus, coupled Hydro-Mechanical (HM), to investigate mechanical property variations when an intact rock experienced step drilling borehole, fluid injection, and fluid pressure acting on the borehole and fracture wall. We tested the peak shear stress of sandstone under different experimental conditions, which showed that drilling borehole, water injection, and increased pore pressure led to the decrease in peak shear stress. Furthermore, as pore pressure increased, peak shear stress dispersion increased due to crack propagation irregularity. Because the peak shear stress changed during the fluid injection steps, we suggest that the risk of triggering geological disaster with injection under shear stress, pore, borehole, and fluid pressure should be considered. PMID:27929142

An Experimental Investigation of the Risk of Triggering Geological Disasters by Injection under Shear Stress.

PubMed

Liu, Yixin; Xu, Jiang; Peng, Shoujian

2016-12-08

Fluid injection has been applied in many fields, such as hazardous waste deep well injection, forced circulation in geothermal fields, hydraulic fracturing, and CO 2 geological storage. However, current research mainly focuses on geological data statistics and the dominating effects of pore pressure. There are only a few laboratory-conditioned studies on the role of drilling boreholes and the effect of injection pressure on the borehole wall. Through experimental phenomenology, this study examines the risk of triggering geological disasters by fluid injection under shear stress. We developed a new direct shear test apparatus, coupled Hydro-Mechanical (HM), to investigate mechanical property variations when an intact rock experienced step drilling borehole, fluid injection, and fluid pressure acting on the borehole and fracture wall. We tested the peak shear stress of sandstone under different experimental conditions, which showed that drilling borehole, water injection, and increased pore pressure led to the decrease in peak shear stress. Furthermore, as pore pressure increased, peak shear stress dispersion increased due to crack propagation irregularity. Because the peak shear stress changed during the fluid injection steps, we suggest that the risk of triggering geological disaster with injection under shear stress, pore, borehole, and fluid pressure should be considered.

a Fractal Permeability Model Coupling Boundary-Layer Effect for Tight Oil Reservoirs

NASA Astrophysics Data System (ADS)

Wang, Fuyong; Liu, Zhichao; Jiao, Liang; Wang, Congle; Guo, Hu

A fractal permeability model coupling non-flowing boundary-layer effect for tight oil reservoirs was proposed. Firstly, pore structures of tight formations were characterized with fractal theory. Then, with the empirical equation of boundary-layer thickness, Hagen-Poiseuille equation and fractal theory, a fractal torturous capillary tube model coupled with boundary-layer effect was developed, and verified with experimental data. Finally, the parameters influencing effective liquid permeability were quantitatively investigated. The research results show that effective liquid permeability of tight formations is not only decided by pore structures, but also affected by boundary-layer distributions, and effective liquid permeability is the function of fluid type, fluid viscosity, pressure gradient, fractal dimension, tortuosity fractal dimension, minimum pore radius and maximum pore radius. For the tight formations dominated with nanoscale pores, boundary-layer effect can significantly reduce effective liquid permeability, especially under low pressure gradient.

Nanometer-Scale Pore Characteristics of Lacustrine Shale, Songliao Basin, NE China

PubMed Central

Wang, Min; Yang, Jinxiu; Wang, Zhiwei; Lu, Shuangfang

2015-01-01

In shale, liquid hydrocarbons are accumulated mainly in nanometer-scale pores or fractures, so the pore types and PSDs (pore size distributions) play a major role in the shale oil occurrence (free or absorbed state), amount of oil, and flow features. The pore types and PSDs of marine shale have been well studied; however, research on lacustrine shale is rare, especially for shale in the oil generation window, although lacustrine shale is deposited widely around the world. To investigate the relationship between nanometer-scale pores and oil occurrence in the lacustrine shale, 10 lacustrine shale core samples from Songliao Basin, NE China were analyzed. Analyses of these samples included geochemical measurements, SEM (scanning electron microscope) observations, low pressure CO2 and N2 adsorption, and high-pressure mercury injection experiments. Analysis results indicate that: (1) Pore types in the lacustrine shale include inter-matrix pores, intergranular pores, organic matter pores, and dissolution pores, and these pores are dominated by mesopores and micropores; (2) There is no apparent correlation between pore volumes and clay content, however, a weak negative correlation is present between total pore volume and carbonate content; (3) Pores in lacustrine shale are well developed when the organic matter maturity (Ro) is >1.0% and the pore volume is positively correlated with the TOC (total organic carbon) content. The statistical results suggest that oil in lacustrine shale mainly occurs in pores with diameters larger than 40 nm. However, more research is needed to determine whether this minimum pore diameter for oil occurrence in lacustrine shale is widely applicable. PMID:26285123

Nanometer-Scale Pore Characteristics of Lacustrine Shale, Songliao Basin, NE China.

PubMed

Wang, Min; Yang, Jinxiu; Wang, Zhiwei; Lu, Shuangfang

2015-01-01

In shale, liquid hydrocarbons are accumulated mainly in nanometer-scale pores or fractures, so the pore types and PSDs (pore size distributions) play a major role in the shale oil occurrence (free or absorbed state), amount of oil, and flow features. The pore types and PSDs of marine shale have been well studied; however, research on lacustrine shale is rare, especially for shale in the oil generation window, although lacustrine shale is deposited widely around the world. To investigate the relationship between nanometer-scale pores and oil occurrence in the lacustrine shale, 10 lacustrine shale core samples from Songliao Basin, NE China were analyzed. Analyses of these samples included geochemical measurements, SEM (scanning electron microscope) observations, low pressure CO2 and N2 adsorption, and high-pressure mercury injection experiments. Analysis results indicate that: (1) Pore types in the lacustrine shale include inter-matrix pores, intergranular pores, organic matter pores, and dissolution pores, and these pores are dominated by mesopores and micropores; (2) There is no apparent correlation between pore volumes and clay content, however, a weak negative correlation is present between total pore volume and carbonate content; (3) Pores in lacustrine shale are well developed when the organic matter maturity (Ro) is >1.0% and the pore volume is positively correlated with the TOC (total organic carbon) content. The statistical results suggest that oil in lacustrine shale mainly occurs in pores with diameters larger than 40 nm. However, more research is needed to determine whether this minimum pore diameter for oil occurrence in lacustrine shale is widely applicable.

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

Effects of intermediate wettability on entry capillary pressure in angular pores.

PubMed

Rabbani, Harris Sajjad; Joekar-Niasar, Vahid; Shokri, Nima

2016-07-01

Entry capillary pressure is one of the most important factors controlling drainage and remobilization of the capillary-trapped phases as it is the limiting factor against the two-phase displacement. It is known that the entry capillary pressure is rate dependent such that the inertia forces would enhance entry of the non-wetting phase into the pores. More importantly the entry capillary pressure is wettability dependent. However, while the movement of a meniscus into a strongly water-wet pore is well-defined, the invasion of a meniscus into a weak or intermediate water-wet pore especially in the case of angular pores is ambiguous. In this study using OpenFOAM software, high-resolution direct two-phase flow simulations of movement of a meniscus in a single capillary channel are performed. Interface dynamics in angular pores under drainage conditions have been simulated under constant flow rate boundary condition at different wettability conditions. Our results shows that the relation between the half corner angle of pores and contact angle controls the temporal evolution of capillary pressure during the invasion of a pore. By deviating from pure water-wet conditions, a dip in the temporal evolution of capillary pressure can be observed which will be pronounced in irregular angular cross sections. That enhances the pore invasion with a smaller differential pressure. The interplay between the contact angle and pore geometry can have significant implications for enhanced remobilization of ganglia in intermediate contact angles in real porous media morphologies, where pores are very heterogeneous with small shape factors. Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.

How Pore-Fluid Pressure due to Heavy Rainfall Influences Volcanic Eruptions, Example of 1998 and 2008 Eruptions of Cerro Azul (Galapagos)

NASA Astrophysics Data System (ADS)

Albino, F.; Amelung, F.; Gregg, P. M.

2016-12-01

About 30 worldwide seismic studies have shown a strong correlation between rainfall and earthquakes in the past 22 years (e.g. Costain and Bollinger, 2010). Such correlation has been explained by the phenomenon of hydro-seismicity via pore pressure diffusion: an increase of pore-fluid in the upper crust reduces the normal stress on faults, which can trigger shear failure. Although this pore pressure effect is widely known for earthquakes, this phenomenon and more broadly poro-elasticity process are not widely studied on volcanoes. However, we know from our previous works that tensile failures that open to propagate magma through the surface are also pore pressure dependent. We have demonstrated that an increase of pore pressure largely reduces the overpressure required to rupture the magma reservoir. We have shown that the pore pressure has more influence on reservoir stability than other parameters such as the reservoir depth or the edifice loading. Here, we investigate how small pore-fluid changes due to hydrothermal or aquifer refill during heavy rainfall may perturb the conditions of failure around magma reservoirs and, what is more, if these perturbations are enough to trigger magma intrusions. We quantify the pore pressure effect on magmatic system by combining 1) 1D pore pressure diffusion model to quantify how pore pressure changes from surface to depth after heavy rainfall events and 2) 2D poro-elastic numerical model to provide the evolution of failure conditions of the reservoir as a consequence of these pore pressure changes. Sensitivity analysis is also performed to characterize the influence on our results of the poro-elastic parameters (hydraulic diffusivity, permeability and porosity) and the geometry of the magma reservoir and the aquifer (depth, size, shape). Finally, we apply our methodology to Cerro Azul volcano (Galapagos) where both last eruptions (1998 and 2008) occurred just after heavy rainfall events, without any pre-eruptive inflation. In the event eruptions can be triggered by pore pressure changes in the crust, meteorological records and water table measurements are required to improve eruption's forecast, especially in regions where heavy rainfall events (typhoons, monsoons or hurricanes) frequently occur, such as in Indonesia, Philippines or Central America.

Fluid Pressure in the Shallow Plate Interface at the Nankai Trough Subduction Zone

NASA Astrophysics Data System (ADS)

Tobin, H. J.; Saffer, D.

2003-12-01

The factors controlling the occurrence, magnitude, and other characteristics of great earthquakes is a fundamental outstanding question in fault physics. Pore fluid pressure is perhaps the most critical yet poorly known parameter governing the strength and seismogenic character of plate boundary faults, but unfortunately cannot be directly inferred through available geophysical sensing methods. Moreover, true in situ fluid pressure has proven difficult to measure even in boreholes. At the Nankai Trough, several hundred meters of sediment are subducted beneath the frontal portion of the accretionary prism. The up-dip portion of the plate interface is therefore hosted in these fine-grained marine sedimentary rocks. ODP Leg 190 and 196 showed that these rapidly-loaded underthrust sediments are significantly overpressured near the deformation front. Here, we attempt to quantitatively infer porosity, pore pressure, and effective normal stress at the plate interface at depths currently inaccessible to drilling. Using seismic reflection interval velocity calibrated at the boreholes to porosity, we quantitatively infer pore pressure to ˜ 20 km down-dip of the deformation front, to a plate interface depth of ˜ 6 km. We have developed a Nankai-specific velocity-porosity transform using ODP cores and logs. We use this function to derive a porosity profile for each of two down-dip seismic sections extracted from a 3-D dataset from the Cape Muroto region. We then calculate pore fluid pressure and effective vertical (fault-normal) stress for the underthrust sediment section using a compaction disequilibrium approach and core-based consolidation test data. Because the pore fluid pressure at the fault interface is likely controlled by that of the top of the underthrust section, this calculation represents a quantitative profile of effective stress and pore pressure at the plate interface. Results show that seismic velocity and porosity increase systematically downdip in the underthrust section, but the increase is suppressed relative to that expected from normally consolidating sediments. The computed pore pressure increases landward from an overpressure ratio (λ * = hydrostatic pressure divided by the lithostatic overburden) of ˜ 0.6 at the deformation front to ˜ 0.77 where sediments have been subducted 15 km. The results of this preliminary analysis suggest that a 3-dimensional mapping of predicted effective normal stress in the seismic data volume is possible.

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.

Optimization of a high-pressure pore water extraction device.

PubMed

Cyr, Martin; Daidié, Alain

2007-02-01

High-pressure squeezing is a technique used for the extraction of the pore water of porous materials such as sediments, soils, rocks, and concrete. The efficiency of extraction depends on the maximum pressures on the materials. This article presents the design of a high-pressure device reaching an axial pressure of 1000 MPa which has been developed to improve the efficiency of extraction. The increase in squeezing pressure implies high stresses inside the chamber, so specialized expertise was required to design a safe, functional device that could withstand pressures significantly higher than common laboratory equipment. The design includes finite element calculations, selection of appropriate materials, and descriptive construction details for the apparatus. It also includes an experimental study of the performance of the apparatus in terms of extraction efficiency.

New general pore size distribution model by classical thermodynamics application: Activated carbon

USGS Publications Warehouse

Lordgooei, M.; Rood, M.J.; Rostam-Abadi, M.

2001-01-01

A model is developed using classical thermodynamics to characterize pore size distributions (PSDs) of materials containing micropores and mesopores. The thermal equation of equilibrium adsorption (TEEA) is used to provide thermodynamic properties and relate the relative pore filling pressure of vapors to the characteristic pore energies of the adsorbent/adsorbate system for micropore sizes. Pore characteristic energies are calculated by averaging of interaction energies between adsorbate molecules and adsorbent pore walls as well as considering adsorbate-adsorbate interactions. A modified Kelvin equation is used to characterize mesopore sizes by considering variation of the adsorbate surface tension and by excluding the adsorbed film layer for the pore size. The modified-Kelvin equation provides similar pore filling pressures as predicted by density functional theory. Combination of these models provides a complete PSD of the adsorbent for the micropores and mesopores. The resulting PSD is compared with the PSDs from Jaroniec and Choma and Horvath and Kawazoe models as well as a first-order approximation model using Polanyi theory. The major importance of this model is its basis on classical thermodynamic properties, less simplifying assumptions in its derivation compared to other methods, and ease of use.

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.

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

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.

Effect of Stepwise Pressure Change on Porosity Evolution during Directional Solidification in Small Cylindrical Channels

NASA Technical Reports Server (NTRS)

Grugel, R.N.; Lee, C.P.; Cox, M.C.; Blandford, B.T.; Anilkumar, A.V.

2008-01-01

Controlled directional solidification experiments were performed in capillary channels, using nitrogen-saturated succinonitrile, to examine the effect of an in-situ stepwise processing pressure increase on an isolated pore evolution. Two experiments were performed using different processing pressure input profiles. The results indicate that a processing pressure increase has a transient effect on pore growth geometry characterized by an initial phase of decreasing pore diameter, followed by a recovery phase of increasing pore diameter. The experimental results also show that processing pressure can be used as a control parameter to either increase or terminate porosity formation. A theoretical model is introduced which indicates that the pore formation process is limited by the diffusion of solute-gas through the melt, and that the observed response toa pressure increase is attributed to the re-equilibration of solute concentration in the melt associated with the increased melt pressure.

Using Advanced Tensiometers to Monitor Temporal Variations in Pore Pressure

NASA Astrophysics Data System (ADS)

Nichols, R. L.; Young, M. H.; Dixon, K. L.; Rossabi, J.; Hyde, W. K.; Holmes-Burns, H.

2002-12-01

The Savannah River Site has installed a comprehensive vadose zone monitoring system (VZMS) at it's low level radioactive waste disposal facility to collect the necessary information to calculate contaminant flux. The VZMS includes water content reflectometers, suction lysimeters, advanced tensiometers (ATs), water flux meters, access ports for neutron probes, and a tipping bucket rain gauge. Forty one ATs were installed from 1999 to 2001 at depths ranging from 2 to 60 feet and have been operated continuously. The installation depths were based on a hydrostatigraphic model developed from core logs, cone penetrometer logs, moisture content profiles, water retention curves model that were obtained during the phased installation of the VZMS. An AT consists of a porous cup installed at a prescribed depth with casing back to the surface and a pressure transducer that is lowered into the casing and connects with the porous cup. The pressure transducer transmits it's signal to a datalogger where the data is stored for future retrieval using a cellular phone communications package. Results from the 2 year operating period show that the AT calibrations are stable and t ATs are capable of extended monitoring of pore pressures in the 0 to 300 cm H2 O range. The ATs had sufficient resolution to detect the naturally occurring fluctuations in pore pressure (1 to 100 cm H2 O over 1 to 72 hours) that resulted from infiltration events at the site. The stable performance of the ATs combined with their ability to detect naturally occurring fluctuations in pore pressure make the ATs a useful tool in measuring temporal pore pressure variations for use in calibrating numerical models of fluid flow in variably saturated porous media.

Full paleostress tensor reconstruction: case study of the Panasqueira Mine, Portugal.

NASA Astrophysics Data System (ADS)

Pascal, C.; Jaques Ribeiro, L. M.

2017-12-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 3D exposures of mineralised quartz veins can be found. To carry out full paleostress restoration we needed to determine (1) pore (paleo)pressure and (2) vein attitudes. To these aims we conducted an extensive fluid inclusion study to derive fluid isochores from the quartz of the studied veins. To further 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. As a second step, we measured 600 subhorizontal quartz veins in all the levels of the mine. The inversion of the attitudes of the veins allowed for reconstructing the orientations of the principal axes of stress, the unscaled Mohr circle and the relative pore pressure. After merging these results with the previously obtained absolute pore pressure we reconstructed the six parameters of the paleostress tensor.

Effect of Processing Pressure on Isolated Pore Formation during Controlled Directional Solidification in Small Channels

NASA Technical Reports Server (NTRS)

Cox, Matthew C.; Anilkumar, Amrutur V.; Grugel, RIchard N.; Lee, Chun P.

2008-01-01

Directional solidification experiments were performed, using succinonitrile saturated with nitrogen gas, to examine the effects of in-situ processing pressure changes on the formation growth, and evolution of an isolated, cylindrical gaseous pore. A novel solidification facility, capable of processing thin cylindrical samples (I.D. < 1.0 mm), under controlled pressure conditions, was used for the experiments. A new experimental method for growing the isolated pore from a seed bubble is introduced. The experimental results indicate that an in-situ processing pressure change will result in either a transient change in pore diameter or a complete termination of pore growth, indicating that pressure changes can be used as a control parameter to terminate bubble growth. A simple analytical model has been introduced to explain the experimental observations.

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.

Micromechanics of cataclastic pore collapse in limestone

NASA Astrophysics Data System (ADS)

Zhu, Wei; Baud, Patrick; Wong, Teng-Fong

2010-04-01

The analysis of compactant failure in carbonate formations hinges upon a fundamental understanding of the mechanics of inelastic compaction. Microstructural observations indicate that pore collapse in a limestone initiates at the larger pores, and microcracking dominates the deformation in the periphery of a collapsed pore. To capture these micromechanical processes, we developed a model treating the limestone as a dual porosity medium, with the total porosity partitioned between macroporosity and microporosity. The representative volume element is made up of a large pore which is surrounded by an effective medium containing the microporosity. Cataclastic yielding of this effective medium obeys the Mohr-Coulomb or Drucker-Prager criterion, with failure parameters dependent on porosity and pore size. An analytic approximation was derived for the unconfined compressive strength associated with failure due to the propagation and coalescence of pore-emanated cracks. For hydrostatic loading, identical theoretical results for the pore collapse pressure were obtained using the Mohr-Coulomb or Drucker-Prager criterion. For nonhydrostatic loading, the stress state at the onset of shear-enhanced compaction was predicted to fall on a linear cap according to the Mohr-Coulomb criterion. In contrast, nonlinear caps in qualitative agreement with laboratory data were predicted using the Drucker-Prager criterion. Our micromechanical model implies that the effective medium is significantly stronger and relatively pressure-insensitive in comparison to the bulk sample.

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.

Constitutive model for porous materials

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

Weston, A.M.; Lee, E.L.

1982-01-01

A simple pressure versus porosity compaction model is developed to calculate the response of granular porous bed materials to shock impact. The model provides a scheme for calculating compaction behavior when relatively limited material data are available. While the model was developed to study porous explosives and propellants, it has been applied to a much wider range of materials. The early development of porous material models, such as that of Hermann, required empirical dynamic compaction data. Erkman and Edwards successfully applied the early theory to unreacted porous high explosives using a Gruneisen equation of state without yield behavior and withoutmore » trapped gas in the pores. Butcher included viscoelastic rate dependance in pore collapse. The theoretical treatment of Carroll and Holt is centered on the collapse of a circular pore and includes radial inertia terms and a complex set of stress, strain and strain rate constitutive parameters. Unfortunately data required for these parameters are generally not available. The model described here is also centered on the collapse of a circular pore, but utilizes a simpler elastic-plastic static equilibrium pore collapse mechanism without strain rate dependence, or radial inertia terms. It does include trapped gas inside the pore, a solid material flow stress that creates both a yield point and a variation in solid material pressure with radius. The solid is described by a Mie-Gruneisen type EOS. Comparisons show that this model will accurately estimate major mechanical features which have been observed in compaction experiments.« less

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.

High fluid pressure and triggered earthquakes in the enhanced geothermal system in Basel, Switzerland

NASA Astrophysics Data System (ADS)

Terakawa, T.; Miller, S. A.; Deichmann, N.

2011-12-01

We estimate the pore fluid pressure field of the stimulated region during the fluid injection experiment in Basel, Switzerland by analyzing 118 well-constrained focal mechanisms. This technique, termed focal mechanism tomography (FMT), uses the orientations of the slip planes within the prevailing regional stress field as indicator of the fluid pressure along the plane at the time of slip. Elevated pore fluid pressures were concentrated within 500 m of the open hole section, and we find average earthquake triggering excess pressures of about 10MPa, with a peak value of 19.3 MPa, consistent with the known wellhead pressure applied at the borehole. Our results demonstrate that FMT is a robust approach, being validated at the macroscopic scale of the Basel stimulation experiment. Over-pressurized fluids induced many small events (M < 3) along faults unfavourably-oriented relative to the tectonic stress pattern, while larger events tended to occur along optimally-oriented faults. This suggests that small-scale hydraulic networks, developed from the high pressure stimulation, interact to load (hydraulically isolated) high strength bridges that produce the larger events. The triggering pore fluid pressures are substantially higher than that predicted from a linear pressure diffusion process from the source boundary, showing that the system is highly permeable along flow paths, allowing fast pressure diffusion to the boundaries of the stimulated region.

Developing an Effective Model for Shale Gas Flow in Nano-scale Pore Clusters based on FIB-SEM Images

NASA Astrophysics Data System (ADS)

Jiang, W. B.; Lin, M.; Yi, Z. X.; Li, H. S.

2016-12-01

Nano-scale pores existed in the form of clusters are the controlling void space in shale gas reservoir. Gas transport in nanopores which has a significant influence on shale gas' recoverability displays multiple transport regimes, including viscous, slippage flow and Knudsen diffusion. In addition, it is also influenced by pore space characteristics. For convenience and efficiency consideration, it is necessary to develop an upscaling model from nano pore to pore cluster scale. Existing models are more like framework functions that provide a format, because the parameters that represent pore space characteristics are underdetermined and may have multiple possibilities. Therefore, it is urgent to make them clear and obtained a model that is closer to reality. FIB-SEM imaging technology is able to acquire three dimensional images with nanometer resolution that nano pores can be visible. Based on the images of two shale samples, we used a high-precision pore network extraction algorithm to generate equivalent pore networks and simulate multiple regime (non-Darcy) flow in it. Several structural parameters can be obtained through pore network modelling. It is found that although the throat-radius distributions are very close, throat flux-radius distributions of different samples can be divided into two categories. The variation of tortuosity with pressure and the overall trend of throat-flux distribution changes with pressure are disclosed. A deeper understanding of shale gas flow in nano-scale pore clusters is obtained. After all, an upscaling model that connects absolute permeability, apparent permeability and other characteristic parameters is proposed, and the best parameter scheme considering throat number-radius distribution and flowing porosity for this model is selected out of three schemes based on pore scale results, and it can avoid multiple-solution problem and is useful in reservoir modelling and experiment result analysis, etc. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the National Natural Science Foundation of China (41574129), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the 973 Program (2014CB239004).

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.

Pore Pressure Distribution and Flank Instability in Hydrothermally Altered Stratovolcanoes

NASA Astrophysics Data System (ADS)

Ball, J. L.; Taron, J.; Hurwitz, S.; Reid, M. E.

2015-12-01

Field and geophysical investigations of stratovolcanoes with long-lived hydrothermal systems commonly reveal that initially permeable regions (such as brecciated layers of pyroclastic material) can become both altered and water-bearing. Hydrothermal alteration in these regions, including clay formation, can turn them into low-permeability barriers to fluid flow, which could increase pore fluid pressures resulting in flank slope instability. We examined elevated pore pressure conditions using numerical models of hydrothermal flow in stratovolcanoes, informed by geophysical data about internal structures and deposits. Idealized radially symmetric meshes were developed based on cross-sectional profiles and alteration/permeability structures of Cascade Range stratovolcanoes. We used the OpenGeoSys model to simulate variably saturated conditions in volcanoes heated only by regional heat fluxes, as well as 650°C intrusions at two km depth below the surface. Meteoric recharge was estimated from precipitation rates in the Cascade Range. Preliminary results indicate zones of elevated pore pressures form: 1) where slopes are underlain by continuous low-permeability altered layers, or 2) when the edifice has an altered core with saturated, less permeable limbs. The first scenario might control shallow collapses on the slopes above the altered layers. The second could promote deeper flank collapses that are initially limited to the summit and upper slopes, but could progress to the core of an edifice. In both scenarios, pore pressures can be further elevated by shallow intrusions, or evolve over longer time scales under forcing from regional heat flux. Geometries without confining low-permeability layers do not show these pressure effects. Our initial scenarios use radially symmetric models, but we are also simulating hydrothermal flow under real 3D geometries with asymmetric subsurface structures (Mount Adams). Simulation results will be used to inform 3D slope-stability models.

A novel method to study single-particle dynamics by the resistive pulse technique

NASA Astrophysics Data System (ADS)

Berge, L. I.; Feder, J.; Jøssang, T.

1989-08-01

We have developed a new method, a pressure-reversal technique, which extends the uses of the resistive pulse (Coulter counter) technique to single-particle dynamics. The resistive pulse technique measures the increase in resistance when particles suspended in an electrolyte are transported through a current-carrying aperture. By the new method, the pressure is reversed when a particle exits the pore. A trigger signal, derived from the particle pulses, is used to activate two miniature solenoid valves which serve as pressure switches. In this way, the particle reenters the pore. A single particle flowing back and forth may be studied over a long period of time. The time the particle spends outside the pore between reversals is variable from a few milliseconds to several seconds. We have so far used pore diameters in the range of 3-30 μm. The new technique enables us to study single-particle dissolution and single-particle flow dynamics. The experimental arrangement and the details of the new method are described together with some illustrative measurements.

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.

Numerical models of pore pressure and stress changes along basement faults due to wastewater injection: Applications to the 2014 Milan, Kansas Earthquake

USGS Publications Warehouse

Hearn, Elizabeth H.; Koltermann, Christine; Rubinstein, Justin R.

2018-01-01

We have developed groundwater flow models to explore the possible relationship between wastewater injection and the 12 November 2014 Mw 4.8 Milan, Kansas earthquake. We calculate pore pressure increases in the uppermost crust using a suite of models in which hydraulic properties of the Arbuckle Formation and the Milan earthquake fault zone, the Milan earthquake hypocenter depth, and fault zone geometry are varied. Given pre‐earthquake injection volumes and reasonable hydrogeologic properties, significantly increasing pore pressure at the Milan hypocenter requires that most flow occur through a conductive channel (i.e., the lower Arbuckle and the fault zone) rather than a conductive 3‐D volume. For a range of reasonable lower Arbuckle and fault zone hydraulic parameters, the modeled pore pressure increase at the Milan hypocenter exceeds a minimum triggering threshold of 0.01 MPa at the time of the earthquake. Critical factors include injection into the base of the Arbuckle Formation and proximity of the injection point to a narrow fault damage zone or conductive fracture in the pre‐Cambrian basement with a hydraulic diffusivity of about 3–30 m2/s. The maximum pore pressure increase we obtain at the Milan hypocenter before the earthquake is 0.06 MPa. This suggests that the Milan earthquake occurred on a fault segment that was critically stressed prior to significant wastewater injection in the area. Given continued wastewater injection into the upper Arbuckle in the Milan region, assessment of the middle Arbuckle as a hydraulic barrier remains an important research priority.

Predicting Stress vs. Strain Behaviors of Thin-Walled High Pressure Die Cast Magnesium Alloy with Actual Pore Distribution

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

Choi, Kyoo Sil; Barker, Erin; Cheng, Guang

2016-01-06

In this paper, a three-dimensional (3D) microstructure-based finite element modeling method (i.e., extrinsic modeling method) is developed, which can be used in examining the effects of porosity on the ductility/fracture of Mg castings. For this purpose, AM60 Mg tensile samples were generated under high-pressure die-casting in a specially-designed mold. Before the tensile test, the samples were CT-scanned to obtain the pore distributions within the samples. 3D microstructure-based finite element models were then developed based on the obtained actual pore distributions of the gauge area. The input properties for the matrix material were determined by fitting the simulation result to themore » experimental result of a selected sample, and then used for all the other samples’ simulation. The results show that the ductility and fracture locations predicted from simulations agree well with the experimental results. This indicates that the developed 3D extrinsic modeling method may be used to examine the influence of various aspects of pore sizes/distributions as well as intrinsic properties (i.e., matrix properties) on the ductility/fracture of Mg castings.« less

Heterogeneity, pore pressure, and injectate chemistry: Control measures for geologic carbon storage

DOE PAGES

Dewers, Thomas; Eichhubl, Peter; Ganis, Ben; ...

2017-11-28

Desirable outcomes for geologic carbon storage include maximizing storage efficiency, preserving injectivity, and avoiding unwanted consequences such as caprock or wellbore leakage or induced seismicity during and post injection. Here, to achieve these outcomes, three control measures are evident including pore pressure, injectate chemistry, and knowledge and prudent use of geologic heterogeneity. Field, experimental, and modeling examples are presented that demonstrate controllable GCS via these three measures. Observed changes in reservoir response accompanying CO 2 injection at the Cranfield (Mississippi, USA) site, along with lab testing, show potential for use of injectate chemistry as a means to alter fracture permeabilitymore » (with concomitant improvements for sweep and storage efficiency). Further control of reservoir sweep attends brine extraction from reservoirs, with benefit for pressure control, mitigation of reservoir and wellbore damage, and water use. State-of-the-art validated models predict the extent of damage and deformation associated with pore pressure hazards in reservoirs, timing and location of networks of fractures, and development of localized leakage pathways. Experimentally validated geomechanics models show where wellbore failure is likely to occur during injection, and efficiency of repair methods. Use of heterogeneity as a control measure includes where best to inject, and where to avoid attempts at storage. Lastly, an example is use of waste zones or leaky seals to both reduce pore pressure hazards and enhance residual CO 2 trapping.« less

Heterogeneity, pore pressure, and injectate chemistry: Control measures for geologic carbon storage

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

Dewers, Thomas; Eichhubl, Peter; Ganis, Ben

Desirable outcomes for geologic carbon storage include maximizing storage efficiency, preserving injectivity, and avoiding unwanted consequences such as caprock or wellbore leakage or induced seismicity during and post injection. Here, to achieve these outcomes, three control measures are evident including pore pressure, injectate chemistry, and knowledge and prudent use of geologic heterogeneity. Field, experimental, and modeling examples are presented that demonstrate controllable GCS via these three measures. Observed changes in reservoir response accompanying CO 2 injection at the Cranfield (Mississippi, USA) site, along with lab testing, show potential for use of injectate chemistry as a means to alter fracture permeabilitymore » (with concomitant improvements for sweep and storage efficiency). Further control of reservoir sweep attends brine extraction from reservoirs, with benefit for pressure control, mitigation of reservoir and wellbore damage, and water use. State-of-the-art validated models predict the extent of damage and deformation associated with pore pressure hazards in reservoirs, timing and location of networks of fractures, and development of localized leakage pathways. Experimentally validated geomechanics models show where wellbore failure is likely to occur during injection, and efficiency of repair methods. Use of heterogeneity as a control measure includes where best to inject, and where to avoid attempts at storage. Lastly, an example is use of waste zones or leaky seals to both reduce pore pressure hazards and enhance residual CO 2 trapping.« less

Effects of high shock pressures and pore morphology on hot spot mechanisms in HMX

NASA Astrophysics Data System (ADS)

Springer, H. K.; Tarver, C. M.; Bastea, S.

2017-01-01

The shock initiation and detonation behavior of heterogeneous solid explosives is governed by its microstructure and reactive properties. New additive manufacturing techniques offer unprecedented control of explosive microstructures previously impossible, enabling us to develop novel explosives with tailored shock sensitivity and detonation properties. Since microstructure-performance relationships are not well established for explosives, there is little material design guidance for these manufacturing techniques. In this study, we explore the effects of high shock pressures (15-38 GPa) with long shock durations and different pore morphologies on hot spot mechanisms in HMX. HMX is chosen as the model material because we have experimental data on many of the chemical-thermal-mechanical properties required for pore collapse simulations. Our simulations are performed using the multi-physics arbitrary Lagrangian Eulerian finite element hydrocode, ALE3D, with Cheetah-based models for the unreacted and the product equation-of-states. We use a temperature-dependent specific heat with the unreacted equation-of-state and a temperature-dependent viscosity model to ensure accurate shock temperatures for subsequent chemistry. The Lindemann Law model is used for shock melting in HMX. In contrast to previous pore collapse studies at lower shock pressures (≤10 GPa) in HMX and shorter post-collapse burning times, our calculations show that shock melting occurs above 15 GPa due to higher bulk heating and a prominent elongated ("jet-like") hot spot region forms at later times. The combination of the elongated, post-collapse hot spot region and the higher bulk heating with increasing pressure dramatically increases the growth rate of reaction. Our calculations show that the reaction rate, dF/dt, increases with increasing shock pressure. We decompose the reaction rate into ignition ((dF/dt)ig) and growth ((dF/dt)gr) phases to better analyze our results. We define the ignition phase to primarily include pore collapse and growth phase to primarily include post-collapse grain burning. We are able to track late-time, post-collapse burning due to the unique loading conditions employed in these calculations. We find that (dF/dt)gr > (dF/dt)ig for all pressures considered. (dF/dt)gr changes more significantly from 25 to 38 GPa (from 0.05/µs to >10-100/µs) than from 15 to 25 GPa (from 0.005/µs to 0.05/µs). There is a three order-of-magnitude difference in the reaction from 15 to 38 GPa just after pore collapse. This is qualitatively consistent with fitting the (macroscopic) Ignition and Growth model to high pressure shock initiation data, where much larger reaction fractions are needed to capture the early stages of reaction. Calculated burn rates demonstrate better agreement with data at intermediate times in the growth phase for 15 to 25 GPa and late times for 30 GPa then at any time in the growth phase for 38 GPa. Our calculations are much higher than burn rate data at the earliest times in the growth phase for all pressures, which may reflect the higher localized pressures and temperatures just after pore collapse in the ignition phase. Our calculations with spherical, conical, and elliptical pores show that the influence of morphology on reaction rate is pressure dependent and the most influential pore shapes at lower pressures aren't the same at higher pressures in the regime studied. Altogether these studies provide the basis for developing microstructure-aware models that can be used to design new explosives with optimal performance-safety characteristics. Such models can be used to guide additive manufacturing of explosives and fully exploit their disruptive nature.

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

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

NASA Astrophysics Data System (ADS)

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

2016-12-01

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

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

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.

Elements of an improved model of debris-flow motion

USGS Publications Warehouse

Iverson, R.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. ?? 2009 American Institute of Physics.

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

Explicit Pore Pressure Material Model in Carbon-Cloth Phenolic

NASA Technical Reports Server (NTRS)

Gutierrez-Lemini, Danton; Ehle, Curt

2003-01-01

An explicit material model that uses predicted pressure in the pores of a carbon-cloth phenolic (CCP) composite has been developed. This model is intended to be used within a finite-element model to predict phenomena specific to CCP components of solid-fuel-rocket nozzles subjected to high operating temperatures and to mechanical stresses that can be great enough to cause structural failures. Phenomena that can be predicted with the help of this model include failures of specimens in restrained-thermal-growth (RTG) tests, pocketing erosion, and ply lifting

Particle Deformation and Concentration Polarization in Electroosmotic Transport of Hydrogels through Pores

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

Vlassiouk, Ivan V

2013-01-01

In this article, we report detection of deformable, hydrogel particles by the resistive-pulse technique using single pores in a polymer film. The hydrogels pass through the pores by electroosmosis and cause formation of a characteristic shape of resistive pulses indicating the particles underwent dehydration and deformation. These effects were explained via a non-homogeneous pressure distribution along the pore axis modeled by the coupled Poisson-Nernst-Planck and Navier Stokes equations. The local pressure drops are induced by the electroosmotic fluid flow. Our experiments also revealed the importance of concentration polarization in the detection of hydrogels. Due to the negative charges as wellmore » as branched, low density structure of the hydrogel particles, concentration of ions in the particles is significantly higher than in the bulk. As a result, when electric field is applied across the membrane, a depletion zone can be created in the vicinity of the particle observed as a transient drop of the current. Our experiments using pores with openings between 200 and 1600 nm indicated the concentration polarization dominated the hydrogels detection for pores wider than 450 nm. The results are of importance for all studies that involve transport of molecules, particles and cells through pores with charged walls. The developed inhomogeneous pressure distribution can potentially influence the shape of the transported species. The concentration polarization changes the interpretation of the resistive pulses; the observed current change does not necessarily reflect only the particle size but also the size of the depletion zone that is formed in the particle vicinity.« less

Numerical Simulation of Multiphase Flow in Nanoporous Organic Matter With Application to Coal and Gas Shale Systems

NASA Astrophysics Data System (ADS)

Song, Wenhui; Yao, Jun; Ma, Jingsheng; Sun, Hai; Li, Yang; Yang, Yongfei; Zhang, Lei

2018-02-01

Fluid flow in nanoscale organic pores is known to be affected by fluid transport mechanisms and properties within confined pore space. The flow of gas and water shows notably different characteristics compared with conventional continuum modeling approach. A pore network flow model is developed and implemented in this work. A 3-D organic pore network model is constructed from 3-D image that is reconstructed from 2-D shale SEM image of organic-rich sample. The 3-D pore network model is assumed to be gas-wet and to contain initially gas-filled pores only, and the flow model is concerned with drainage process. Gas flow considers a full range of gas transport mechanisms, including viscous flow, Knudsen diffusion, surface diffusion, ad/desorption, and gas PVT and viscosity using a modified van der Waals' EoS and a correlation for natural gas, respectively. The influences of slip length, contact angle, and gas adsorption layer on water flow are considered. Surface tension considers the pore size and temperature effects. Invasion percolation is applied to calculate gas-water relative permeability. The results indicate that the influences of pore pressure and temperature on water phase relative permeabilities are negligible while gas phase relative permeabilities are relatively larger in higher temperatures and lower pore pressures. Gas phase relative permeability increases while water phase relative permeability decreases with the shrinkage of pore size. This can be attributed to the fact that gas adsorption layer decreases the effective flow area of the water phase and surface diffusion capacity for adsorbed gas is enhanced in small pore size.

Experimental Microfracture Permeability Development in Crystalline Rocks Under Different Tectonic Stress Regimes

NASA Astrophysics Data System (ADS)

Faulkner, D. R.; Armitage, P. J.

2011-12-01

Geothermal fields rely on permeable fracture networks that can act for significant periods of time. In crystalline rocks, permeability may be stimulated by injections of fluid pressure at depth. We show how high-pressure laboratory experiments can be used to quantify the effects of different stress states on the permeability of two rocks; Darley Dale sandstone (~10-16 m2 permeability) and Westerly granite (~10-20 m2 permeability). It is well known that microfractures start to grow at stresses around one half of the failure stress. Failure in the experiments was reproduced in several ways: (1) by fixing σ3 and increasing σ1 - equivalent to a compressive or strike-slip tectonic regime (2) by fixing σ1 and decreasing σ3 - equivalent to an extensional tectonic regime (3) by increasing the pore fluid pressure at a fixed differential stress to simulate high pore fluid pressure failure, and (4) by fixing the mean stress while increasing σ1 and decreasing σ3 in sympathy. Permeability was monitored during all of these tests. From these tests we are able to quantify the relative contributions of mean stress, differential stress and pore fluid pressure on the permeability in the pre-failure region. This provides key data on the development of microfracture permeability that might be produced during the stimulation of geothermal fields during injection within different tectonic environments.

Laboratory Evaluation of EGS Shear Stimulation-Test 001

DOE Data Explorer

Bauer, Steve

2014-07-29

this is the results of an initial setup-shakedon test in order to develop the plumbing system for this test design. a cylinder of granite with offset holes was jacketed and subjected to confining pressure and low temperature (85C) and pore water pressure. flow through the sample was developed at different test stages.

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

New Laboratory Observations of Thermal Pressurization Weakening

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

From Red Cells to Soft Porous Lubrication

NASA Astrophysics Data System (ADS)

Gacka, T.; Nathan, R.; Wu, L.; Wu, Q.; Cbmss Laboratory Team; Chinese Academy Of Sci. Team

2011-11-01

Feng and Weinbaum (J. Fluid. Mech., 422, 282, 2000), inspired by the enhanced lift phenomena in downhill skiing, developed a new lubrication theory for highly compressible porous media where significantly increased lifting force was predicted as a planing surface glided over a soft porous layer; suggesting superior potential use of porous media for soft lubrication. In this study, we experimentally examine the lift generation phenomena by developing a novel soft porous bearing that consists of a running conveyer belt covered with a soft, 100% polyester, porous sheet, and a stationary, fully instrumented, inclined, planar, upper board. Pore pressure was generated as the upper boundary glides over the soft porous bearing and was measured by pressure sensors. One observed that the pore pressure distribution is consistent with predictions by Feng and Weinbaum (2000), and is a function of the relative velocity between the planing surface and the running belt, the mechanical properties (e.g. porosity, permeability and stiffness) and thickness of the porous layer, as well as the compression ratios at the leading and trailing edges. A load cell is used to characterize the performance of the porous bearing, by comparing pore pressure to total lifting forces. The study presented herein significantly improves our understanding of the behavior of highly compressible porous media under fast compression.

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.

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.

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.

Micro X-ray CT Imaging of Sediments under Confining Pressure

NASA Astrophysics Data System (ADS)

Schindler, M.; Prasad, M.

2016-12-01

We developed a pressure and temperature control system for use inside the micro X-ray CT scanner Xradia 400. We succeeded in building a pressure vessel that can be pressurized to 34.5 MPa (5000 psi) while being transparent to X-rays. The setup can currently be cooled to -5°C and heated to 40°C. We were able to observe grain damage and porosity reduction due to applied confining pressure in clean quartz sand samples and quartz sand and bentonite samples. By comparing micro CT images at atmospheric pressure and 13.8 MPa (2000 psi) confining pressure, we observed compaction of the samples resulting in grain damage and fracturing of sediment grains (Figure 1). When the confining pressure was decreased some grains experienced further fracturing. The grain damage appears irreversible. Further fracturing of grains in pre-compacted sediment was observed upon repeated confining pressure cycling. We are currently working on feed-throughs for fluid lines and electric wiring to use ultrasonic transducers and pressure control in combination. Further we plan to include pore pressure in addition to confining pressure into the system. The pressure control system in combination with ultrasonic transducers will allow us to visually observe pore scale changes in rock samples while simultaneously identifying their influence on ultrasonic velocities. Such pore-scale changes are usually not taken into account by rock physics models and could help to identify why laboratory data diverges from theoretical models. Further, it is possible to compute compressibility from mCT images at different stress states by image correlation

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.

Sediment compaction and pore pressure prediction in deepwater basin of the South China Sea: Estimation from ODP and IODP drilling well data

NASA Astrophysics Data System (ADS)

Xie, Yangbing; Wu, Tuoyu; Sun, Jin; Zhang, Hanyu; Wang, Jiliang; Gao, Jinwei; Chen, Chuanxu

2018-02-01

Overpressure in deepwater basins not only causes serious soft sediment deformation, but also significantly affects the safety of drilling operations. Therefore, prediction of overpressure in sediments has become an important task in deepwater oil exploration and development. In this study, we analyze the drilling data from ODP Leg 184 Sites 1144, 1146, and 1148, and IODP Leg 349 Sites U1431, U1432, U1433, and U1435 to study the sediment compaction and controls in the northern South China Sea. Sedimentation rate, sediment content, distribution area, and buried depth are the factors that influence sediment compaction in the deepwater basin of the South China Sea. Among these factors, the sediment content is the most important. The fitted normal compacted coefficients and mudline porosity for an interval of 50 m shows disciplinary variation versus depth. The pore pressure predicted from different fitted results shows varying overpressure situations. The normal compaction trend from Site 1144 reflects the porosity variation trend in stable deposition basins in the northern South China Sea. The predicted pore pressure shows overpressure at Site 1144, which is attributed to compaction disequilibrium. Nevertheless, the mixed lithology column may influence the predicted over-pressure at Site 1148, which is responsible for the confusing result. Above all, we find that sediment compaction should serve as a proxy for pore pressure in the deepwater basin of the South China Sea.

Comparison of the Modified Biot-Gassmann Theory and the Kuster-Toksoz Theory in Predicting Elastic Velocities of Sediments

USGS Publications Warehouse

Lee, Myung W.

2008-01-01

Elastic velocities of water-saturated sandstones depend primarily on porosity, effective pressure, and the degree of consolidation. If the dry-frame moduli are known, from either measurements or theoretical calculations, the effect of pore water on velocities can be modeled using the Gassmann theory. Kuster and Toksoz developed a theory based on wave-scattering theory for a variety of inclusion shapes, which provides a means for calculating dry- or wet-frame moduli. In the Kuster-Toksoz theory, elastic wave velocities through different sediments can be predicted by using different aspect ratios of the sediment's pore space. Elastic velocities increase as the pore aspect ratio increases (larger pore aspect ratio describes a more spherical pore). On the basis of the velocity ratio, which is assumed to be a function of (1-0)n, and the Biot-Gassmann theory, Lee developed a semi-empirical equation for predicting elastic velocities, which is referred to as the modified Biot-Gassmann theory of Lee. In this formulation, the exponent n, which depends on the effective pressure and the degree of consolidation, controls elastic velocities; as n increases, elastic velocities decrease. Computationally, the role of exponent n in the modified Biot-Gassmann theory by Lee is similar to the role of pore aspect ratios in the Kuster-Toksoz theory. For consolidated sediments, either theory predicts accurate velocities. However, for unconsolidated sediments, the modified Biot-Gassmann theory by Lee performs better than the Kuster-Toksoz theory, particularly in predicting S-wave velocities.

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.

Induced groundwater flux by increases in the aquifer's total stress.

PubMed

Chang, Ching-Min; Yeh, Hund-Der

2015-01-01

Fluid-filled granular soils experience changes in total stress because of earth and oceanic tides, earthquakes, erosion, sedimentation, and changes in atmospheric pressure. The pore volume may deform in response to the changes in stress and this may lead to changes in pore fluid pressure. The transient fluid flow can therefore be induced by the gradient in excess pressure in a fluid-saturated porous medium. This work demonstrates the use of stochastic methodology in prediction of induced one-dimensional field-scale groundwater flow through a heterogeneous aquifer. A closed-form of mean groundwater flux is developed to quantify the induced field-scale mean behavior of groundwater flow and analyze the impacts of the spatial correlation length scale of log hydraulic conductivity and the pore compressibility. The findings provided here could be useful for the rational planning and management of groundwater resources in aquifers that contain lenses with large vertical aquifer matrix compressibility values. © 2014, National Ground Water Association.

On the development and benchmarking of an approach to model gas transport in fractured media with immobile water storage

NASA Astrophysics Data System (ADS)

Harp, D. R.; Ortiz, J. P.; Pandey, S.; Karra, S.; Viswanathan, H. S.; Stauffer, P. H.; Anderson, D. N.; Bradley, C. R.

2017-12-01

In unsaturated fractured media, the rate of gas transport is much greater than liquid transport in many applications (e.g., soil vapor extraction operations, methane leaks from hydraulic fracking, shallow CO2 transport from geologic sequestration operations, and later-time radionuclide gas transport from underground nuclear explosions). However, the relatively immobile pore water can inhibit or promote gas transport for soluble constituents by providing storage. In scenarios with constant pressure gradients, the gas transport will be retarded. In scenarios with reversing pressure gradients (i.e. barometric pressure variations) pore water storage can enhance gas transport by providing a ratcheting mechanism. Recognizing the computational efficiency that can be gained using a single-phase model and the necessity of considering pore water storage, we develop a Richard's solution approach that includes kinetic dissolution/volatilization of constituents. Henry's Law governs the equilibrium gaseous/aqueous phase partitioning in the approach. The approach is implemented in a development branch of the PFLOTRAN simulator. We verify the approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) sinusoidal barometric pumping of a fracture, and (4) gas transport along a fracture with uniform flow and diffusive walls. We demonstrate the retardation of gas transport in cases with constant pressure gradients and the enhancement of gas transport with reversing pressure gradients. The figure presents the verification of our approach to the analytical solution of barometric pumping of a fracture from Nilson et al (1991) where the x-axis "Horizontal axis" is the distance into the matrix block from the fracture.

On the development and benchmarking of an approach to model gas transport in fractured media with immobile water storage

NASA Astrophysics Data System (ADS)

Harp, D. R.; Ortiz, J. P.; Pandey, S.; Karra, S.; Viswanathan, H. S.; Stauffer, P. H.; Anderson, D. N.; Bradley, C. R.

2016-12-01

In unsaturated fractured media, the rate of gas transport is much greater than liquid transport in many applications (e.g., soil vapor extraction operations, methane leaks from hydraulic fracking, shallow CO2 transport from geologic sequestration operations, and later-time radionuclide gas transport from underground nuclear explosions). However, the relatively immobile pore water can inhibit or promote gas transport for soluble constituents by providing storage. In scenarios with constant pressure gradients, the gas transport will be retarded. In scenarios with reversing pressure gradients (i.e. barometric pressure variations) pore water storage can enhance gas transport by providing a ratcheting mechanism. Recognizing the computational efficiency that can be gained using a single-phase model and the necessity of considering pore water storage, we develop a Richard's solution approach that includes kinetic dissolution/volatilization of constituents. Henry's Law governs the equilibrium gaseous/aqueous phase partitioning in the approach. The approach is implemented in a development branch of the PFLOTRAN simulator. We verify the approach with analytical solutions of: (1) 1D gas diffusion, (2) 1D gas advection, (3) sinusoidal barometric pumping of a fracture, and (4) gas transport along a fracture with uniform flow and diffusive walls. We demonstrate the retardation of gas transport in cases with constant pressure gradients and the enhancement of gas transport with reversing pressure gradients. The figure presents the verification of our approach to the analytical solution of barometric pumping of a fracture from Nilson et al (1991) where the x-axis "Horizontal axis" is the distance into the matrix block from the fracture.

A two-phase debris-flow model that includes coupled evolution of volume fractions, granular dilatancy, and pore-fluid pressure

USGS Publications Warehouse

George, David L.; Iverson, Richard M.

2011-01-01

Pore-fluid pressure plays a crucial role in debris flows because it counteracts normal stresses at grain contacts and thereby reduces intergranular friction. Pore-pressure feedback accompanying debris deformation is particularly important during the onset of debrisflow motion, when it can dramatically influence the balance of forces governing downslope acceleration. We consider further effects of this feedback by formulating a new, depth-averaged mathematical model that simulates coupled evolution of granular dilatancy, solid and fluid volume fractions, pore-fluid pressure, and flow depth and velocity during all stages of debris-flow motion. To illustrate implications of the model, we use a finite-volume method to compute one-dimensional motion of a debris flow descending a rigid, uniformly inclined slope, and we compare model predictions with data obtained in large-scale experiments at the USGS debris-flow flume. Predictions for the first 1 s of motion show that increasing pore pressures (due to debris contraction) cause liquefaction that enhances flow acceleration. As acceleration continues, however, debris dilation causes dissipation of pore pressures, and this dissipation helps stabilize debris-flow motion. Our numerical predictions of this process match experimental data reasonably well, but predictions might be improved by accounting for the effects of grain-size segregation.

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.

Can the collapse of a fly ash heap develop into an air-fluidized flow? - Reanalysis of the Jupille accident (1961)

NASA Astrophysics Data System (ADS)

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

2015-01-01

A fly ash heap collapse occurred in Jupille (Liege, Belgium) in 1961. The subsequent flow of fly ash reached a surprisingly long runout and had catastrophic consequences. Its unprecedented degree of fluidization attracted scientific attention. As drillings and direct observations revealed no water-saturated zone at the base of the deposits, scientists assumed an air-fluidization mechanism, which appeared consistent with the properties of the material. In this paper, the air-fluidization assumption is tested based on two-dimensional numerical simulations. The numerical model has been developed so as to focus on the most prominent processes governing the flow, with parameters constrained by their physical interpretation. Results are compared to accurate field observations and are presented for different stages in the model enhancement, so as to provide a base for a discussion of the relative influence of pore pressure dissipation and pore pressure generation. These results show that the apparently high diffusion coefficient that characterizes the dissipation of air pore pressures is in fact sufficiently low for an important degree of fluidization to be maintained during a flow of hundreds of meters.

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.

Impact of pore size variability and network coupling on electrokinetic transport in porous media

NASA Astrophysics Data System (ADS)

Alizadeh, Shima; Bazant, Martin Z.; Mani, Ali

2016-11-01

We have developed and validated an efficient and robust computational model to study the coupled fluid and ion transport through electrokinetic porous media, which are exposed to external gradients of pressure, electric potential, and concentration. In our approach a porous media is modeled as a network of many pores through which the transport is described by the coupled Poisson-Nernst-Planck-Stokes equations. When the pore sizes are random, the interactions between various modes of transport may provoke complexities such as concentration polarization shocks and internal flow circulations. These phenomena impact mixing and transport in various systems including deionization and filtration systems, supercapacitors, and lab-on-a-chip devices. In this work, we present simulations of massive networks of pores and we demonstrate the impact of pore size variation, and pore-pore coupling on the overall electrokinetic transport in porous media.

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

Geomechanical modelling of induced seismicity using Coulomb stress and pore pressure changes

NASA Astrophysics Data System (ADS)

Zhao, B.; Shcherbakov, R.

2016-12-01

In recent years, there has been a dramatic increase in seismicity (earthquakes) due to anthropogenic activities related to the unconventional oil and gas exploration in the Western Canada Sedimentary Basin (WCSB). There are compelling evidences that hydraulic fracturing and wastewater injection operations play a key role in induced seismicity in the WCSB; however, their physical mechanisms are still not fully understood. Therefore, this study focuses on exploring the physical mechanisms of induced seismicity and developing a realistic geomechanical model by incorporating the past seismicity and well production data. In this work, we model the Coulomb stress changes due to past moderate (magnitude greater than 3 with known fault plane solutions) induced earthquakes and pore pressure changes due to wastewater injection in Alberta, specifically in Fox Creek and Fort St. John areas. Relationships between Coulombs stress changes, fault geometry and orientation and subsequent earthquake locations are tested. Subsurface flow due to injection well operations is studied to model the pore pressure changes in time and space, using known well production data, which include well types, well locations and water extraction and injection rates. By modelling the changes in pore pressure and Coulomb stress, we aim at constraining the time scale of occurrence of possible future earthquakes. The anticipating results can help to control the parameters of anthropogenic energy related operations such as hydraulic fracturing and wastewater injection in mitigating the risk due to induced seismicity.

Investigation of Thermophysical Parameters Properties for Enhancing Overpressure Mechanism Estimation. Case Study: Miri Area, West Baram Delta

NASA Astrophysics Data System (ADS)

Adha, Kurniawan; Yusoff, Wan Ismail Wan; Almanna Lubis, Luluan

2017-10-01

Determining the pore pressure data and overpressure zone is a compulsory part of oil and gas exploration in which the data can enhance the safety with profit and preventing the drilling hazards. Investigation of thermophysical parameters such as temperature and thermal conductivity can enhance the pore pressure estimation for overpressure mechanism determination. Since those parameters are dependent on rock properties, it may reflect the changes on the column of thermophysical parameters when there is abnormally in pore pressure. The study was conducted in “MRI 1” well offshore Sarawak, where a new approach method designed to determine the overpressure generation. The study was insisted the contribution of thermophysical parameters for supporting the velocity analysis method, petrophysical analysis were done in these studies. Four thermal facies were identified along the well. The overpressure developed below the thermal facies 4, where the pressure reached 38 Mpa and temperature was increasing significantly. The velocity and the thermal conductivity cross plots shows a linear relationship since the both parameters mainly are the function of the rock compaction. When the rock more compact, the particles were brought closer into contact and making the sound wave going faster while the thermal conductivity were increasing. In addition, the increment of temperature and high heat flow indicated the presence of fluid expansion mechanism. Since the shale sonic velocity and density analysis were the common methods in overpressure mechanism and pore pressure estimation. As the addition parameters for determining overpressure zone, the presence of thermophysical analysis was enhancing the current method, where the current method was the single function of velocity analysis. The presence of thermophysical analysis will improve the understanding in overpressure mechanism determination as the new input parameters. Thus, integrated of thermophysical technique and velocity analysis are important parameters in investigating the overpressure mechanisms and pore pressure estimation during oil and gas exploitation in the future.

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.

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 size engineering applied to the design of separators for nickel-hydrogen cells and batteries

NASA Technical Reports Server (NTRS)

Abbey, K. M.; Britton, D. L.

1983-01-01

Pore size engineering in starved alkaline multiplate cells involves adopting techniques to widen the volume tolerance of individual cells. Separators with appropriate pore size distributions and wettability characteristics (capillary pressure considerations) to have wider volume tolerances and an ability to resist dimensional changes in the electrodes were designed. The separators studied for potential use in nickel-hydrogen cells consist of polymeric membranes as well as inorganic microporous mats. In addition to standard measurements, the resistance and distribution of electrolyte as a function of total cell electrolyte content were determined. New composite separators consisting of fibers, particles and/or binders deposited on Zircar cloth were developed in order to engineer the proper capillary pressure characteristics in the separator. These asymmetric separators were prepared from a variety of fibers, particles and binders.

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.

CONCEPTUAL MODEL FOR ORIGIN OF ABNORMALLY PRESSURED GAS ACCUMULATIONS IN LOW-PERMEABILITY RESERVOIRS.

USGS Publications Warehouse

Law, B.E.; Dickinson, W.W.

1985-01-01

The paper suggests that overpressured and underpressured gas accumulations of this type have a common origin. In basins containing overpressured gas accumulations, rates of thermogenic gas accumulation exceed gas loss, causing fluid (gas) pressure to rise above the regional hydrostatic pressure. Free water in the larger pores is forced out of the gas generation zone into overlying and updip, normally pressured, water-bearing rocks. While other diagenetic processes continue, a pore network with very low permeability develops. As a result, gas accumulates in these low-permeability reservoirs at rates higher than it is lost. In basins containing underpressured gas accumulations, rates of gas generation and accumulation are less than gas loss. The basin-center gas accumulation persists, but because of changes in the basin dynamics, the overpressured accumulation evolves into an underpressured system.

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.

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.

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.

Generation of nanobubbles by ceramic membrane filters: The dependence of bubble size and zeta potential on surface coating, pore size and injected gas pressure.

PubMed

Ahmed, Ahmed Khaled Abdella; Sun, Cuizhen; Hua, Likun; Zhang, Zhibin; Zhang, Yanhao; Zhang, Wen; Marhaba, Taha

2018-07-01

Generation of gaseous nanobubbles (NBs) by simple, efficient, and scalable methods is critical for industrialization and applications of nanobubbles. Traditional generation methods mainly rely on hydrodynamic, acoustic, particle, and optical cavitation. These generation processes render issues such as high energy consumption, non-flexibility, and complexity. This research investigated the use of tubular ceramic nanofiltration membranes to generate NBs in water with air, nitrogen and oxygen gases. This system injects pressurized gases through a tubular ceramic membrane with nanopores to create NBs. The effects of membrane pores size, surface energy, and the injected gas pressures on the bubble size and zeta potential were examined. The results show that the gas injection pressure had considerable effects on the bubble size, zeta potential, pH, and dissolved oxygen of the produced NBs. For example, increasing the injection air pressure from 69 kPa to 414 kPa, the air bubble size was reduced from 600 to 340 nm respectively. Membrane pores size and surface energy also had significant effects on sizes and zeta potentials of NBs. The results presented here aim to fill out the gaps of fundamental knowledge about NBs and development of efficient generation methods. Copyright © 2018 Elsevier Ltd. All rights reserved.

Permeability changes induced by microfissure closure and opening in tectonized materials. Effect on slope pore pressure regime.

NASA Astrophysics Data System (ADS)

De la Fuente, Maria; Vaunat, Jean; Pedone, Giuseppe; Cotecchia, Federica; Sollecito, Francesca; Casini, Francesca

2015-04-01

Tectonized clays are complex materials characterized by several levels of structures that may evolve during load and wetting/drying processes. Some microstructural patterns, as microfissures, have a particular influence on the value of permeability which is one of the main factors controlling pore pressure regime in slopes. In this work, the pore pressure regime measured in a real slope of tectonized clay in Southern Italy is analyzed by a numerical model that considers changes in permeability induced by microfissure closure and opening during the wetting and drying processes resulting from climatic actions. Permeability model accounts for the changes in Pore Size Distribution observed by Microscopy Intrusion Porosimetry. MIP tests are performed on representative samples of ground in initial conditions ("in situ" conditions) and final conditions (deformed sample after applying a wetting path that aims to reproduce the saturation of the soil under heavy rains). The resulting measurements allow for the characterization at microstructural level of the soil, identifying the distribution of dominant families pores in the sample and its evolution under external actions. Moreover, comparison of pore size density functions allows defining a microstructural parameter that depends on void ratio and degree of saturation and controls the variation of permeability. Model has been implemented in a thermo-hydro-mechanical code provided with a special boundary condition for climatic actions. Tool is used to analyze pore pressure measurements obtained in the tectonized clay slope. Results are analyzed at the light of the effect that permeability changes during wetting and drying have on the pore pressure regime.

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.

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.

Silica incorporated membrane for wastewater based filtration

NASA Astrophysics Data System (ADS)

Fernandes, C. S.; Bilad, M. R.; Nordin, N. A. H. M.

2017-10-01

Membrane technology has long been applied for waste water treatment industries due to its numerous advantages compared to other conventional processes. However, the biggest challenge in pressure driven membrane process is membrane fouling. Fouling decreases the productivity and efficiency of the filtration, reduces the lifespan of the membrane and reduces the overall efficiency of water treatment processes. In this study, a novel membrane material is developed for water filtration. The developed membrane incorporates silica nanoparticles mainly to improve its structural properties. Membranes with different loadings of silica nanoparticles were applied in this study. The result shows an increase in clean water permeability and filterability of the membrane for treating activated sludge, microalgae solution, secondary effluent and raw sewage as feed. Adding silica into the membrane matrix does not significantly alter contact angle and membrane pore size. We believe that silica acts as an effective pore forming agent that increases the number of pores without significantly altering the pore sizes. A higher number of small pores on the surface of the membrane could reduce membrane fouling because of a low specific loading imposed to individual pores.

Physically based estimation of soil water retention from textural data: General framework, new models, and streamlined existing models

USGS Publications Warehouse

Nimmo, J.R.; Herkelrath, W.N.; Laguna, Luna A.M.

2007-01-01

Numerous models are in widespread use for the estimation of soil water retention from more easily measured textural data. Improved models are needed for better prediction and wider applicability. We developed a basic framework from which new and existing models can be derived to facilitate improvements. Starting from the assumption that every particle has a characteristic dimension R associated uniquely with a matric pressure ?? and that the form of the ??-R relation is the defining characteristic of each model, this framework leads to particular models by specification of geometric relationships between pores and particles. Typical assumptions are that particles are spheres, pores are cylinders with volume equal to the associated particle volume times the void ratio, and that the capillary inverse proportionality between radius and matric pressure is valid. Examples include fixed-pore-shape and fixed-pore-length models. We also developed alternative versions of the model of Arya and Paris that eliminate its interval-size dependence and other problems. The alternative models are calculable by direct application of algebraic formulas rather than manipulation of data tables and intermediate results, and they easily combine with other models (e.g., incorporating structural effects) that are formulated on a continuous basis. Additionally, we developed a family of models based on the same pore geometry as the widely used unsaturated hydraulic conductivity model of Mualem. Predictions of measurements for different suitable media show that some of the models provide consistently good results and can be chosen based on ease of calculations and other factors. ?? Soil Science Society of America. All rights reserved.

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.

Kinetic models of controllable pore growth of anodic aluminum oxide membrane

NASA Astrophysics Data System (ADS)

Huang, Yan; Zeng, Hong-yan; Zhao, Ce; Qu, Ye-qing; Zhang, Pin

2012-06-01

An anodized Al2O3 (AAO) membrane with apertures about 72 nm in diameter was prepared by two-step anodic oxidation. The appearance and pore arrangement of the AAO membrane were characterized by energy dispersive x-ray spectroscopy and scanning electron microscopy. It was confirmed that the pores with high pore aspect ratio were parallel, well-ordered, and uniform. The kinetics of pores growth in the AAO membrane was derived, and the kinetic models showed that pores stopped developing when the pressure ( σ) trended to equal the surface tension at the end of anodic oxidation. During pore expansion, the effects of the oxalic acid concentration and expansion time on the pore size were investigated, and the kinetic behaviors were explained with two kinetic models derived in this study. They showed that the pore size increased with extended time ( r= G· t+ G'), but decreased with increased concentration ( r = - K·ln c- K') through the derived mathematic formula. Also, the values of G, G', K, and K' were derived from our experimental data.

Pore size engineering applied to the design of separators for nickel-hydrogen cells and batteries

NASA Technical Reports Server (NTRS)

Abbey, K. M.; Britton, D. L.

1983-01-01

Pore size engineering in starved alkaline multiplate cells involves adopting techniques to widen the volume tolerance of individual cells. Separators with appropriate pore size distributions and wettability characteristics (capillary pressure considerations) to have wider volume tolerances and an ability to resist dimensional changes in the electrodes were designed. The separators studied for potential use in nickel-hydrogen cells consist of polymeric membranes as well as inorganic microporous mats. In addition to standard measurements, the resistance and distribution of electrolyte as a function of total cell electrolyte content were determined. New composite separators consisting of fibers, particles and/or binders deposited on Zircar cloth were developed in order to engineer the proper capillary pressure characteristics in the separator. These asymmetric separators were prepared from a variety of fibers, particles and binders. Previously announced in STAR as N83-24571

Pre-breakdown processes in a dielectric fluid in inhomogeneous pulsed electric fields

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

Shneider, Mikhail N., E-mail: m.n.shneider@gmail.com; Pekker, Mikhail

2015-06-14

We consider the development of pre-breakdown cavitation nanopores appearing in the dielectric fluid under the influence of the electrostrictive stresses in the inhomogeneous pulsed electric field. It is shown that three characteristic regions can be distinguished near the needle electrode. In the first region, where the electric field gradient is greatest, the cavitation nanopores, occurring during the voltage nanosecond pulse, may grow to the size at which an electron accelerated by the field inside the pores can acquire enough energy for excitation and ionization of the liquid on the opposite pore wall, i.e., the breakdown conditions are satisfied. In themore » second region, the negative pressure caused by the electrostriction is large enough for the cavitation initiation (which can be registered by optical methods), but, during the voltage pulse, the pores do not reach the size at which the potential difference across their borders becomes sufficient for ionization or excitation of water molecules. And, in the third, the development of cavitation is impossible, due to an insufficient level of the negative pressure: in this area, the spontaneously occurring micropores do not grow and collapse under the influence of surface tension forces. This paper discusses the expansion dynamics of the cavitation pores and their most probable shape.« less

Investigation of Coupled model of Pore network and Continuum in shale gas

NASA Astrophysics Data System (ADS)

Cao, G.; Lin, M.

2016-12-01

Flow in shale spanning over many scales, makes the majority of conventional treatment methods disabled. For effectively simulating, a coupled model of pore-scale and continuum-scale was proposed in this paper. Based on the SEM image, we decompose organic-rich-shale into two subdomains: kerogen and inorganic matrix. In kerogen, the nanoscale pore-network is the main storage space and migration pathway so that the molecular phenomena (slip and diffusive transport) is significant. Whereas, inorganic matrix, with relatively large pores and micro fractures, the flow is approximate to Darcy. We use pore-scale network models (PNM) to represent kerogen and continuum-scale models (FVM or FEM) to represent matrix. Finite element mortars are employed to couple pore- and continuum-scale models by enforcing continuity of pressures and fluxes at shared boundary interfaces. In our method, the process in the coupled model is described by pressure square equation, and uses Dirichlet boundary conditions. We discuss several problems: the optimal element number of mortar faces, two categories boundary faces of pore network, the difference between 2D and 3D models, and the difference between continuum models FVM and FEM in mortars. We conclude that: (1) too coarse mesh in mortars will decrease the accuracy, while too fine mesh will lead to an ill-condition even singular system, the optimal element number is depended on boundary pores and nodes number. (2) pore network models are adjacent to two different mortar faces (PNM to PNM, PNM to continuum model), incidental repeated mortar nodes must be deleted. (3) 3D models can be replaced by 2D models under certain condition. (4) FVM is more convenient than FEM, for its simplicity in assigning interface nodes pressure and calculating interface fluxes. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the 973 Program (2014CB239004), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the National Natural Science Foundation of China (41574129).

The influence of pore textures on the permeability of volcanic rocks

NASA Astrophysics Data System (ADS)

Mueller, S.; Spieler, O.; Scheu, B.; Dingwell, D.

2006-12-01

The permeability of a porous medium is strongly dependent on its porosity, as a higher proportion of pore volume is generally expected to lead to a greater probability of pore interconnectedness and the formation of a fluid-flow providing pathway. However, the relationship between permeability and porosity is not a unique one, as many other textural parameters may play an important role and substantially affect gas flow properties. Among these parameters are (a) the connection geometry (i.e. intergranular pore spaces in clastic sediments vs. bubble interconnections), (b) the pore sizes, (c) pore shape and (d) pore size distribution. The gas permeability of volcanic rocks may influence various eruptive processes. The transition from a quiescent degassing dome to rock failure (fragmentation) may, for example, be controlled by the rock's permeability, in as much as it affects the speed by which a gas overpressure in vesicles is reduced in response to decompression. It is therefore essential to understand and quantify influences of different pore textures on the degassing properties of volcanic rocks, as well as investigate the effects of permeability on eruptive processes. Using a modified shock-tube-based fragmentation apparatus, we have measured unsteady-state permeability at a high initial pressure differential. Following sudden decompression above the rock cylinder, pressurized gas flows through the sample in a steel autoclave. A transient 1D filtration code has been developed to calculate permeability using the experimental pressure decay curve within a defined volume below the sample. An external furnace around the autoclave and the use of compressed salt as sealant allows also measurements at high temperatures up to 800 °C. Over 130 permeability measurements have been performed on samples of different volcanic settings, covering a wide range of porosity. The results show a general positive relationship between porosity and permeability with a high data scatter. Analysis of the samples eruptive origin as well as the pore sizes, shapes and size distribution allow an estimation of the contribution of various textural effects to the overall permeability.

Analysis of Effective Interconnectivity of DegraPol-foams Designed for Negative Pressure Wound Therapy

PubMed Central

Milleret, Vincent; Bittermann, Anne Greet; Mayer, Dieter; Hall, Heike

2009-01-01

Many wounds heal slowly and are difficult to manage. Therefore Negative Pressure Wound Therapy (NPWT) was developed where polymer foams are applied and a defined negative pressure removes wound fluid, reduces bacterial burden and increases the formation of granulation tissue. Although NPWT is used successfully, its mechanisms are not well understood. In particular, different NPWT dressings were never compared. Here a poly-ester urethane Degrapol® (DP)-foam was produced and compared with commercially available dressings (polyurethane-based and polyvinyl-alcohol-based) in terms of apparent pore sizes, swelling and effective interconnectivity of foam pores. DP-foams contain relatively small interconnected pores; PU-foams showed large pore size and interconnectivity; whereas PVA-foams displayed heterogeneous and poorly interconnected pores. PVA-foams swelled by 40 %, whereas DP- and PU-foams remained almost without swelling. Effective interconnectivity was investigated by submitting fluorescent beads of 3, 20 and 45 μm diameter through the foams. DP- and PU-foams removed 70-90 % of all beads within 4 h, independent of the bead diameter or bead pre-adsorption with serum albumin. For PVA-foams albumin pre-adsorbed beads circulated longer, where 20 % of 3 μm and 10 % of 20 μm diameter beads circulated after 96 h. The studies indicate that efficient bead perfusion does not only depend on pore size and swelling capacity, but effective interconnectivity might also depend on chemical composition of the foam itself. In addition due to the efficient sieve-effect of the foams uptake of wound components in vivo might occur only for short time suggesting other mechanisms being decisive for success of NPWT.

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.

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.

Toward direct pore-scale modeling of three-phase displacements

NASA Astrophysics Data System (ADS)

Mohammadmoradi, Peyman; Kantzas, Apostolos

2017-12-01

A stable spreading film between water and gas can extract a significant amount of bypassed non-aqueous phase liquid (NAPL) through immiscible three-phase gas/water injection cycles. In this study, the pore-scale displacement mechanisms by which NAPL is mobilized are incorporated into a three-dimensional pore morphology-based model under water-wet and capillary equilibrium conditions. The approach is pixel-based and the sequence of invasions is determined by the fluids' connectivity and the threshold capillary pressure of the advancing interfaces. In addition to the determination of three-phase spatial saturation profiles, residuals, and capillary pressure curves, dynamic finite element simulations are utilized to predict the effective permeabilities of the rock microtomographic images as reasonable representations of the geological formations under study. All the influential features during immiscible fluid flow in pore-level domains including wetting and spreading films, saturation hysteresis, capillary trapping, connectivity, and interface development strategies are taken into account. The capabilities of the model are demonstrated by the successful prediction of saturation functions for Berea sandstone and the accurate reconstruction of three-phase fluid occupancies through a micromodel.

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.