Investigation of ultrasonic wave interactions with fluid-saturated porous rocks
Adler, L.
1990-01-01
Investigation of ultrasonic waves, especially the slow compressional wave, with fluid-saturated porous solids, especially rocks. This research effort should fine applications in the geophysical evaluation of fluid-bearing porous rocks where parameters such as tortuosity, permeability, saturation level, and internal impurities are difficult to measure by conventional techniques. The proposed investigation may be divided into three major subtasks: Experimental study of surface wave propagation on fluid-saturated porous materials. A new, so-called direct excitation technique will be used on both air- and water-saturated samples; further development of the Lamb wave technique recently introduced to study guided wave propagation in thin fluid-saturated porous plates. The analytical treatment will be extended to account for viscous losses and scattering inhomogeneities; and theoretical and experimental study of slow wave propagation in fluid-saturated natural rocks. A new technique based on the transmission of airborne ultrasound through air-saturated porous plates will be used to determine properties such as tortuosity, permeability, etc. 22 refs., 37 figs., 3 tabs.
Elastic and inelastic deformation of fluid-saturated rock.
Makhnenko, Roman Y; Labuz, Joseph F
2016-10-13
In situ rock is often saturated with fluid, the presence of which affects both elastic parameters and inelastic deformation processes. Techniques were developed for testing fluid-saturated porous rock under the limiting conditions of drained (long-term), undrained (short-term) and unjacketed (solid matrix) response in hydrostatic, axisymmetric and plane-strain compression. Drained and undrained poroelastic parameters, including bulk modulus, Biot and Skempton coefficients, of Berea sandstone were found to be stress dependent up to 35 MPa mean stress, and approximately constant at higher levels of loading. The unjacketed bulk modulus was measured to be constant for pressure up to 60 MPa, and it appears to be larger than the unjacketed pore bulk modulus. An elasto-plastic constitutive model calibrated with parameters from drained tests provided a first-order approximation of undrained inelastic deformation: dilatant hardening was observed due to pore pressure decrease during inelastic deformation of rock specimens with constant fluid content.This article is part of the themed issue 'Energy and the subsurface'. PMID:27597783
Standing Torsional Waves in Fluid-Saturated Porous Circular Cylinder
NASA Astrophysics Data System (ADS)
Solorza, S.; Sahay, P. N.
2002-12-01
For dynamic measurement of elastic constants of a porous material saturated with viscous fluid when resonance-bar technique is applied, one also observes attenuation of the wave field. The current practice is to interpret it in terms of solid-viscosity by assuming a viscoelastic rheology for porous material. The likely mechanisms of attenuation in a fluid saturated porous material are: 1) motion of the fluid with respect to the solid frame and 2) viscous loss within the pore fluid. Therefore, it is appropriate to assume a poroelastic rheology and link the observed attenuation value to fluid properties and permeability. In the framework of poroelastic theory, the explicit formula linking attenuation to the properties of solid and fluid constituents and permeability are not worked out yet. In order to established such a link one has to workout solutions of appropriate boundary value problems in such a framework. Here, we have carried out the solution of boundary value problem associated with torsional oscillation of a finite poroelastic circular cylinder, casted in the framework of volume-averaged theory of poroelasticity. Analysing this solution by a perturbative approach we are able to develop explicit expressions for resonance frequency and attenuation for this mode of vibration. It shows how the attenuation is controlled by the permeability and the fluid properties, and how the resonance frequency drops over its value for the dry porous frame due to the effect of the fluid-mass.
A coupled deformation-diffusion theory for fluid-saturated porous solids
NASA Astrophysics Data System (ADS)
Henann, David; Kamrin, Ken; Anand, Lallit
2012-02-01
Fluid-saturated porous materials are important in several familiar applications, such as the response of soils in geomechanics, food processing, pharmaceuticals, and the biomechanics of living bone tissue. An appropriate constitutive theory describing the coupling of the mechanical behavior of the porous solid with the transport of the fluid is a crucial ingredient towards understanding the material behavior in these varied applications. In this work, we formulate and numerically implement in a finite-element framework a large-deformation theory for coupled deformation-diffusion in isotropic, fluid-saturated porous solids. The theory synthesizes the classical Biot theory of linear poroelasticity and the more-recent Coussy theory of poroplasticity in a large deformation framework. In this talk, we highlight several salient features of our theory and discuss representative examples of the application of our numerical simulation capability to problems of consolidation as well as deformation localization in granular materials.
Seismic attenuation: Laboratory measurements in fluid saturated rocks
NASA Astrophysics Data System (ADS)
Subramaniyan, Shankar; Madonna, Claudio; Tisato, Nicola; Saenger, Erik; Quintal, Beatriz
2014-05-01
Seismic wave attenuation could be used as an indicator of reservoir fluids due to its dependence on rock and fluid properties. Over the past 30 years, many laboratory methodologies to study attenuation in rocks have been employed, such as ultrasonic (MHz), resonant bar (kHz) and forced oscillation methods in the low frequency range (0.01-100Hz) (Tisato & Madonna 2012; Madonna & Tisato 2013). Forced oscillation methods have gained prominence over time as the frequency range of measurements correspond to that of field seismic data acquired for oil/gas exploration. These experiments measure attenuation as the phase shift between the applied stress (sinusoidal) and measured strain. Since the magnitudes of measured phase shifts are quite low (Q-1 ~0.01-0.1) and the amplitudes of strain applied to the rock samples are of the order ~10-6 (i.e., similar orders of magnitude to seismic waves), it is challenging. A comparison of such forced oscillation setups will be presented to provide an overview of the various possibilities of design and implementation for future setups. In general, there is a lack of laboratory data and most of the published data are for sandstones. Currently, attenuation measurements are being carried out on carbonate and sandstone samples. We employ the Seismic Wave Attenuation Module (SWAM, Madonna & Tisato 2013) to measure seismic attenuation in these samples for different saturation degrees (90% and 100% water) and under three different confining pressures (5, 10 and 15MPa). Preliminary results from these investigations will be discussed. REFERENCES Madonna, C. & Tisato, N. 2013: A new seismic wave attenuation module to experimentally measure low-frequency attenuation in extensional mode. Geophysical Prospecting, doi: 10.1111/1365-2478.12015. Tisato, N. & Madonna, C. 2012: Attenuation at low seismic frequencies in partially saturated rocks: Measurements and description of a new apparatus. Journal of Applied Geophysics, 86, 44-53.
Nagy, P.B.; Nayfeh, A.H.
1995-09-25
The surface stiffness of a fluid-saturated porous solid is defined as the ratio between a small change in capillary pressure and the average displacement of the boundary due to the resulting rise or fall of the fluid level in the pore channels. When the surface pores are structurally open, the surface stiffness is entirely due to the stiffness of the microscopic fluid membranes extended by capillary forces over the surface pores. Due to interfacial tension between the immiscible wetting fluid in the pores and nonwetting fluid (air) above the surface, essentially closed-pore boundary conditions can prevail at the interface. It has recently been shown that the surface stiffness of a porous material containing cylindrical pores can be calculated simply as the surface tension of the saturating fluid divided by the static permeability of the porous solid [P. B. Nagy, Appl. Phys. Lett. {bold 60}, 2735 (1992)]. In this letter, we show that the same simple relationship can be generalized for the surface stiffness of fluid-saturated porous media containing parallel prismatic pore channels of any number, size, or shape. {copyright} {ital 1995} {ital American} {ital Institute} {ital of} {ital Physics}.
Effect of thermal stratification on free convection in a fluid-saturated porous enclosure
Kumar, B.V.R.; Singh, P.
1998-08-28
Natural convection heat transfer from an isothermal vertical surface in a fluid-saturated porous medium under the influence of thermal stratification has been numerically analyzed using the finite element method. The combined effect of Rayleigh number and thermal stratification on the global heat flux in porous cavities of varying aspect ratios has been analyzed. It is observed that the global heat flux (1) decreases with increasing values of the thermal stratification b and (2) increases with increasing values of Rayleigh number Ra. However, the cumulative heat flux along the midsegment of the vertical wall increases with increasing b and decreases with Ra. With increasing b, a reduction in the convection zone with a downward drift is noticed. Global heat flux can be boosted by reducing the width L of the cavity by a factor of {approximately}8 or more. The flow field and the heat transfer results are presented through the streamlines, isotherms, and global heat flux plots.
On wave propagation characteristics in fluid saturated porous materials by a nonlocal Biot theory
NASA Astrophysics Data System (ADS)
Tong, Lihong; Yu, Yang; Hu, Wentao; Shi, Yufeng; Xu, Changjie
2016-09-01
A nonlocal Biot theory is developed by combing Biot theory and nonlocal elasticity theory for fluid saturated porous material. The nonlocal parameter is introduced as an independent variable for describing wave propagation characteristics in poroelastic material. A physical insight on nonlocal term demonstrates that the nonlocal term is a superposition of two effects, one is inertia force effect generated by fluctuation of porosity and the other is pore size effect inherited from nonlocal constitutive relation. Models for situations of excluding fluid nonlocal effect and including fluid nonlocal effect are proposed. Comparison with experiment confirms that model without fluid nonlocal effect is more reasonable for predicting wave characteristics in saturated porous materials. The negative dispersion is observed theoretically which agrees well with the published experimental data. Both wave velocities and quality factors as functions of frequency and nonlocal parameter are examined in practical cases. A few new physical phenomena such as backward propagation and disappearance of slow wave when exceeding critical frequency and disappearing shear wave in high frequency range, which were not predicted by Biot theory, are demonstrated.
Tao, Chao; Jiang, Jack J.; Czerwonka, Lukasz
2011-01-01
The human vocal fold is treated as a continuous, transversally isotropic, porous solid saturated with liquid. A set of mathematical equations, based on the theory of fluid-saturated porous solids, is developed to formulate the vibration of the vocal fold tissue. As the fluid-saturated porous tissue model degenerates to the continuous elastic tissue model when the relative movement of liquid in the porous tissue is ignored, it can be considered a more general description of vocal fold tissue than the continuous, elastic model. Using the fluid-saturated porous tissue model, the vibration of a bunch of one-dimensional fibers in the vocal fold is analytically solved based on the small amplitude assumption. It is found that the vibration of the tissue will lead to the accumulation of excess liquid in the midmembranous vocal fold. The degree of liquid accumulation is positively proportional to the vibratory amplitude and frequency. The correspondence between the liquid distribution predicted by the porous tissue theory and the location of vocal nodules observed in clinical practice, provides theoretical evidence for the liquid accumulation hypothesis of vocal nodule formation (Jiang, Ph. D., dissertation, 1991, University of Iowa). PMID:19660905
Nelson, J.T. . Dept. of Mechanical Engineering Lawrence Berkeley Lab., CA )
1988-11-01
A theoretical model for acoustic emission in a vertically heterogeneous porous layer bounded by semi-infinite solid regions is developed using linearized equations of motion for a fluid/solid mixture and a reflectivity method. Green's functions are derived for both point loads and moments. Numerically integrated propagators represent solutions for intermediate heterogeneous layers in the porous region. These are substituted into a global matrix for solution by Gaussian elimination and back-substitution. Fluid partial stress and seismic responses to dislocations associated with fracturing of a layer of rock with a hydraulically conductive fracture network are computed with the model. A constitutive model is developed for representing the fractured rock layer as a porous material, using commonly accepted relationships for moduli. Derivations of density, tortuosity, and sinuosity are provided. The main results of the model application are the prediction of a substantial fluid partial stress response related to a second mode wave for the porous material. The response is observable for relatively large distances, on the order of several tens of meters. The visco-dynamic transition frequency associated with parabolic versus planar fluid velocity distributions across micro-crack apertures is in the low audio or seismic range, in contrast to materials with small pore size, such as porous rocks, for which the transition frequency is ultrasonic. Seismic responses are predicted for receiver locations both in the layer and in the outlying solid regions. In the porous region, the seismic response includes both shear and dilatational wave arrivals and a second-mode arrival. The second-mode arrival is not observable outside of the layer because of its low velocity relative to the dilatational and shear wave propagation velocities of the solid region.
NASA Astrophysics Data System (ADS)
Schmidt, Max W.; Ulmer, Peter
2004-04-01
Fluid saturated high-pressure experiments often result in strongly zoned experimental charges, this hinders experimentation in chemically homogeneous systems which in turn has serious consequences on equilibration, reaction progress, and (apparent) phase stabilities. In order to overcome these problems, a 600-ton press accommodating either a multianvil or end-loaded piston cylinder module has been mounted in such a way that it can be turned by 180°, thus inverting its position in the gravity field. During turning, hydraulic pressure, heating power, and cooling water remain connected allowing fully controlled pressures and temperatures during experiments. A series of experiments at 13 GPa, 950°C, on a serpentine bulk composition in the MgO-SiO 2-H 2O system demonstrates that continuous turning at a rate of 2 turns/min results in a nearly homogeneous charge composed of phase E + enstatite. The same experiment at static conditions resulted in four mineral zones: quench phase E, enstatite, enstatite + phase E, and phase E + phase A. Phase A disappears in experiments at a turning rate ≥1 turn/min. A static 15-min experiment shows that zonation already forms within this short time span. Placing two short capsules within a single static experiment reveals that the fluid migrates to the hot spot in each capsule and is not gravitationally driven toward the top. The zonation pattern follows isotherms within the capsule, and the degree of zonation increases with temperature gradient (measured as 10 °C within a capsule) and run time. Our preferred interpretation is that Soret diffusion causes a density-stratified fluid within the capsule that does not convect in a static experiment and results in temperature dependant chemical zonation. The aggravation of zonation and appearance of additional phases with run time can be explained with a dissolution-reprecipitation process where the cold spot of the capsule is relatively MgO enriched and the hot spot relatively SiO 2 and H
Yang, Zhixin; Wang, Shaowei; Zhao, Moli; Li, Shucai; Zhang, Qiangyong
2013-01-01
The onset of double diffusive convection in a viscoelastic fluid-saturated porous layer is studied when the fluid and solid phase are not in local thermal equilibrium. The modified Darcy model is used for the momentum equation and a two-field model is used for energy equation each representing the fluid and solid phases separately. The effect of thermal non-equilibrium on the onset of double diffusive convection is discussed. The critical Rayleigh number and the corresponding wave number for the exchange of stability and over-stability are obtained, and the onset criterion for stationary and oscillatory convection is derived analytically and discussed numerically. PMID:24312193
NASA Astrophysics Data System (ADS)
Song, Yongjia; Hu, Hengshan; Rudnicki, John W.
2016-07-01
Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori-Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.
Thermal instability of a fluid-saturated porous medium bounded by thin fluid layers
Pillatsis, G.; Taslim, M.E.; Narusawa, U. )
1987-08-01
A linear stability analysis is performed for a horizontal Darcy porous layer of depth 2d{sub m} sandwiched between two fluid layers of depth d (each) with the top and bottom boundaries being dynamically free and kept at fixed temperatures. The Beavers-Joseph condition is employed as one of the interfacial boundary conditions between the fluid and the porous layer. The critical Rayleigh number and the horizontal wave number for the onset of convective motion depend on the following four dimensional parameters: {cflx d} (= d{sub m}/d, the depth ratio), {delta} (= {radical}K/d{sub m} with K being the permeability of the porous medium) {alpha} (the proportionality constant in the Beavers-Joseph condition), and k/k{sub m} (the thermal conductivity ratio). In order to analyze the effect of these parameters on the stability condition, a set of numerical solutions is obtained in terms of a convergent series for the respective layers, for the case in which the thickness of the porous layer is much greater than that of the fluid layer. A comparison of this study with the previously obtained exact solution for the case of constant heat flux boundaries is made to illustrate quantitative effects of the interfacial and the top/bottom boundaries on the thermal instability of a combined system of porous and fluid layers.
Imaging of multiphase fluid saturation within a porous material via sodium NMR.
Washburn, Kathryn E; Madelin, Guillaume
2010-01-01
We present in this paper a method to monitor multiphase fluid core saturation through measurement of the sodium NMR signal. In a rock core saturated with water and oil, sodium will be present only in the water phase, and therefore can be used to separate the two fluids. Two dimensional sodium images were taken to monitor the movement of brine into oil saturated rock cores. The measured fluid exchange agrees well with expected behavior from traditional core analysis methods. Indications of damage to the rock structure can be seen from the patterns of fluid imbibition. PMID:19864169
Turbulent flow over a channel with fluid-saturated porous bed
Technology Transfer Automated Retrieval System (TEKTRAN)
The characteristics of fully developed turbulent flow in a hybrid domain channel, which consists of a clear fluid region and a porous bed, are examined numerically using a model based on the macroscopic Reynolds-averaged Navier–Stokes equations. By adopting the classical continuity interface conditi...
Frequency-Dependent Seismic Waves in Fluid-Saturated Fractured Rock
NASA Astrophysics Data System (ADS)
Korneev, V. A.; Goloshubin, G.
2015-12-01
Fractures are the natural and essential elements of rock. Fracture systems are the most important features that define rock permeability and strength, as well as their anisotropy properties. Recent advancement in induced fracturing is a core part of the gas/oil shale technology, where fracture monitoring and control became a special topic of interest. Krauklis wave (K-wave) is the result of interaction between a fluid mass and elasticity of fracture walls, and it propagates primarily along the fracture systems in the fluid. At the fracture tips and fracture intersections it partially converts into the body waves. It is quite clear that incorporation of K-waves in a theory of wave propagation in fractured rock is one of the most important problems to solve for understanding of their seismic properties. One of the most fundamental properties of fractured rock is a fractal fracture distribution and it is rarely, if ever, taken into account in existing wave propagation theories. However, this property exists on a widest variety of scales and in particular reveals itself in a form of Gutenberg-Richter Law experimentally proven, starting from laboratory measurements and up to the global seismicity. We computed P and S-wave velocities of the rock containing fluid (and proppant) filled fractures, considering the effect of extremely slow and dispersive wave propagation within individual fractures. This was made possible by introducing the concept of "effective fracture-wave volume," and by evaluating the elastic constants of rock containing a complex, fractal network of fractures. These velocities were used to compute seismic waves reflected normally from a fractured reservoir. We demonstrate that by taking into account the Krauklis wave phenomenon for the fractally distributed fluid-filled fractures, it is possible to explain the observed low-frequency anomalies above the underground natural reservoirs. These anomalies include increase of amplitude and a phase delay of
NASA Astrophysics Data System (ADS)
Taleb, A.; BenHamed, H.; Ouarzazi, M. N.; Beji, H.
2016-05-01
We report theoretical and numerical results on bifurcations in thermal instability for a viscoelastic fluid saturating a porous square cavity heated from below. The modified Darcy law based on the Oldroyd-B model was used for modeling the momentum equation. In addition to Rayleigh number ℜ, two more dimensionless parameters are introduced, namely, the relaxation time λ1 and the retardation time λ2. Temporal stability analysis showed that the first bifurcation from the conductive state may be either oscillatory for sufficiently elastic fluids or stationary for weakly elastic fluids. The dynamics associated with the nonlinear interaction between the two kinds of instabilities is first analyzed in the framework of a weakly nonlinear theory. For sufficiently elastic fluids, analytical expressions of the nonlinear threshold above which a second hysteretic bifurcation from oscillatory to stationary convective pattern are derived and found to agree with two-dimensional numerical simulations of the full equations. Computations performed with high Rayleigh number indicated that the system exhibits a third transition from steady single-cell convection to oscillatory multi-cellular flows. Moreover, we found that an intermittent oscillation regime may exist with steady state before the emergence of the secondary Hopf bifurcation. For weakly elastic fluids, we determined a second critical value ℜ2 Osc ( λ 1 , λ 2 ) above which a Hopf bifurcation from steady convective pattern to oscillatory convection occurs. The well known limit of ℜ2 Osc ( λ 1 = 0 , λ 2 = 0 ) = 390 for Newtonian fluids is recovered, while the fluid elasticity is found to delay the onset of the Hopf bifurcation. The major new findings were presented in the form of bifurcation diagrams as functions of viscoelastic parameters for ℜ up to 420.
NASA Astrophysics Data System (ADS)
RamReddy, Chetteti; Pradeepa, Teegala
2016-05-01
Based on the nonlinear variation of density with temperature (NDT) in the buoyancy term, the mixed convection flow along a vertical plate of a micropolar fluid saturated porous medium is considered. In addition, the effect of homogeneous-heterogeneous reaction and convective boundary condition has been taken into account. Using lie scaling group transformations, the similarity representation is attained for the system of partial differential equations, prior to being solved by a spectral quasilinearization method. The results show that in the presence of aiding and opposing flow situations, both the species concentration and mass transfer rate decreases when the strength of homogeneous and heterogeneous reaction parameters are enhanced.
Zhao, C.; Muehlaus, H.B.; Hobbs, B.E.
1998-03-01
A parametric study is carried out to investigate how geological inhomogeneity affects the pore-fluid convective flow field, the temperature distribution, and the mass concentration distribution in a fluid-saturated porous medium. The related numerical results have demonstrated that (1) the effects of both medium permeability inhomogeneity and medium thermal conductivity inhomogeneity are significant on the pore-fluid convective flow and the species concentration distribution in the porous medium; (2) the effect of medium thermal conductivity inhomogeneity is dramatic on the temperature distribution in the porous medium, but the effect of medium permeability inhomogeneity on the temperature distribution may be considerable, depending on the Rayleigh number involved in the analysis; (3) if the coupling effect between pore-fluid flow and mass transport is weak, the effect of the Lewis number is negligible on the pore-fluid convective flow and temperature distribution, but it is significant on the species concentration distribution in the medium.
Dilatant hardening of fluid-saturated sandstone
NASA Astrophysics Data System (ADS)
Makhnenko, Roman Y.; Labuz, Joseph F.
2015-02-01
The presence of pore fluid in rock affects both the elastic and inelastic deformation processes, yet laboratory testing is typically performed on dry material even though in situ the rock is often saturated. Techniques were developed for testing fluid-saturated porous rock under the limiting conditions of drained, undrained, and unjacketed response. Confined compression experiments, both conventional triaxial and plane strain, were performed on water-saturated Berea sandstone to investigate poroelastic and inelastic behavior. Measured drained response was used to calibrate an elasto-plastic constitutive model that predicts undrained inelastic deformation. The experimental data show good agreement with the model: dilatant hardening in undrained triaxial and plane strain compression tests under constant mean stress was predicted and observed.
NASA Astrophysics Data System (ADS)
Kubik, J.; Cieszko, M.
2005-12-01
The compatibility conditions matching macroscopic mechanical fields at the contact surface between a fluid-saturated porous solid and an adjacent bulk fluid are considered. The general form of balance equations at that discontinuity surface are analyzed to obtain the compatibility conditions for the tangent and normal components of the velocity and the stress vector fields. Considerations are based on the procedure similar to that used in the phenomenological thermodynamics for derivation of constitutive relations, where the entropy inequality and the concept of Lagrange multipliers are applied. This procedure made possible to derive the compatibility conditions for the viscous fluid flowing tangentially and perpendicularly to the boundary surface of the porous solid and to formulate the generalized form of the so called slip condition for the fluid velocity field, postulated earlier by Beavers and Joseph, J. Fluid. Mech. 30, 197-207 (1967).
NASA Astrophysics Data System (ADS)
Labrousse, L.; Duretz, T.; Gerya, T.
2015-10-01
We present two-dimensional numerical models of plate subduction and collision inspired by the Scandinavian Caledonian orogeny to investigate the possible impact of continental crust partial melting on the exhumation of ultra-high pressure metamorphic rocks. Three possible reactions were tested: low temperature solidus representing H2O-fluid-saturated partial melting, and two end-member reaction curves for dehydration melting. Thermo-mechanical effects of partial melting were implemented as (1) a viscosity decrease as a determined rheologically critical melt percentage was reached (here 0.1), (2) a change in effective heat capacity and adiabatic heating/cooling accounting for a latent heat term in the heat equation. Among the 3 tested reactions, only H2O-fluid-saturated partial melting drastically modifies the collision dynamics from the non-melting reference model holding all other parameters constant. A substantially low general viscosity truncation (here 1017 Pa s) is needed to properly resolve the effect of partial melting on deep collision processes. Low temperature melting indeed induces the development of a low viscosity buoyant plume prior to slab detachment, where migmatites exhume from UHP conditions at rates and with pressure-temperature paths similar to the natural values acknowledged for the Norwegian Caledonides. High temperature melting has no drastic influence on early collision dynamics. While positive buoyancy remains the first order driver for the exhumation of buried continental rocks, exhumation initiates in these cases with eduction subsequent to slab detachment. Melting and formation of a migmatite plume can later occur along decompression path while continental crust undergoes thermal reequilibration at temperatures above 900 °C. Some of the partially molten material can also relaminate in the overriding plate rather than exhume within the collision zone. Even if minor in terms of amount of magma produced, H2O-fluid-saturated partial melting
NASA Astrophysics Data System (ADS)
Mehta, C. B.; Singh, M.; Kumar, S.
2016-02-01
An investigation is made on the effect of Hall currents on thermal instability of a compressible couple-stress fluid in the presence of a horizontal magnetic field saturated in a porous medium. The analysis is carried out within the framework of the linear stability theory and normal mode technique. A dispersion relation governing the effects of viscoelasticity, Hall currents, compressibility, magnetic field and porous medium is derived. For the stationary convection a couple-stress fluid behaves like an ordinary Newtonian fluid due to the vanishing of the viscoelastic parameter. Compressibility, the magnetic filed and couple-stress parameter have stabilizing effects on the system whereas Hall currents and medium permeability have a destabilizing effect on the system, but in the absence of Hall current couple-stress has a destabilizing effect on the system. It has been observed that oscillatory modes are introduced due to the presence of viscoelasticity, magnetic field porous medium and Hall currents which were non-existent in their absence.
Biogenic Cracks in Porous Rock
NASA Astrophysics Data System (ADS)
Hemmerle, A.; Hartung, J.; Hallatschek, O.; Goehring, L.; Herminghaus, S.
2014-12-01
Microorganisms growing on and inside porous rock may fracture it by various processes. Some of the mechanisms of biofouling and bioweathering are today identified and partially understood but most emphasis is on chemical weathering, while mechanical contributions have been neglected. However, as demonstrated by the perseverance of a seed germinating and cracking up a concrete block, the turgor pressure of living organisms can be very significant. Here, we present results of a systematic study of the effects of the mechanical forces of growing microbial populations on the weathering of porous media. We designed a model porous medium made of glass beads held together by polydimethylsiloxane (PDMS), a curable polymer. The rheological properties of the porous medium, whose shape and size are tunable, can be controlled by the ratio of crosslinker to base used in the PDMS (see Fig. 1). Glass and PDMS being inert to most chemicals, we are able to focus on the mechanical processes of biodeterioration, excluding any chemical weathering. Inspired by recent measurements of the high pressure (~0.5 Mpa) exerted by a growing population of yeasts trapped in a microfluidic device, we show that yeast cells can be cultured homogeneously within porous medium until saturation of the porous space. We investigate then the effects of such an inner pressure on the mechanical properties of the sample. Using the same model system, we study also the complex interplay between biofilms and porous media. We focus in particular on the effects of pore size on the penetration of the biofilm within the porous sample, and on the resulting deformations of the matrix, opening new perspectives into the understanding of life in complex geometry. Figure 1. Left : cell culture growing in a model porous medium. The white spheres represent the grains, bonds are displayed in grey, and microbes in green. Right: microscopy picture of glass beads linked by PDMS bridges, scale bar: 100 μm.
NASA Astrophysics Data System (ADS)
Chand, Ramesh
2015-12-01
Thermal instability in a horizontal layer of Oldroydian visco-elastic fluid in a porous medium is investigated. For porous medium the Brinkman-Darcy model is considered. A linear stability analysis based upon perturbation method and normal mode technique is used to find solution of the fluid layer confined between two free-free boundaries. The onset criterion for stationary and oscillatory convection is derived analytically. The influence of the Brinkman-Darcy, Prandtl-Darcy number, stress relaxation parameter on the stationary and oscillatory convection is studied both analytically and graphically. The sufficient condition for the validity of PES has also been derived.
NASA Astrophysics Data System (ADS)
Pal, Dulal; Chatterjee, Sewli
2011-03-01
A numerical model is developed to examine the combined effects of Soret and Dufour on mixed convection magnetohydrodynamic heat and mass transfer in micropolar fluid-saturated Darcian porous medium in the presence of thermal radiation, non-uniform heat source/sink and Ohmic dissipation. The governing boundary layer equations for momentum, angular momentum (microrotation), energy and species transfer are transformed to a set of non-linear ordinary differential equations by using similarity solutions which are then solved numerically based on shooting algorithm with Runge-Kutta-Fehlberg integration scheme over the entire range of physical parameters with appropriate boundary conditions. The influence of Darcy number, Prandtl number, Schmidt number, Soret number and Dufour number, magnetic parameter, local thermal Grashof number and local solutal Grashof number on velocity, temperature and concentration fields are studied graphically. Finally, the effects of related physical parameters on local Skin-friction, local Nusselt number and local Sherwood number are also studied. Results showed that the fields were influenced appreciably by the Soret and Dufour effects, thermal radiation and magnetic field, etc.
NASA Astrophysics Data System (ADS)
Hsu, S. Y.; Tsai, J. P.; Chang, L. C.
2014-12-01
The flow of three immiscible fluids - water, NAPL, air - in porous media is important in many subsurface processes. To model the three-fluid flow, the relation of relative permeability-saturation-capillary pressure (k-S-P) of three fluids is of central importance. In this experimental study, we directly measure the k-S-P of the water (wetting phase) when three fluids are coexist in a micromodel during the water drainage and imbibition. The results show that the sequence of the non-wetting fluids (air and NAPL) entering into the micromodel affects the fluid distributions as well as the relative permeability of water. During the drainage process, the relative permeability of water dropped drastically when the pathway of water from inlet to outlet of the micromodel was visually blocked by the non-wetting fluids. At this stage, the relative permeability of water was low but not down to zero. The water was still able to move via corner flows or thin-film flows. During the imbibition process, the water displaced two non-wetting liquids via both "snap-off" and "piston-type" motions. The relative permeability of water jumped when the water pathway was formed again. In addition, we found that the well-known scaling format proposed by Parker et al. [1] might fail when the interfaces between the most non-wetting (air) and the most wetting (water) fluids occurs in the three-fluids system. References[1] J. C. Parker, R. J. Lenhard, and T. Kuppusamy, Water Resources Research, 23, 4, 618-624 (1987)
Experimental studies of electrokinetic conversions in fluid-saturated borehole models
Zhu, Z.; Haartsen, M.W.; Toksoez, M.N.
1999-10-01
Experimental and theoretical studies show that there are electromagnetic (EM) fields generated by seismic waves with two kinds of conversion mechanisms in a fluid-saturated, porous medium. Within a homogeneous formation, the seismic wave generates a seismoelectric field that exists only in the area disturbed by the seismic wave and whose apparent velocity is that of the seismic wave. At an interface between differing formation properties, the generated seismoelectric wave is a propagating EM wave that can be detected everywhere, An electrode, used as a receiver on the ground surface, can detect the propagating EM wave generated at an interface, but cannot detect the seismoelectric field generated in a homogeneous formation. When the electrode is in a borehole and close to a porous formation, it can detect both the EM waves and the seismoelectric field. In this paper, electrokinetic measurements are performed with borehole models made of natural rocks or artificial materials. Experimental results show that the Stoneley wave and other acoustic modes, excited by a monopole source in the borehole models, generate seismoelectric fields in fluid-saturated formations. The electric components of the seismoelectric fields can be detected by an electrode in the borehole or on the borehole wall. The amplitude and frequency of the seismoelectric fields are related not only to the seismic wave, but also to formation properties such as permeability, conductivity, etc. Comparison between the waveforms of the seismoelectric signals and acoustic logging waves suggests that seismoelectric well logging may explore the different properties of the formation. Electroseismic measurements are also performed with these borehole models. The electric pulse through the electrode in the borehole or on the borehole wall induces Stoneley waves in fluid-saturated models that can be received by a monopole transducer in the same borehole. These measurement methods (seismoelectric logging or
Dispersion of elastic moduli in a porous-cracked rock: Theoretical predictions for squirt-flow
NASA Astrophysics Data System (ADS)
Adelinet, M.; Fortin, J.; Guéguen, Y.
2011-04-01
Crustal rocks contain variable amount of both cracks and equant pores depending on tectonic and thermal stresses but also on their geological origin. Crack damage and porosity change result in effects on elastic waves velocities. When rocks are fluid saturated, dispersion of the P- and S-waves should be taken into account. This paper deals with frequency dispersion of elastic moduli in a fluid saturated porous and cracked rock with the assumption that squirt-flow is the dominant process. We develop a theoretical approach to calculate both high (HF) and low (LF) frequency bulk and shear moduli. The HF moduli are derived from a new effective medium model, called CPEM, with an isotropic distribution of pores or cracks with idealized geometry, respectively spheres and ellipsoids. LF moduli are obtained by taking HF dry moduli from the CPEM and substituting into Gassmann's equations. In the case of a porosity only supported by equant pores, the calculated dispersion in elastic moduli is equal to zero. In the case of a crack porosity, no bulk dispersion is predicted but a shear dispersion appears. Finally in the general case of a mixed porosity (pores and cracks), dispersion in bulk and in shear is predicted. Our results show that the maximum dispersion is predicted for a mixture of pores and spheroidal cracks with a very small aspect ratio (≤ 10 - 3 ). Our theoretical predictions are compared to experimental data obtained during hydrostatic experiment performed on a basaltic rock and a good agreement is observed. We also used our theoretical model to predict elastic waves velocities and Vp/Vs ratio dispersion. We show that the P-waves dispersion can reach almost 20% and the Vp/Vs dispersion a maximum value of 9% for a crack porosity of about 1%. Since laboratory data are ultrasonic measurements and field data are obtained at much lower frequencies, these results are useful for geophysicists to interpret seismic data in terms of fluid and rock interactions.
Theoretical and numerical aspects of fluid-saturated elasto-plastic soils
Ehlers, W.
1995-12-31
The theoretical and numerical treatment of fluid-saturated porous solid materials generally falls into the category of porous media models, which are described within the framework of the classical theory of mixtures extended by the concept of volume fractions (porous media theories). In particular, this concept allows for the description of saturated, unsaturated and empty porous matrix materials, thus offering a well-founded theoretical background for a lot of engineering problems occurring, for instance, in the fields of geomechanics (soil and rock mechanics as well as glacier and rock ice mechanics), oil producing industries, sintering technologies, biomechanics, etc. In the present contribution, theoretical and numerical studies are outlined to describe a two-phase material composed of an incompressible elasto-plastic soil matrix saturated by an incompressible viscous pore fluid. In this context, the phenomenon of phase incompressibility is well known as a microscopic effect not implying bulk incompressibility in the macro regime. This is seen from the fact that even if the material density functions of the individual constituents are constant during deformation, the corresponding bulk densities can still change through changes in the volume fractions. Within the framework of a pure mechanical theory, constitutive equations are given for both the solid and the fluid partial stress tensors and for the interaction force acting between the two materials. Concerning the porous soil matrix, the elastic properties are described by an elasticity law of Hookean type, while the plastic range is governed by a {open_quote}single surface{close_quote} yield function exhibiting a smooth and closed shape in the principal stress space together with a non-associated flow rule. The viscosity effects of the pore fluid are included in the fluid stress tensor and in the drag force.
Pressure diffusion waves in porous media
Silin, Dmitry; Korneev, Valeri; Goloshubin, Gennady
2003-04-08
Pressure diffusion wave in porous rocks are under consideration. The pressure diffusion mechanism can provide an explanation of the high attenuation of low-frequency signals in fluid-saturated rocks. Both single and dual porosity models are considered. In either case, the attenuation coefficient is a function of the frequency.
Computational Modeling of Seismic Wave Propagation Velocity-Saturation Effects in Porous Rocks
NASA Astrophysics Data System (ADS)
Deeks, J.; Lumley, D. E.
2011-12-01
Compressional and shear velocities of seismic waves propagating in porous rocks vary as a function of the fluid mixture and its distribution in pore space. Although it has been possible to place theoretical upper and lower bounds on the velocity variation with fluid saturation, predicting the actual velocity response of a given rock with fluid type and saturation remains an unsolved problem. In particular, we are interested in predicting the velocity-saturation response to various mixtures of fluids with pressure and temperature, as a function of the spatial distribution of the fluid mixture and the seismic wavelength. This effect is often termed "patchy saturation' in the rock physics community. The ability to accurately predict seismic velocities for various fluid mixtures and spatial distributions in the pore space of a rock is useful for fluid detection, hydrocarbon exploration and recovery, CO2 sequestration and monitoring of many subsurface fluid-flow processes. We create digital rock models with various fluid mixtures, saturations and spatial distributions. We use finite difference modeling to propagate elastic waves of varying frequency content through these digital rock and fluid models to simulate a given lab or field experiment. The resulting waveforms can be analyzed to determine seismic traveltimes, velocities, amplitudes, attenuation and other wave phenomena for variable rock models of fluid saturation and spatial fluid distribution, and variable wavefield spectral content. We show that we can reproduce most of the published effects of velocity-saturation variation, including validating the Voigt and Reuss theoretical bounds, as well as the Hill "patchy saturation" curve. We also reproduce what has been previously identified as Biot dispersion, but in fact in our models is often seen to be wave multi-pathing and broadband spectral effects. Furthermore, we find that in addition to the dominant seismic wavelength and average fluid patch size, the
Not Available
1994-06-01
Calculations showed that capillary forces can easily produce closed- pore boundary conditions at interface between nonwetting fluid (air) and a porous solid saturated by a wetting fluid (water). The direct excitation technique was used to measure surface wave velocity and attenuation on both wet and dry rocks. The strong correlation between the observed surface wave velocity change caused by water saturation and the formation permeability can be used for ultrasonic assessment of the dynamic permeability. The experimental system was improved further by introducing laser interferometric detection, which was adapted to surface wave inspection of fluid-saturated permeable materials. In a separate effort, the surface stiffness of different water-saturated porous solids was studied by a novel acoustical method. Areas for further study are described.
Goloshubin, Gennady M.; Korneev, Valeri A.
2005-09-06
A method for identifying, imaging and monitoring dry or fluid-saturated underground reservoirs using seismic waves reflected from target porous or fractured layers is set forth. Seismic imaging the porous or fractured layer occurs by low pass filtering of the windowed reflections from the target porous or fractured layers leaving frequencies below low-most corner (or full width at half maximum) of a recorded frequency spectra. Additionally, the ratio of image amplitudes is shown to be approximately proportional to reservoir permeability, viscosity of fluid, and the fluid saturation of the porous or fractured layers.
Goloshubin, Gennady M.; Korneev, Valeri A.
2006-11-14
A method for identifying, imaging and monitoring dry or fluid-saturated underground reservoirs using seismic waves reflected from target porous or fractured layers is set forth. Seismic imaging the porous or fractured layer occurs by low pass filtering of the windowed reflections from the target porous or fractured layers leaving frequencies below low-most corner (or full width at half maximum) of a recorded frequency spectra. Additionally, the ratio of image amplitudes is shown to be approximately proportional to reservoir permeability, viscosity of fluid, and the fluid saturation of the porous or fractured layers.
Dynamic flow localization in porous rocks under combined pressure and shear loading
NASA Astrophysics Data System (ADS)
Yarushina, Viktoriya; Podladchikov, Yuri; Simon, Nina
2015-04-01
Flow localization occurs in deforming porous fluid saturated rocks. It exhibits itself as veins, pockmarks on the ocean floor or gas chimneys visible on seismic images from several chalk fields of the Central North Sea and from the Utsira formation at Sleipner in the Norwegian North Sea, which is one of the best documented CO2 storage sites. Porosity waves were repeatedly shown to be a viable mechanism of flow self-localization that does not require the pre-existence of a connected fracture network. Porosity waves result from an instability of the Darcy flow that occurs in porous rocks with time-dependent viscous or viscoelastoplastic rheology. Local fluid overpressure generated by fluid injection or chemical reactions aided by buoyancy force drives upward fluid migration. Viscous deformation delays pressure diffusion thus maintaining local overpressure for considerable periods of time. Development of an under-pressured region just below the over-pressured domain leads to separation of the fluid-filled high-porosity blob from the source and the background flow. The instability organizes the flow into separate vertical channels. Pressure distribution, shape and scaling of these channels are highly sensitive to the rheology of the porous rock. In this contribution, based on a micromechanical approach, we consider the complex rheology of brittle, ductile and transitional regimes of deformation of porous rocks in the presence of combined pressure and shear loading. Accurate description of transitional brittle-ductile deformation is a challenging task due to a large number of microscopic processes involved. We use elastoplastic and viscoplastic analytical solutions for the non-hydrostatic deformation of a singular cavity in the representative volume element in order to deduce expected behavior of the porous rock. The model provides micro-mechanisms for various failure modes (localized and homogeneous) and dilatancy onset. In particular, the model predicts that dilatancy
Reaction-driven fracturing of porous rock
NASA Astrophysics Data System (ADS)
Ulven, O. I.; Jamtveit, B.; Malthe-Sørenssen, A.
2014-10-01
A 2-D computer model has been developed to investigate fluid-mediated transformation processes such as chemical weathering, mineral carbonation, and serpentinization that require transport of H2O and/or CO2 into reacting rock volumes. Hydration and carbonation cause local volume expansion, and the resulting nonuniform stresses may drive fracturing, which increases both the rate of transport and the accessible reactive surface area in the system and thus accelerates the rate of the transformation process. The model couples reactions, fracturing, and fluid transport for systems with a range of initial porosities, assumed constant throughout the process. With low initial porosity, a sharp reaction front between completely reacted material and unreacted material propagates into unaltered rock, while for high porosities, diffuse reaction fronts are formed in which a large fraction of the initial volume is partly reacted. When diffusive transport is rate limiting, the total reaction rate depends on porosity to a power N, where N is in the range 0.45-2. The exponent N increases as the reaction-generated expansion decreases. In high-porosity rocks, the total reaction rate is limited by reaction kinetics, and it is thus insensitive to porosity variations. As the volume increasing reaction proceeds, fracturing divides the unreacted porous material into subdomains, which may undergo further subdivision as they are consumed by the reaction. The total reaction rate and progress depend on the initial geometry of a reacting domain, and this significantly affects the weathering profiles for systems that evolve from an initial assembly of blocks with different sizes and shapes.
NASA Astrophysics Data System (ADS)
Ding, Pinbo; Di, Bangrang; Wei, Jianxin; Li, Xiangyang; Deng, Yinghua
2014-02-01
In this study, we analyse the influence of fluid on P- and S-wave anisotropy in a fractured medium. Equivalent medium theories are used to describe the relationship between the fluid properties and the rock physics characteristics in fractured rocks, and P-wave and S-wave velocities and anisotropy are considered to be influenced by fluid saturation. However, these theoretical predictions require experimental measurement results for calibration. A new construction method was used to create synthetic rock samples with controlled fracture parameters. The new construction process provides synthetic rocks that have a more realistic mineral composition, porous structure, cementation and pressure sensitivity than samples used in previous research on fractured media. The synthetic rock samples contain fractures which have a controlled distribution, diameter, thickness and fracture density. In this study, the fracture diameter was about 4 mm, the thickness of fractures was about 0.06 mm, and the fracture density in the two fractured rock samples was about 3.45%. SEM images show well-defined penny-shaped fractures of 4 mm in length and 0.06 mm in width. The rock samples were saturated with air, water and oil, and P- and S-wave velocities were measured in an ultrasonic measurement system. The laboratory measurement results show that the P-wave anisotropy is strongly influenced by saturated fluid, and the P-wave anisotropy parameter, ɛ, has a much larger value in air saturation than in water and oil saturations. The S-wave anisotropy decreases when the samples are saturated with oil, which can be caused by high fluid viscosity. In the direction perpendicular to the fractures (the 0° direction), shear-wave splitting is negligible, and is similar to the blank sample without fractures, as expected. In the direction parallel to the fractures (the 90° direction) shear-wave splitting is significant. The fractured rock samples show significant P- and S-wave anisotropy caused by
Magma degassing during eruption through water-saturated porous rocks
NASA Astrophysics Data System (ADS)
Melnik, O. E.; Afanasyev, A. A.; Zarin, G. A.
2016-05-01
In the case of extrusive eruption, we consider the problem on magma degassing which rises in a volcano conduit crossing porous water-saturated rocks. We show that the intensity of outflow of volcanic gases into the rocks is comparable to the intensity of their transport with the rising magma. The magma degassing in the rocks substantially affects the eruption dynamics, in particular, the duration of the periods of eruptive activity.
Seismic attenuation and velocity dispersion in heterogeneous partially saturated porous rocks
NASA Astrophysics Data System (ADS)
Rubino, J. Germán.; Holliger, Klaus
2012-03-01
Using a numerical approach, we explore wave-induced fluid flow effects in partially saturated porous rocks in which the gas-water saturation patterns are governed by mesoscopic heterogeneities associated with the dry frame properties. The link between the dry frame properties and the gas saturation is defined by the assumption of capillary pressure equilibrium, which in the presence of heterogeneity implies that neighbouring regions can exhibit different levels of saturation. To determine the equivalent attenuation and phase velocity of the synthetic rock samples considered in this study, we apply a numerical upscaling procedure, which permits to take into account mesoscopic heterogeneities associated with the dry frame properties as well as spatially continuous variations of the pore fluid properties. The multiscale nature of the fluid saturation is taken into account by locally computing the physical properties of an effective fluid, which are then used for the larger-scale simulations. We consider two sets of numerical experiments to analyse such effects in heterogeneous partially saturated porous media, where the saturation field is determined by variations in porosity and clay content, respectively. In both cases we also evaluate the seismic responses of corresponding binary, patchy-type saturation patterns. Our results indicate that significant attenuation and modest velocity dispersion effects take place in this kind of media for both binary patchy-type and spatially continuous gas saturation patterns and in particular in the presence of relatively small amounts of gas. The numerical experiments also show that the nature of the gas distribution patterns is a critical parameter controlling the seismic responses of these environments, since attenuation and velocity dispersion effects are much more significant and occur over a broader saturation range for binary patchy-type gas-water distributions. This analysis therefore suggests that the physical mechanisms
Blast energy mitigation in porous rocks
NASA Astrophysics Data System (ADS)
Essink, Brittany C.
Geo-materials are commonly used and sought after for blast mitigation applications due to their wide availability and low cost compared to industry trademarked materials. Characterization of these natural geo-materials such as volcanic rocks is of paramount importance in determining their blast mitigation capabilities. While there is a large amount of information available for materials such as concrete or sand blasts, information on the properties of volcanic rocks is far more scarce. This lack of data is due to the wide range of existing natural volcanic rocks and the variation in the minerals and pore structures of the rocks. In this thesis, silicate volcanic rock samples are characterized both through static and dynamic experimental methods. Initial X-ray powder diffraction scans have been conducted and analyzed to obtain the mineral composition information of the rock samples. Additional tomographic scans under quasi-static loading have been recorded to better understand the internal composition of the material pore structure and the material fracture. For this study, standard compression experiments were conducted at two separate strain rates for ten samples each on a UTM test frame to characterize the behavior of the rock under quasi-static conditions. High strain rate uniaxial compression tests were conducted for three strain rates using a split-Hopkinson pressure bar with pulse shaping to determine the dynamic response of the material. The stress-strain data from the experiments was used to determine the modulus of toughness of the material. Due to the high porosity and heterogeneity of the material, 25 samples were used for dynamic experimentation to attempt to capture and minimize the effects of scatter in the natural material. High speed photography was used to capture the sample deformation during two separate strain rates and to visualize crack propagation and strain rate in the samples. It was found that after an initial yielding, the material is
Seismic stress mobilization of natural colloids in a porous rock
Roberts, Peter M; Abdel-fattah, Amr I
2008-01-01
Stress oscillations at 26 Hz enhanced the release of natural micro-particles (colloids) in a porous rock sample. Micron-scale effects were induced by meter-scale wavelengths. The results are attributed to altering the release rate coefficient for colloids trapped in pores. The rate change did not depend on colloid size and thus is not due to altering colloid-pore-wall interactions. Enhanced colloid detachment from pore walls and flushing from dead-end pores are likely mechanisms. This phenomenon could impact a broad range of physical sciences involving colloid dynamics and porous transport.
NASA Astrophysics Data System (ADS)
Ulven, Ole Ivar; Sun, WaiChing
2016-04-01
Fluid transport in a porous medium has important implications for understanding natural geological processes. At a sufficiently large scale, a fluid-saturated porous medium can be regarded as a two-phase continuum, with the fluid constituent flowing in the Darcian regime. Nevertheless, a fluid mediated chemical reaction can in some cases change the permeability of the rock locally: Mineral dissolution can cause increased permeability, whereas mineral precipitation can reduce the permeability. This might trigger a complicated hydro-chemo-mechanical coupling effect that causes channeling of fluids or clogging of the system. If the fluid is injected or produced at a sufficiently high rate, the pressure might increase enough to cause the onset and propagation of fractures. Fractures in return create preferential flow paths that enhance permeability, localize fluid flow and chemical reaction, prevent build-up of pore pressure and cause anisotropy of the hydro-mechanical responses of the effective medium. This leads to a complex coupled process of solid deformation, chemical reaction and fluid transport enhanced by the fracture formation. In this work, we develop a new coupled numerical model to study the complexities of feedback among fluid pressure evolution, fracture formation and permeability changes due to a chemical process in a 2D system. We combine a discrete element model (DEM) previously used to study a volume expanding process[1, 2] with a new fluid transport model based on poroelasticity[3] and a fluid-mediated chemical reaction that changes the permeability of the medium. This provides new insights into the hydro-chemo-mechanical process of a transforming porous medium. References [1] Ulven, O. I., Storheim, H., Austrheim, H., and Malthe-Sørenssen, A. "Fracture Initiation During Volume Increasing Reactions in Rocks and Applications for CO2 Sequestration", Earth Planet. Sc. Lett. 389C, 2014a, pp. 132 - 142, doi:10.1016/j.epsl.2013.12.039. [2] Ulven, O. I
Quantifying Biofilm in Porous Media Using Rock Physics Models
NASA Astrophysics Data System (ADS)
Alhadhrami, F. M.; Jaiswal, P.; Atekwana, E. A.
2012-12-01
Biofilm formation and growth in porous rocks can change their material properties such as porosity, permeability which in turn will impact fluid flow. Finding a non-intrusive method to quantify biofilms and their byproducts in rocks is a key to understanding and modeling bioclogging in porous media. Previous geophysical investigations have documented that seismic techniques are sensitive to biofilm growth. These studies pointed to the fact that microbial growth and biofilm formation induces heterogeneity in the seismic properties. Currently there are no rock physics models to explain these observations and to provide quantitative interpretation of the seismic data. Our objectives are to develop a new class of rock physics model that incorporate microbial processes and their effect on seismic properties. Using the assumption that biofilms can grow within pore-spaces or as a layer coating the mineral grains, P-wave velocity (Vp) and S-wave (Vs) velocity models were constructed using travel-time and waveform tomography technique. We used generic rock physics schematics to represent our rock system numerically. We simulated the arrival times as well as waveforms by treating biofilms either as fluid (filling pore spaces) or as part of matrix (coating sand grains). The preliminary results showed that there is a 1% change in Vp and 3% change in Vs when biofilms are represented discrete structures in pore spaces. On the other hand, a 30% change in Vp and 100% change in Vs was observed when biofilm was represented as part of matrix coating sand grains. Therefore, Vp and Vs changes are more rapid when biofilm grows as grain-coating phase. The significant change in Vs associated with biofilms suggests that shear velocity can be used as a diagnostic tool for imaging zones of bioclogging in the subsurface. The results obtained from this study have significant implications for the study of the rheological properties of biofilms in geological media. Other applications include
Computed microtomography and numerical study of porous rock samples
NASA Astrophysics Data System (ADS)
Bielecki, J.; Jarzyna, J.; Bożek, S.; Lekki, J.; Stachura, Z.; Kwiatek, W. M.
2013-12-01
Measurement of physical properties of porous geological materials is a crucial issue in oil and gas recovery industry. A conventional experimental way to obtain information on porosity, pore size distribution, specific surface area and permeability are the intrusion porosimetry and permeameter measurement. However, in this approach Washburn's equation is usually used, thus approximation of cylindrical pore shapes is made. However, in recent times the computed microtomography (CMT) technique is more widely used in geoscience (Appoloni et al., 2007). We have already reported preliminary results of investigation of elemental content, microporosity, and specific surface area of porous rocks by means of the CMT technique based on a laboratory source (Bielecki et al., 2012). In this paper, results of complex study of porous rock samples with the use of X-ray CMT (laboratory-source-based facility and synchrotron radiation source) combined with permeability tensor computation by means of the Lattice Boltzmann Method (LBM) are presented. Moreover, the proton induced X-ray emission (PIXE) method was additionally used for elemental content determination of pores-filling substance.
Method of determining interwell oil field fluid saturation distribution
Donaldson, Erle C.; Sutterfield, F. Dexter
1981-01-01
A method of determining the oil and brine saturation distribution in an oil field by taking electrical current and potential measurements among a plurality of open-hole wells geometrically distributed throughout the oil field. Poisson's equation is utilized to develop fluid saturation distributions from the electrical current and potential measurement. Both signal generating equipment and chemical means are used to develop current flow among the several open-hole wells.
An Experimental Study and Constitutive Modeling of Saturated Porous Rocks
NASA Astrophysics Data System (ADS)
Xie, S. Y.; Shao, J. F.
2015-01-01
This paper is devoted to the experimental characterization and constitutive modeling of saturated porous rocks. A typical porous chalk is investigated. Drained hydrostatic and triaxial compression tests are first performed to characterize the basic mechanical behavior of chalk. Drained triaxial tests with constant interstitial pressure are then carried out to study the effects of interstitial pressure on the plastic deformation and failure criterion. Finally, undrained triaxial compression tests are performed to investigate poromechanical coupling in saturated conditions. Based on the experimental data and some relevant micromechanical considerations, a micromechanics-based plastic model is proposed and extended to poroplastic coupling using the effective stress concept. The proposed model is verified through comparisons between the numerical results and experimental data for both drained and undrained tests.
Harmut Spetzler
2005-11-28
This paper describes the culmination of a research project in which we investigated the complex modulus change in partially fluid saturated porous rocks. The investigation started with simple flow experiments over ''clean'' and ''contaminated'' surfaces, progressed to moduli measurements on partially filled single cracks, to measurements in ''clean'' and ''contaminated'' porous rocks and finally to a feasibility study in the field. For the experiments with the simple geometries we were able to measure fundamental physical properties such as contact angles of the meniscus and time dependent forces required to get the meniscus moving and to keep it moving at various velocities. From the data thus gathered we were able to interpret the complex elastic moduli data we measured in the partially saturated single cracks. While the geometry in real rocks is too complex to make precise calculations we determined that we had indeed identified the mechanisms responsible for the changes in the moduli we had measured. Thus encouraged by the laboratory studies we embarked on a field experiment in the desert of Arizona. The field site allowed for controlled irrigation. Instrumentation for fluid sampling and water penetration were already in place. The porous loosely consolidated rocks at the site were not ideal for finding the effects of the attenuation mechanism we had identified in the lab, but for logistic and cost constraint reasons we chose to field test the idea at that site. Tiltmeters and seismometers were installed and operated nearly continuously for almost 3 years. The field was irrigated with water in the fall of 2003 and with water containing a biosurfactant in the fall of 2004. We have indications that the biosurfactant irrigation has had a notable effect on the tilt data.
Identification and isolation of closed pore in porous rock using digital rock physics approach
NASA Astrophysics Data System (ADS)
Latief, Fourier Dzar Eljabbar
2015-09-01
The presence of closed pore in porous rock provide various effect with regard to its structural, elastic and flow properties. Physical based approach to measure porosity such as mercury porosimetry injection is unable to locate closed pore inside porous rock even though it is still possible to quantify the closed porosity. Digital data of porous rock in the form of three dimensional image can now be obtained by means of several methods such as micro-CT scan. Using the digital data, closed pore can be identified and isolated using digital rock physics approach. We first construct a synthetic three dimensional porous sample which consist of two simple side-to-side connected pore (cylinder and box shaped) and two spherical isolated pore which has closed porosity of 1.41 %. The digital image analysis which implemented in software CTAn (Bruker Micro-CT) still produce error of 0.04% which is very low. However, analysis using Lattice Boltzmann Method based simulation of fluid flow provide exact match to the closed porosity of the synthetic sample. Nevertheless, there are two disadvantages of this method, i.e., the simulation could take hours compared to the digital image analysis which only took several minutes and the limitation of numerical definition of zero velocity. Thus we apply both methods in order to overcome the drawbacks of each methods to analyze the digital sample of Fontainebleau sandstone. Using CTAn, we obtained the closed porosity of 0.02891845 % and using the LBM based fluid flow simulation of 0.028948346 %. The closed pore can then be isolated to further calculate the surface area. The result also confirmed that pore space of Fontainebleau sandstone is well connected.
Modelling of reactive fluid transport in deformable porous rocks
NASA Astrophysics Data System (ADS)
Yarushina, V. M.; Podladchikov, Y. Y.
2009-04-01
One outstanding challenge in geology today is the formulation of an understanding of the interaction between rocks and fluids. Advances in such knowledge are important for a broad range of geologic settings including partial melting and subsequent migration and emplacement of a melt into upper levels of the crust, or fluid flow during regional metamorphism and metasomatism. Rock-fluid interaction involves heat and mass transfer, deformation, hydrodynamic flow, and chemical reactions, thereby necessitating its consideration as a complex process coupling several simultaneous mechanisms. Deformation, chemical reactions, and fluid flow are coupled processes. Each affects the others. Special effort is required for accurate modelling of the porosity field through time. Mechanical compaction of porous rocks is usually treated under isothermal or isoentropic simplifying assumptions. However, joint consideration of both mechanical compaction and reactive porosity alteration requires somewhat greater than usual care about thermodynamic consistency. Here we consider the modelling of multi-component, multi-phase systems, which is fundamental to the study of fluid-rock interaction. Based on the conservation laws for mass, momentum, and energy in the form adopted in the theory of mixtures, we derive a thermodynamically admissible closed system of equations describing the coupling of heat and mass transfer, chemical reactions, and fluid flow in a deformable solid matrix. Geological environments where reactive transport is important are located at different depths and accordingly have different rheologies. In the near surface, elastic or elastoplastic properties would dominate, whereas viscoplasticity would have a profound effect deeper in the lithosphere. Poorly understood rheologies of heterogeneous porous rocks are derived from well understood processes (i.e., elasticity, viscosity, plastic flow, fracturing, and their combinations) on the microscale by considering a
Rock Physics Models of Biofilm Growth in Porous Media
NASA Astrophysics Data System (ADS)
Jaiswal, P.; alhadhrami, F. M.; Atekwana, E. A.
2013-12-01
Recent studies suggest the potential to use acoustic techniques to image biofilm growth in porous media. Nonetheless the interpretation of the seismic response to biofilm growth and development remains speculative because of the lack of quantitative petrophysical models that can relate changes in biofilm saturation to changes in seismic attributes. Here, we report our efforts in developing quantitative rock physics models to biofilm saturation with increasing and decreasing P-wave velocity (VP) and amplitudes recorded in the Davis et al. [2010] physical scale experiment. We adapted rock physics models developed for modeling gas hydrates in unconsolidated sediments. Two distinct growth models, which appear to be a function of pore throat size, are needed to explain the experimental data. First, introduction of biofilm as an additional mineral grain in the sediment matrix (load-bearing mode) is needed to explain the increasing time-lapse VP. Second, introduction of biofilm as part of the pore fluid (pore-filling mode) is required to explain the decreasing time-lapse VP. To explain the time-lapse VP, up to 15% of the pore volume was required to be saturated with biofilm. The recorded seismic amplitudes, which can be expressed as a function of porosity, permeability and grain size, showed a monotonic time-lapse decay except on Day 3 at a few selected locations, where it increased. Since porosity changes are constrained by VP, amplitude increase could be modeled by increasing hydraulic conductivity. Time lapse VP at locations with increasing amplitudes suggest that these locations have a load-bearing growth style. We conclude that permeability can increase by up to 10% at low (~2%) biofilm saturation in load-bearing growth style due to the development of channels within the biofilm structure. Developing a rock physics model for the biofilm growth in general may help create a field guide for interpreting porosity and permeability changes in bioremediation, MEOR and
Hoversten, G.M.; Gritto, Roland; Washbourne, John; Daley, Tom
2002-06-10
This paper presents a method for combining seismic and electromagnetic measurements to predict changes in water saturation, pressure, and CO{sub 2} gas/oil ratio in a reservoir undergoing CO{sub 2} flood. Crosswell seismic and electromagnetic data sets taken before and during CO{sub 2} flooding of an oil reservoir are inverted to produce crosswell images of the change in compressional velocity, shear velocity, and electrical conductivity during a CO{sub 2} injection pilot study. A rock properties model is developed using measured log porosity, fluid saturations, pressure, temperature, bulk density, sonic velocity, and electrical conductivity. The parameters of the rock properties model are found by an L1-norm simplex minimization of predicted and observed differences in compressional velocity and density. A separate minimization, using Archie's law, provides parameters for modeling the relations between water saturation, porosity, and the electrical conductivity. The rock-properties model is used to generate relationships between changes in geophysical parameters and changes in reservoir parameters. Electrical conductivity changes are directly mapped to changes in water saturation; estimated changes in water saturation are used along with the observed changes in shear wave velocity to predict changes in reservoir pressure. The estimation of the spatial extent and amount of CO{sub 2} relies on first removing the effects of the water saturation and pressure changes from the observed compressional velocity changes, producing a residual compressional velocity change. This velocity change is then interpreted in terms of increases in the CO{sub 2}/oil ratio. Resulting images of the CO{sub 2}/oil ratio show CO{sub 2}-rich zones that are well correlated to the location of injection perforations, with the size of these zones also correlating to the amount of injected CO{sub 2}. The images produced by this process are better correlated to the location and amount of injected
Effective Hydro-Mechanical Properties of Fluid-Saturated Fracture Networks
NASA Astrophysics Data System (ADS)
Pollmann, N.; Vinci, C.; Renner, J.; Steeb, H.
2015-12-01
Consideration of hydro-mechanical processes is essential for the characterization of liquid-resources as well as for many engineering applications. Furthermore, the modeling of seismic waves in fractured porous media finds application not only in geophysical exploration but also reservoir management. Fractures exhibit high-aspect-ratio geometries, i.e. they constitute thin and long hydraulic conduits. Motivated by this peculiar geometry, the investigation of the hydro-mechanically coupled processes is performed by means of a hybrid-dimensional modeling approach. The effective material behavior of domains including complex fracture patterns in a porous rock is assessed by investigating the fluid pressure and the solid displacement of the skeleton saturated by compressible fluids. Classical balance equations are combined with a Poiseuille-type flow in the dimensionally reduced fracture. In the porous surrounding rock, the classical Biot-theory is applied. For simple geometries, our findings show that two main fluid-flow processes occur, leak-off from fractures to the surrounding rock and fracture flow within and between the connected fractures. The separation of critical frequencies of the two flow processes is not straightforward, in particular for systems containing a large number of fractures. Our aim is to model three dimensional hydro-mechanically coupled processes within complex fracture patterns and in particular determine the frequency-dependent attenuation characteristics. Furthermore, the effect of asperities of the fracture surfaces on the fracture stiffness and on the hydraulic conductivity will be added to the approach.
Unsaturated hydraulic properties of porous sedimentary rocks explained by mercury porosimetry
NASA Astrophysics Data System (ADS)
Clementina Caputo, Maria; Turturro, Celeste; Gerke, Horst H.
2016-04-01
The understanding of hydraulic properties is essential in the modeling of flow and solute transport including contaminants through the vadose zone, which consists of the soil as well as of the underlying porous sediments or rocks. The aim of this work is to study the relationships between unsaturated hydraulic properties of porous rocks and their pore size distribution. For this purpose, two different lithotypes belonging to Calcarenite di Gravina Formation, a Plio-Pleistocene sedimentary rock of marine origin, were investigated. The two lithotypes differ mainly in texture and came from two distinct quarry districts, Canosa di Puglia (C) and Massafra (M) in southern Italy, respectively. This relatively porous rock formation (porosities range between 43% for C and 41% for M) often constitutes a thick layer of vadose zone in several places of Mediterranean basin. The water retention curves (WRCs) and the unsaturated hydraulic conductivity functions were determined using four different experimental methods that cover the full range from low to high water contents: the WP4 psychrometer test, the Wind's evaporation method, the Stackman's method and the Quasi-steady centrifuge method. Pore size estimation by means of mercury intrusion porosimetry (MIP) was performed. WRCs were compared with the pore size distributions to understand the influence of fabric, in terms of texture and porosity, features of pores and pore size distribution on the hydraulic behavior of rocks. The preliminary results show that the pore size distributions obtained by MIP do not cover the entire pore size range of the investigated Calcarenite. In fact, some pores in the rock samples of both lithotypes were larger than the maximum size that could be investigated by MIP. This implies that for explaining the unsaturated hydraulic properties over the full moisture range MIP results need to be combined with results obtained by other methods such as image analysis and SEM.
Modeling added compressibility of porosity and the thermomechanical response of wet porous rock
Rubin, M.B.; Elata, D.; Attia, A.V.
1995-06-01
This paper concerned with modeling the response of a porous brittle solid whose pores may be dry or partially filled with fluid. A form for the Helmholtz free energy is proposed which incorporated known Mie-Grueneisen constitutive equations for the nonporous solid and for the fluid, and which uses an Eilnstein formulation with variable specific heat. In addition, a functional form for porosity is postulated which porous rock. Restrictions on constitutive assumptions for the composite of porous solid ad fluid are obtained which ensure thermodynamic consistency. Examples show that although the added compressibility of porosity is determined by fitting data for dry Mt. Helen Tuff, the predicted responses of saturated and partially saturated tuff agree well with experimental data.
On the Relationship between Stress and Elastic Strain for Porous and Fractured Rock
Liu, Hui-Hai; Rutqvist, Jonny; Berryman, James G.
2008-02-25
Modeling the mechanical deformations of porous and fractured rocks requires a stress-strain relationship. Experience with inherently heterogeneous earth materials suggests that different varieties of Hook's law should be applied within regions of the rock having significantly different stress-strain behavior, e.g., such as solid phase and various void geometries. We apply this idea by dividing a rock body conceptually into two distinct parts. The natural strain (volume change divided by rock volume at the current stress state), rather than the engineering strain (volume change divided by the unstressed rock volume), should be used in Hooke's law for accurate modeling of the elastic deformation of that part of the pore volume subject to a relatively large degree of relative deformation (i.e., cracks or fractures). This approach permits the derivation of constitutive relations between stress and a variety of mechanical and/or hydraulic rock properties. We show that the theoretical predictions of this method are generally consistent with empirical expressions (from field data) and also laboratory rock experimental data.
NASA Astrophysics Data System (ADS)
Sharifi Haddad, Amin
Fractured porous media are important structures in petroleum engineering and geohydrology. The accelerating global demand for energy has turned the focus to fractured formations. The fractured porous media are also found in conventional naturally fractured reservoirs and the water supply from karst (carbonate) aquifers. Studying mass transfer processes allows us to explore the complexities and uncertainties encountered with fractured rocks. This dissertation is developing an analytical methodology for the study of mass transfer in fractured reservoirs. The dissertation begins with two cases that demonstrate the importance of the rock matrix block size distribution and dispersivity through a transient mass exchange mechanism between rock matrix blocks and fractures. The first case assumes a medium with no surface adsorption, and the second case includes the surface adsorption variable. One of the main focuses of this work is the characterization of the rock matrix block size distribution in fractured porous media. Seismic surveying, well test analysis, well logging, and geomechanical tools are currently used to characterize this property, based on measurements of different variables. This study explores an innovative method of using solute transport to determine the fracture intensity. This methodology is applied to slab-shaped rock matrix blocks and can easily be extended to other geometries. Another focus of this dissertation is the characterization of dispersivity in field scale studies. Improving our knowledge of dispersivity will enable more accurate mass transfer predictions and advance the study of transport processes. Field tracer tests demonstrated that dispersivity is scale-dependent. Proposed functions for the increasing trend of dispersivity include linear and asymptotic scale-dependence. This study investigated the linear dispersivity trend around the injection wellbore. An analysis of the tracer concentration in a monitoring well was used to
Anisotropy of a synthetic layered porous medium
Furre, A.K.; Holt, R.M.
1995-12-31
Artificial, periodically laminated porous sandstones are manufactured in order to study acoustic behavior compared to permeability of porous rocks. Cores are made by mixing sand and epoxy. The laminated media are built up by two alternating types of layers, one highly permeable layer with low velocity, the other with lower permeability but higher velocity (the permeability of each layer is, however, of the order of 1-50 Darcy). Permeability and acoustic measurements are performed at different incidence angles to the layers. The wave velocities in dry samples show anisotropy largely confirming the Backus theory. Fluid saturation increases the velocities and introduces an attenuation anisotropy which may be linked to fluid flow anisotropy. The spectral contents of the waves traces do not depend on the angle of incidence, in contrast to the case when anisotropy is induced by fractures.
Use of ``rock-typing`` to characterize carbonate reservoir heterogeneity. Final report
Ikwuakor, K.C.
1994-03-01
The objective of the project was to apply techniques of ``rock-typing`` and quantitative formation evaluation to borehole measurements in order to identify reservoir and non-reservoir rock-types and their properties within the ``C`` zone of the Ordovician Red River carbonates in the northeast Montana and northwest North Dakota areas of the Williston Basin. Rock-typing discriminates rock units according to their pore-size distribution. Formation evaluation estimates porosities and pore fluid saturation. Rock-types were discriminated using crossplots involving three rock-typing criteria: (1) linear relationship between bulk density and porosity, (2) linear relationship between acoustic interval transit-time and porosity, and (3) linear relationship between acoustic interval transit-time and bulk density. Each rock-type was quantitatively characterized by the slopes and intercepts established for different crossplots involving the above variables, as well as porosities and fluid saturations associated with the rock-types. All the existing production was confirmed through quantitative formation evaluation. Highly porous dolomites and anhydritic dolomites contribute most of the production, and constitute the best reservoir rock-types. The results of this study can be applied in field development and in-fill drilling. Potential targets would be areas of porosity pinchouts and those areas where highly porous zones are downdip from non-porous and tight dolomites. Such areas are abundant. In order to model reservoirs for enhanced oil recovery (EOR) operations, a more localized (e.g. field scale) study, expanded to involve other rock-typing criteria, is necessary.
A visco-poroelastic damage model for modelling compaction and brittle failure of porous rocks
NASA Astrophysics Data System (ADS)
Jacquey, Antoine B.; Cacace, Mauro; Blöcher, Guido; Milsch, Harald; Scheck-Wenderoth, Magdalena
2016-04-01
Hydraulic stimulation of geothermal wells is often used to increase heat extraction from deep geothermal reservoirs. Initiation and propagation of fractures due to pore pressure build-up increase the effective permeability of the porous medium. Understanding the processes controlling the initiation of fractures, the evolution of their geometries and the hydro-mechanical impact on transport properties of the porous medium is therefore of great interest for geothermal energy production. In this contribution, we will present a thermodynamically consistent visco-poroelastic damage model which can deal with the multi-scale and multi-physics nature of the physical processes occurring during deformation of a porous rock. Deformation of a porous medium is crucially influenced by the changes in the effective stress. Considering a strain-formulated yield cap and the compaction-dilation transition, three different regimes can be identified: quasi-elastic deformation, cataclastic compaction with microcracking (damage accumulation) and macroscopic brittle failure with dilation. The governing equations for deformation, damage accumulation/healing and fluid flow have been implemented in a fully-coupled finite-element-method based framework (MOOSE). The MOOSE framework provides a powerful and flexible platform to solve multiphysics problems implicitly and in a tightly coupled manner on unstructured meshes which is of interest for such non-linear context. To validate and illustrate the model, simulations of the deformation behaviour of cylindrical porous Bentheimer sandstone samples under different confining pressures are compared to experiments. The first experiment under low confining pressure leads to shear failure, the second for high confining pressure leads to cataclastic compaction and the third one with intermediate confining pressure correspond to a transitional regime between the two firsts. Finally, we will demonstrate that this approach can also be used at the field
Water saturation effects on elastic wave attenuation in porous rocks with aligned fractures
NASA Astrophysics Data System (ADS)
Amalokwu, Kelvin; Best, Angus I.; Sothcott, Jeremy; Chapman, Mark; Minshull, Tim; Li, Xiang-Yang
2014-05-01
Elastic wave attenuation anisotropy in porous rocks with aligned fractures is of interest to seismic remote sensing of the Earth's structure and to hydrocarbon reservoir characterization in particular. We investigated the effect of partial water saturation on attenuation in fractured rocks in the laboratory by conducting ultrasonic pulse-echo measurements on synthetic, silica-cemented, sandstones with aligned penny-shaped voids (fracture density of 0.0298 ± 0.0077), chosen to simulate the effect of natural fractures in the Earth according to theoretical models. Our results show, for the first time, contrasting variations in the attenuation (Q-1) of P and S waves with water saturation in samples with and without fractures. The observed Qs/Qp ratios are indicative of saturation state and the presence or absence of fractures, offering an important new possibility for remote fluid detection and characterization.
Multinuclear NMR microscopy of two-phase fluid systems in porous rock.
Doughty, D A; Tomutsa, L
1996-01-01
The high-field magnetic resonance (MR) characteristics of fluids in porous reservoir rock exhibit short T2 relaxation times and broad natural line widths. These characteristics severely restrict which MR imaging (MRI) methodology can be used to obtain high-resolution porescale images of fluids in porous rock. An MR microscopy protocol based on 3D backprojection using strong imaging gradients was developed to overcome many of these constraints. To improve the image quality of two-phase systems, multinuclear MRI using proton MR to image the brine phase and 19F MR of a fluorinated hydrocarbon to image the oil phase was used. Resolution as high as 25 microns per pixel has been obtained for fluid systems in Bentheim and Fontainebleau sandstones. Separate proton and 19F images of brine and oil phases show good agreement with total saturation images. Software has been developed to perform 3D erosion/dilations and to extract the pore size distribution from binarized 3D images of fluid filled porosity. Results from pore size measurements show significant differences in the nature of the pore network in Fontainebleau and Bentheim sandstones. PMID:8970097
A damage mechanics approach for quantifying stress changes due to brittle failure of porous rocks
NASA Astrophysics Data System (ADS)
Jacquey, Antoine B.; Cacace, Mauro; Blöcher, Guido; Milsch, Harald; Scheck-Wenderoth, Magdalena
2016-04-01
Natural fault zones or man-made injection or production of fluid impact the regional stress distribution in Earth's crust and can be responsible for localized stress discontinuities. Understanding the processes controlling fracturing of the porous rocks and mechanical behaviour of fault zones is therefore of interest for several applications including geothermal energy production. In this contribution, we will present a thermodynamically consistent visco-poroelastic damage model which can deal with the multi-scale and multi-physics nature of the physical processes controlling the deformation of porous rocks during and after brittle failure. Deformation of a porous medium is crucially influenced by the changes in the effective stress. Considering a strain-formulated yield cap and the compaction-dilation transition, three different regimes can be identified: quasi-elastic deformation, cataclastic compaction with microcracking (damage accumulation) and macroscopic brittle failure with dilation. The governing equations for deformation, damage accumulation/healing and fluid flow have been implemented in a fully-coupled finite-element-method based framework (MOOSE). The MOOSE framework provides a powerful and flexible platform to solve multiphysics problems implicitly and in a tightly coupled manner on unstructured meshes which is of interest for such non-linear context. To illustrate the model, simulation of a compaction experiment of a sandstone leading to shear failure will be presented which allows to quantify the stress drop accompanying the failure. Finally, we will demonstrate that this approach can also be used at the field scale to simulate hydraulic fracturing and assess the resulting changes in the stress field.
An H2O-CO2 mixed fluid saturation model compatible with rhyolite-MELTS
NASA Astrophysics Data System (ADS)
Ghiorso, Mark S.; Gualda, Guilherme A. R.
2015-06-01
A thermodynamic model for estimating the saturation conditions of H2O-CO2 mixed fluids in multicomponent silicate liquids is described. The model extends the capabilities of rhyolite-MELTS (Gualda et al. in J Petrol 53:875-890, 2012a) and augments the water saturation model in MELTS (Ghiorso and Sack in Contrib Mineral Petrol 119:197-212, 1995). The model is internally consistent with the fluid-phase thermodynamic model of Duan and Zhang (Geochim Cosmochim Acta 70:2311-2324, 2006). It may be used independently of rhyolite-MELTS to estimate intensive variables and fluid saturation conditions from glass inclusions trapped in phenocrysts. The model is calibrated from published experimental data on water and carbon dioxide solubility, and mixed fluid saturation in silicate liquids. The model is constructed on the assumption that water dissolves to form a hydroxyl melt species, and that carbon dioxide both a molecular species and a carbonate ion, the latter complexed with calcium. Excess enthalpy interaction terms in part compensate for these simplistic assumptions regarding speciation. The model is restricted to natural composition liquids over the pressure range 0-3 GPa. One characteristic of the model is that fluid saturation isobars at pressures greater than ~100 MPa always display a maximum in melt CO2 at nonzero H2O melt concentrations, regardless of bulk composition. This feature is universal and can be attributed to the dominance of hydroxyl speciation at low water concentrations. The model is applied to four examples. The first involves estimation of pressures from H2O-CO2-bearing glass inclusions found in quartz phenocrysts of the Bishop Tuff. The second illustrates H2O and CO2 partitioning between melt and fluid during fluid-saturated equilibrium and fractional crystallization of MORB. The third example demonstrates that the position of the quartz-feldspar cotectic surface is insensitive to melt CO2 contents, which facilitates geobarometry using phase
Initial conditions or emergence; what determines dissolution patterns in heterogeneous porous rocks?
NASA Astrophysics Data System (ADS)
Szymczak, Piotr; Upadhyay, Virat; Ladd, Anthony
2016-04-01
Dissolution of fractured or porous rocks is often accompanied by the formation of highly localized flow paths. Dissolution, in general, does not proceed uniformly, as it is influenced both by the heterogeneities in the rock matrix and by the instabilities associated with the positive feedback loops between the flow, dispersion, and chemical reactions. As a result, distinct channels or "wormholes" develop within the rocks in which both the flow and dissolution focus. In this communication, we aim to investigate how these emerging flow paths are influenced by the initial local inhomogeneities of the porosity field. Our results indicate a surprising insensitivity of the evolving dissolution patterns and flow rates to the amplitude and correlation length characterizing the inhomogeneities. At long times wormhole competition overwhelms the initial variations in aperture distribution, resulting in a universal relation between the separation of the wormholes and their length. This hierarchy of scales even persists in the presence of relatively large inhomogeneities (vugs), which focus the flow at the beginning of the dissolution process, but - if the sample is large enough - with time tend to be overwhelmed by the spontaneous growth of instabilities. A natural consequence of wormhole competition is that the separation between growing wormholes corresponds roughly to their length, something that is borne out by field observations.
NASA Astrophysics Data System (ADS)
Guillon, S.; Vu, M. T.; Pili, E.; Adler, P. M.
2013-05-01
Air permeability is measured in the fractured crystalline rocks of the Roselend Natural Laboratory (France). Single-hole pneumatic injection tests as well as differential barometric pressure monitoring are conducted on scales ranging from 1 to 50 m, in both shallow and deep boreholes, as well as in an isolated 60 m3 chamber at 55 m depth. The field experiments are interpreted using numerical simulations in equivalent homogeneous porous media with their real 3-D geometry in order to estimate pneumatic parameters. For pneumatic injection tests, steady-state data first allow to estimate air permeability. Then, pressure recovery after a pneumatic injection test allows to estimate the air-filled porosity. Comparison between the various studied cases clarifies the influence of the boundary conditions on the accuracy of the often used 1-D estimate of air permeability. It also shows that permeabilities correlate slightly with fracture density. In the chamber, a 1 order-of-magnitude difference is found between the air permeabilities obtained from pneumatic injection tests and from differential barometric pressure monitoring. This discrepancy is interpreted as a scale effect resulting from the approximation of the heterogeneous fractured rock by a homogeneous numerical model. The difference between the rock volumes investigated by pneumatic injection tests and by differential barometric pressure monitoring may also play a role. No clear dependence of air permeability on saturation has been found so far.
Quasistatic Shock Waves: A Mechanism for Nonuniform Compaction in Porous Rock
OLSSON,WILLIAM A.
2000-09-08
Recent studies have observed compaction zones pass through porous rock under axisymmetric compression. An initially thin, compacted layer appears at the yield point of the stress-strain curve and then grows by thickening in the direction of maximum compression at constant stress. Strain localization theory has been applied to compaction to explain the formation of these features. This paper describes the growth of the compaction zones, that is, the propagation of their boundaries, in terms of shock wave analysis. The ratio of the applied shortening rate to the velocity of the boundary is related to the porosity change across the boundary. Certain features of the stress-strain curve are explained by the model.
Jang, Long-Kuan; Chang, Philip W.; Findley, John E.; Yen, Teh Fu
1983-01-01
This paper presents a bench-scale study on the transport in highly permeable porous rock of three bacterial species—Bacillus subtilis, Pseudomonas putida, and Clostridium acetobutylicum—potentially applicable in microbial-enhanced oil recovery processes. The transport of cells during the injection of bacterial suspension and nutrient medium was simulated by a deep bed filtration model. Deep bed filtration coefficients and the maximum capacity of cells in porous rock were measured. Low to intermediate (∼106/ml) injection concentrations of cellular suspensions are recommended because plugging of inlet surface is less likely to occur. In addition to their resistance to adverse environments, spores of clostridia are strongly recommended for use in microbial-enhanced oil recovery processes since they are easiest among the species tested to push through porous rock. After injection, further transport of bacteria during incubation can occur by growth and mobility through the stagnant nutrient medium which fills the porous rock. We have developed an apparatus to study the migration of bacteria through a Berea sandstone core containing nutrient medium. PMID:16346414
Jang, L K; Chang, P W; Findley, J E; Yen, T F
1983-11-01
This paper presents a bench-scale study on the transport in highly permeable porous rock of three bacterial species-Bacillus subtilis, Pseudomonas putida, and Clostridium acetobutylicum-potentially applicable in microbial-enhanced oil recovery processes. The transport of cells during the injection of bacterial suspension and nutrient medium was simulated by a deep bed filtration model. Deep bed filtration coefficients and the maximum capacity of cells in porous rock were measured. Low to intermediate ( approximately 10/ml) injection concentrations of cellular suspensions are recommended because plugging of inlet surface is less likely to occur. In addition to their resistance to adverse environments, spores of clostridia are strongly recommended for use in microbial-enhanced oil recovery processes since they are easiest among the species tested to push through porous rock. After injection, further transport of bacteria during incubation can occur by growth and mobility through the stagnant nutrient medium which fills the porous rock. We have developed an apparatus to study the migration of bacteria through a Berea sandstone core containing nutrient medium. PMID:16346414
Use of GPU Computing to Study Coupled Deformation and Fluid Flow in Porous Rocks
NASA Astrophysics Data System (ADS)
Räss, Ludovic; Omlin, Samuel; Simon, Nina; Podladchikov, Yuri
2015-04-01
Actual challenges in computational geodynamics put high requirements for the development of new coupled models. These need to solve accurate physics, on high resolution and in reasonable computation time. Multi-scale problems such as deformation of porous rocks triggered by fluid flow require both high temporal and spatial resolution. The resulting preferential flow paths involve complex physics and a strong coupling between deformation and fluid flow processes. Shortcuts such as sequential or iterative coupling of two existing solvers will not be sufficient in these difficult cases to localize the deformation and flow. We base our numerical implementation on the physically and thermodynamically consistent mathematical model for fluid flow in porous rocks, taking nonlinear stress dependent visco-elasto-plastic rheology into account. The effective permeability used for the Darcy flow is obtained through the nonlinear Karman-Cozeny relation. The model is not restricted by the lithostatic stress assumption, allowing for background stress regime as it occurs in natural conditions. We have developed a fully three-dimensional numerical application based on an iterative finite difference scheme. The application is written in C-CUDA, is enabled for GPU accelerators and is parallelized with MPI to run on multi-GPU clusters. The parallelization on a rectangular grid is straightforward (at each iteration, the boundaries of the local problem are updated by the neighboring processes) and requires no MPI global operations, only MPI point-to-point communication between neighboring processes. This parallelization method should allow by construction for linear weak scaling on any number of processors. Our linearly scaling numerical application predicts the formation of dynamically evolving fluid pathways. These supercomuting applications are vital for resolving actual challenging high-resolution three-dimensional models.
A fully coupled porous flow and geomechanics model for fluid driven cracks: a peridynamics approach
NASA Astrophysics Data System (ADS)
Ouchi, Hisanao; Katiyar, Amit; York, Jason; Foster, John T.; Sharma, Mukul M.
2015-03-01
A state-based non-local peridynamic formulation is presented for simulating fluid driven fractures in an arbitrary heterogeneous poroelastic medium. A recently developed peridynamic formulation of porous flow has been coupled with the existing peridynamic formulation of solid and fracture mechanics resulting in a peridynamic model that for the first time simulates poroelasticity and fluid-driven fracture propagation. This coupling is achieved by modeling the role of pore pressure on the deformation of porous media and vice versa through porosity variation with medium deformation, pore pressure and total mean stress. The poroelastic model is verified by simulating the one-dimensional consolidation of fluid saturated rock. An additional porous flow equation with material permeability dependent on fracture width is solved to simulate fluid flow in the fractured region. Finally, single fluid-driven fracture propagation with a two-dimensional plane strain assumption is simulated and verified against the corresponding classical analytical solution.
Houseworth, J.E.
2004-09-16
Exact analytical solutions are presented for solute transport in an unsaturated fracture and porous rock matrix. The problem includes advective transport in the fracture and rock matrix as well as advective and diffusive fracture-matrix exchange. Linear sorption in the fracture and matrix and radioactive decay are also treated. The solution is for steady, uniform transport velocities within the fracture and matrix, but allows for independent specification of each of the velocities. The problem is first solved in terms of the solute concentrations that result from an instantaneous point source. Superposition integrals are then used to derive the solute mass flux at a fixed downstream position from an instantaneous point source and for the solute concentrations that result from a continuous point source. Solutions are derived for cases with the solute source in the fracture and the solute source in the matrix. The analytical solutions are closed-form and are expressed in terms of algebraic functions, exponentials, and error functions. Comparisons between the analytical solutions and numerical simulations, as well as sensitivity studies, are presented. Increased sensitivity to cross-flow and solute source location is found for increasing Peclet number. The numerical solutions are found to compare well with the analytical solutions at lower Peclet numbers ,but show greater deviation at higher Peclet numbers.
Impact of hygrometry changes on creep behaviour of a porous rock and associated acoustic emission.
NASA Astrophysics Data System (ADS)
Grgic, D.; Amitrano, D.
2009-04-01
Static fatigue of a polycrystalline porous rock (iron ore) was studied by performing multi step uniaxial creep tests under partially saturated conditions, and the impact of water saturation was analyzed. The samples were, in a first step set to a partial water saturation of 90%. In a second step, the samples were saturated completely in order to simulate the impact of flooding corresponding to the conditions of abandoned iron mines. We recorded axial and transversal strain and acoustic emission (EA). The experimental results show that the water saturation induces a strong increase in AE activity and dilatant inelastic volumetric strain. This is associated with a notable decrease in Young's modulus and in the b-value of the Gutenberg-Richter law (i.e., the relative number of large-amplitude events increases) as the rock approaches failure, indicating that microfracturing plays an important role in the creep process. Water saturation accelerates static fatigue through hydro-mechanical coupling and subcritical stress corrosion cracking. The chemical reactions involved in the corrosion of iron ore and leading to a decrease in its intrinsic mechanical properties are described. These reactions play a major role in the static fatigue of iron ore, which on a large scale is probably the main mechanism explaining certain collapses in underground iron mines. It is also shown that creep straining of iron ore is partially reversible after stress removal, indicating that it results also from time-dependent viscoplastic mechanism (i.e., dislocation creep).
Numerical simulation of salt cementation in the porous rocks adjacent to salt diapirs
NASA Astrophysics Data System (ADS)
Allstadt, Raphael; Li, Shiyuan; Marquart, Gabriele; Reuning, Lars; Niederau, Jan
2015-04-01
Porosity and permeability are among the most important petrophysical properties of reservoirs rocks in oil systems. Observations during exploration indicate that in the vicinity of salt domes the porosity of reservoir rocks is often reduced by halite cementation. In this study we present results of simulating the process of salt precipitation near salt diapirs by using a schematic model of a Zechstein diapir in the North Sea basin. The numerical simulation is based on solving the transport equations for heat, porous flow and dispersive and reactive chemical species. Chemical reaction and equilibrium is based on the PHREEQC computer code. In our model over-pressured brine is entering from below and is deflected towards the diapir due to an intermediate layer of low permeability. The high thermal conductivity of salt yields a lateral temperature gradient starting from the diapir. Due to this effect the simulated temperature profile shows lower temperatures close to the salt dome than in comparable depths further away. Caused by the temperature-controlled solubility of NaCl in the brine and supplied ions by the diapir, halite first precipitates near the salt diapir by cementing the pore spaces and thus reducing the porosity. Salt-precipitation in the simulation starts after 840 000 years and reduces the porosity from 10 % to 5.5 % after 19 Mill. years. The permanent influx of brine causes growth of the cementation area and the related reduction of porosity in the reservoir.
Multiphase fluid simulations through porous rock using Shan-Chen type lattice Boltzmann method
NASA Astrophysics Data System (ADS)
Qohar, Ulin Nuha Abdul; Latief, Fourier Dzar Eljabbar; Fauzi, Umar
2015-09-01
Fluid flow with multiphase components is one of the daily problems that interesting to learn and widely used in various fields, one of them in rock physics. Euler approach and Lagrange approach, the two models are applied to study the fluid flow, which often known as the microscopic and macroscopic approaches. Lattice Boltzmann Method combines the advantages and appropriate of both approaches, that is used as a more efficient model approach, also known as mesoscopic. The LB method has been used to model the two-phase fluid flow with different viscosities using the Shan-Chen model, this model explain the interaction between two different fluid phases. Validated numerical models compiled using two ways, namely analytical models and physical models. Based on this research, numerical models are designed to meet the analytical model with an error on the lattice at the edge of the lattice. When compared with physical models, numerical models have qualitatively suitability. Based on the results of the validation of numerical models, modeling is done to the image of porous rock which gives the results of multiphase fluid flow profile inside the pore. Moreover, the results obtained indicate that there are effective pores that fluid can not be drained.
Pore scale heterogeneity in the mineral distribution and surface area of porous rocks
NASA Astrophysics Data System (ADS)
Lai, Peter; Moulton, Kevin; Krevor, Samuel
2014-05-01
There are long-standing challenges in characterizing reactive transport in porous media at scales larger than individual pores. This hampers the prediction of the field-scale impact of geochemical processes on fluid flow [1]. This is a source of uncertainty for carbon dioxide injection, which results in a reactive fluid-rock system, particularly in carbonate rock reservoirs. A potential cause is the inability of the continuum approach to incorporate the impact of heterogeneity in pore-scale reaction rates. This results in part from pore-scale heterogeneities in surface area of reactive minerals [2,3]. The objective of this study was to quantify heterogeneity in reactive surface and observe the extent of its non-normal character. In this study we describe our work in using micron-scale x-ray imaging and other spectroscopic techniques for the purpose of describing the statistical distribution of reactive surface area within a porous medium, and identifying specific mineral phases and their distribution in 3-dimensions. Using in-house image processing techniques and auxilary charactersation with thin section, electron microscope and spectroscopic techniques we quantified the surface area of each mineral phase in the x-ray CT images. This quantification was validated against nitrogen BET surface area and backscattered electron imaging measurements of the CT-imaged samples. Distributions in reactive surface area for each mineral phase were constructed by calculating surface areas in thousands of randomly selected subvolume images of the total sample, each normalized to the pore volume in that image. In all samples, there is little correlation between the reactive surface area fraction and the volumetric fraction of a mineral in a bulk rock. Berea sandstone was far less heterogeneous and has a characteristic pore size at which a surface area distribution may be used to quantify heterogeneity. In carbonates, heterogeneity is more complex and surface area must be
Wave propagation in fluid-saturated media: waveform and spectral analysis
NASA Astrophysics Data System (ADS)
Boadu, F. K.
2000-04-01
A new numerical approach to the solution of waves propagating in a fluid-saturated medium, using Biot's theory as a foundation, has important implications for oil reservoir management and earthquake prediction. A numerical scheme is developed using an exponential transformation that explicitly treats the petrophysical and fluid properties of the medium within the framework of a generalized model. The scheme accounts for wave dissipation and velocity modifications. The numerical solution is used to perform numerical experiments to study the dynamic behaviour of waves in a fluid-saturated medium at well-logging frequencies (15kHz). The results from the numerical experiments indicate that the degree of saturation by a high-viscosity fluid (HVF) such as oil, the temperature and the porosity of a medium strongly influence the spectral power distribution, frequency content and the velocity of waves propagating through the medium. An increase in HVF saturation causes enhanced attenuation of the low-frequency components, and increases the seismic velocity. An increase in porosity, however, enriches the low-frequency components and decreases the seismic velocity. A spectral quantification procedure is suggested and used to obtain information about the petrophysical and fluid properties of the medium from the spectral characteristics of the transmitted waveform. The procedure involves segmentation of the energy or power distribution of the transmitted waveforms into specified energy bands. The energy or power in these bands is then estimated. The extracted quantification variables are found to have strong correlations with the degree of HVF saturation, and the temperature and the porosity of the medium.
Pure water injection into porous rock with superheated steam and salt in a solid state
NASA Astrophysics Data System (ADS)
Montegrossi, G.; Tsypkin, G.; Calore, C.
2012-04-01
Most of geothermal fields require injection of fluid into the hot rock to maintain pressure and productivity. The presence of solid salt in porous space may cause an unexpected change in the characteristics of the reservoir and produced fluids, and dramatically affect the profitability of the project. We consider an injection problem of pure water into high temperature geothermal reservoir, saturated with superheated vapour and solid salt. Pure water moves away from injection point and dissolves solid salt. When salty water reaches the low-pressure hot domain, water evaporation occurs and, consequently, salt precipitates. We develop a simplified analytical model of the process and derive the similarity solutions for a 1-D semi-infinite reservoir. These solutions are multi-valued and describe the reduction in permeability and porosity due to salt precipitation at the leading boiling front. If the parameters of the system exceed critical values, then similarity solution ceases to exist. We identify this mathematical behaviour with reservoir sealing in the physical system. The TOUGH2-EWASG code has been used to verify this hypothesis and investigate the precipitate formation for an idealized bounded 1-D geothermal system of a length of 500 m with water injection at one extreme and fluid extraction at the other one. Both boundaries are kept at constant pressure and temperature. The result for the semi-infinite numerical model show that the monotonic grow of the solid salt saturation to reach asymptotic similarity solution generally occurs over a very large length starting from the injection point. Reservoir sealing occurs if solid salt at the initial state occupies a considerable part of the porous space. Numerical experiments for the bounded 500 m system demonstrate that a small amount of salt is enough to get reservoir sealing. Generally, salt tend to accumulate near the production well, and salt plug forms at the elements adjacent to the extraction point. This type
Small watershed response to porous rock check dams in a semiarid watershed
NASA Astrophysics Data System (ADS)
Nichols, Mary; Polyakov, Viktor; Nearing, Mark
2016-04-01
Rock check dams are used throughout the world as technique for mitigating erosion problems on degraded lands. Increasingly, they are being used in restoration efforts on rangelands in the southwestern US, however, their impact on watershed response and channel morphology is not well quantified. In 2008, 37 porous rock structures were built on two small (4.0 and 3.1 ha) instrumented watersheds on an alluvial fan at the base of the Santa Rita Mountains in southern Arizona, USA. 35 years of historical rainfall and runoff, and sediment data are available to compare with 7 years of data collected after check dam construction. In addition, post construction measurements of channel geometry and longitudinal channel profiles were compared with pre-construction measurements to characterize the impact of check dams on sediment retention and channel morphology. The primary impact of the check dams is was retention of channel sediment and reduction in channel gradient; however response varied between the proximal watersheds with 80% of the check dams on one of the watersheds filled to 100% of their capacity after 7 runoff seasons. In addition, initial impact on precipitation runoff ratios is was not persistent. The contrasting watershed experiences lower sediment yields and only 20% of the check dams on this watershed are were filled to capacity and continue to influence runoff during small events. Within the watersheds the mean gradient of the channel reach immediately upstream of the structures has been reduced by 35% (from 0.061 to 0.039) and 34% on (from 0.071 to 0.047).
Visco-poroelastic damage model for brittle-ductile failure of porous rocks
NASA Astrophysics Data System (ADS)
Lyakhovsky, Vladimir; Zhu, Wenlu; Shalev, Eyal
2015-04-01
The coupling between damage accumulation, dilation, and compaction during loading of sandstones is responsible for different structural features such as localized deformation bands and homogeneous inelastic deformation. We distinguish and quantify the role of each deformation mechanism using new mathematical model and its numerical implementation. Formulation includes three different deformation regimes: (I) quasi-elastic deformation characterized by material strengthening and compaction; (II) cataclastic flow characterized by damage increase and compaction; and (III) brittle failure characterized by damage increase, dilation, and shear localization. Using a three-dimensional numerical model, we simulate the deformation behavior of cylindrical porous Berea sandstone samples under different confining pressures. The obtained stress, strain, porosity changes and macroscopic deformation features well reproduce the laboratory results. The model predicts different rock behavior as a function of confining pressures. The quasi-elastic and brittle regimes associated with formation of shear and/or dilatant bands occur at low effective pressures. The model also successfully reproduces cataclastic flow and homogeneous compaction under high pressures. Complex behavior with overlap of common features of all regimes is simulated under intermediate pressures, resulting with localized compaction or shear enhanced compaction bands. Numerical results elucidate three steps in the formation of compaction bands: (1) dilation and subsequent shear localization, (2) formation of shear enhanced compaction band, and (3) formation of pure compaction band.
Andrade, JosÃÂ© E; Rudnicki, John W
2012-12-14
In this project, a predictive multiscale framework will be developed to simulate the strong coupling between solid deformations and fluid diffusion in porous rocks. We intend to improve macroscale modeling by incorporating fundamental physical modeling at the microscale in a computationally efficient way. This is an essential step toward further developments in multiphysics modeling, linking hydraulic, thermal, chemical, and geomechanical processes. This research will focus on areas where severe deformations are observed, such as deformation bands, where classical phenomenology breaks down. Multiscale geometric complexities and key geomechanical and hydraulic attributes of deformation bands (e.g., grain sliding and crushing, and pore collapse, causing interstitial fluid expulsion under saturated conditions), can significantly affect the constitutive response of the skeleton and the intrinsic permeability. Discrete mechanics (DEM) and the lattice Boltzmann method (LBM) will be used to probe the microstructure---under the current state---to extract the evolution of macroscopic constitutive parameters and the permeability tensor. These evolving macroscopic constitutive parameters are then directly used in continuum scale predictions using the finite element method (FEM) accounting for the coupled solid deformation and fluid diffusion. A particularly valuable aspect of this research is the thorough quantitative verification and validation program at different scales. The multiscale homogenization framework will be validated using X-ray computed tomography and 3D digital image correlation in situ at the Advanced Photon Source in Argonne National Laboratories. Also, the hierarchical computations at the specimen level will be validated using the aforementioned techniques in samples of sandstone undergoing deformation bands.
Pore Scale Heterogeneity in the Mineral Distribution, Surface Area and Adsorption in Porous Rocks
NASA Astrophysics Data System (ADS)
Lai, P. E. P.; Krevor, S. C.
2014-12-01
The impact of heterogeneity in chemical transport and reaction is not understood in continuum (Darcy/Fickian) models of reactive transport. This is manifested in well-known problems such as scale dependent dispersion and discrepancies in reaction rate observations made at laboratory and field scales [1]. Additionally, this is a source of uncertainty for carbon dioxide injection, which produces a reactive fluid-rock system particularly in carbonate rock reservoirs. A potential cause is the inability of the continuum approach to incorporate the impact of heterogeneity in pore-scale reaction rates. This results in part from pore-scale heterogeneities in surface area of reactive minerals [2, 3]. We use x-ray micro tomography to describe the non-normal 3-dimensional distribution of reactive surface area within a porous medium according to distinct mineral groups. Using in-house image processing techniques, thin sections, nitrogen BET surface area, backscattered electron imaging and energy dispersive spectroscopy, we compare the surface area of each mineral phase to those obtained from x-ray CT imagery. In all samples, there is little correlation between the reactive surface area fraction and the volumetric fraction of a mineral in a bulk rock. Berea sandstone was far less heterogeneous and has a characteristic pore size at which a surface area distribution may be used to quantify heterogeneity. In carbonates, heterogeneity is more complex and surface area must be characterized at multiple length scales for an accurate description of reactive transport. We combine the mineral specific surface area characterisation to dynamic tomography, imaging the flow of water and solutes, to observe flow dependent and mineral specific adsorption. The observations may contribute to the incorporation of experimentally based statistical descriptions of pore scale heterogeneity in reactive transport into upscaled models, moving it closer to predictive capabilities for field scale
Approach to failure in a discrete element model of the compressive failure of porous rocks (Invited)
NASA Astrophysics Data System (ADS)
Kun, F.; Varga, I.; Lennartz-Sassinek, S.; Main, I. G.
2013-12-01
We investigate how a porous rock sample approaches failure under uniaxial compression. Computer simulations are carried out in the framework of a discrete element model (DEM) which takes into account both the microstructure of the material and the dynamics of local fracturing, revealing much more detail and observation bandwidth in then granular mechanics than possible during standard laboratory tests. The synthetic sample is generated by sedimentation of randomly-sized spherical particles with a log-normal size distribution inside a cylindrical container. The cohesive interaction of particles is represented by beam elements that break when overstressed. The breaking rule takes into account both stretching and shear of particle contacts. When particles not connected by a beam come into contact their interaction is described by the Hertz contact law. The time evolution of the system is generated by molecular dynamics simulations in three dimensions. Computer simulations showed that under strain controlled uniaxial loading of the system micro-cracks initially nucleate in an uncorrelated way all over the sample. As loading proceeds localization occurs, i.e. the damage concentrates into a narrow damage band. Inside the damage band the material is crushed, into a poorly sorted mixture of fine powder and larger fragments with a power-law mass distribution, as observed in fault wear products (gouge) in natural and laboratory faults. Dynamic bursts of radiated energy, analogous to acoustic emissions observed in laboratory experiments, are identified as correlated trails of local fracture emerging as the consequence of stress redistribution. Characteristic quantities of burst such as size/rupture area, released elastic energy, and duration proved to have power law probability-size distributions over a broad range. The energy and duration of bursts have power law dependence on the rupture area created. As the system approaches macroscopic failure consecutive bursts become
Constitutive Theory for Velocity Dispersion in Rock with Dual Porosity
Wang, H F; Berryman, J G
2002-03-28
The high frequency behavior of the bulk modulus of fluid-saturated rock can be obtained from a double-porosity constitutive model, which is a direct conceptual extension of Biot's (1941) constitutive equations and which provides additional stiffening due to unrelaxed induced pore pressures in the soft porosity phase. Modeling the stiffening of the shear modulus at high frequency requires an effective medium average over the unequal induced pore pressures in cracks of different orientations. The implicit assumptions are that pore fluid equilibration does not occur between cracks of different orientations and between cracks and porous matrix. The correspondence between the constitutive equations of Berryman and Wang (1995) and Mavko and Jizba (1991) is explicitly noted.
NASA Astrophysics Data System (ADS)
Fernø, M. A.; Gauteplass, J.; Hauge, L. P.; Abell, G. E.; Adamsen, T. C. H.; Graue, A.
2015-09-01
Here we show for the first time the combined positron emission tomography (PET) and computed tomography (CT) imaging of flow processes within porous rocks to quantify the development in local fluid saturations. The coupling between local rock structure and displacement fronts is demonstrated in exploratory experiments using this novel approach. We also compare quantification of 3-D temporal and spatial water saturations in two similar CO2 storage tests in sandstone imaged separately with PET and CT. The applicability of each visualization technique is evaluated for a range of displacement processes, and the favorable implementation of combining PET/CT for laboratory core analysis is discussed. We learn that the signal-to-noise ratio (SNR) is over an order of magnitude higher for PET compared with CT for the studied processes.
NASA Astrophysics Data System (ADS)
Han, Tongcheng; Gurevich, Boris; Pervukhina, Marina; Clennell, Michael Ben; Zhang, Junfang
2016-04-01
Knowledge about the pressure dependency of elastic and electrical properties is important for a variety of geophysical applications. We present a technique to invert for the stiff and compliant porosity from velocity measurements made as a function of differential pressure on saturated sandstones. A dual porosity concept is used for dry rock compressibility and a squirt model is employed for the pressure and frequency dependent elastic properties of the rocks when saturated. The total porosity obtained from inversion shows satisfactory agreement with experimental results. The electrical cementation factor was determined using the inverted porosity in combination with measured electrical conductivity. It was found that cementation factor increased exponentially with increasing differential pressure during isostatic loading. Elastic compressibility, electrical cementation factor and electrical conductivity of the saturated rocks correlate linearly with compliant porosity, and electrical cementation factor and electrical conductivity exhibit linear correlations with elastic compressibility of the saturated rocks under loading. The results show that the dual porosity concept is sufficient to explain the pressure dependency of elastic, electrical and joint elastic-electrical properties of saturated porous sandstones.
NASA Astrophysics Data System (ADS)
Schulze-Makuch, Dirk; Cherkauer, Douglas S.
Previous studies have shown that hydraulic conductivity of an aquifer seems to increase as the portion of the aquifer tested increases. To date, such studies have all relied on different methods to determine hydraulic conductivity at each scale of interest, which raises the possibility that the observed increase in hydraulic conductivity is due to the measurement method, not to the scale. This study analyzes hydraulic conductivity with respect to scale during individual aquifer tests in porous, heterogeneous carbonate rocks in southeastern Wisconsin, USA. Results from this study indicate that hydraulic conductivity generally increases during an individual test as the volume of aquifer impacted increases, and the rate of this increase is the same as the rate of increase determined by using different measurement methods. Thus, scale dependence of hydraulic conductivity during single tests does not depend on the method of measurement. This conclusion is supported by 22 of 26 aquifer tests conducted in porous-flow-dominated carbonate units within the aquifer. Instead, scale dependency is probably caused by heterogeneities within the aquifer, a conclusion supported by digital simulation. All of the observed types of hydraulic-conductivity variations with scale during individual aquifer tests can be explained by a conceptual model of a simple heterogeneous aquifer composed of high-conductivity zones within a low-conductivity matrix. Résumé Certaines études ont montré que la conductivité hydraulique d'un aquifère semble augmenter en même temps que la partie testée de l'aquifère s'étend. Jusqu'à présent, ces études ont toutes reposé sur des méthodes de détermination de la conductivité hydraulique différentes pour chaque niveau d'échelle, ce qui a conduit à penser que l'augmentation observée de la conductivité hydraulique pouvait être due aux méthodes de mesure et non à l'effet d'échelle. Cette étude analyse la conductivité hydraulique par
NASA Astrophysics Data System (ADS)
Wang, W.; Regueiro, R. A.
2014-12-01
The coupling between multiphase flow, heat transfer, and poromechanics in fractured geomaterials has aroused great interest in the areas of geomechanics, geoenvironmental engineering, and petroleum engineering. Relevant applications include nuclear waste repositories, geological sequestration of CO2, geothermal systems, and exploitation of shale gas reservoirs. The paper presents a fully coupled thermo-poro-mechanical (TPM) cohesive interface element (CIE) model, which can represent fluid and heat flow along and across the fracture, and shear/normal deformation of the fracture surfaces. The proposed model is then applied to analyze two popular geological engineering problems using the finite element method (FEM) with a small strain formulation. The first application is the fracturing process in organic-rich shale due to heating. In the finite element analysis, multiple horizontal microcracks parallel to the bedding plane are assumed to preexist in the porous source rock, and are represented by coupled TPM cohesive interface elements. The porous bulk rock is assumed to be homogeneous, isotropic (for the time being, with transverse isotropy a natural extension), and linearly elastic. The excess pore fluid pressure, which mainly causes the development of the fractures, is actually induced by the rapid decomposition of organic matter during heating according to the literature. However, the involved complex chemical reaction process is beyond the scope of the paper, and is therefore substituted by a fluid injection process within the cracks under room temperature (25C) and high temperature (400C) in the paper. We investigate the fracture propagation due to pore fluid pressure increase and the development of fracture-induced permeability. The second application is a nuclear waste repository in a partially saturated fractured rock. Multiphase transport of moisture and heat, thermally-induced stress, as well as the change of fracture apertures are investigated due to short
ERIC Educational Resources Information Center
Lee, Alice
This science unit is designed for limited- and non-English speaking students in a Chinese bilingual education program. The unit covers rock material, classification, characteristics of types of rocks, and rock cycles. It is written in Chinese and simple English. At the end of the unit there is a list of main terms in both English and Chinese, and…
Micro-poromechanics model of fluid-saturated chemically active fibrous media
Misra, Anil; Parthasarathy, Ranganathan; Singh, Viraj; Spencer, Paulette
2014-01-01
We have developed a micromechanics based model for chemically active saturated fibrous media that incorporates fiber network microstructure, chemical potential driven fluid flow, and micro-poromechanics. The stress-strain relationship of the dry fibrous media is first obtained by considering the fiber behavior. The constitutive relationships applicable to saturated media are then derived in the poromechanics framework using Hill’s volume averaging. The advantage of this approach is that the resultant continuum model accounts for the discrete nature of the individual fibers while retaining a form suitable for porous materials. As a result, the model is able to predict the influence of micro-scale phenomena, such as the fiber pre-strain caused by osmotic effects and evolution of fiber network structure with loading, on the overall behavior and in particular, on the poromechanics parameters. Additionally, the model can describe fluid-flow related rate-dependent behavior under confined and unconfined conditions and varying chemical environments. The significance of the approach is demonstrated by simulating unconfined drained monotonic uniaxial compression under different surrounding fluid bath molarity, and fluid-flow related creep and relaxation at different loading-levels and different surrounding fluid bath molarity. The model predictions conform to the experimental observations for saturated soft fibrous materials. The method can potentially be extended to other porous materials such as bone, clays, foams and concrete. PMID:25755301
Misra, Anil; Parthasarathy, Ranganathan; Singh, Viraj; Spencer, Paulette
2013-01-01
The authors have derived macroscale poromechanics parameters for chemically active saturated fibrous media by combining microstructure-based homogenization with Hill's volume averaging. The stress-strain relationship of the dry fibrous media is first obtained by considering the fiber behavior. The constitutive relationships applicable to saturated media are then derived in the poromechanics framework using Hill's Lemmas. The advantage of this approach is that the resultant continuum model assumes a form suited to study porous materials, while retaining the effect of discrete fiber deformation. As a result, the model is able to predict the influence of microscale phenomena such as fiber buckling on the overall behavior, and in particular, on the poromechanics constants. The significance of the approach is demonstrated using the effect of drainage and fiber nonlinearity on monotonic compressive stress-strain behavior. The model predictions conform to the experimental observations for articular cartilage. The method can potentially be extended to other porous materials such as bone, clays, foams, and concrete. PMID:25419475
Parthasarathy, Ranganathan; Singh, Viraj; Spencer, Paulette
2014-01-01
The authors have derived macroscale poromechanics parameters for chemically active saturated fibrous media by combining microstructure-based homogenization with Hill's volume averaging. The stress-strain relationship of the dry fibrous media is first obtained by considering the fiber behavior. The constitutive relationships applicable to saturated media are then derived in the poromechanics framework using Hill's Lemmas. The advantage of this approach is that the resultant continuum model assumes a form suited to study porous materials, while retaining the effect of discrete fiber deformation. As a result, the model is able to predict the influence of microscale phenomena such as fiber buckling on the overall behavior, and in particular, on the poromechanics constants. The significance of the approach is demonstrated using the effect of drainage and fiber nonlinearity on monotonic compressive stress-strain behavior. The model predictions conform to the experimental observations for articular cartilage. The method can potentially be extended to other porous materials such as bone, clays, foams, and concrete. PMID:25419475
NASA Astrophysics Data System (ADS)
Nurhandoko, Bagus Endar B.; Wardaya, Pongga Dikdya; Adler, John; Siahaan, Kisko R.
2012-06-01
Seismic wave parameter plays very important role to characterize reservoir properties whereas pore parameter is one of the most important parameter of reservoir. Therefore, wave propagation phenomena in pore media is important to be studied. By referring this study, in-direct pore measurement method based on seismic wave propagation can be developed. Porosity play important role in reservoir, because the porosity can be as compartment of fluid. Many type of porosity like primary as well as secondary porosity. Carbonate rock consist many type of porosity, i.e.: inter granular porosity, moldic porosity and also fracture porosity. The complexity of pore type in carbonate rocks make the wave propagation in these rocks is more complex than sand reservoir. We have studied numerically wave propagation in carbonate rock by finite difference modeling in time-space domain. The medium of wave propagation was modeled by base on the result of pattern recognition using artificial neural network. The image of thin slice of carbonate rock is then translated into the velocity matrix. Each mineral contents including pore of thin slice image are translated to velocity since mineral has unique velocity. After matrix velocity model has been developed, the seismic wave is propagated numerically in this model. The phenomena diffraction is clearly shown while wave propagates in this complex carbonate medium. The seismic wave is modeled in various frequencies. The result shows dispersive phenomena where high frequency wave tends to propagate in matrix instead pores. In the other hand, the low frequency waves tend to propagate through pore space even though the velocity of pore is very low. Therefore, this dispersive phenomena of seismic wave propagation can be the future indirect measurement technology for predicting the existence or intensity of pore space in reservoir rock. It will be very useful for the future reservoir characterization.
Prediction of crack density in porous-cracked rocks from elastic wave velocities
NASA Astrophysics Data System (ADS)
Byun, Ji-Hwan; Lee, Jong-Sub; Park, Keunbo; Yoon, Hyung-Koo
2015-04-01
The stability of structures that are built over rock is affected by cracks in the rock that result from weathering, thawing and freezing processes. This study investigates a new method for determining rock crack densities using elastic wave velocities. The Biot-Gassmann model, which consists of several elastic moduli and Poisson's ratio, was used to determine a theoretical equation to predict the crack density of rocks. Ten representative specimens were extracted from ten boreholes to highlight the spatial variability. Each specimen was characterized using X-Ray Diffraction (XRD) analysis. The specimens were carved into cylinders measuring 50 mm in diameter and 30 mm in height using an abrasion process. A laboratory test was performed to obtain the elastic wave velocity using transducers that can transmit and receive compressional and shear waves. The measured compressional wave and shear wave velocities were approximately 2955 m/s-5209 m/s and 1652 m/s-2845 m/s, respectively. From the measured elastic wave velocities, the analyzed crack density and crack porosity were approximately 0.051-0.185 and 0.03%-0.14%, respectively. The calculated values were compared with the results of previous studies, and they exhibit similar values and trends. The sensitivity of the suggested theoretical equation was analyzed using the error norm technique. The results show that the compressional wave velocity and the shear modulus of a particle are the most influential factors in this equation. The study demonstrates that rock crack density can be estimated using the elastic wave velocities, which may be useful for investigating the stability of structures that are built over rock.
Semiarid watershed response to low-tech porous rock check dams
Technology Transfer Automated Retrieval System (TEKTRAN)
Rock check dams are used throughout the world as a technique for mitigating erosion problems on degraded lands. Increasingly, they are being used in restoration efforts on rangelands in the southwestern US; however, their impact on watershed response and channel morphology is not well quantified. In...
The effects of rock heterogeneity on compaction localization in porous carbonates
NASA Astrophysics Data System (ADS)
Cilona, Antonino; Faulkner, Daniel Roy; Tondi, Emanuele; Agosta, Fabrizio; Mancini, Lucia; Rustichelli, Andrea; Baud, Patrick; Vinciguerra, Sergio
2014-10-01
Recent field-based studies document the presence of bed-parallel compaction bands within the Oligocene-Miocene carbonates of Bolognano Formation exposed at the Majella Mountain of central Italy. These compaction bands are interpreted as burial-related structures, which accommodate volumetric strain by means of grain rotation/sliding, grain crushing, intergranular pressure solution and pore collapse. In order to constrain the pressure conditions at which these compaction bands formed, and investigate the role exerted by rock heterogeneity (grain and pore size and cement amount) on compaction localization, we carried out a suite of triaxial compression experiments, under dry conditions and room temperature on representative host rock samples of the Bolognano Formation. The experiments were performed at confining pressures that are proxy of those experienced by the rock during burial (5-35 MPa). Cylinders were cored out from a sample of the carbonate lithofacies most commonly affected by natural compaction bands. Natural structures were sampled and compared to the laboratory ones. During the experiments, the samples displayed shear-enhanced compaction and strain hardening associated with various patterns of strain localization. The brittle-ductile transition occurred at 12.5 MPa whereas compaction bands nucleated at 25 MPa confining pressure. A positive correlation between confining pressure and the angle formed by the deformation bands and the major principal stress axis was documented. Additional experiments were performed at 25 MPa on specimens cored oblique (parallel and at 45°) to the bedding. Detailed microstructural analyses, performed on pristine and deformed rocks by using optical microscopy, scanning electron microscopy and X-ray computed microtomography techniques, showed that grain crushing and mechanical twinning are the dominant deformation processes in the laboratory structures. Conversely, pressure solution appears to be dominant in the natural
NASA Astrophysics Data System (ADS)
Li, Y.; Jackson, I.; David, E.; Schmitt, D. R.
2013-12-01
The stiffness of rocks is significantly affected by the presence of cracks as well as pore fluids, the latter potentially increasing the effective stiffness of cracks. Reversible pore-fluid flow within the crack network, occurring during seismic wave propagation, may result in strongly frequency dependent seismic properties. Theoretical models for fluid flow induced seismic wave dispersion have been proposed but have so far not been subject to thorough experimental testing. Soda-lime-silica glass beads, of ~300 μm diameter were sintered near the glass transition temperature to produce a synthetic analogue for sedimentary rock with low porosity (~2%) and a simpler microstructure. Widely distributed cracks with uniformly low aspect ratio (~0.0007) and crack porosity ~0.2% were introduced by quenching heated cylindrical samples into liquid water at room temperature. Combined use of low-frequency (mHz-Hz) forced oscillation techniques at the Australian National University with ultrasonic pulse transmission methods (MHz) at the University of Alberta, is allowing a broadband measurement of seismic velocities and attenuation on a thermally cracked glass-bead sample. A recent upgrade of the data acquisition system on the apparatus for forced oscillation measurements is providing improved precision in determining shear and Young's moduli, measured at seismic frequencies, reveal a strong systematic variation with effective pressure (Peff=Pc-Pf) and some relaxation at longer oscillation periods tentatively attributed to fluid flow. Under water-saturated conditions, at low frequencies, both shear and Young's moduli are noticeably higher than under dry or argon-saturated conditions, possibly attributed to spatial restricted flow of water during forced-oscillation tests. Ongoing measurement of ultrasonic velocities should thus provide the 'intermediate' to 'high' frequency bounds on elastic moduli.
Seismic attenuation in fractured porous media: insights from a hybrid numerical and analytical model
NASA Astrophysics Data System (ADS)
Ekanem, A. M.; Li, X. Y.; Chapman, M.; Main, I. G.
2015-04-01
Seismic attenuation in fluid-saturated porous rocks can occur by geometric spreading, wave scattering or the internal dissipation of energy, most likely due to the squirt-flow mechanism. In principle, the pattern of seismic attenuation recorded on an array of sensors contains information about the medium, in terms of material heterogeneity and anisotropy, as well as material properties such as porosity, crack density, and pore-fluid composition and mobility. In practice, this inverse problem is challenging. Here we provide some insights into the effects of internal dissipation by analysing synthetic data produced by a hybrid numerical and analytical model for seismic wave propagation in a fractured medium embedded within a layered geological structure. The model is made up of one anisotropic and three isotropic horizontal layers. The anisotropic layer consists of a porous, fluid-saturated material containing vertically aligned inclusions representing a set of fractures. This combination allows squirt-flow to occur between the pores in the matrix and the model fractures. Our results show that the fluid mobility and the associated relaxation time of the fluid-pressure gradient control the frequency range over which attenuation occurs. This induced attenuation increases with incidence angle and azimuth away from the fracture strike-direction. Azimuthal variations in the induced attenuation are elliptical allowing the fracture orientations to be obtained from the axes of the ellipse. These observations hold out the potential of using seismic attenuation as an additional diagnostic in the characterisation of rock formations for a variety of applications including hydrocarbon exploration and production, subsurface storage of CO2, and geothermal energy extraction.
NASA Astrophysics Data System (ADS)
Pan, Li-Hua; Hou, Peng-Fei; Chen, Jia-Yun
2016-08-01
The 2D steady-state solutions regarding the expressions of stress and strain for fluid-saturated, orthotropic, poroelastic plane are derived in this paper. For this object, the general solutions of the corresponding governing equation are first obtained and expressed in harmonic functions. Based on these compact general solutions, the suitable harmonic functions with undetermined constants for line fluid source in the interior of infinite poroelastic body and a line fluid source on the surface of semi-infinite poroelastic body are presented, respectively. The fundamental solutions can be obtained by substituting these functions into the general solution, and the undetermined constants can be obtained by the continuous conditions, equilibrium conditions and boundary conditions.
NASA Astrophysics Data System (ADS)
Pollet-Villard, M.; Daval, D.; Saldi, G.; Knauss, K.; Wild, B.; Guyot, F. J.; Cabié, M.; Fritz, B.
2014-12-01
Prediction of reaction kinetics of fluid/rock interactions represents a critical issue for several geological and engineering concerns. In the specific context of geothermal energy, the relative intensities of primary mineral leaching and secondary phase formation significantly affect porosity and permeability of the reservoir, thereby influencing its hydraulic performance and the efficiency of a geothermal site. Moreover, it is noteworthy that in general, the circulation of aqueous fluids induces only modest modifications of their chemical composition, which slightly deviate from an equilibrium state. Therefore, fluid rock interactions take place at close-to-equilibrium conditions, where the rate-affinity relations are poorly known and intensively debated [1]. To unravel these points, in the context of the geothermal power station of Soultz-sous-Forêts (Alsace, France), our strategy consists in (1) investigating the dissolution of the main cleavages of K-spar, one of the prevalent primary minerals in the reservoir, in order to decipher the impact of crystallographic orientation and microstructural surface modifications on the dissolution kinetics and (2) proposing a relation between K-spar dissolution rate and the Gibbs free energy of reaction (∆G) over a wide range of ∆G conditions. Our experimental work relies on a coupled approach which combines classical experiments of K-spar dissolution monitored by aqueous chemical analyses (ICPAES) and innovative techniques of nm to μm-scale characterization of solid surface (FIB-TEM, AFM, VSI) [2]. Our ongoing experiments evidence that K-spar dissolution is an anisotropic process, with faces (-1-1 1) dissolving up to ten times faster than the slowest (001) faces. The complex evolution of surface topography with the occurrence of etch pits is at odds with the shrinking core model implemented in most of reactive-transport codes, representing a possible cause of an apparent modification of silicate dissolution rate over
Fluid-Rock Characterization and Interactions in NMR Well Logging
Hirasaki, George J.; Mohanty, Kishore, K.
2001-07-13
The objective of this project is to characterize the fluid properties and fluid-rock interactions that are needed for formation evaluation by NMR well logging. This is the first annual progress report submitted to the DOE. It reports on the work completed during the reporting period even if it may have started before this period. This project is a partnership between Professor George J. Hirasaki at Rice University and Professor Kishore Mohanty at University of Houston. In addition to the DOE, this project is supported by a consortium of oil companies and service companies. The fluid properties characterization has emphasized the departure of live oils from correlations based on dead oils. Also, asphaltic components can result in a difference between the T1 and T2 relaxation time distributions as well as reduce the hydrogen index. The fluid rock characterizations that are reported here are the effects of wettability and internal magnetic field gradients. A pore reconstruction method ha s been developed to recreate three-dimensional porous media from two-dimensional images that reproduce some of their key statistical properties. A Monte Carlo simulation technique has been developed to calculate the magnetization decay in fluid saturated porous media given their pore structure.
Effect of isolated fractures on accelerated flow in unsaturated porous rock
Su, G.W.; Nimmo, J.R.; Dragila, M.I.
2003-01-01
Fractures that begin and end in the unsaturated zone, or isolated fractures, have been ignored in previous studies because they were generally assumed to behave as capillary barriers and remain nonconductive. We conducted a series of experiments using Berea sandstone samples to examine the physical mechanisms controlling flow in a rock containing a single isolated fracture. The input fluxes and fracture orientation were varied in these experiments. Visualization experiments using dyed water in a thin vertical slab of rock were conducted to identify flow mechanisms occurring due to the presence of the isolated fracture. Two mechanisms occurred: (1) localized flow through the rock matrix in the vicinity of the isolated fracture and (2) pooling of water at the bottom of the fracture, indicating the occurrence of film flow along the isolated fracture wall. These mechanisms were observed at fracture angles of 20 and 60 degrees from the horizontal, but not at 90 degrees. Pooling along the bottom of the fracture was observed over a wider range of input fluxes for low-angled isolated fractures compared to high-angled ones. Measurements of matrix water pressures in the samples with the 20 and 60 degree fractures also demonstrated that preferential flow occurred through the matrix in the fracture vicinity, where higher pressures occurred in the regions where faster flow was observed in the visualization experiments. The pooling length at the terminus of a 20 degree isolated fracture was measured as a function of input flux. Calculations of the film flow rate along the fracture were made using these measurements and indicated that up to 22% of the flow occurred as film flow. These experiments, apparently the first to consider isolated fractures, demonstrate that such features can accelerate flow through the unsaturated zone and should be considered when developing conceptual models.
NASA Astrophysics Data System (ADS)
Stefansson, A.
2013-12-01
Reaction and reactive transport modeling is becoming an increasingly popular method to study fluid-rock interaction and fluid transport on small to large scales. In this study, fluid-rock experiments were carried out and the observations compared with the results of reaction and reactive transport models. The systems studied included fluid-rock interaction of olivine on one hand and basaltic glass on the other hand with dilute aqueous solutions containing CO2 at acid to neutral pH and temperatures from ambient to 250 °C. The experiments were conducted using batch type experiments in closed reactors and 1-D plug experiments in flow-through reactors and the solution chemistry, the reaction progress, secondary mineralization and porosity changes analyzed as a function of time. The reaction and 1-D reactive transport simulations were conducted with the aid of the PHREEQC program. For the simulations the thermodynamic database for mineral reactions was largely updated and the kinetics of mineral dissolution as well as mineral nucleation and crystal growth was incorporated. According to the experimental results and the reactive transport simulations, olivine and basaltic glass progressively dissolves forming secondary minerals and solutes that are partially transported out of them column (system). The exact reaction path was found to depend on solution composition and pH and reaction progress (time). The mass movement of the system at a particular steady state as well as porosity changes may be divided into three stages. Stage I is characterized by initial olivine or basaltic glass leaching, stage II is characterized by progressive mineral formation and decrease in porosity and stage III is characterized by remobilization of the previously formed secondary minerals and eventual increase in porosity. The reaction and reactive transport modeling was found to simulate reasonable the reaction path as a function of reaction time. However, exact mass movement and time
Gutierrez, Marte
2013-05-31
Colorado School of Mines conducted research and training in the development and validation of an advanced CO{sub 2} GS (Geological Sequestration) probabilistic simulation and risk assessment model. CO{sub 2} GS simulation and risk assessment is used to develop advanced numerical simulation models of the subsurface to forecast CO2 behavior and transport; optimize site operational practices; ensure site safety; and refine site monitoring, verification, and accounting efforts. As simulation models are refined with new data, the uncertainty surrounding the identified risks decrease, thereby providing more accurate risk assessment. The models considered the full coupling of multiple physical processes (geomechanical and fluid flow) and describe the effects of stochastic hydro-mechanical (H-M) parameters on the modeling of CO{sub 2} flow and transport in fractured porous rocks. Graduate students were involved in the development and validation of the model that can be used to predict the fate, movement, and storage of CO{sub 2} in subsurface formations, and to evaluate the risk of potential leakage to the atmosphere and underground aquifers. The main major contributions from the project include the development of: 1) an improved procedure to rigorously couple the simulations of hydro-thermomechanical (H-M) processes involved in CO{sub 2} GS; 2) models for the hydro-mechanical behavior of fractured porous rocks with random fracture patterns; and 3) probabilistic methods to account for the effects of stochastic fluid flow and geomechanical properties on flow, transport, storage and leakage associated with CO{sub 2} GS. The research project provided the means to educate and train graduate students in the science and technology of CO{sub 2} GS, with a focus on geologic storage. Specifically, the training included the investigation of an advanced CO{sub 2} GS simulation and risk assessment model that can be used to predict the fate, movement, and storage of CO{sub 2} in
Textural evidence for jamming and dewatering of a sub-surface, fluid-saturated granular flow
NASA Astrophysics Data System (ADS)
Sherry, T. J.; Rowe, C. D.; Kirkpatrick, J. D.; Brodsky, E. E.
2011-12-01
Sand injectites are spectacular examples of large-scale granular flows involving migration of hundreds of cubic meters of sand slurry over hundreds of meters to kilometers in the sub-surface. By studying the macro- and microstructural textures of a kilometer-scale sand injectite, we interpret the fluid flow regimes during emplacement and define the timing of formation of specific textures in the injected material. Fluidized sand sourced from the Santa Margarita Fm., was injected upward into the Santa Cruz Mudstone, Santa Cruz County, California. The sand injectite exposed at Yellow Bank Beach records emplacement of both hydrocarbon and aqueous sand slurries. Elongate, angular mudstone clasts were ripped from the wall rock during sand migration, providing evidence for high velocity, turbid flow. However, clast long axis orientations are consistently sub-horizontal suggesting the slurry transitioned to a laminar flow as the flow velocity decreased in the sill-like intrusion. Millimeter to centimeter scale laminations are ubiquitous throughout the sand body and are locally parallel to the mudstone clast long axes. The laminations are distinct in exposure because alternating layers are preferentially cemented with limonite sourced from later groundwater infiltration. Quantitative microstructural analyses show that the laminations are defined by subtle oscillations in grain alignment between limonite and non-limonite stained layers. Grain packing, size and shape distributions do not vary. The presence of limonite in alternating layers results from differential infiltration of groundwater, indicating permeability changes between the layers despite minimal grain scale differences. Convolute dewatering structures deform the laminations. Dolomite-cemented sand, a signature of hydrocarbon saturation, forms irregular bodies that cross-cut the laminations and dewatering structures. Laminations are not formed in the dolomite-cemented sand. The relative viscosity difference
Stability of sharp reaction fronts in porous rocks and implications for non-sharp reaction zones
NASA Astrophysics Data System (ADS)
Wangen, Magnus
2014-05-01
The flow of reactive fluids in the subsurface, like for instance acids, may create reaction fronts. A sharp reaction front is an idealization of the narrow zone where the reaction takes place. Narrow reaction zones are studied with a one-component reaction transport model, where a first order reaction changes the porosity. The porosity field is coupled to the permeability field, where an increasing porosity leads to an increasing permeability. Therefore, the reaction has a feed-back on the flow field. We have derived 1D approximate solutions for the change in concentration and porosity across the reaction zone. These solutions are used to derive a condition for reaction fronts to be narrow. The condition gives a minimum reaction rate necessary for 90% of the reaction to be restricted to the given area. Sharp fronts are idealizations of narrow fronts that are more amendable for analytical treatment. A condition has recently been derived for the stability of sharp reaction fronts in homogeneous porous medium using linear stability analysis. The condition gives that a perturbation of a flat reaction front of any wave-length becomes unstable if the permeability behind the front increases. The front instability grows faster for short wave lengths than for long wave lengths. Similarly, the perturbations of the front will die out if the permeability behind the front decreases, and short wave length perturbations will die out faster than long wave length perturbations. It is a condition that applies for both 2D and 3D porous media. Numerical experiments are shown that demonstrate the front stability criterion, when the fronts are narrow, but not sharp. The sharp front approximation turns out to be useful for the interpretation of reactions that are not sufficiently fast to give narrow reaction zones, when the reaction alters the porosity- and the permeability fields. Dissolution is an important example of reactions that increase the porosity and therefore the permeability
Rupture cascades in a discrete element model of a porous sedimentary rock.
Kun, Ferenc; Varga, Imre; Lennartz-Sassinek, Sabine; Main, Ian G
2014-02-14
We investigate the scaling properties of the sources of crackling noise in a fully dynamic numerical model of sedimentary rocks subject to uniaxial compression. The model is initiated by filling a cylindrical container with randomly sized spherical particles that are then connected by breakable beams. Loading at a constant strain rate the cohesive elements fail, and the resulting stress transfer produces sudden bursts of correlated failures, directly analogous to the sources of acoustic emissions in real experiments. The source size, energy, and duration can all be quantified for an individual event, and the population can be analyzed for its scaling properties, including the distribution of waiting times between consecutive events. Despite the nonstationary loading, the results are all characterized by power-law distributions over a broad range of scales in agreement with experiments. As failure is approached, temporal correlation of events emerges accompanied by spatial clustering. PMID:24580692
Quantitative model of vapor dominated geothermal reservoirs as heat pipes in fractured porous rock
Pruess, K.
1985-03-01
We present a numerical model of vapor-dominated reservoirs which is based on the well-known conceptual model of White, Muffler, and Truesdell. Computer simulations show that upon heat recharge at the base, a single phase liquid-dominated geothermal reservoir in fractured rock with low matrix permeability will evolve into a two-phase reservoir with B.P.D. (boiling point-for-depth) pressure and temperature profiles. A rather limited discharge event through cracks in the caprock, involving loss of only a few percent of fluids in place, is sufficient to set the system off to evolve a vapor-dominated state. The attributes of this state are discussed, and some features requiring further clarification are identified. 26 refs., 5 figs.
Spetzler, H.
1991-01-01
Here we report on the present status of the development of an apparatus for attenuation measurements. At the recent (April 91) DOE Meeting on Review of Underground Imaging the following paper was published. We quote it here for convenience since it contains rather recent information describing our efforts. Attenuation (Q{sup {minus}1}) measurements both in the field and the laboratory may prove to be powerful tools in geophysics in the future. Some seismic signatures based on Q (e.g. bright spots) only provide a qualitative indication of gas saturation because similar amplitude anomalies occur for small amounts of gas saturation as for large amounts. A better understanding of Q may provide additional information about rock properties required in finding petroleum reserves. Field Q measurements are also being utilized for engineering purposes (Murphy and Rosenbaum, 1989, Hatherly, 1986, and Young and Hill, 1986). Q measurements may even be used to identify toxic plumes. 18 refs., 10 figs.
Trapping non-wettable fluid in porous rock: Implication to CO2 sequestration
NASA Astrophysics Data System (ADS)
Kim, Y.; Yun, T.
2013-12-01
The residual saturation of CO2 mainly determines the effective storage capacity in geological formation whereas its transport and fate are dominated by fluid properties and pore characteristics. This experimental study evaluates the relative permeability of brine and non-wettable fluids in Berea sandstone. The surrogate fluids representing CO2 are continuously injected into the brine-saturated sandstone and the effluent is simultaneously separated to measure the residual volume. The variables under consideration include the viscosity and surface tension of injected fluids, porosity, anisotropy of rock, and injection pressure and the residual saturation of non-wettable fluids is quantified based on the proposed variables. Results highlight that the storage capacity can be readily modulated and maximized by controlling the cyclic injection, initial saturation of non-wettable fluids, and injection pressure.
NASA Astrophysics Data System (ADS)
Walker, E.; Glover, P. W. J.; Ruel, J.
2014-02-01
High-quality streaming potential coupling coefficient measurements have been carried out using a newly designed cell with both a steady state methodology and a new pressure transient approach. The pressure transient approach has shown itself to be particularly good at providing high-quality streaming potential coefficient measurements as each transient increase or decrease allows thousands of measurements to be made at different pressures to which a good linear regression can be fitted. Nevertheless, the transient method can be up to 5 times as fast as the conventional measurement approaches because data from all flow rates are taken in the same transient measurement rather than separately. Test measurements have been made on samples of Berea and Boise sandstone as a function of salinity (approximately 18 salinities between 10-5 mol/dm3 and 2 mol/dm3). The data have also been inverted to obtain the zeta potential. The streaming potential coefficient becomes greater (more negative) for fluids with lower salinities, which is consistent with existing measurements. Our measurements are also consistent with the high-salinity streaming potential coefficient measurements made by Vinogradov et al. (2010). Both the streaming potential coefficient and the zeta potential have also been modeled using the theoretical approach of Glover (2012). This modeling allows the microstructural, electrochemical, and fluid properties of the saturated rock to be taken into account in order to provide a relationship that is unique to each particular rock sample. In all cases, we found that the experimental data were a good match to the theoretical model.
Permeability Estimation of Porous Rock by Means of Fluid Flow Simulation and Digital Image Analysis
NASA Astrophysics Data System (ADS)
Winardhi, C. W.; Maulana, F. I.; Latief, F. D. E.
2016-01-01
Permeability plays an important role to determine the characteristics of how fluids flow through a porous medium which can be estimated using various methods. Darcy's law and the Kozeny-Carman equation are two of the most utilized methods in estimating permeability. In Darcy's law, permeability can be calculated by applying a pressure gradient between opposing sides of inlet-outlet of a certain direction. The permeability then depends on the fluid viscosity and the flowrate. The Kozeny-Carman equation is an empirical equation which depends on several parameters such as shape factor of the pore, tortuosity, specific surface area, and porosity to determine the permeability. For both methods, digital image obtained by means of Micro CT-Scan is used. In this research, the permeability estimation using the Darcy's law was conducted by simulating fluid flow through the digital image using Lattice Boltzmann Method (LBM). As for the Kozeny-Carman equation, digital image analysis was used to obtain the required parameters. Two Kozeny-Carman equations were used to calculate the permeability of the samples. The first equation (KC1) depends on pore shape factor, porosity, tortuosity, and specific surface area while the second equation (KC2) only depends on pore radius, porosity, and tortuosity. We investigate the methods by first testing on three simple pipe models which vary in the radii. By using the result from Darcy's law as a reference, we compare the results from the Kozeny-Carman equations. From the calculation, KC2 yield smaller difference to the reference. The three methods were then applied to the Fontainebleau sandstone to verify the previous result.
Mechanistic models of biofilm growth in porous media
NASA Astrophysics Data System (ADS)
Jaiswal, Priyank; Al-Hadrami, Fathiya; Atekwana, Estella A.; Atekwana, Eliot A.
2014-07-01
Nondestructive acoustics methods can be used to monitor in situ biofilm growth in porous media. In practice, however, acoustic methods remain underutilized due to the lack of models that can translate acoustic data into rock properties in the context of biofilm. In this paper we present mechanistic models of biofilm growth in porous media. The models are used to quantitatively interpret arrival times and amplitudes recorded in the 29 day long Davis et al. (2010) physical scale biostimulation experiment in terms of biofilm morphologies and saturation. The model pivots on addressing the sediment elastic behavior using the lower Hashin-Shtrikman bounds for grain mixing and Gassmann substitution for fluid saturation. The time-lapse P wave velocity (VP; a function of arrival times) is explained by a combination of two rock models (morphologies); "load bearing" which assumes the biofilm as an additional mineral in the rock matrix and "pore filling" which assumes the biofilm as an additional fluid phase in the pores. The time-lapse attenuation (QP-1; a function of amplitudes), on the other hand, can be explained adequately in two ways; first, through squirt flow where energy is lost from relative motion between rock matrix and pore fluid, and second, through an empirical function of porosity (φ), permeability (κ), and grain size. The squirt flow model-fitting results in higher internal φ (7% versus 5%) and more oblate pores (0.33 versus 0.67 aspect ratio) for the load-bearing morphology versus the pore-filling morphology. The empirical model-fitting results in up to 10% increase in κ at the initial stages of the load-bearing morphology. The two morphologies which exhibit distinct mechanical and hydraulic behavior could be a function of pore throat size. The biofilm mechanistic models developed in this study can be used for the interpretation of seismic data critical for the evaluation of biobarriers in bioremediation, microbial enhanced oil recovery, and CO2
Does wettability influence seismic wave propagation in liquid-saturated porous rocks?
NASA Astrophysics Data System (ADS)
Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe
2015-12-01
In order to evaluate the effects of different fluids, we measured the P and S wave speeds through a series of synthetic epoxy-bonded carbonate composites with porosities from 21.9 to 31.5 per cent and saturated with air, with kerosene and with brine. These observed speeds were converted to the saturated bulk and shear moduli and compared to predictions made using Biot and Gassmann formulations. The observed bulk moduli agreed with those calculated for both kerosene and brine saturation as did the high frequency shear modulus under kerosene saturation. The observed shear modulus under water saturation, however, was significantly lower than the prediction. After excluding the currently known mechanisms of shear weakening, we suspect this disparity may be due to variations in the wetting of the epoxy that coats the pores surfaces. The kerosene completely wets this surface while the brine is only weakly hydrophilic with a wetting angle of 73.6°. At the molecular scale, this means that Stoke's no-slip boundary condition may not always apply as has been more recently demonstrated by many researchers in other disciplines such as microfluidic engineering but the implications for wave propagation in liquid-saturated rocks have not been considered.
A simple, effective-medium model for water saturation in porous rocks
Berg, C.R.
1995-07-01
A general equation for water saturation is derived from effective-medium theory. A simpler low-frequency equation is also derived. Both equations are directly solvable for water saturation. The model should be applicable to any granular, water-wet formation. Additional relationships are derived specifically for application to shaly sandstones, but the model should be applicable to a wide range of rock types, water conductivities, and tool frequencies. In the derivation, hydrocarbons are included in the matrix component of the equation and the combined ``matrix`` elements are treated as resistors in parallel. The low-frequency equation is compared to various approaches to calculation of water saturation, such as Dual-Water, Waxman-Smits, and Bussian. The general equation is compared to three-component effective medium, porosity index, and complex refractive index models (CRIM). The model is proven to work on experimental data under a wide range of frequencies and water conductivities. It is recommended that the new saturation model be used for nearly all types of electrical saturation calculation, whether the measurements are from standard resistivity tools or from dielectric tools.
Time-dependent Deformation in Porous Rocks Driven by Chemo-mechanical Coupling
NASA Astrophysics Data System (ADS)
Meredith, P. G.; Brantut, N.; Heap, M. J.; Baud, P.
2015-12-01
We report results from triaxial deformation of porous sandstone and limestone conducted to determine the time-dependency of deformation. Experiments were run on water-saturated samples under constant differential stress (creep) conditions. In sandstone, the deformation is entirely brittle for all levels of stress and for all resulting strain rates. The strain rate during creep is very stress sensitive, with a change of only 20 MPa in differential stress producing three orders of magnitude change in strain rate. Failure occurs by localized shear faulting after an extended period of dilatant microcracking, as evidenced by the output of acoustic emissions. By contrast, the behaviour of limestone is more complex. At low effective pressure, the creep behavior is brittle and characterised by the same features as observed for sandstone; a decelerating phase of creep, followed by an inflection and then an accelerating creep phase leading to macroscopic failure. Similarly, only a small amount of inelastic strain is accommodated before failure, and P wave speeds measured throughout deformation decrease continuously, indicating a continuous increase in dilatant crack damage. At higher effective pressure, brittle creep still occurs, but the details of the time-dependent deformation behavior are quite different. First, the total amount of accumulated creep strain increases dramatically with decreasing strain rate, and no localized failure occurs even at these higher strains. Second, the rate of decrease in P wave speeds during deformation decreases with decreasing strain rate; indicating that less damage is accumulated per unit strain when the strain rate is lower. Third, complementary strain rate stepping experiments indicate that the deformation becomes more compactant at lower strain rates. Taken together, these observations suggest that rate-dependent compactive deformation mechanisms compete with dilatant subcritical crack growth during creep in limestone under low
NASA Astrophysics Data System (ADS)
Yurie Khachay, Professor; Mindubaev, Mansur
2016-04-01
One of the main problems of the study of the intrusion thermal effects on the maturation of the organic matter is to estimate the volume, intensity, thermal effects of the intrusion and its redistribution in porous media by convection. A numerical algorithm for solving the problem of the developed convection in two-dimensional and three-dimensional models of the porous medium depending on the incline angle is developed. It is defined that the convective stability in the medium decreases with increasing incline angle. It was found that depending on the incline angle the structure of convection from many cells for a flat horizontal layer changes and it transfers to more elongated structures along the layer. It is shown that depending on the incline angles, invading sill and imbedding volume of the porous medium it can be realized either stationary or non-stationary convection that provides a principal different thermal conditions of hydrocarbons maturation in the motherboard porous medium. We give numerical examples of the influence of the incline angle on the flow structure inside the porous inclusion. By the stationary convection the volume of the boundary layers between the convective sells increases. That can lead to increasing of the part of motherboard rocks that are outer the temperature conditions of oil catalysis and as a consequence to the overestimation of the deposits.
NASA Astrophysics Data System (ADS)
Varsakelis, Christos; Papalexandris, Miltiadis
2012-11-01
In this talk we present results from a numerical study of unsteady, shear flows of fluid-saturated granular materials in the presence of gravity. In our study, we employ a two-pressure, two-velocity continuum model for the mixtures of interest. The governing equations are integrated via a predictor-corrector algorithm that combines a projection method for the pressure of each phase and an interface-tracking scheme. Initially, a high particle concentration ball is placed between two parallel plates while the rest of the domain is filled with a carrier fluid. The mixture is set in motion by the horizontal movement of the upper plate with constant speed. Because of the developing shear stresses and the onset of the Rayleigh-Taylor instability, the ball deforms to a wavy finger-like shape whose length increases with time. Further, fluid entrainment produces a mushroom pattern in its frontal part. At the same time, this granular finger descends due to gravity and once it reaches the bottom plate it forms an asymmetric granular pile. This talk concludes with results from a parametric study with respect to the shear rate and the diameter of the particles.
Spetzler, Hartmut
2006-05-01
We have been exploring a new technology that is based on using low-frequency seismic attenuation data to monitor changes in fluid saturation conditions in two-fluid phase porous materials. The seismic attenuation mechanism is related to the loss of energy due to the hysteresis of resistance to meniscus movement (changes in surface tension, wettability) when a pore containing two fluids is stressed at very low frequencies (< 10 Hz). This technology has potential applications to monitoring changes in (1) leakage at buried waste sites, (2) contaminant remediation, and (3) flooding during enhanced petroleum recovery. We have concluded a three year field study at the Maricopa Agricultural Center site of the University of Arizona. Three sets of instruments were installed along an East-West line perpendicular to the 50m by 50m inigation site. Each set of instruments consisted of one three component seismometer and one tiltmeter. Microseisms and solid Earth-tides served as strain sources. The former have a power peak at a period of about 6 seconds and the tides have about two cycles per day. Installation of instruments commenced in late summer of 2002. The instruments operated nearly continuously until April 2005. During the fall of 2003 the site was irrigated with water and one year later with water containing 150 ppm of a biosurfactant additive. This biodegradable additive served to mimic a class of contaminants that change the surface tension of the inigation fluid. Tilt data clearly show tidal tilts superimposed on local tilts due to agricultural irrigation and field work. When the observed signals were correlated with site specific theoretical tilt signals we saw no anomalies for the water irrigation in 2003, but large anomalies on two stations for the surfactant irrigation in 2004. Occasional failures of seismometers as well as data acquisition systems contributed to less than continuous coverage. These data are noisier than the tilt data, but do also show possible
The Effect of Micrite and Micro-porous Skeletal Grains on the NMR Response of Carbonate Rocks
NASA Astrophysics Data System (ADS)
El Husseiny, A.; Knight, R. J.
2015-12-01
Nuclear magnetic resonance (NMR) measurements are commonly used to obtain an estimate of the pore size distribution in a saturated material. The NMR relaxation time T2 is inversely proportional to the surface-area-to-volume (S/V) ratio of the pore space. An estimate of pore size is obtained by defining a geometrical parameter (α) linking pore radius to (S/V) -1. For example, for a spherical pore α =3; i.e. the radius is equal to 3(S/V)-1. Carbonate rocks have a complex pore structure, being composed of many different materials. In this study we assessed the use of NMR to differentiate between the water held in the pore space in micrite (carbonate mud) and in the pore space of micro-porous skeletal grains, two major constituents in carbonates. We expected to find significant differences in NMR T2 measurements on these two materials, given that the pore and particle size of the micrite is ~100 times smaller than that of the skeletal grains. However, we found that the T2 distributions for the two materials overlap. Measurements of porosity and measurements of surface area with a nitrogen BET method showed that S/V of the two materials is very similar. The skeletal grains have very high surface area due to surface roughness and micro-porosity. The geometric factor (α), relating (S/V) -1 to pore radius, was found to be two orders of magnitude larger for the skeletal grains than for micrite. We concluded that knowledge of the geometric value is needed in order to accurately estimate pore size from NMR T2 data.
NASA Astrophysics Data System (ADS)
De Paoli, Marco; Zonta, Francesco; Soldati, Alfredo
2016-05-01
Solute convection in porous media at high Rayleigh-Darcy numbers has important fundamental features and may also bear implications for geological CO2 sequestration processes. With the aid of direct numerical simulations, we examine the role of anisotropic permeability on the distribution of solutal concentration in fluid saturated porous medium. Our computational analyses span over few decades of Rayleigh-Darcy number and confirm the linear scaling of Nusselt number that was previously found in the literature. In addition, we find that anisotropic permeability γ < 1, i.e., with vertical permeability smaller than horizontal permeability, effectively increases the Nusselt number compared with isotropic conditions. We link this seemingly counterintuitive effect with the occurring modifications to the flow topology in the anisotropic conditions. Finally, we use our data computed for the two-sided configuration (i.e., Dirichlet conditions on upper and lower boundaries) to examine the time evolution of solutal dynamics in the one-sided configuration, and we demonstrate that the finite-time (short-term) amount of CO2 that can be dissolved in anisotropic sedimentary rocks is much larger than in isotropic rocks.
NASA Astrophysics Data System (ADS)
Almrabat, Abdulhadi M.
The thesis presents the results of a study of the characterization and modeling of the stress and pore-fluid dependent acoustic properties of fractured porous rocks. A new laboratory High Pressure and High Temperature (HPHT) triaxial testing system was developed to characterize the seismic properties of sandstone under different levels of effective stress confinement and changes in pore-fluid composition. An intact and fractured of Berea sandstones core samples were used in the experimental studies. The laboratory test results were used to develop analytical models for stress-level and pore-fluid dependent seismic velocity of sandstones. Models for stress-dependent P and S-wave seismic velocities of sandstone were then developed based on the assumption that stress-dependencies come from the nonlinear elastic response of micro-fractures contained in the sample under normal and shear loading. The contact shear stiffness was assumed to increase linearly with the normal stress across a micro-fracture, while the contact normal stiffness was assumed to vary as a power law with the micro-fracture normal stress. Both nonlinear fracture normal and shear contact models were validated by experimental data available in the literature. To test the dependency of seismic velocity of sandstone on changes in pore-fluid composition, another series of tests were conducted where P and S-wave velocities were monitored during injection of supercritical CO 2 in samples of Berea sandstone initially saturated with saline water and under constant confining stress. Changes in seismic wave velocity were measured at different levels of supercritical CO2 saturation as the initial saline water as pore-fluid was displaced by supercritical CO 2. It was found that the P- iv wave velocity significantly decreased while the S-wave velocity remained almost constant as the sample supercritical CO2 saturation increased. The dependency of the seismic velocity on changes on pore fluid composition during
Mixed convection opposing flow in porous annulus
NASA Astrophysics Data System (ADS)
Salman, Ahmed N. J.; Kamangar, Sarfaraz; Al-Rashed, Abdullah A. A. A.; Khan, T. M. Yunus; Khaleed, H. M. T.
2016-06-01
The current work investigates the mixed convection flow in a vertical porous annulus embedded with fluid saturated porous medium. The annulus is isothermally heated discretely at 20%, 35% and 50% of the height of cylinder at the center of annulus. Darcy law with thermal non-equilibrium approach is considered. The governing partial differential equations are solved using Finite Element Method (FEM). The effects of Peclet number Pe and conductivity ratio Kr on heat transfer and fluid flow is discussed It is found that the applied velocity in the downward direction, in case of an opposing flow, does not allow the thermal energy to reach from a hot to a cold surface.
Bounds on Transport Coefficients of Porous Media
Berryman, J G
2005-03-21
An analytical formulation of conductivity bounds by Bergman and Milton is used in a different way to obtain rigorous bounds on the real transport coefficients (electrical conductivity, thermal conductivity, and/or fluid permeability) of a fluid-saturated porous medium. These bounds do not depend explicitly on the porosity, but rather on two formation factors--one associated with the pore space and the other with the solid frame. Hashin-Shtrikman bounds for transport in random polycrystals of porous-material laminates will also be discussed.
NASA Astrophysics Data System (ADS)
Huang, Yi-Jyun; Wu, Cheng-Yueh; Hsieh, Bieng-Zih
2016-04-01
Gas hydrates are crystalline compounds in which guest gas molecules are trapped in host lattices of ice crystals. In Taiwan, the significant efforts have recently begun to evaluate the reserves of hydrate because the vast accumulations of gas hydrates had been recognized in southwestern offshore Taiwan. Class-3 type hydrate accumulations are referred to an isolated hydrate layer without an underlying zone of mobile fluids, and the entire hydrate layer may be well within the hydrate stability zone. The depressurization method is a useful dissociation method for gas production from Class-3 hydrate accumulations. The dissociation efficiency is controlled by the responses of hydrate to the propagating pressure disturbance, and the pressure propagation is relating to the amount (or saturation) of the mobile fluid in pore space of the hydrate layer. The purpose of this study is to study the effects of fluid saturation on the gas recovery from a class-3 hydrate accumulation using depressurization method. The case of a class-3 hydrate deposit of Yuan-An Ridge in southwestern offshore Taiwan is studied. The numerical method was used in this study. The reservoir simulator we used to study the dissociation of hydrate and the production of gas was the STARS simulator developed by CMG, which coupled heat transfer, geo-chemical, geo-mechanical, and multiphase fluid flow mechanisms. The study case of Yuan-An Ridge is located in southwestern offshore Taiwan. The hydrate deposit was found by the bottom simulating reflectors (BSRs). The geological structure of the studied hydrate deposit was digitized to build the geological model (grids) of the case. The formation parameters, phase behavior data, rock and fluid properties, and formation's initial conditions were assigned sequentially to grid blocks, and the completion and operation conditions were designed to wellbore blocks to finish the numerical model. The changes of reservoir pressure, temperature, saturation due to the hydrate
NASA Astrophysics Data System (ADS)
Kung, J.; Chien, Y. V.; Wu, W.; Dong, J.; Chang, Y.; Tsai, C.; Yang, M.; Wang, K.
2012-12-01
Previous studies showed that the voids and their geometry in the sedimentary rocks have great influence on the compressibility of rock, which reflects on its elastic velocities. Some models were developed to discuss the relations among velocity, porosity and void geometry. Therefore, the information of porosity, and void geometry and its distribution in rock is essential for understanding how the elastic properties of porous rocks affected by their poregeometry. In this study, we revisited a well-studied porous rock, Darley Dale sandstone, which has been studied by different groups with different purposes. Most of them are the deformation experiments. Different from previous studies, we measured the sound velocity of Darley dale sandstone under hydrostatic conditions. Also, we employed different techniques to investigate the pore geometry and porosity of Darley Dale sandstone to gain the insight of velocity changing behavior under the crustal conditions. Here, we measured a fully-dense copper block for a comparison. We performed X-ray CT scanning (XCT) to image the pore space of sandstone to construct the 3-D image of pore geometry, distribution and the pore size. The CT image data are allowed us to estimate the porosity of sandstone, too. One the other hand, the porosity of sample was measured using imbibitions method at ambient conditions and helium porosimeter at high pressure (up to 150 MPa). A set of specimens were cored from Darley Dale sandstone block. P and S wave velocities of specimens were measured at ambient conditions. We also performed high pressure velocity measurements on a selected rock specimen and a copper block up to 150 MPa under dry condition. Porosity of a set of rock specimens measured by imbibitions method was spanned from 6% to 15%, largely distributed within a range of 8%-11%. Compared the porosity obtained from three different techniques, imbibitions method, helium porosimeter and XCT, values from those measurements are in good agreement
Transient electrokinetic response of finely layered, fluid-filled porous media
NASA Astrophysics Data System (ADS)
Dietrich, M.; Delprat-Jannaud, F.; Garambois, S.
2012-04-01
Transient electrokinetic coupling phenomena created at the microscopic scale by the passage of seismic waves through fluid-saturated porous media generate conversions between seismic and electromagnetic (EM) energy which can be observed at the macroscopic scale. Far from being a mere scientific curiosity, transient seismoelectric or electroseismic phenomena are especially appealing to oil and gas exploration and hydrogeology as they open up the (fairly unique) possibility to characterize fluid-bearing geological formations with the resolution of seismic methods. Indeed, electrokinetic effects are likely to reconcile the sensitivity of electromagnetic exploration methods to fluids with the high resolving power of seismic prospecting techniques for structural imaging, thus naturally bridging the gap between these two important geophysical investigation means. Accounting for the electromagnetic dimension of the seismic wave propagation, or conversely, accounting for the seismic dimension of electromagnetic wave propagation gives new insights into the microstructure and physico-chemistry of fluid-filled porous or fractured media. We present full-waveform simulations of the coupled seismoelectromagnetic wave propagation in fluid-saturated, finely stratified porous media of interest to oil and gas exploration. Our simulation code uses the macroscopic governing equations derived by Pride [1994], which couple Biot's theory and Maxwell equations via flux/force transport equations. The synthetic seismoelectrograms and seismomagnetrograms are computed by extending the generalized reflection and transmission matrix method and by using a discrete wave number integration of the global reflectivity obtained in the frequency wave number domain. The theoretical signals clearly display the coseismic electric and magnetic fields travelling with the seismic disturbances as well as the seismic-to-electromagnetic conversions taking place at contrasts in solid and fluid properties. Our
NASA Astrophysics Data System (ADS)
Almrabat, Abdulhadi M.
The thesis presents the results of a study of the characterization and modeling of the stress and pore-fluid dependent acoustic properties of fractured porous rocks. A new laboratory High Pressure and High Temperature (HPHT) triaxial testing system was developed to characterize the seismic properties of sandstone under different levels of effective stress confinement and changes in pore-fluid composition. An intact and fractured of Berea sandstones core samples were used in the experimental studies. The laboratory test results were used to develop analytical models for stress-level and pore-fluid dependent seismic velocity of sandstones. Models for stress-dependent P and S-wave seismic velocities of sandstone were then developed based on the assumption that stress-dependencies come from the nonlinear elastic response of micro-fractures contained in the sample under normal and shear loading. The contact shear stiffness was assumed to increase linearly with the normal stress across a micro-fracture, while the contact normal stiffness was assumed to vary as a power law with the micro-fracture normal stress. Both nonlinear fracture normal and shear contact models were validated by experimental data available in the literature. To test the dependency of seismic velocity of sandstone on changes in pore-fluid composition, another series of tests were conducted where P and S-wave velocities were monitored during injection of supercritical CO 2 in samples of Berea sandstone initially saturated with saline water and under constant confining stress. Changes in seismic wave velocity were measured at different levels of supercritical CO2 saturation as the initial saline water as pore-fluid was displaced by supercritical CO 2. It was found that the P- iv wave velocity significantly decreased while the S-wave velocity remained almost constant as the sample supercritical CO2 saturation increased. The dependency of the seismic velocity on changes on pore fluid composition during
Characterization of Porous Medium Properties Using 2D NMR
NASA Astrophysics Data System (ADS)
Sun, Boqin; Dunn, Keh-Jim
2003-03-01
We have successfully applied the concept of 2D NMR to the characterization of properties of fluid-saturated porous medium. Using a two-windowed modified CPMG pulse sequence, we were able to explore the magnetic internal filed gradient distribution within the pore space of a fluid-saturated porous medium due to magnetic susceptibility contrast between the solid matrix and pore fluid. Similar scheme is used to identify and quantify different types of pore fluids, such as oil, water, and gas, based on the contrast in their diffusion coefficients. The magic angle spinning technique (MAS) can also be applied in the 2D NMR framework for delineating the chemical shift spectra of the pore fluids in a porous medium at different T1 or T2 relaxation times. The results can be displayed in a two-dimensional plot, with one axis being the T1 or T2 relaxation times, the other axis being the internal field gradient, diffusion coefficient, or chemical shift, and the third axis being the proton population. Our preliminary laboratory work indicates that the 2D NMR approach can be a powerful tool for the characterization of properties of fluid-saturated porous medium, such as fluid typing, oil viscosity determination, surface wettability, etc.
BISQ model based on a Kelvin-Voigt viscoelastic frame in a partially saturated porous medium
NASA Astrophysics Data System (ADS)
Nie, Jian-Xin; Ba, Jing; Yang, Ding-Hui; Yan, Xin-Fei; Yuan, Zhen-Yu; Qiao, Hai-Peng
2012-06-01
Taking into account three important porous media mechanisms during wave propagation (the Biot-flow, squirt-flow, and solid-skeleton viscoelastic mechanisms), we introduce water saturation into the dynamic governing equations of wave propagation by analyzing the effective medium theory and then providing a viscoelastic Biot/squirt (BISQ) model which can analyze the wave propagation problems in a partially viscous pore fluid saturated porous media. In this model, the effects of pore fluid distribution patterns on the effective bulk modulus at different frequencies are considered. Then we derive the wave dynamic equations in the time-space domain. The phase velocity and the attenuation coefficient equations of the viscoelatic BISQ model in the frequency-wavenumber domain are deduced through a set of plane harmonic solution assumptions. Finally, by means of numerical simulations, we investigate the effects of water saturation, permeability, and frequency on compressional wave velocity and attenuation. Based on tight sandstone and carbonate experimental observed data, the compressional wave velocities of partially saturated reservoir rocks are calculated. The compressional wave velocity in carbonate reservoirs is more sensitive to gas saturation than in sandstone reservoirs.
Tomutsa, L.; Brinkmeyer, A.; Doughty, D.
1993-04-01
A synergistic rock characterization methodology has been developed. It derives reservoir engineering parameters from X-ray tomography (CT) scanning, computer assisted petrographic image analysis, minipermeameter measurements, and nuclear magnetic resonance imaging (NMRI). This rock characterization methodology is used to investigate the effect of small-scale rock heterogeneity on oil distribution and recovery. It is also used to investigate the applicability of imaging technologies to the development of scaleup procedures from core plug to whole core, by comparing the results of detailed simulations with the images ofthe fluid distributions observed by CT scanning. By using the rock and fluid detailed data generated by imaging technology describe, one can verify directly, in the laboratory, various scaling up techniques. Asan example, realizations of rock properties statistically and spatially compatible with the observed values are generated by one of the various stochastic methods available (fuming bands) and are used as simulator input. The simulation results were compared with both the simulation results using the true rock properties and the fluid distributions observed by CT. Conclusions regarding the effect of the various permeability models on waterflood oil recovery were formulated.
Zhang, Yan-Hong; Ye, Shu-Jun; Wu, Ji-Chun
2014-06-01
Based on light transmission method in quantification of liquid saturation and its application in two-phase flow system, two groups of sandbox experiments were set up to study the migration of gas or Dense Non-Aqueous Phase Liquids (DNAPLs) in water saturated porous media. The migration of gas or DNAPL was monitored in the study. Two modified Light Intensity-Saturation (LIS) models for water/gas two-phase system were applied and verified by the experiment data. Moreover two new LIS models for NAPL/water system were developed and applied to simulate the DNAPL infiltration experiment data. The gas injection experiment showed that gas moved upward to the top of the sandbox in the form of 'fingering' and finally formed continuous distribution. The results of DNAPL infiltration experiment showed that TCE mainly moved downward as the result of its gravity, eventually formed irregular plume and accumulated at the bottom of the sandbox. The outcomes of two LIS models for water/gas system (WG-A and WG-B) were consistent to the measured data. The results of two LIS models for NAPL/water system (NW-A and NW-B) fit well with the observations, and Model NW-A based on assumption of individual drainage gave better results. It could be a useful reference for quantification of NAPL/water saturation in porous media system. PMID:25158486
NASA Astrophysics Data System (ADS)
Liu, Zhongxian; Liang, Jianwen; Wu, Chengqing
2016-06-01
Two dimensional diffraction of Rayleigh waves by a fluid-saturated poroelastic alluvial valley of arbitrary shape in a poroelastic half-space is investigated using the method of fundamental solutions (MFS). To satisfy the free surface boundary conditions exactly, Green's functions of compressional (PI and PII) and shear (SV) wave sources buried in a fluid-saturated poroelastic half-space are adopted. Next, the procedure for solving the scattering wave field is presented. It is verified that the MFS is of excellent accuracy and numerical stability. Numerical results illustrate that the dynamic response strongly depends on such factors as the incident frequency, the porosity of alluvium, the boundary drainage condition, and the valley shape. There is a significant difference between the diffraction of Rayleigh waves for the saturated soil case and for the corresponding dry soil case. The wave focusing effect both on the displacement and pore pressure can be observed inside the alluvial valley and the amplification effect seems most obvious in the case of higher porosity and lower frequency. Additionally, special attention should also be paid to the concentration of pore pressure, which is closely related to the site liquefaction in earthquakes.
NASA Astrophysics Data System (ADS)
Tian, Z.; Liu, J.
2015-12-01
Abstract: Although various factors have impact on the resistivity of subsurface rock formation, in depth range of general electrical prospecting, the conductive actions of rocks are basically realized relying on the aqueous solutions filled in the pores. Therefore, quantitatively studying the impact of the water level on rock resistivity is important to analyze and classify strata, investigate the underground structures. In this research, we proposed a feasible research on building electric property rock formation models with different porosity and water saturation based on theories of two-phase media. The propagation of audio-magnetotelluric (AMT) waves is simulated by using finite-difference (FD) scheme, and theoretic resistivity distribution is calculated on account of the response of AMT. According to a sequence of synthetic examples, through comparing and analyzing the simulated results with various porosity and water saturation respectively, we discuss the impact on layers resistivity while porosity and water saturation of rock stratum are changing. The results shows the extent that the mentioned factors can have impact on the propagation of AMT waves. Key words: audio-magnetotelluric modeling, two-phase media, porosity, water saturation, finite-difference
NASA Astrophysics Data System (ADS)
Webster, J. D.; Rebbert, C. R.
2001-01-01
Silicate glasses in crystal-free to crystal-poor melt inclusions from two coarse-grained granite xenoliths in alkaline volcanic rocks of Ascension Island were analyzed for 29 major, minor, and trace elements. For most constituents, the glass compositions are similar to those of the volcanic whole rocks and the xenoliths; we interpret the glasses to be chemically representative of granite melt. The melt inclusions are silicic; alkaline; contain low S and P abundances; and are enriched in H 2O, Cl, F, and Na relative to K. Inclusions from one xenolith contain 1.3 wt.% F, on average, whereas those from the other xenolith contain half that amount. The melt inclusion compositions allow investigation of the means and extent of granite magma evolution. The presence of magmatic fluid inclusions in close proximity to melt inclusions in the phenocrysts indicates that the granite melt was saturated in one or more volatile phases (Roedder and Coombs, 1967). The Cl and H 2O contents of the melt inclusions are consistent with the exsolution of volatile phase(s) at pressures of 3000 to 2000 bars. The glasses also show trends involving the (Cl/H 2O) ratio and the incompatible trace elements in melt that are indicative of melt evolution via crystal fractionation of volatile phase-saturated melt. These trends should be useful for recognizing volatile phase saturation in other Cl-enriched melts.
NASA Astrophysics Data System (ADS)
Kitamura, Keigo; Takahashi, Miki; Mizoguchi, Kazuo; Masuda, Koji; Ito, Hisao; Song, Sheng-Rong
2010-09-01
Changes in Vp/Vs (Poisson's ratio) around a fault are related to changes in the fluid transport properties of rocks, which play a significant role in seismogenic processes. Here we report a notable relationship between Vp/Vs and the permeability of porous fault-related rocks (Chelungpu fault, Taiwan) by direct and simultaneous measurement of elastic wave velocities (Vp and Vs) and permeability under increasing effective confining pressure (Peff) up to 25 MPa. Vp and Vs for all samples increased with Peff in the range up to 20 MPa, then were nearly constant as Peff increased to 25 MPa. Most silty sandstones with large proportions of fine-grained material showed positive correlations between Vp/Vs and permeability with rising pressure. On the other hand, well-sorted sandstones showed only slight changes in permeability with respect to Vp/Vs with rising pressure. We infer that grain size distributions, in particular the amount of silt- and clay-size grains, are responsible for the change in permeability with pressure as small particles clog pore networks with increasing Peff, causing the decrease in permeability. These findings may be useful to explain changes in permeability and pore pressure in the deep crust.
Cook, N.G.W.
1993-05-01
The principal objective of research on the seismic properties of reservoir rocks is to develop a basic understanding of the effects of rock microstructure and its contained pore fluids on seismic velocities and attenuation. Ultimately, this knowledge would be used to extract reservoir properties information such as the porosity, permeability, clay content, fluid saturation, and fluid type from borehole, cross-borehole, and surface seismic measurements to improve the planning and control of oil and gas recovery. This thesis presents laboratory ultrasonic measurements for three granular materials and attempts to relate the microstructural properties and the properties of the pore fluids to P- and S-wave velocities and attenuation. These experimental results show that artificial porous materials with sintered grains and a sandstone with partially cemented grains exhibit complexities in P- and S-wave attenuation that cannot be adequately explained by existing micromechanical theories. It is likely that some of the complexity observed in the seismic attenuation is controlled by details of the rock microstructure, such as the grain contact area and grain shape, and by the arrangement of the grain packing. To examine these effects, a numerical method was developed for analyzing wave propagation in a grain packing. The method is based on a dynamic boundary integral equation and incorporates generalized stiffness boundary conditions between individual grains to account for viscous losses and grain contact scattering.
NASA Astrophysics Data System (ADS)
Zaitsev, V. Yu.; Matveev, L. A.
2012-05-01
Presently, experimental evidence for extremely high strain-sensitivity of dissipation in rocks and similar microstructured materials is obtained both in laboratory and field conditions, in particular observations of pronounced amplitude modulation of the radiation of high-stability seismo-acoustic sources by tidal deformations of rocks with typical strains ~ 10-8. Such data indicate the presence of some thresholdless in amplitude and very efficient mechanism of strain-dependent dissipation. Conventionally, its origin is discussed in the context of frictional or adhesion-hysteretic loss at cracks in rocks. However, such dissipation mechanisms are not relevant to weak perturbations with displacements smaller than atomic size. Here, we revise thresholdless thermoelastic loss in dry cracks and viscous loss in saturated cracks taking into account wavy asperities typical of real cracks, which can create elongated (strip-like) contacts or almost closed "waists" in cracks. Thermoelastic loss at these contacts can be very efficient. Besides, the state of such contacts can already be strongly perturbed by the average strain which yet practically does not change the mean opening of the entire crack. Thus the dissipation localized at such contacts can be significantly affected by quite small average strain (e.g., 10-8), which is usually believed to be unable to produce any appreciable effect on the dissipation. Next, for liquid-saturated cracks, the presence of inner elongated asperities also drastically changes the character of squirt-type viscous dissipation. Velocity gradients and consequently the dissipation are localized in the vicinity of the nearly-closed waists which almost harness the liquid flow in the crack. This dissipation can be comparable in magnitude with viscous dissipation at the entire crack with smooth interface, but the decrement maximum is strongly shifted downwards on the frequency axis. Since near the waist the gap is much smaller than the average crack
NASA Astrophysics Data System (ADS)
Schultz, R. A.; Soliva, R.; Fossen, H.
2013-12-01
Deformation bands in porous rocks tend to develop into spatially organized arrays that display a variety of lengths and thicknesses, and their geometries and arrangements are of interest with respect to fluid flow in reservoirs. Field examples of deformation band arrays in layered clastic sequences suggest that the development of classic deformation band arrays, such as ladders and conjugate sets, and the secondary formation of through-going faults appear to be related to the physical properties of the host rock, the orientation of stratigraphic layers relative to the far-field stress state, and the evolution of the local stress state within the developing array. We have identified several field examples that demonstrate changes in band properties, such as type and orientation, as a function of one or more of these three main factors. Normal-sense deformation-band arrays such as those near the San Rafael Swell (Utah) develop three-dimensional ladder-style arrays at a high angle to the maximum compression direction; these cataclastic shear bands form at acute angles to the maximum compression not very different from that of the optimum frictional sliding plane, thus facilitating the eventual nucleation of a through-going fault. At Orange quarry (France), geometrically conjugate sets of reverse-sense compactional shear bands form with angles to the maximum compression direction that inhibit fault nucleation within them; the bands in this case also form at steep enough angles to bedding that stratigraphic heterogeneities within the deforming formation were apparently not important. Two exposures of thrust-sense ladders at Buckskin Gulch (Utah) demonstrate the importance of host-rock properties, bedding-plane involvement, and local stress perturbations on band-array growth. In one ladder, thrust-sense shear deformation bands nucleated along suitably oriented bedding planes, creating overprinting sets of compaction bands that can be attributed to layer properties and
NASA Astrophysics Data System (ADS)
Hallbauer-Zadorozhnaya, Valeriya; Santarato, Giovanni; Abu Zeid, Nasser
2015-08-01
In this paper, two separate but related goals are tackled. The first one is to demonstrate that in some saturated rock textures the non-linear behaviour of induced polarization (IP) and the violation of Ohm's law not only are real phenomena, but they can also be satisfactorily predicted by a suitable physical-mathematical model, which is our second goal. This model is based on Fick's second law. As the model links the specific dependence of resistivity and chargeability of a laboratory sample to the injected current and this in turn to its pore size distribution, it is able to predict pore size distribution from laboratory measurements, in good agreement with mercury injection capillary pressure test results. This fact opens up the possibility for hydrogeophysical applications on a macro scale. Mathematical modelling shows that the chargeability acquired in the field under normal conditions, that is at low current, will always be very small and approximately proportional to the applied current. A suitable field test site for demonstrating the possible reliance of both resistivity and chargeability on current was selected and a specific measuring strategy was established. Two data sets were acquired using different injected current strengths, while keeping the charging time constant. Observed variations of resistivity and chargeability are in agreement with those predicted by the mathematical model. These field test data should however be considered preliminary. If confirmed by further evidence, these facts may lead to changing the procedure of acquiring field measurements in future, and perhaps may encourage the design and building of a new specific geo-resistivity meter. This paper also shows that the well-known Marshall and Madden's equations based on Fick's law cannot be solved without specific boundary conditions.
Numerical investigation of the hydro-mechanical contribution to seismic attenuation in damaged rocks
NASA Astrophysics Data System (ADS)
Pollmann, Nele; Jänicke, Ralf; Renner, Jörg; Steeb, Holger
2016-04-01
The investigation of hydro-mechanical processes, in particular the modeling of seismic waves in fractured porous media, is essential for the physical interpretation of data obtained from seismic exploration. Here, we specifically investigate attenuation processes in fluid-saturated porous rock containing fracture networks to identify effective hydro-mechanical properties by numerical simulation. The main purpose of this work is the characterization of the overall hydro-mechanical properties by computational homogenization. We determine an effective Skempton coefficient by investigating the fluid pressure and the solid displacement of the skeleton saturated by compressible fluids. Fracture networks are stochastically generated to mimic geological in-situ situations. The fractures are approximated as ellipses with aspect ratios up to 1/100, i.e. they constitute thin and long hydraulic conduits with high permeabilities. Simulations are designed on the material scale with and without conservation of fluid mass in the control volume. Using computational homogenization approaches, we define an effective Skempton coefficient. A range of fracture networks with different characteristic properties is studied for different varieties of fractures. On the material scale we find strongly heterogeneous pressure propagation in the fracture network and the surrounding rock, respectively. The pressure diffusion is much faster in the fracture network than in the matrix, rendering the macroscopic hydro-mechanical behavior strongly time dependent. The effective Skempton coefficient converges to an ensemble-specific instantaneous value and to 1 for long-time studies. The ultimate objective of our study is to evaluate whether constraints on the structure of fracture networks can be deduced from observations of attenuation and its frequency dependence.
NASA Astrophysics Data System (ADS)
Besedina, A. N.; Vinogradov, E. A.; Gorbunova, E. M.; Kabychenko, N. V.; Svintsov, I. S.; Pigulevskiy, P. I.; Svistun, V. K.; Shcherbina, S. V.
2015-01-01
The first part of this work is dedicated to the response of different-age structures to lunisolar tides, which can be considered as a sounding signal for monitoring the state of fluid-saturated reservoirs. The complex approach to processing the data obtained at the testing sites of the Institute of Geosphere Dynamics of the Russian Academy of Sciences, Institute of Geophysics of the National Academy of Sciences of Ukraine, and KIEV station of the IRIS seismic network is applied for recognizing the tides against the hydrogeological, barometric, and seismic series. The comparative analysis of the experimental and theoretical values of the diurnal and semidiurnal tidal components in the time series of ground displacements is carried out. The tidal variations in the groundwater level are compared with the tidal components revealed in the ground displacement of the different-age structure of the Moscow Basin and Ukrainian Shield, which are parts of the East European artesian region. The differences in the tidal responses of the groundwater level and ground displacement probably suggest that the state of the massif is affected by certain additional factors associated, e.g., with the passage of earthquake-induced seismic waves and the changes in the hydrogeodynamic environment.
Capillary-Driven Solute Transport and Precipitation in Porous Media during Dry-Out
NASA Astrophysics Data System (ADS)
Ott, Holger; Andrew, Matthew; Blunt, Martin; Snippe, Jeroen
2014-05-01
The injection of dry or under-saturated gases or supercritical (SC) fluids into water bearing formations might lead to a formation dry-out in the vicinity of the injection well. The dry-out is caused by the evaporation/dissolution of formation water into the injected fluid and the subsequent transport of dissolved water in the injected fluid away from the injection well. Dry-out results in precipitation from solutes of the formation brine and consequently leads to a reduction of the rock's pore space (porosity) and eventually to a reduction of permeability near the injection well, or even to the loss of injectivity. Recently evidence has been found that the complexity of the pore space and the respective capillary driven solute transport plays a key role. While no effective-permeability (Keff) reduction was observed in a single-porosity sandstone, multi porosity carbonate rocks responded to precipitation with a strong reduction of Keff. The reason for the different response of Keff to salt precipitation is suspected to be in the exact location of the precipitate (solid salt) in the pore space. In this study, we investigate dry-out and salt precipitation due to supercritical CO2 injection in single and multi-porosity systems under near well-bore conditions. We image fluid saturation changes by means of μCT scanning during desaturation. We are able to observe capillary driven transport of the brine phase and the respective transport of solutes on the rock's pore scale. Finally we have access to the precipitated solid-salt phase and their distribution. The results can proof the thought models behind permeability porosity relationships K(φ) for injectivity modeling. The topic and the mechanisms we show are of general interest for drying processes in porous material such as soils and paper.
Homogenization of two fluid flow in porous media
Daly, K. R.; Roose, T.
2015-01-01
The macroscopic behaviour of air and water in porous media is often approximated using Richards' equation for the fluid saturation and pressure. This equation is parametrized by the hydraulic conductivity and water release curve. In this paper, we use homogenization to derive a general model for saturation and pressure in porous media based on an underlying periodic porous structure. Under an appropriate set of assumptions, i.e. constant gas pressure, this model is shown to reduce to the simpler form of Richards' equation. The starting point for this derivation is the Cahn–Hilliard phase field equation coupled with Stokes equations for fluid flow. This approach allows us, for the first time, to rigorously derive the water release curve and hydraulic conductivities through a series of cell problems. The method captures the hysteresis in the water release curve and ties the macroscopic properties of the porous media with the underlying geometrical and material properties.
NMR imaging of immiscible displacements in porous media
Majors, P.D.; Li, P.; Peters, E.J.
1995-12-31
We introduce a rapid, quantitative nuclear magnetic resonance imaging (NMRI) technique to resolve and measure multiple fluid phases in porous media. Liquids are resolved on the basis of their NMR spin-spin (T{sub 2}) relaxation times, and their intensities are corrected via attenuation analysis. The spatially resolved and corrected NMRI intensities are normalized to yield fluid saturations. In-situ saturation measurements are presented for three immiscible (oil and water) displacements in the same Berea sandstone core. NMRI and effluent recovery methods were compared. T{sub 2} of the displacement fluids were observed to be sensitive to displacement conditions.
Non-Darcian forced convection analysis in an annulus partially filled with a porous material
Chikh, S.; Boumedien, A.; Bouhadef, K.; Lauriat, G.
1995-12-01
Numerical solutions are presented for fully developed forced convection in concentric annuli partially filled with a porous medium. The porous medium is attached at the inner cylinder, which is maintained at uniform heat flux or at uniform wall temperature while the outer cylinder is adiabatic. The Brinkman-Forchheimer-extended Darcy model was used to model the flow inside the porous medium. The dependence of the fluid flow and heat transfer on several parameters of the problem is thoroughly documented. The inertia coefficient at which the inertial effects reduce the flow rate by 5% is determined as a function of the Darcy number for various thicknesses of the porous substrate. It is also shown that a critical thickness at which the value of the Nusselt number reaches a minimum does not exist if the effective thermal conductivity of the fluid-saturated porous medium is much higher than the fluid conductivity.
Hayat, Tasawar; Rafiq, Maimona; Ahmad, Bashir
2016-01-01
The objective of present paper is to examine the peristaltic flow of magnetohydrodynamic (MHD) Jeffrey fluid saturating porous space in a channel through rotating frame. Unlike the previous attempts, the flow formulation is based upon modified Darcy's law porous medium effect in Jeffrey fluid situation. In addition the impacts due to Soret and Dufour effects in the radiative peristaltic flow are accounted. Rosseland's approximation has been utilized for the thermal radiative heat flux. Lubrication approach is implemented for the simplification. Resulting problems are solved for the stream function, temperature and concentration. Graphical results are prepared and analyzed for different parameters of interest entering into the problems. PMID:26808387
Hayat, Tasawar; Rafiq, Maimona; Ahmad, Bashir
2016-01-01
The objective of present paper is to examine the peristaltic flow of magnetohydrodynamic (MHD) Jeffrey fluid saturating porous space in a channel through rotating frame. Unlike the previous attempts, the flow formulation is based upon modified Darcy's law porous medium effect in Jeffrey fluid situation. In addition the impacts due to Soret and Dufour effects in the radiative peristaltic flow are accounted. Rosseland’s approximation has been utilized for the thermal radiative heat flux. Lubrication approach is implemented for the simplification. Resulting problems are solved for the stream function, temperature and concentration. Graphical results are prepared and analyzed for different parameters of interest entering into the problems. PMID:26808387
NASA Astrophysics Data System (ADS)
Lebedev, M.; Clennell, B.; Pervukhina, M.; Shulakova, V.; Mueller, T.; Gurevich, B.
2009-04-01
Porous rocks in hydrocarbon reservoirs are often saturated with a mixture of two or more fluids. Interpretation of exploration seismograms requires understanding of the relationship between distribution of the fluids patches and acoustic properties of rocks. The sizes of patches as well as their distribution affect significantly the seismic response. If the size of the fluid patch is smaller than the diffusion wavelength then pressure equilibration is achieved and the bulk modulus of the rock saturated with a mixture is defined by the Gassmann equations (Gassmann, 1951) with the saturation-weighted average of the fluid bulk modulus given by Wood's law (Wood, 1955, Mavko et al., 1998). If the fluid patch size is much larger than the diffusion wavelength then there is no pressure communication between different patches. In this case, fluid-flow effects can be neglected and the overall rock may be considered equivalent to an elastic composite material consisting of homogeneous parts whose properties are given by Gassmann theory with Hill's equation for the bulk modulus (Hill, 1963, Mavko et al., 1998). At intermediate values of fluid saturation the velocity-saturation relationship is significantly affected by the fluid patch distribution. In order to get an improved understanding of factors influencing the patch distribution and the resulting seismic wave response we performed simultaneous measurements of P-wave velocities and rock sample CT imaging. The CT imaging allows us to map the fluid distribution inside rock sample during saturation (water imbibition). We compare the experimental results with theoretical predictions. In this paper we will present results of simultaneous measurements of longitudinal wave velocities and imaging mapping of fluid distribution inside rock sample during sample saturation. We will report results of two kinds of experiments: "dynamic" and "quasi static" saturation. In both experiments Casino Cores Otway Basin sandstone, Australia core
Study on Two-Phase Flow in Heterogeneous Porous Media by Light Transmission Method
NASA Astrophysics Data System (ADS)
Qiao, W.
2015-12-01
The non-aqueous phase liquid (NAPL) released to the subsurface can form residual ganglia and globules occupying pores and also accumulate and form pools, in which multiphase system forms. Determining transient fluid saturations in a multiphase system is essential to understand the flow characteristics of systems and to perform effective remediation strategies. As a non-destructive and non-invasive laboratory technique utilized for the measurement of liquid saturation in porous media, light transmission is of the lowest cost and safe. Utilization of Coupled Charge Device camera in light transmission systems provides a nearly instantaneous high-density array of spatial measurements over a very large dynamic range. The migration of NAPL and air spariging technique applied to remove NAPL in aquifer systems are typically two-phase flow problem. Because of the natural aquifer normally being heterogeneous, two 2-D sandboxes (Length55cm×width1.3cm×hight45cm) are set up to study the migration of gas and DNAPL in heterogeneous porous media based on light transmission method and its application in two-phase flow. Model D for water/gas system developed by Niemet and Selker (2001) and Model NW-A for water/NAPL system developed by Zhang et al. (2014) are applied for the calculation of fluid saturation in the two experiments, respectively. The gas injection experiments show that the gas moves upward in the irregular channels, piling up beneath the low permeability lenses and starting lateral movement. Bypassing the lenses, the gas moves upward and forms continuous distribution in the top of the sandbox. The faster of gas injects, the wider of gas migration will be. The DNAPL infiltration experiment shows that TCE mainly moves downward as the influence of gravity, stopping vertical infiltration when reaching the low permeability lenses because of its failure to overcome the capillary pressure. Then, TCE accumulates on the surface and starts transverse movement. Bypassing the
Thermal convection at low Rayleigh number from concentrated sources in porous media
Hickox, C. E.
1980-01-01
A simple mathematical theory is proposed for the analysis of natural convective motion, at low Rayleigh number, from a concentrated source of heat in a fluid-saturated porous medium. The theory consists of retaining only the leading terms of series expansions of the dependent variables in terms of the Rayleigh number, is thus linear, and is valid only in the limit of small Rayleigh number. Based on fundamental results for a variety of isolated sources, superposition is used to provide solutions for situations of practical interest. Special emphasis is given to the analysis of sub-seabed disposal of nuclear waste. 8 figures.
Polgari, Marta; Szabo, Zoltan; Szabo-Drubina, Magda; Hein, James R.; Yeh, Hsueh-Wen
2005-01-01
The mineralogical, chemical, and isotopic compositions were determined for a white tripoli from the footwall of the Jurassic Úrkút Mn-oxide ore deposit in the Bakony Mountains, Hungary. The tripoli consists of quartz and chalcedony, with SiO2 contents up to 100 wt.%; consequently, trace-element contents are very low. Oxygen isotopes and quartz crystallinity indicate a low-temperature diagenetic origin for this deposit. The tripoli was formed by dissolution of the carbonate portion of the siliceous (sponge spicules) Isztimér Limestone. Dissolution of the carbonate was promoted by inorganic and organic acids generated during diagensis and left a framework composed of diagenetic silica that preserved the original volume of the limestone layer. The relative enrichment of silica and high porosity is the result of that carbonate dissolution. The silty texture of this highly friable rock is due to the structurally weak silica framework.
Geophysical and transport properties of reservoir rocks. Summary annual report
Cook, N.G.W.
1990-04-29
Definition of petrophysical properties, such as porosity, permeability and fluid saturation, on the scale of meters, is the key to planning and control of successful Enhanced Oil Recovery techniques for domestic reservoirs. Macroscopic transport properties in reservoir rocks depend critically upon processes at the pore level involving interactions between the pore topology and the physical and chemical properties of the rock minerals and interstitial fluids. Similar interactions at the pore level determine also the macroscopic electrical and seismic properties of reservoir rocks. The objective of this research is to understand, using analysis and experiment, how fluids in pores affect the geophysical and sport properties of reservoir rocks. The goal is to develop equations-relating seismic and electrical properties of rock to the porosity, permeability and fluid saturations so as to invert geophysical images for improved reservoir management. Results from seismic measurements performed so far in this study suggest that even subtle changes in fluid contacts and the in-situ state of effective stress can be detected using geophysical imaging techniques. The experiments using Wood`s metal and wax are revealing the topology and sport properties of the pore space in clastic sedimentary rocks. A deeper understanding of these properties is considered-to be the key to the recovery of much of the mobile oil left in domestic reservoirs and to the effective management of enhanced oil recovery techniques. The results of Wood`s metal percolation tests indicate that most of the permeability of Berea sandstone resides in the critical percolating paths and these paths occupy only a small fraction of the total porosity. This result may have important implications for flooding in terms of override and efficiency as a function of saturation.
NASA Astrophysics Data System (ADS)
Baccheschi, P.; Pastori, M.; Margheriti, L.; Piccinini, D.
2016-03-01
The Abruzzi region is located in the Central Apennines Neogene fold-and-thrust belt and has one of the highest seismogenic potential in Italy, with high and diffuse crustal seismicity related to NE-SW oriented extension. In this study, we investigate the detailed spatial variation in shear wave splitting providing high-resolution anisotropic structure beneath the L'Aquila region. To accomplish this, we performed a systematic analysis of crustal anisotropic parameters: fast polarization direction (ϕ) and delay time (δt). We benefit from the dense coverage of seismic stations operating in the area and from a catalogue of several accurate earthquake locations of the 2009 L'Aquila seismic sequence, related to the Mw 6.1 2009 L'Aquila main shock, to describe in detail the geometry of the anisotropic volume around the active faults that ruptured. The spatial variations both in ϕ and δt suggest a complex anisotropic structure beneath the region caused by a combination of both structural- and stress-induced mechanisms. The average ϕ is NNW-SSE oriented (N141°), showing clear similarity both with the local fault strike and the SHmax. In the central part of the study area fast axes are oriented NW-SE, while moving towards the northeastern and northwestern sectors the fast directions clearly diverge from the general trend of NW-SE and rotate accordingly to the local fault strikes. The above-mentioned fault-parallel ϕ distribution suggests that the observed anisotropy is mostly controlled by the local fault-related structure. Toward the southeast fast directions become orthogonal both to strike of the local mapped faults and to the SHmax. Here, ϕ are predominantly oriented NE-SW; we interpret this orientation as due to the presence of a highly fractured and overpressurized rock volume which should be responsible of the 90° flips in ϕ and the increase in δt. Another possible mechanism for NE-SW orientation of ϕ in the southeastern sector could be ascribed to the
High-pressure mechanical instability in rocks
Byerlee, J.D.; Brace, W.F.
1969-01-01
At a confining pressure of a few kilobars, deformation of many sedimentary rocks, altered mafic rocks, porous volcanic rocks, and sand is ductile, in that instabilities leading to audible elastic shocks are absent. At pressures of 7 to 10 kilobars, however, unstable faulting and stick-slip in certain of these rocks was observed. This high pressure-low temperature instability might be responsible for earthquakes in deeply buried sedimentary or volcanic sequences.
NASA Astrophysics Data System (ADS)
Vlahou, Ioanna; Grae Worster, M.
2009-11-01
Frost damage, caused by the freezing of water-saturated media, affects plant roots, pavements and the foundations of buildings, and is a major erosional force in rocks. The process has been studied extensively in the case of soils, and mechanisms such as the formation of ice lenses have been identified. Here, we consider the freezing of water in a three-dimensional cavity in a water-saturated, porous, elastic rock. Initially, the expansion of water as it freezes causes flow away from the solidification front, into the porous rock. The Darcy flow in the porous medium controls the pressure field and therefore the freezing temperature. At later times, disjoining thermomolecular forces create a pre-melted film of water between the ice and the rock and cause flow of pore water from the surrounding rock into the cavity. We find that the disjoining forces between the ice and the rock have the dominant effect, so we focus on those later times when the cavity is ice-filled. We solve the coupled set of integro-differential equations governing the elastic stress in the rock and the flow through its pores to determine the evolution of the shape and extent of the ice-filled cavity.
Khan, Sami Ullah; Ali, Nasir; Abbas, Zaheer
2015-01-01
An analysis is carried out to study the heat transfer in unsteady two-dimensional boundary layer flow of a magnetohydrodynamics (MHD) second grade fluid over a porous oscillating stretching surface embedded in porous medium. The flow is induced due to infinite elastic sheet which is stretched periodically. With the help of dimensionless variables, the governing flow equations are reduced to a system of non-linear partial differential equations. This system has been solved numerically using the finite difference scheme, in which a coordinate transformation is used to transform the semi-infinite physical space to a bounded computational domain. The influence of the involved parameters on the flow, the temperature distribution, the skin-friction coefficient and the local Nusselt number is shown and discussed in detail. The study reveals that an oscillatory sheet embedded in a fluid-saturated porous medium generates oscillatory motion in the fluid. The amplitude and phase of oscillations depends on the rheology of the fluid as well as on the other parameters coming through imposed boundary conditions, inclusion of body force term and permeability of the porous medium. It is found that amplitude of flow velocity increases with increasing viscoelastic and mass suction/injection parameters. However, it decreases with increasing the strength of the applied magnetic field. Moreover, the temperature of fluid is a decreasing function of viscoelastic parameter, mass suction/injection parameter and Prandtl number. PMID:26657931
Khan, Sami Ullah; Ali, Nasir; Abbas, Zaheer
2015-01-01
An analysis is carried out to study the heat transfer in unsteady two-dimensional boundary layer flow of a magnetohydrodynamics (MHD) second grade fluid over a porous oscillating stretching surface embedded in porous medium. The flow is induced due to infinite elastic sheet which is stretched periodically. With the help of dimensionless variables, the governing flow equations are reduced to a system of non-linear partial differential equations. This system has been solved numerically using the finite difference scheme, in which a coordinate transformation is used to transform the semi-infinite physical space to a bounded computational domain. The influence of the involved parameters on the flow, the temperature distribution, the skin-friction coefficient and the local Nusselt number is shown and discussed in detail. The study reveals that an oscillatory sheet embedded in a fluid-saturated porous medium generates oscillatory motion in the fluid. The amplitude and phase of oscillations depends on the rheology of the fluid as well as on the other parameters coming through imposed boundary conditions, inclusion of body force term and permeability of the porous medium. It is found that amplitude of flow velocity increases with increasing viscoelastic and mass suction/injection parameters. However, it decreases with increasing the strength of the applied magnetic field. Moreover, the temperature of fluid is a decreasing function of viscoelastic parameter, mass suction/injection parameter and Prandtl number. PMID:26657931
NASA Astrophysics Data System (ADS)
Jaunat, J.; Dupuy, A.; Huneau, F.; Celle-Jeanton, H.; Le Coustumer, P.
2016-04-01
A numerical groundwater model of the weathered crystalline aquifer of Ursuya (a major water source for the north-western Pyrenees region, south-western France) has been computed based on monitoring of hydrological, hydrodynamic and meteorological parameters over 3 years. The equivalent porous media model was used to simulate groundwater flow in the different layers of the weathered profile: from surface to depth, the weathered layer (5 · 10-8 ≤ K ≤ 5 · 10-7 m s-1), the transition layer (7 · 10-8 ≤ K ≤ 1 · 10-5 m s-1, the highest values being along major discontinuities), two fissured layers (3.5 · 10-8 ≤ K ≤ 5 · 10-4 m s-1, depending on weathering profile conditions and on the existence of active fractures), and the hard-rock basement simulated with a negligible hydraulic conductivity (K = 1 10 -9 ). Hydrodynamic properties of these five calculation layers demonstrate both the impact of the weathering degree and of the discontinuities on the groundwater flow. The great agreement between simulated and observed hydraulic conditions allowed for validation of the methodology and its proposed use for application on analogous aquifers. With the aim of long-term management of this strategic aquifer, the model was then used to evaluate the impact of climate change on the groundwater resource. The simulations performed according to the most pessimistic climatic scenario until 2050 show a low sensitivity of the aquifer. The decreasing trend of the natural discharge is estimated at about -360 m3 y-1 for recharge decreasing at about -5.6 mm y-1 (0.8 % of annual recharge).
Weak Elastic Anisotropy in a Cracked Rock
NASA Astrophysics Data System (ADS)
Zhu, W.; Wong, T.
2006-12-01
Crack and textural fabrics have significant control over the development of mechanical anisotropy in a rock. Bedding in sedimentary rocks, cleavage in slates, preferred orientation of anisotropic minerals and anisotropic distribution of microcracks can all contribute to elastic anisotropy. Using Kachanov's (1992, 1993) formulation we analyzed the effects of an axisymmetric system of microcracks on seismic anisotropy. The elastic behavior of such a cracked rock is transversely isotropic, and its seismic properties can be characterized by the three Thomsen parameters. In this study we calculated the parameters ɛ, δ and γ under dry and saturated conditions. We derived analytic expressions for the model proposed by Sayers & Kachanov (1995), which assumes that the contribution from the fourth rank crack density tensor is negligible. This model predicts that the elliptic anisotropy condition ɛ=δ is obeyed in a dry rock. Guided by microstructural observations we adopted a two-parameter axisymmetric distribution to characterize the crack density, which predicts that δ and γ in a fluid saturated rock are related to ɛ in a nonlinear manner. All three Thomsen parameters are sensitively dependent on the crack density difference. While our model shows basic agreement with some of the laboratory data on seismic anisotropy in saturated shale, there are discrepancies which suggest that the petrofabric associated with preferred orientation of clay minerals and elastic anisotropy of the rock matrix may have considerable influence which should not be neglected in model. Preliminary comparison with borehole log data suggests rock physics tests which may be useful for interpreting the shear wave anisotropy observations.
Seismoelectric Phenomena in Fluid-Saturated Sediments
Block, G I; Harris, J G
2005-04-22
Seismoelectric phenomena in sediments arise from acoustic wave-induced fluid motion in the pore space, which perturbs the electrostatic equilibrium of the electric double layer on the grain surfaces. Experimental techniques and the apparatus built to study this electrokinetic (EK) effect are described and outcomes for studies of seismoelectric phenomena in loose glass microspheres and medium-grain sand are presented. By varying the NaCl concentration in the pore fluid, we measured the conductivity dependence of two kinds of EK behavior: (1) the electric fields generated within the samples by the passage of transmitted acoustic waves, and (2) the electromagnetic wave produced at the fluid-sediment interface by the incident acoustic wave. Both phenomena are caused by relative fluid motion in the sediment pores--this feature is characteristic of poroelastic (Biot) media, but not predicted by either viscoelastic fluid or solid models. A model of plane-wave reflection from a fluid-sediment interface using EK-Biot theory leads to theoretical predictions that compare well to the experimental data for both sand and glass microspheres.
ERIC Educational Resources Information Center
Rice, Dale; Corley, Brenda
1987-01-01
Discusses some of the ways that rocks can be used to enhance children's creativity and their interest in science. Suggests the creation of a dramatic production involving rocks. Includes basic information on sedimentary, igneous, and metamorphic rocks. (TW)
Exact effective-stress rules in rock mechanics
Berryman, J.G. )
1992-09-15
The standard paradigm for analysis of rock deformation arises from postulating the existence of an equivalent homogeneous porous rock.'' However, data on the pore-pressure dependence of fluid permeability for some rocks cannot be explained using any equivalent homogeneous porous medium. In contrast, a positive result shows that deformation measurements on both high-porosity sandstones and low-porosity granites can be explained adequately in terms of an equivalent two-constituent model of porous rocks, for which exact results have recently been discovered.
Thermo-hydro-mechanical processes in fractured rock formations during a glacial advance
NASA Astrophysics Data System (ADS)
Selvadurai, A. P. S.; Suvorov, A. P.; Selvadurai, P. A.
2015-07-01
The paper examines the coupled thermo-hydro-mechanical (THM) processes that develop in a fractured rock region within a fluid-saturated rock mass due to loads imposed by an advancing glacier. This scenario needs to be examined in order to assess the suitability of potential sites for the location of deep geologic repositories for the storage of high-level nuclear waste. The THM processes are examined using a computational multiphysics approach that takes into account thermo-poroelasticity of the intact geological formation and the presence of a system of sessile but hydraulically interacting fractures (fracture zones). The modelling considers coupled thermo-hydro-mechanical effects in both the intact rock and the fracture zones due to contact normal stresses and fluid pressure at the base of the advancing glacier. Computational modelling provides an assessment of the role of fractures in modifying the pore pressure generation within the entire rock mass.
Thermo-hydro-mechanical processes in fractured rock formations during glacial advance
NASA Astrophysics Data System (ADS)
Selvadurai, A. P. S.; Suvorov, A. P.; Selvadurai, P. A.
2014-11-01
The paper examines the coupled thermo-hydro-mechanical (THM) processes that develop in a fractured rock region within a fluid-saturated rock mass due to loads imposed by an advancing glacier. This scenario needs to be examined in order to assess the suitability of potential sites for the location of deep geologic repositories for the storage of high-level nuclear waste. The THM processes are examined using a computational multiphysics approach that takes into account thermo-poroelasticity of the intact geological formation and the presence of a system of sessile but hydraulically interacting fractures (fracture zones). The modeling considers coupled thermo-hydro-mechanical effects in both the intact rock and the fracture zones due to contact normal stresses and fluid pressure at the base of the advancing glacier. Computational modelling provides an assessment of the role of fractures that can modify the pore pressure generation within the entire rock mass.
Effect of segmental heating on mixed convection aiding flow in a vertical porous annulus
NASA Astrophysics Data System (ADS)
Salman, Ahmed N. J.; Al-Rashed, Abdullah A. A. A.; Kamangar, Sarfaraz; Khan, T. M. Yunus; Khaleed, H. M. T.
2016-06-01
Mixed convection flow in a vertical porous annulus embedded with fluid saturated porous medium for aiding is investigated. The annulus is imposed by 20%, 35% and 50% heater length at the bottom, mid and top sections of the annulus respectively. Darcy law with thermal non-equilibrium approach is considered. The governing partial differential equations are converted to simple algebraic equations using Finite Element Method (FEM). The effects of Peclet number Pe and conductivity ratio Kr on heat transfer and fluid flow behaviour are examined and it is found that for lower conductivity ratio, the heat transfer rate was higher with the increase in the Peclet number Pe, whereas this trend reversed when thermal conductivity ratio Kr is increased.
Natural convection flow in porous enclosure with localized heating from below with heat flux
NASA Astrophysics Data System (ADS)
Siddiki, Md. Noor-A.-Alam; Molla, Md. Mamun; Saha, Suvash C.
2016-07-01
Unsteady natural convection flow in a two dimensional fluid saturated porous enclosure with localized heating from below with heat flux, symmetrical cooling from the sides and the insulated top wall has been investigated numerically. The governing equations are the Darcy's law for the porous media and the energy equation for the temperature field has been considered. The non-dimensional Darcy's law in terms of the stream function is solved by finite difference method using the successive over-relaxation (SOR) scheme and the energy equation is solved by Alternative Direction Alternative (ADI) scheme. The uniform heat flux source is located centrally at the bottom wall. The numerical results are presented in terms of the streamlines and isotherms, as well as the local and average rate of heat transfer for the wide range of the Darcy's Rayleigh number and the length of the heat flux source at the bottom wall.
Dispersivity as an oil reservoir rock characteristic
Menzie, D.E.; Dutta, S.
1989-12-01
The main objective of this research project is to establish dispersivity, {alpha}{sub d}, as an oil reservoir rock characteristic and to use this reservoir rock property to enhance crude oil recovery. A second objective is to compare the dispersion coefficient and the dispersivity of various reservoir rocks with other rock characteristics such as: porosity, permeability, capillary pressure, and relative permeability. The dispersivity of a rock was identified by measuring the physical mixing of two miscible fluids, one displacing the other in a porous medium. 119 refs., 27 figs., 12 tabs.
NASA Astrophysics Data System (ADS)
Song, Yongjia; Hu, Hengshan; Rudnicki, John W.; Duan, Yunda
2016-09-01
An exact analytical solution is presented for the effective dynamic transverse shear modulus in a heterogeneous fluid-filled porous solid containing cylindrical inclusions. The complex and frequency-dependent properties of the dynamic shear modulus are caused by the physical mechanism of mesoscopic-scale wave-induced fluid flow whose scale is smaller than wavelength but larger than the size of pores. Our model consists of three phases: a long cylindrical inclusion, a cylindrical shell of poroelastic matrix material with different mechanical and/or hydraulic properties than the inclusion and an outer region of effective homogeneous medium of laterally infinite extent. The behavior of both the inclusion and the matrix is described by Biot's consolidation equations, whereas the surrounding effective medium which is used to describe the effective transverse shear properties of the inner poroelastic composite is assumed to be a viscoelastic solid whose complex transverse shear modulus needs to be determined. The determined effective transverse shear modulus is used to quantify the S-wave attenuation and velocity dispersion in heterogeneous fluid-filled poroelastic rocks. The calculation shows the relaxation frequency and relative position of various fluid saturation dispersion curves predicted by this study exhibit very good agreement with those of a previous 2-D finite-element simulation. For the double-porosity model (inclusions having a different solid frame than the matrix but the same pore fluid as the matrix) the effective shear modulus also exhibits a size-dependent characteristic that the relaxation frequency moves to lower frequencies by two orders of magnitude if the radius of the cylindrical poroelastic composite increases by one order of magnitude. For the patchy-saturation model (inclusions having the same solid frame as the matrix but with a different pore fluid from the matrix), the heterogeneity in pore fluid cannot cause any attenuation in the
NASA Astrophysics Data System (ADS)
Song, Yongjia; Hu, Hengshan; Rudnicki, John W.; Duan, Yunda
2016-06-01
An exact analytical solution is presented for the effective dynamic transverse shear modulus in a heterogeneous fluid-filled porous solid containing cylindrical inclusions. The complex and frequency-dependent properties of the dynamic shear modulus are caused by the physical mechanism of mesoscopic-scale wave-induced fluid flow whose scale is smaller than wavelength but larger than the size of pores. Our model consists of three phases: a long cylindrical inclusion, a cylindrical shell of poroelastic matrix material with different mechanical and/or hydraulic properties than the inclusion and an outer region of effective homogeneous medium of laterally infinite extent. The behavior of both the inclusion and the matrix is described by Biot's consolidation equations, whereas the surrounding effective medium which is used to describe the effective transverse shear properties of the inner poroelastic composite is assumed to be a viscoelastic solid whose complex transverse shear modulus needs to be determined. The determined effective transverse shear modulus is used to quantify the S-wave attenuation and velocity dispersion in heterogeneous fluid-filled poroelastic rocks. The calculation shows the relaxation frequency and relative position of various fluid saturation dispersion curves predicted by this study exhibit very good agreement with those of a previous two-dimensional finite-element simulation. For the double-porosity model (inclusions having a different solid frame than the matrix but the same pore fluid as the matrix) the effective shear modulus also exhibits a size-dependent characteristic that the relaxation frequency moves to lower frequencies by two orders of magnitude if the radius of the cylindrical poroelastic composite increases by one order of magnitude. For the patchy-saturation model (inclusions having the same solid frame as the matrix but with a different pore fluid from the matrix), the heterogeneity in pore fluid cannot cause any attenuation in
Underground Research Laboratories for Crystalline Rock and Sedimentary Rock in Japan
Shigeta, N.; Takeda, S.; Matsui, H.; Yamasaki, S.
2003-02-27
The Japan Nuclear Cycle Development Institute (JNC) has started two off-site (generic) underground research laboratory (URL) projects, one for crystalline rock as a fractured media and the other for sedimentary rock as a porous media. This paper introduces an overview and current status of these projects.
ERIC Educational Resources Information Center
Henn, Cynthia A.
2004-01-01
There are many interpretations for the symbols that are seen in rock art, but no decoding key has ever been discovered. This article describes one classroom's experiences with a lesson on rock art--making their rock art and developing their own personal symbols. This lesson allowed for creativity, while giving an opportunity for integration…
ERIC Educational Resources Information Center
Rommel-Esham, Katie; Constable, Susan D.
2006-01-01
In this article, the authors discuss a literature-based activity that helps students discover the importance of making detailed observations. In an inspiring children's classic book, "Everybody Needs a Rock" by Byrd Baylor (1974), the author invites readers to go "rock finding," laying out 10 rules for finding a "perfect" rock. In this way, the…
ERIC Educational Resources Information Center
Barker, Rachel M.
One of a series of general interest publications on science topics, the booklet provides those interested in rock collecting with a nontechnical introduction to the subject. Following a section examining the nature and formation of igneous, sedimentary, and metamorphic rocks, the booklet gives suggestions for starting a rock collection and using…
Brusseau, M.L.; Narter, M.; Schnaar, G.; Marble, J.
2009-06-01
The objective of this study was to quantitatively characterize the impact of porous-medium texture on interfacial area between immiscible organic liquid and water residing within natural porous media. Synchrotron X-ray microtomography was used to obtain high-resolution, three-dimensional images of solid and liquid phases in packed columns. The image data were processed to generate quantitative measurements of organic-liquid/water interfacial area and of organic-liquid blob sizes. Ten porous media, comprising a range of median grain sizes, grain-size distributions, and geochemical properties, were used to evaluate the impact of porous-medium texture on interfacial area. The results show that fluid-normalized specific interfacial area (A{sub f}) and maximum specific interfacial area (A{sub m}) correlate very well to inverse median grain diameter. These functionalities were shown to result from a linear relationship between effective organic-liquid blob diameter and median grain diameter. These results provide the basis for a simple method for estimating specific organic-liquid/water interfacial area as a function of fluid saturation for a given porous medium. The availability of a method for which the only parameter needed is the simple-to-measure median grain diameter should be of great utility for a variety of applications.
Chamkha, A.J.
1997-10-01
There has been considerable interest in studying natural or buoyancy-induced flows in fluid-saturated porous media adjacent to surfaces in recent years. This interest stems from numerous possible industrial and technological applications. Examples of some applications include geothermal reservoirs, drying of porous solids, heat exchanger design, petroleum production, filtration, chemical catalytic reactor, nuclear waste repositories, and geophysical flows. Here, continuum equations governing steady, laminar, buoyancy-induced flow and heat transfer of a power-law fluid over a horizontal surface immersed in a uniform porosity and permeability porous medium are developed. These partial differential equations are transformed into ordinary differential equations by using a general similarity transformation for variable surface temperature and constant heat flux cases. The resulting equations are solved numerically by an implicit finite-difference method. Numerical results for typical velocity and temperature profiles are presented and discussed.
NASA Astrophysics Data System (ADS)
Singh, M.
2015-12-01
The instability of plane interface between two superposed Rivlin-Ericksen elastico-viscous fluids saturated through a porous medium has been studied to include the suspended (dust) particles effect. Following the linearized stability theory and normal mode analysis the dispersion relation is obtained. For stationary convection, the Rivlin-Ericksen elastico-viscous fluid behaves like Newtonian fluids. It found that for a potentially stable arrangement the Rivlin-Ericksen elastico-viscous fluid of different permeabilities in the presence of suspended particles in a porous medium is stable, whereas in a potentially unstable case instability of the system occurs. In the presence of a magnetic field for a potentially stable arrangement the system is always stable and for the potentially unstable arrangement, the magnetic field succeeds in stabilizing certain wave-number band which was unstable in the absence of the magnetic field.
NASA Astrophysics Data System (ADS)
Purcell, C. C.; Mur, A. J.; Delany, D.; Haljasmaa, I. V.; Soong, Y.; Harbert, W.
2011-12-01
The exploration of velocity differences in various fluid saturated rock types under reservoir conditions should prove to be useful in seismic monitoring of sequestration and hydrothermal sites. Different saturation values, along with mixtures of other common pore fluids could help delineate various areas of a CO2 flood or enhanced geothermal pressurization, in addition to estimating a minimum saturation amount needed to be seen in seismic surveys. We also explore the effects of varying parameters on the saturated velocities, including porosity, bulk frame composition, pressure, temperature, different pore filling phases, fluid mixtures, and compliant porosity. A software toolkit is currently in development that would allow exploration of these parameters to be easily achieved and visualized. Fluid substitution using Gassmann's equation (Gassmann [1]) is an important tool in the analysis of velocity dispersion in saturated rocks. Mavko and Jizba [2] created a model of squirt dispersion for elastic wave velocities at ultrasonic frequencies that predicts total dispersion for fluid filled rocks. Gurevich et al. [3] extend the Mavko-Jizba expressions to low fluid bulk modulus situations, such as gas filled rocks. These equations are typically used to calculate velocities of rocks filled with typical pore filling phases such as brine or gas. Purcell et al. [4] compared these equations to CO2 saturated limestone samples at reservoir pressures and temperatures. This paper compares the accuracy of these equations over various pressures and temperature ranges for a variety of rock types. Dry rock ultrasonic lab measurements of velocity have been made for carbonate, sandstone, rhyolite and coal and incorporated into a rock physics database. In addition, waveforms for each measurement have been used to estimate Q. Measurements were made between 2.3 and 50 MPa with generally a minimum of 40 measurements per sample completed. Various saturating phases, including supercritical CO
NASA Astrophysics Data System (ADS)
Beckermann, C.; Ramadhyani, S.; Viskanta, R.
1987-05-01
A numerical and experimental study is performed to analyze the steady-state natural convection fluid flow and heat transfer in a vertical rectangular enclosure that is partially filled with a vertical layer of a fluid-saturated porous medium. The flow in the porous layer is modeled utilizing the Brinkman-Forchheimer-extended Darcy equations. The numerical model is verified by conducting a number of experiments, with spherical glass beads as the porous medium and water and glycerin as the fluids, in rectangular test cells. The agreement between the flow visualization results and temperature measurements and the numerical model is, in general, good. It is found that the amount of fluid penetrating from the fluid region into the porous layer depends strongly on the Darcy (Da) and Rayleigh (Ra) numbers. For a relatively low product of Ra x Da, the flow takes place primarily in the fluid layer, and heat transfer in the porous layer is by conduction only. On other hand, fluid penetration into a relatively highly permeable porous layer has a significant impact on the natural convection flow patterns in the entire enclosure.
NASA Technical Reports Server (NTRS)
2004-01-01
This false-color image taken by NASA's Mars Exploration Rover Opportunity shows a rock informally named 'Earhart' on the lower slopes of 'Endurance Crater.' The rock was named after the pilot Amelia Earhart. Like 'Escher' and other rocks dotting the bottom of Endurance, scientists believe fractures in Earhart could have been formed by one of several processes. They may have been caused by the impact that created Endurance Crater, or they might have arisen when water leftover from the rock's formation dried up. A third possibility is that much later, after the rock was formed, and after the crater was created, the rock became wet once again, then dried up and developed cracks. Rover team members do not have plans to investigate Earhart in detail because it is located across potentially hazardous sandy terrain. This image was taken on sol 219 (Sept. 4) by the rover's panoramic camera, using its 750-, 530- and 430-nanometer filters.
Pore Scale Simulations of Rock Deformation, Fracture, and Fluid Flow in Three Dimensions
Wang, Herbert F.
2005-04-01
The pore-scale examination of rock deformation and fluid flow consisted of three separate tasks. (1) New laboratory measurements were made of poroelastic properties of Berea sandstone and a new method was developed to measure both the poroelastic constants and the hydraulic conductivity on the same sample of rock in a single test. (2) The second task was to develop constitutive theories of elastic and poroelastic properties of dual-porosity rocks and rocks with cracks. The new constitutive relations explain wave-velocity dispersion in fluid-saturated rock and the stiffening of shear modulus when dry rock is saturated. (3) The third task involved pore-scale percolation modeling of two-phase fluid flow in granular media. The model properly simulates fractal geometries of nonwetting clusters and saturations for flow in unstable gradients. The percolation model was coupled with a water-vapor diffusion model to produce saturation maps in a rock core during evaporative drying. The realistic patchy saturation was used in a heuristic model for predicting elastic properties of partially-saturated rock, which mimicked laboratory results.
ERIC Educational Resources Information Center
Prestwich, Dorothy; Sumrall, Joseph; Chessin, Debby A.
2010-01-01
It all began one Monday morning. Raymond could not wait to come to large group. In his hand, he held a chunk of white granite he had found. "Look at my beautiful rock!" he cried. The rock was passed around and examined by each student. "I wonder how rocks are made?" wondered one student. "Where do they come from?" asked another. At this moment, a…
NASA Technical Reports Server (NTRS)
Matveyev, S. N.
1986-01-01
Rock flows are defined as forms of spontaneous mass movements, commonly found in mountainous countries, which have been studied very little. The article considers formations known as rock rivers, rock flows, boulder flows, boulder stria, gravel flows, rock seas, and rubble seas. It describes their genesis as seen from their morphological characteristics and presents a classification of these forms. This classification is based on the difference in the genesis of the rubbly matter and characterizes these forms of mass movement according to their source, drainage, and deposit areas.
Rock Pore Structure as Main Reason of Rock Deterioration
NASA Astrophysics Data System (ADS)
Ondrášik, Martin; Kopecký, Miloslav
2014-03-01
Crashed or dimensional rocks have been used as natural construction material, decoration stone or as material for artistic sculptures. Especially old historical towns not only in Slovakia have had experiences with use of stones for construction purposes for centuries. The whole buildings were made from dimensional stone, like sandstone, limestone or rhyolite. Pavements were made especially from basalt, andesite, rhyolite or granite. Also the most common modern construction material - concrete includes large amounts of crashed rock, especially limestone, dolostone and andesite. However, rock as any other material if exposed to exogenous processes starts to deteriorate. Especially mechanical weathering can be very intensive if rock with unsuitable rock properties is used. For long it had been believed that repeated freezing and thawing in relation to high absorption is the main reason of the rock deterioration. In Slovakia for many years the high water absorption was set as exclusion criterion for use of rocks and stones in building industry. Only after 1989 the absorption was accepted as merely informational rock property and not exclusion. The reason of the change was not the understanding of the relationship between the porosity and rock deterioration, but more or less good experiences with some high porous rocks used in constructions exposed to severe weather conditions and proving a lack of relationship between rock freeze-thaw resistivity and water absorption. Results of the recent worldwide research suggest that understanding a resistivity of rocks against deterioration is hidden not in the absorption but in the structure of rock pores in relation to thermodynamic properties of pore water and tensile strength of rocks and rock minerals. Also this article presents some results of research on rock deterioration and pore structure performed on 88 rock samples. The results divide the rocks tested into two groups - group N in which the pore water does not freeze
NASA Technical Reports Server (NTRS)
2004-01-01
This image, taken by the panoramic camera on NASA's Mars Exploration Rover Spirit during the rover's trek through the 'Columbia Hills' at 'Gusev Crater,' shows the horizontally layered rock dubbed 'Tetl.' Scientists hope to investigate this rock in more detail, aiming to determine whether the rock's layering is volcanic or sedimentary in origin. If for some reason this particular rock is not favorably positioned for grinding and examination by the toolbox of instruments on the rover's robotic arm, Spirit will be within short reach of another similar rock, dubbed 'Coba.' Spirit took this image on its 264th martian day, or sol (Sept. 29, 2004). This is a false-color composite image generated from the panoramic camera's 750-, 530-, and 430-nanometer filters.
A versatile facility for laboratory studies of viscoelastic and poroelastic behaviour of rocks
Jackson, Ian; Schijns, Heather; Schmitt, Douglas R.; Mu Junjie; Delmenico, Alison
2011-06-15
Novel laboratory equipment has been modified to allow both torsional and flexural oscillation measurements at sub-microstrain amplitudes, thereby providing seismic-frequency constraints on both the shear and compressional wave properties of cylindrical rock specimens within the linear regime. The new flexural mode capability has been tested on experimental assemblies containing fused silica control specimens. Close consistency between the experimental data and the results of numerical modelling with both finite-difference and finite-element methods demonstrates the viability of the new technique. The capability to perform such measurements under conditions of independently controlled confining and pore-fluid pressure, with emerging strategies for distinguishing between local (squirt) and global (specimen-wide) fluid flow, will have particular application to the study of frequency-dependent seismic properties expected of cracked and fluid-saturated rocks of the Earth's upper crust.
A versatile facility for laboratory studies of viscoelastic and poroelastic behaviour of rocks
NASA Astrophysics Data System (ADS)
Jackson, Ian; Schijns, Heather; Schmitt, Douglas R.; Mu, Junjie; Delmenico, Alison
2011-06-01
Novel laboratory equipment has been modified to allow both torsional and flexural oscillation measurements at sub-microstrain amplitudes, thereby providing seismic-frequency constraints on both the shear and compressional wave properties of cylindrical rock specimens within the linear regime. The new flexural mode capability has been tested on experimental assemblies containing fused silica control specimens. Close consistency between the experimental data and the results of numerical modelling with both finite-difference and finite-element methods demonstrates the viability of the new technique. The capability to perform such measurements under conditions of independently controlled confining and pore-fluid pressure, with emerging strategies for distinguishing between local (squirt) and global (specimen-wide) fluid flow, will have particular application to the study of frequency-dependent seismic properties expected of cracked and fluid-saturated rocks of the Earth's upper crust.
A versatile facility for laboratory studies of viscoelastic and poroelastic behaviour of rocks.
Jackson, Ian; Schijns, Heather; Schmitt, Douglas R; Mu, Junjie; Delmenico, Alison
2011-06-01
Novel laboratory equipment has been modified to allow both torsional and flexural oscillation measurements at sub-microstrain amplitudes, thereby providing seismic-frequency constraints on both the shear and compressional wave properties of cylindrical rock specimens within the linear regime. The new flexural mode capability has been tested on experimental assemblies containing fused silica control specimens. Close consistency between the experimental data and the results of numerical modelling with both finite-difference and finite-element methods demonstrates the viability of the new technique. The capability to perform such measurements under conditions of independently controlled confining and pore-fluid pressure, with emerging strategies for distinguishing between local (squirt) and global (specimen-wide) fluid flow, will have particular application to the study of frequency-dependent seismic properties expected of cracked and fluid-saturated rocks of the Earth's upper crust. PMID:21721712
Thermo-mechanical pressurization of experimental faults in cohesive rocks during seismic slip
NASA Astrophysics Data System (ADS)
Violay, M.; Di Toro, G.; Nielsen, S.; Spagnuolo, E.; Burg, J. P.
2015-11-01
Earthquakes occur because fault friction weakens with increasing slip and slip rates. Since the slipping zones of faults are often fluid-saturated, thermo-mechanical pressurization of pore fluids has been invoked as a mechanism responsible for frictional dynamic weakening, but experimental evidence is lacking. We performed friction experiments (normal stress 25 MPa, maximal slip-rate ∼3 ms-1) on cohesive basalt and marble under (1) room-humidity and (2) immersed in liquid water (drained and undrained) conditions. In both rock types and independently of the presence of fluids, up to 80% of frictional weakening was measured in the first 5 cm of slip. Modest pressurization-related weakening appears only at later stages of slip. Thermo-mechanical pressurization weakening of cohesive rocks can be negligible during earthquakes due to the triggering of more efficient fault lubrication mechanisms (flash heating, frictional melting, etc.).
Fluid effects on seismic waves in hard rocks with fractures and in soft granular media
Berryman, James G.
2009-03-01
When fractures in otherwise hard rocks are filled with fluids (oil, gas, water, CO{sub 2}), the type and physical state of the fluid (liquid or gas) can make a large difference in the wave speeds and attenuation properties of seismic waves. The present work summarizes methods of deconstructing theses effects of fractures, together with any fluids contained within them, on wave propagation as observed in reflection seismic data. Additional studies of waves in fluid-saturated granular media show that the behavior can be quite different from that for fractured media, since these materials are typically much softer mechanically than are the fractured rocks (i.e., having a very small drained moduli). Important fluid effects in such media are often governed as much by fluid viscosity as by fluid bulk modulus.
Liquefaction of sedimentary rocks during impact crater development
NASA Astrophysics Data System (ADS)
Hippertt, J. P.; Lana, C.; Weinberg, R. F.; Tohver, E.; Schmieder, M.; Scholz, R.; Gonçalves, L.; Hippertt, J. F.
2014-12-01
Impact crater development on every planetary body requires catastrophic movement of large volumes of crustal rocks. The process produces well-known features such as brecciation and frictional melting, but a mechanism that explains how rocks accommodate the strain during the cratering flow remains unclear. Here, we investigate target rocks from the Araguainha impact crater (central Brazil) that typify what happens to a consolidated, fluid-saturated sedimentary rock at ˜ 2 km below the surface prior to the impact event. Sandstone units record a pattern of chaotic large-scale folds and pervasive microscopic (grain-to-grain) brecciation that result from rock strength degradation triggered by the impact. Field mapping and extensive textural observations indicate that these sandstones experienced initial microstructural damage from the shock wave and that this process may have weakened grain-to-grain bonds and started the process of pervasive microbrecciation. Accompanying heating and decompression lead to vaporization and expansion of fluids in the sandstone pores, magnifying the process of brecciation by effectively liquefying the rock mass and allowing for chaotic folding (at a range of scales up to blocks 100 m in length) in the central uplift. This is a vaporization-assisted microbrecciation, and it may have inhibited the formation of pseudotachylites, because energy was dissipated by pervasive microcracking, vaporization of pore fluids, and large scale chaotic folding, rather than localized displacement on brittle faults and frictional heating. We suggest that impact liquefaction of sedimentary rocks depends on whether the presence of pore-fluids and related micro-brecciation are sufficient to dissipate most of the impact energy.
ERIC Educational Resources Information Center
Bickett, Marianne
2011-01-01
This article discusses rock art which was the very first "art." Rock art, such as the images created on the stone surfaces of the caves of Lascaux and Altimira, is the true origin of the canvas, paintbrush, and painting media. For there, within caverns deep in the earth, the first artists mixed animal fat, urine, and saliva with powdered minerals…
NASA Technical Reports Server (NTRS)
2006-01-01
27 January 2006 This Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) image shows some of the light-toned, layered, sedimentary rock outcrops in northern Terby Crater. Terby is located along the north edge of Hellas Planitia. The sedimentary rocks might have been deposited in a greater, Hellas-filling sea -- or not. Today, the rocks are partly covered by dark-toned sediment and debris.
Location near: 27.2oS, 285.3oW Image width: 3 km (1.9 mi) Illumination from: upper left Season: Southern Summer
Freeze fracturing of elastic porous media: a mathematical model
Vlahou, I.; Worster, M. G.
2015-01-01
We present a mathematical model of the fracturing of water-saturated rocks and other porous materials in cold climates. Ice growing inside porous rocks causes large pressures to develop that can significantly damage the rock. We study the growth of ice inside a penny-shaped cavity in a water-saturated porous rock and the consequent fracturing of the medium. Premelting of the ice against the rock, which results in thin films of unfrozen water forming between the ice and the rock, is one of the dominant processes of rock fracturing. We find that the fracture toughness of the rock, the size of pre-existing faults and the undercooling of the environment are the main parameters determining the susceptibility of a medium to fracturing. We also explore the dependence of the growth rates on the permeability and elasticity of the medium. Thin and fast-fracturing cracks are found for many types of rocks. We consider how the growth rate can be limited by the existence of pore ice, which decreases the permeability of a medium, and propose an expression for the effective ‘frozen’ permeability. PMID:25792954
NASA Technical Reports Server (NTRS)
1997-01-01
This false color composite image of the Rock Garden shows the rocks 'Shark' and 'Half Dome' at upper left and middle, respectively. Between these two large rocks is a smaller rock (about 0.20 m wide, 0.10 m high, and 6.33 m from the Lander) that was observed close-up with the Sojourner rover (see PIA00989).
Mars Pathfinder is the second in NASA's Discovery program of low-cost spacecraft with highly focused science goals. The Jet Propulsion Laboratory, Pasadena, CA, developed and manages the Mars Pathfinder mission for NASA's Office of Space Science, Washington, D.C. JPL is a division of the California Institute of Technology (Caltech). The Imager for Mars Pathfinder (IMP) was developed by the University of Arizona Lunar and Planetary Laboratory under contract to JPL. Peter Smith is the Principal Investigator.
Method and apparatus for determining two-phase flow in rock fracture
Persoff, Peter; Pruess, Karsten; Myer, Larry
1994-01-01
An improved method and apparatus as disclosed for measuring the permeability of multiple phases through a rock fracture. The improvement in the method comprises delivering the respective phases through manifolds to uniformly deliver and collect the respective phases to and from opposite edges of the rock fracture in a distributed manner across the edge of the fracture. The improved apparatus comprises first and second manifolds comprising bores extending within porous blocks parallel to the rock fracture for distributing and collecting the wetting phase to and from surfaces of the porous blocks, which respectively face the opposite edges of the rock fracture. The improved apparatus further comprises other manifolds in the form of plenums located adjacent the respective porous blocks for uniform delivery of the non-wetting phase to parallel grooves disposed on the respective surfaces of the porous blocks facing the opposite edges of the rock fracture and generally perpendicular to the rock fracture.
Instability of fluid flow over saturated porous medium
NASA Astrophysics Data System (ADS)
Lyubimova, Tatyana; Kolchanova, Ekaterina; Lyubimov, Dmitry
2013-04-01
We investigate the stability of a fluid flow over a saturated porous medium. The problem is of importance due to the applications to washing out of contaminants from the bottom layer of vegetation, whose properties are similar to the properties of porous medium. In the case of porous medium with the relatively high permeability and porosity the flow involves a part of the fluid saturating the porous medium, with the tangential fluid velocity drop occurring because of the resistance of the solid matrix. The drop leads to the instability analogous to Kelvin-Helmholtz one accompanied by the formation of travelling waves. In the present paper we consider a two-layer system consisting of a pure fluid layer and a porous layer saturated by the fluid located underneath. The system is bounded by a rigid surface at the bottom and a non-deformable free surface at the top. It is under the gravity and inclined at a slight angle to the horizontal axis. The boundary conditions at the interface between the fluid and porous layers are the continuity of fluid velocities and the balance of normal and tangential stresses taking into account the resistance of the solid matrix with respect to the fluid flow near the interface [1-2]. The problem is solved in the framework of the Brinkman model applying the classical shooting algorithm with orthogonalization. The stability boundaries of the stationary fluid flow over the saturated porous medium with respect to the small oscillatory perturbations are obtained for the various values of the Darcy number and the ratio of the porous layer thickness to the full thickness of the system d. It was shown that at the d > 0.5 with increasing the porous layer thickness (or with decreasing of the fluid layer thickness) the stability threshold rises. This is because of the fact that the instability is primarily caused by perturbations located in the fluid layer. At the d < 0.5 the reduction of the porous layer thickness leads to the stability threshold
NASA Astrophysics Data System (ADS)
Mahood, Gail A.
Most of the advances in volcanology during the past 20 years have concerned the recognition, interpretation, and mode of emplacement of pyroclastic rocks. The literature on pyroclastic rocks is widely scattered, in part because the field draws from sedimentology, igneous petrology, physics, and fluid mechanics, and there have been few review papers on the topic. Fisher and Schmincke have done the discipline of volcanology and all field-oriented geologists a great service in assembling material from a wide range of sources in this comprehensive treatment of pyroclastic rocks. With its introduction to the petrology of magmas involved in explosive eruptions in chapter 2 and a complete treatment of magma rheology and the behavior of dissolved and exsolving magmatic volatiles in chapter 3, they lay sufficient groundwork that anyone with a rudimentary knowledge of geology can understand the book.
Evidence of Ancient Blisters in Rocks
NASA Technical Reports Server (NTRS)
2004-01-01
This image from the panoramic camera on the Mars Exploration Rover Spirit shows scoriaceous rocks (rocks containing holes or cavities) on the ground, as well as a transition from rocky terrain (foreground) to smoother terrain (background). Spirit is heading toward the smoother terrain on its way to the 'Columbia Hills.' The holes in some of the rocks may have resulted from 'blisters' formed by water vapor as it escaped lava. This indicates that the rocks were chilled atop an ancient lava flow. Porous rocks such as these, now appearing in abundance, have not been seen since early in the mission. Scientists believe they may have been covered by crater ejecta. This image was taken on sol 110 (April 24, 2004) at a region dubbed 'site 35.'
ERIC Educational Resources Information Center
Beem, Edgar Allen
2004-01-01
While "early college" programs designed for high-school-age students are beginning to proliferate nationwide, a small New England school has been successfully educating teens for nearly four decades. In this article, the author features Simon's Rock, a small liberal arts college located in the Great Barrington, Massachusetts, that has been…
NASA Astrophysics Data System (ADS)
Postelnicu, Adrian
2010-10-01
Dufour and Soret effects on flow at a stagnation point in a fluid-saturated porous medium are studied in this paper. A two dimensional stagnation-point flow with suction/injection of a Darcian fluid is considered. By using an appropriate similarity transformation, the boundary layer equations of momentum, energy and concentration are reduced to a set of ordinary differential equations, which are solved numerically using the Keller-box method, which is a very efficient finite differences technique. Nusselt and Sherwood numbers are obtained, together with the velocity, temperature and concentration profiles in the boundary layer. For the large suction case, asymptotic analytical solutions of the problem are obtained, which compare favourably with the numerical solutions. A critical view of the problem is presented finally.
Theory of wave propagation in partially saturated double-porosity rocks: a triple-layer patchy model
NASA Astrophysics Data System (ADS)
Sun, Weitao; Ba, Jing; Carcione, José M.
2016-04-01
Wave-induced local fluid flow is known as a key mechanism to explain the intrinsic wave dissipation in fluid-saturated rocks. Understanding the relationship between the acoustic properties of rocks and fluid patch distributions is important to interpret the observed seismic wave phenomena. A triple-layer patchy (TLP) model is proposed to describe the P-wave dissipation process in a double-porosity media saturated with two immiscible fluids. The double-porosity rock consists of a solid matrix with unique host porosity and inclusions which contain the second type of pores. Two immiscible fluids are considered in concentric spherical patches, where the inner pocket and the outer sphere are saturated with different fluids. The kinetic and dissipation energy functions of local fluid flow (LFF) in the inner pocket are formulated through oscillations in spherical coordinates. The wave propagation equations of the TLP model are based on Biot's theory and the corresponding Lagrangian equations. The P-wave dispersion and attenuation caused by the Biot friction mechanism and the local fluid flow (related to the pore structure and the fluid distribution) are obtained by a plane-wave analysis from the Christoffel equations. Numerical examples and laboratory measurements indicate that P-wave dispersion and attenuation are significantly influenced by the spatial distributions of both, the solid heterogeneity and the fluid saturation distribution. The TLP model is in reasonably good agreement with White's and Johnson's models. However, differences in phase velocity suggest that the heterogeneities associated with double-porosity and dual-fluid distribution should be taken into account when describing the P-wave dispersion and attenuation in partially saturated rocks.
Is biodegradable waste a porous environment? A review.
Agostini, Francesco; Sundberg, Cecilia; Navia, Rodrigo
2012-10-01
This article presents a review of the porous physical characteristics, phenomena and simulation models so far investigated and applied in the management of biodegradable wastes (BW), summarising the main properties of porous media and the dynamics of fluids within its voids. The aim is to highlight how the description of biodegradable wastes as porous media and the use of porous media models can facilitate the development of new sustainable and affordable technologies for BW recycling. However, it is pointed out how the lack of physical experimental data and of tailored modelling tools has so far hampered the use of this approach. Therefore, it is suggested that a simpler way to design and implement modelling tools simulating BW treatment technologies is by modifying available models designed originally for other porous media, such as soil and rock. PMID:22782804
NASA Technical Reports Server (NTRS)
1997-01-01
This image, taken by Sojourner's front right camera, was taken when the rover was next to Poohbear (rock at left) and Piglet (not seen) as it looked out toward Mermaid Dune. The textures differ from the foreground soil containing a sorted mix of small rocks, fines and clods, from the area a bit ahead of the rover where the surface is covered with a bright drift material. Soil experiments where the rover wheels dug in the soil revealed that the cloudy material exists underneath the drift.
Mars Pathfinder is the second in NASA's Discovery program of low-cost spacecraft with highly focused science goals. The Jet Propulsion Laboratory, Pasadena, CA, developed and manages the Mars Pathfinder mission for NASA's Office of Space Science, Washington, D.C. JPL is a division of the California Institute of Technology (Caltech).
NASA Technical Reports Server (NTRS)
2005-01-01
16 September 2005 This Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) image shows the complex surfaces of some of the light- and intermediate-toned sedimentary rock exposed by erosion in eastern Sinus Meridiani. Similar rocks occur at the Mars Exploration Rover, Opportunity, site, but they are largely covered by windblown sand and granules. The dark feature with a rayed pattern is the product of a meteor impact.
Location near: 0.8oN, 355.2oW Image width: width: 3 km (1.9 mi) Illumination from: lower left Season: Northern Autumn
NASA Technical Reports Server (NTRS)
2005-01-01
25 August 2005 This Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) image shows light-toned, layered, sedimentary rock outcrops in the crater, Terby. The crater is located on the north rim of Hellas Basin. If one could visit the rocks in Terby, one might learn from them whether they formed in a body of water. It is possible, for example, that Terby was a bay in a larger, Hellas-wide sea.
Location near: 27.9oS, 285.7oW Image width: width: 3 km (1.9 mi) Illumination from: upper left Season: Southern Winter
NASA Technical Reports Server (NTRS)
2002-01-01
(Released 19 April 2002) The Science 'White Rock' is the unofficial name for this unusual landform which was first observed during the Mariner 9 mission in the early 1970's. As later analysis of additional data sets would show, White Rock is neither white nor dense rock. Its apparent brightness arises from the fact that the material surrounding it is so dark. Images from the Mars Global Surveyor MOC camera revealed dark sand dunes surrounding White Rock and on the floor of the troughs within it. Some of these dunes are just apparent in the THEMIS image. Although there was speculation that the material composing White Rock could be salts from an ancient dry lakebed, spectral data from the MGS TES instrument did not support this claim. Instead, the White Rock deposit may be the erosional remnant of a previously more continuous occurrence of air fall sediments, either volcanic ash or windblown dust. The THEMIS image offers new evidence for the idea that the original deposit covered a larger area. Approximately 10 kilometers to the southeast of the main deposit are some tiny knobs of similarly bright material preserved on the floor of a small crater. Given that the eolian erosion of the main White Rock deposit has produced isolated knobs at its edges, it is reasonable to suspect that the more distant outliers are the remnants of a once continuous deposit that stretched at least to this location. The fact that so little remains of the larger deposit suggests that the material is very easily eroded and simply blows away. The Story Fingers of hard, white rock seem to jut out like icy daggers across a moody Martian surface, but appearances can be deceiving. These bright, jagged features are neither white, nor icy, nor even hard and rocky! So what are they, and why are they so different from the surrounding terrain? Scientists know that you can't always trust what your eyes see alone. You have to use other kinds of science instruments to measure things that our eyes can
NASA Technical Reports Server (NTRS)
2005-01-01
14 November 2005 This Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) image shows a portion of the famous 'White Rock' feature in Pollack Crater in the Sinus Sabaeus region of Mars. The light-toned rock is not really white, but its light tone caught the eye of Mars geologists as far back as 1972, when it was first spotted in images acquired by Mariner 9. The light-toned materials are probably the remains of a suite of layered sediments that once spread completely across the interior of Pollack Crater. Dark materials in this image include sand dunes and large ripples.
Location near: 8.1oS, 335.1oW Image width: width: 3 km (1.9 mi) Illumination from: lower left Season: Southern Summer
Ultrasonic Nondestructive Characterization of Porous Materials
NASA Astrophysics Data System (ADS)
Yang, Ningli
2011-12-01
Wave propagation in porous media is studied in a wide range of technological applications. In the manufacturing industry, determining porosity of materials in the manufacturing process is required for strict quality control. In the oil industry, acoustic signals and seismic surveys are used broadly to determine the physical properties of the reservoir rock which is a porous media filled with oil or gas. In porous noise control materials, a precise prediction of sound absorption with frequency and evaluation of tortuosity are necessary. Ultrasonic nondestructive methods are a very important tool for characterization of porous materials. The dissertation deals with two types of porous media: materials with relatively low and closed porosity and materials with comparatively high and open porosity. Numerical modeling, Finite Element simulations and experimental characterization are all discussed in this dissertation. First, ultrasonic scattering is used to determine the porosity in porous media with closed pores. In order get a relationship between the porosity in porous materials and ultrasonic scattering independently and to increase the sensitivity to obtain scattering information, ultrasonic imaging methods are applied and acoustic waves are focused by an acoustic lens. To verify the technique, engineered porous acrylic plates with varying porosity are measured by ultrasonic scanning and ultrasonic array sensors. Secondly, a laser based ultrasonic technique is explored for predicting the mechanical integrity and durability of cementitious materials. The technique used involves the measurement of the phase velocity of fast and slow longitudinal waves in water saturated cement paste. The slow wave velocity is related to the specimen's tortuosity. The fast wave speed is dependent on the elastic properties of porous solid. Experimental results detailing the generation and detection of fast and slow wave waves in freshly prepared and aged water-saturated cement samples
The relation between seismic P- and S-wave velocity dispersion in saturated rocks
Mavko, G.; Jizba, D.
1994-01-01
Seismic velocity dispersion in fluid-saturated rocks appears to be dominated by two mechanisms: the large scale mechanism modeled by Biot, and the local flow or squirt mechanism. The two mechanisms can be distinguished by the ratio of P- to S-wave dispersions, or more conveniently, by the ratio of dynamic bulk to shear compliance dispersions derived from the wave velocities. The authors` formulation suggests that when local flow dominates, the dispersion of the shear compliance will be approximately 4/15 the dispersion of the compressibility. When the Biot mechanism dominates, the constant of proportionality is much smaller. Their examination of ultrasonic velocities from 40 sandstones and granites shows that most, but not all, of the samples were dominated by local flow dispersion, particularly at effective pressures below 40 MPa.
Hearon, Keith; Singhal, Pooja; Horn, John; Small, Ward; Olsovsky, Cory; Maitland, Kristen C.; Wilson, Thomas S.; Maitland, Duncan J.
2013-01-01
Porous shape memory polymers (SMPs) include foams, scaffolds, meshes, and other polymeric substrates that possess porous three-dimensional macrostructures. Porous SMPs exhibit active structural and volumetric transformations and have driven investigations in fields ranging from biomedical engineering to aerospace engineering to the clothing industry. The present review article examines recent developments in porous SMPs, with focus given to structural and chemical classification, methods of characterization, and applications. We conclude that the current body of literature presents porous SMPs as highly interesting smart materials with potential for industrial use. PMID:23646038
BARTON,THOMAS J.; BULL,LUCY M.; KLEMPERER,WALTER G.; LOY,DOUGLAS A.; MCENANEY,BRIAN; MISONO,MAKOTO; MONSON,PETER A.; PEZ,GUIDO; SCHERER,GEORGE W.; VARTULI,JAMES C.; YAGHI,OMAR M.
1999-11-09
Tailoring of porous materials involves not only chemical synthetic techniques for tailoring microscopic properties such as pore size, pore shape, pore connectivity, and pore surface reactivity, but also materials processing techniques for tailoring the meso- and the macroscopic properties of bulk materials in the form of fibers, thin films and monoliths. These issues are addressed in the context of five specific classes of porous materials: oxide molecular sieves, porous coordination solids, porous carbons, sol-gel derived oxides, and porous heteropolyanion salts. Reviews of these specific areas are preceded by a presentation of background material and review of current theoretical approaches to adsorption phenomena. A concluding section outlines current research needs and opportunities.
Berryman, J.G.; Bonner, B.P.; Chin, R.C.Y.
1983-07-01
The measured amplitude A of ultrasonic pulses in intact and fractured samples of water-saturated gabbro and granite is observed to decrease as the permeability kappa increases according to the proportionality Aproportionalkappa/sup -1/2/. This relation is predicted by Biot's theory of elastic waves in fluid-saturated porous media and, therefore, suggests that Biot's attenuation mechanism may play a significant role in low porosity materials at ultrasonic frequencies. The evidence is not conclusive. The limited data set studied here is also consistent with correlations of the form Aproportionalkappa/sup -Epsilon/ where 0.2
A new understanding of fluid-rock deformation
NASA Astrophysics Data System (ADS)
Crampin, Stuart; Gao, Yuan
2015-04-01
Cracks in the pavement show that rock is weak to shear stress. Consequently we have a conundrum. How does in situ rock accumulate the enormous shear-stress energy necessary for release by a large magnitude earthquake without fracturing in smaller earthquakes? For example: observations of changes in seismic shear-wave splitting (SWS) were observed in Iceland before the 2004 Mw9.2 Sumatra-Andaman Earthquake (SAE) at a distance of ~10,500km (the width of the Eurasian Plate) from Indonesia. Observations of SWS monitor microcrack geometry, and the changes in SWS in Iceland indicated that stress-changes before the Sumatra earthquake modified microcrack geometry the width of Eurasia from Indonesia. What is the mechanism for such widespread accumulation of necessarily weak stress? We show that stress is stored in in situ rock by the stress-controlled geometry of the fluid-saturated stress-aligned microcrack. Microcrack aspect-ratios are aligned by fluid flow or dispersion along pressure-gradients between neighbouring microcracks at different orientations to the stress-field by a mechanism known as Anisotropic Poro-Elasticity or APE. Since the minimum stress is typically horizontal, the microcracks are typically vertically-oriented parallel to the maximum horizontal stress as is confirmed by observations of SWS. Such azimuthally varying shear-wave splitting (SWS) is observed in situ rocks in the upper crust, lower crust, and uppermost ~400km of the mantle. (The 'microcracks' in the mantle are intergranular films of hydrolysed melt.) SWS shows that the microcracks are so closely spaced that they verge on fracturing/earthquakes. Phenomena verging on failure are critical-systems with 'butterfly wings' sensitivity. Critical-systems are very common and it must be expected that the Earth, an archetypal complex heterogeneous interactive phenomena is a critical-system. Monitoring SWS above small earthquakes allows stress-accumulation before earthquakes to be recognised and the time
NASA Technical Reports Server (NTRS)
2005-01-01
6 November 2005 This Mars Global Surveyor (MGS) Mars Orbiter Camera (MOC) image shows outcrops of sedimentary rocks in a crater located just north of the Sinus Meridiani region. Perhaps the crater was once the site of a martian lake.
Location near: 2.9oN, 359.0oW Image width: width: 3 km (1.9 mi) Illumination from: lower left Season: Northern Autumn
Thermal Inertia of Rocks and Rock Populations
NASA Technical Reports Server (NTRS)
Golombek, M. P.; Jakosky, B. M.; Mellon, M. T.
2001-01-01
The effective thermal inertia of rock populations on Mars and Earth is derived from a model of effective inertia versus rock diameter. Results allow a parameterization of the effective rock inertia versus rock abundance and bulk and fine component inertia. Additional information is contained in the original extended abstract.
NASA Technical Reports Server (NTRS)
Peterson, Thomas M.
2001-01-01
The next series of planetary exploration missions require a method of extracting rock and soil core samples. Therefore a prototype ultrasonic core driller (UTCD) was developed to meet the constraints of Small Bodies Exploration and Mars Sample Return Missions. The constraints in the design are size, weight, power, and axial loading. The ultrasonic transducer requires a relatively low axial load, which is one of the reasons this technology was chosen. The ultrasonic generator breadboard section can be contained within the 5x5x3 limits and weighs less than two pounds. Based on results attained the objectives for the first phase were achieved. A number of transducer probes were made and tested. One version only drills, and the other will actually provide a small core from a rock. Because of a more efficient transducer/probe, it will run at very low power (less than 5 Watts) and still drill/core. The prototype generator was built to allow for variation of all the performance-effecting elements of the transducer/probe/end effector, i.e., pulse, duty cycle, frequency, etc. The heart of the circuitry is what will be converted to a surface mounted board for the next phase, after all the parameters have been optimized and the microprocessor feedback can be installed.
Displacement propagators of brine flowing within different types of sedimentary rock.
Verganelakis, Dimitris A; Crawshaw, John; Johns, Michael L; Mantle, Michael D; Scheven, Ulrich; Sederman, Andrew J; Gladden, Lynn F
2005-02-01
This paper explores the correlation between different microstructural characteristics of porous sedimentary rocks and the flow properties of a Newtonian infiltrating fluid. Preliminary results of displacement propagator measurements of brine solution flowing through two types of sedimentary rock cores are reported. The two types of rocks, Bentheimer and Portland, are characterized by different porosities, pore-size distributions and permeabilities. Propagators have been measured for brine flow rates of 1 and 5 ml/min. Significant differences are seen between the propagators recorded for the two rocks, and these are related to the spatial distribution of porosity within these porous media. PMID:15833644
Correlation between roughness and porosity in rocks
NASA Astrophysics Data System (ADS)
Rebollo, M. A.; Hogert, E. N.; Albano, J.; Raffo, C. A.; Gaggioli, N. G.
1996-02-01
The porosity of rocks is a very important parameter in the determination of the performance of oil wells. Optical methods allow us to study surface roughness and different materials that have surface properties with random characteristics. Therefore, we have extended these applications to porosity analysis. In our method, we have used the speckle produced by the scattered light from a porous rock, illuminated by a laser beam, and found a linear relationship between the decorrelation of the speckle intensity distribution and the porosity magnitude. In this paper we present the results for samples extracted from oil wells in Argentina.
Guan, Wei; Hu, Hengshan; He, Xiao
2009-04-01
Monopole acoustic logs in a homogeneous fluid-saturated porous formation can be simulated by the real-axis integration (RAI) method to analytically solve Biot's equations [(1956a) J. Acoust. Soc. Am. 28, 168-178; (1956b) J. Acoust. Soc. Am. 28, 179-191; (1962) J. Appl. Phys. 33, 1482-1498], which govern the wave propagation in poro-elastic media. Such analytical solution generally is impossible for horizontally stratified formations which are common in reality. In this paper, a velocity-stress finite-difference time-domain (FDTD) algorithm is proposed to solve the problem. This algorithm considers both the low-frequency viscous force and the high-frequency inertial force in poro-elastic media, extending its application to a wider frequency range compared to existing algorithms which are only valid in the low-frequency limit. The perfectly matched layer (PML) is applied as an absorbing boundary condition to truncate the computational region. A PML technique without splitting the fields is extended to the poro-elastic wave problem. The FDTD algorithm is validated by comparisons against the RAI method in a variety of formations with different velocities and permeabilities. The acoustic logs in a horizontally stratified porous formation are simulated with the proposed FDTD algorithm. PMID:19354370
Berryman, J.G.
2009-11-20
Poroelastic analysis usually progresses from assumed knowledge of dry or drained porous media to the predicted behavior of fluid-saturated and undrained porous media. Unfortunately, the experimental situation is often incompatible with these assumptions, especially when field data (from hydrological or oil/gas reservoirs) are involved. The present work considers several different experimental scenarios typified by one in which a set of undrained poroelastic (stiffness) constants has been measured using either ultrasound or seismic wave analysis, while some or all of the dry or drained constants are normally unknown. Drained constants for such a poroelastic system can be deduced for isotropic systems from available data if a complete set of undrained compliance data for the principal stresses are available - together with a few other commonly measured quantities such as porosity, fluid bulk modulus, and grain bulk modulus. Similar results are also developed here for anisotropic systems having up to orthotropic symmetry if the system is granular (i.e., composed of solid grains assembled into a solid matrix, either by a cementation process or by applied stress) and the grains are known to be elastically homogeneous. Finally, the analysis is also fully developed for anisotropic systems with nonhomogeneous (more than one mineral type), but still isotropic, grains - as well as for uniform collections of anisotropic grains as long as their axes of symmetry are either perfectly aligned or perfectly random.
Natural convection heat transfer from a horizontal wavy surface in a porous enclosure
Murthy, P.V.S.N.; Kumar, B.V.R.; Singh, P.
1997-02-07
The effect of surface undulations on the natural convection heat transfer from an isothermal surface in a Darcian fluid-saturated porous enclosure has been numerically analyzed using the finite element method on a graded nonuniform mesh system. The flow-driving Rayleigh number Ra together with the geometrical parameters of wave amplitude a, wave phase {phi}, and the number of waves N considered in the horizontal dimension of the cavity are found to influence the flow and heat transfer process in the enclosure. For Ra around 50 and above, the phenomenon of flow separation and reattachment is noticed on the walls of the enclosure. A periodic shift in the reattachment point from the bottom wall to the adjacent walls in the clockwise direction, leading to the manifestation of cycles of unicellular and bicellular clockwise and counterclockwise flows, is observed, with the phase varying between 0{degree} and 350{degree}. The counterflow in the secondary circulation zone is intensified with the increase in the value of Ra. The counterflow on the wavy wall hinders the heat transfer into the system. An increase in either wave amplitude or the number of waves considered per unit length decreases the global heat flux into the system. Only marginal changes in global heat flux are noticed with increasing Ra. On the whole, the comparison of global heat flux results in the wavy wall case with those of the horizontal flat wall case shows that, in a porous enclosure, the wavy wall reduces the heat transfer into the system.
Bacteria transport through porous material: Final technical report
Yen, T.F.
1989-02-13
The injection and penetration of bacteria into a reservoir is the most problematic and crucial of the steps in microbial enhanced recovery (MEOR). In the last phase of our work valuable information on bacterial transport in porous media was obtained. A great deal of progress was made to determine chemical bonding characteristics between adsorbed bacteria and the rock surfaces. In order to further enhance our knowledge of the effects of surface tensions on bacteria transport through porous media, a new approach was taken to illustrate the effect of liquid surface tension on bacterial transport through a sandpack column. Work in surface charge characterization of reservoir rock as a composite oxide system was also accomplished. In the last section of this report a mathematical model to simulate the simultaneous diffusion and growth of bacteria cells in a nutrient-enriched porous media is proposed.
Note on the validity of the Kozeny-Carman formulas for consolidated porous media
Schlueter, E.M.; Witherspoon, P.A.
1995-12-31
The region of validity of the well-known Kozeny-Carmen permeability formulas for consolidated porous media and their relationship to the microscopic spatial variations of channel dimensions are established. For highly inhomogeneous rock-pore-space systems, the validity of the critical path analysis for estimating the permeability of consolidated rocks is examined also.
Rock Degradation by Alkali Metals: A Possible Lunar Erosion Mechanism.
Naughton, J J; Barnes, I L; Hammond, D A
1965-08-01
When rocks melt under ultrahigh-vacuum conditions, their alkali components volatilize as metals. These metal vapors act to comminute polycrystalline rocks to their component minerals. The resultant powder is porous and loosely packed and its characteristics may be compatible with the lunar surface as revealed by the Ranger photographs. If meteorite impact or lunar volcanism has produced vaporization or areas of molten lava, alkali erosion may have given dust of this character in adjacent solid areas. PMID:17747570
Neeper, D.A.
1994-02-22
Methods are presented for distributing gases throughout the interstices of porous materials and removing volatile substances from the interstices of porous materials. Continuous oscillation of pressures and flows results in increased penetration of the interstices by flowing gases and increased transport of gaseous components out of the interstices. The invention is particularly useful in soil vapor extraction. 10 figures.
Neeper, Donald A.
1994-01-01
Methods for distributing gases throughout the interstices of porous materials and removing volatile substances from the interstices of porous materials. Continuous oscillation of pressures and flows results in increased penetration of the interstices by flowing gases and increased transport of gaseous components out of the interstices. The invention is particularly useful in soil vapor extraction.
X-ray scattering by porous silicon modulated structures
Lomov, A. A.; Punegov, V. I.; Karavanskii, V. A.; Vasil'ev, A. L.
2012-03-15
A multilayered porous structure formed as a result of the anodization of a Si(111)(Sb) substrate in an HF:C{sub 2}H{sub 5}OH (1: 2) solution with a periodically alternating current has been investigated by high-resolution X-ray diffraction. It is established that, despite 50% porosity, a thickness of 30 {mu}m, and significant strain (4 Multiplication-Sign 10{sup -3)}, the porous silicon structure consists mainly of coherent crystallites. A model of coherent scattering from a multilayered periodic porous structure is proposed within the dynamic theory of diffraction. It is shown that the presence of gradient strains of 5 Multiplication-Sign 10{sup -4} or higher leads to phase loss upon scattering from porous superlattices and the suppression of characteristic satellites in rocking curves.
Tsuo, Y.S.; Menna, P.; Al-Jassim, M.
1995-08-01
We have studied a novel extrinsic gettering method that utilizes the very large surface areas, produced by porous silicon etch on both front and back surfaces of the silicon wafer, as gettering sites. In this method, a simple and low-cost chemical etching is used to generate the porous silicon layers. Then, a high-flux solar furnace (HFSF) is used to provide high-temperature annealing and the required injection of silicon interstitials. The gettering sites, along with the gettered impurities, can be easily removed at the end the process. The porous silicon removal process consists of oxidizing the porous silicon near the end the gettering process followed by sample immersion in HF acid. Each porous silicon gettering process removes up to about 10 {mu}m of wafer thickness. This gettering process can be repeated so that the desired purity level is obtained.
A Rock Encyclopedia That Includes Rock Samples.
ERIC Educational Resources Information Center
Laznicka, Peter
1981-01-01
Described is a rock encyclopedia combining rock sample sets and encyclopedic word and picture entries which can be used as a realistic information resource for independent study or as a part of a course. (JT)
Thermal conductivity of carbonate rocks
Thomas, J., Jr.; Frost, R.R.; Harvey, R.D.
1973-01-01
The thermal conductivities of several well-defined carbonate rocks were determined near 40??C. Values range from 1.2 W m-1 C-1 for a highly porous chalk to 5.1 W m-1 C-1 for a dolomite. The thermal conductivity of magnesite (5.0) is at the high end of the range, and that for Iceland Spar Calcite (3.2) is near the middle. The values for limestones decrease linearly with increasing porosity. Dolomites of comparable porosity have greater thermal conductivities than limestones. Water-sorbed samples have expected greater thermal conductivities than air-saturated (dry) samples of the same rock. An anomalously large increase in the thermal conductivity of a water-sorbed clayey dolomite over that of the same sample when dry is attributed to the clay fraction, which swells during water inhibition, causing more solid-to-solid contacts within the dolomite framework. Measurements were made with a Colora Thermoconductometer. Chemical and mineralogical analyses were made and tabulated. Porosity of the rocks was determined by mercury porosimetry and also from density measurements. The Iceland Spar Calcite and magnesite were included for reference. ?? 1973.
High=porosity Cenozoic carbonate rocks of South Florida: progressive loss of porosity with depth
Halley, Robert B.; Schmoker, James W.
1983-01-01
Porosity measurements by borehole gravity meter in subsurface Cenozoic carbonates of South Florida reveal an extremely porous mass of limestone and dolomite which is transitional in total pore volume between typical porosity values for modern carbonate sediments and ancient carbonate rocks. A persistent decrease of porosity with depth, similar to that of chalks of the Gulf Coast, occurs in these rocks. Carbonate strata with less than 20% porosity are absent from the rocks studied here. Aquifers and aquicludes cannot be distinguished on the basis of porosity. Aquifers are not exceptionally porous when compared to other Tertiary carbonate rocks in South Florida. Permeability in these strata is governed more by the spacial distribution of pore space and matrix than by total volume of porosity present. Dolomite is as porous as, or slightly less porous than, limestones in these rocks. This observation places limits on any model proposed for dolomitization and suggests that dolomitization does not take place by a simple ion-for-ion replacement of magnesium for calcium. Dolomitization may be selective for less porous limestone, or it may involve the incorporation of significant amounts of carbonate as well as magnesium into the rock. The great volume of pore space in these rocks serves to highlight the inefficiency of early diagenesis in reducing carbonate porosity and to emphasize the importance of later porosity reduction which occurs during the burial or late near-surface history of limestones and dolomites.
Uncertainty quantification for porous media flows
NASA Astrophysics Data System (ADS)
Christie, Mike; Demyanov, Vasily; Erbas, Demet
2006-09-01
Uncertainty quantification is an increasingly important aspect of many areas of computational science, where the challenge is to make reliable predictions about the performance of complex physical systems in the absence of complete or reliable data. Predicting flows of oil and water through oil reservoirs is an example of a complex system where accuracy in prediction is needed primarily for financial reasons. Simulation of fluid flow in oil reservoirs is usually carried out using large commercially written finite difference simulators solving conservation equations describing the multi-phase flow through the porous reservoir rocks. This paper examines a Bayesian Framework for uncertainty quantification in porous media flows that uses a stochastic sampling algorithm to generate models that match observed data. Machine learning algorithms are used to speed up the identification of regions in parameter space where good matches to observed data can be found.
ELECTROKINETIC WAVE PHENOMENA IN FLUID-SATURATED GRANULAR MEDIA
Block, G
2005-03-29
Electrokinetic (EK) phenomena in sediments arise from relative fluid motion in the pore space, which perturbs the electrostatic equilibrium of the double layer at the grain surface. We have developed EK techniques in the laboratory to monitor acoustic wave propagation in electrolyte-saturated, unconsolidated sediments. Our experimental results indicate that as an acoustic wave travels through electrolyte-saturated sand, it can generate electric potentials greater than 1 mV. A careful study of these potentials was performed using medium-grain sand and loose glass microspheres for a range of pore fluid salinities and ultrasonic frequencies. Experimental results are also shown to compare well with numerical and analytical modeling based on the coupled electrokinetic-Biot theory developed by Pride (1994).
Discrimination of porosity and fluid saturation using seismic velocity analysis
Berryman, James G.
2001-01-01
The method of the invention is employed for determining the state of saturation in a subterranean formation using only seismic velocity measurements (e.g., shear and compressional wave velocity data). Seismic velocity data collected from a region of the formation of like solid material properties can provide relatively accurate partial saturation data derived from a well-defined triangle plotted in a (.rho./.mu., .lambda./.mu.)-plane. When the seismic velocity data are collected over a large region of a formation having both like and unlike materials, the method first distinguishes the like materials by initially plotting the seismic velocity data in a (.rho./.lambda., .mu./.lambda.)-plane to determine regions of the formation having like solid material properties and porosity.
Solidus of carbonated fertile peridotite under fluid-saturated conditions
Falloon, T.J.; Green, D.H. )
1990-03-01
The solidus for a fertile peridotite composition (Hawaiian pyrolite) in the presence of a CO{sub 2}-H{sub 2}O fluid phase has been determined from 10 to 35 kbar. The intersection of the decarbonation reaction (olivine + diopside + CO{sub 2} {l reversible} orthopyroxene + dolomite) with the pyrolite solidus defines the point Q{prime}, located at 22 kbar and 940 C. At pressures less than Q{prime}, the solidus passes through a temperature maximum at kbar, 1060 C. The solidus is coincident with amphibole breakdown at pressures less than 16 kbar. At pressures above Q{prime}, the solidus is defined by the dissolution of crystalline carbonate into a sodic, dolomitic carbonatite melt. The solidus is at a temperature of 925 C at {approximately} 28 kbar. The solidus temperature above the point Q{prime} is similar to the solidus determined for Hawaiian pyrolite-H{sub 2}O-CO{sub 2} for small contents of H{sub 2}O (<0.3 wt%) and CO{sub 2} (<5 wt%), thus indicating that the primary sodic dolomitic carbonatite melt at both solidi has a very low and limited H{sub 2}O solubility. The new data clarify the roles of carbonatite melt, carbonated silicate melt, and H{sub 2}O-rich fluid in mantle conditions that are relatively oxidized (f{sub O{sub 2}} {approximately} MW to FMQ). In particular, a carbonatite melt + garnet lherzolite region is intersected by continental shield geothermal gradients, but such geotherms only intersect regions with carbonated silicate melt if perturbed to higher temperatures (kinked geotherm).
Characterization of rock for constraining reservoir scale tomography at the Geysers geothermal field
Boitnott, G.N.; Bonner, B.P.
1994-01-20
A suite of laboratory measurements are being conducted on Geysers graywacke recovered from a drilled depth of 2599 meters in NEGU-17. The tests are being conducted to characterize the effect of pressure and fluid saturation on the seismic properties of the graywacke matrix. The measurements indicate that the graywacke is an unusual rock in many respects. Both compressional and shear velocities exhibit relatively little change with pressure. Water saturation causes a slight increase in the compressional velocity, quantitatively consistent with predictions from the Biot-Gassmann equations. Shear velocity decreases with water saturation by an amount greater than that predicted by the Biot-Gassmann equations. This decrease is attributed to chemomechanical weakening caused by the presence of water. Measurements of Q, from torsion experiments on room dry samples at seismic frequencies indicate unusually high Q, (~500). Water saturation decreases the shear modulus by 12 percent, again indicative of chemomechanical weakening. Q, is lower for the water saturated condition, but still relatively high for rock at low stress. Results of ultrasonic pulse propagation experiments on partially saturated samples are typical of low porosity rocks, being characterized by a monotonic decrease in compressional and shear velocity with decrease in saturation. An increase in shear velocity and low frequency shear modulus after vacuum drying indicates the presence of chemo-mechanical weakening resulting from the presence of small amounts of water.
Impact of fluid-rock chemical interactions on tracer transport in fractured rocks.
Mukhopadhyay, Sumit; Liu, H-H; Spycher, N; Kennedy, B M
2013-11-01
In this paper, we investigate the impact of chemical interactions, in the form of mineral precipitation and dissolution reactions, on tracer transport in fractured rocks. When a tracer is introduced in fractured rocks, it moves through the fracture primarily by advection and it also enters the stagnant water of the surrounding rock matrix through diffusion. Inside the porous rock matrix, the tracer chemically interacts with the solid materials of the rock, where it can precipitate depending on the local equilibrium conditions. Alternatively, it can be dissolved from the solid phase of the rock matrix into the matrix pore water, diffuse into the flowing fluids of the fracture and is advected out of it. We show that such chemical interactions between the fluid and solid phases have significant impact on tracer transport in fractured rocks. We invoke the dual-porosity conceptualization to represent the fractured rocks and develop a semi-analytical solution to describe the transient transport of tracers in interacting fluid-rock systems. To test the accuracy and stability of the semi-analytical solution, we compare it with simulation results obtained with the TOUGHREACT simulator. We observe that, in a chemically interacting system, the tracer breakthrough curve exhibits a pseudo-steady state, where the tracer concentration remains more or less constant over a finite period of time. Such a pseudo-steady condition is not observed in a non-reactive fluid-rock system. We show that the duration of the pseudo-state depends on the physical and chemical parameters of the system, and can be exploited to extract information about the fractured rock system, such as the fracture spacing and fracture-matrix interface area. PMID:24077359
Impact of fluid-rock chemical interactions on tracer transport in fractured rocks
NASA Astrophysics Data System (ADS)
Mukhopadhyay, Sumit; Liu, H.-H.; Spycher, N.; Kennedy, B. M.
2013-11-01
In this paper, we investigate the impact of chemical interactions, in the form of mineral precipitation and dissolution reactions, on tracer transport in fractured rocks. When a tracer is introduced in fractured rocks, it moves through the fracture primarily by advection and it also enters the stagnant water of the surrounding rock matrix through diffusion. Inside the porous rock matrix, the tracer chemically interacts with the solid materials of the rock, where it can precipitate depending on the local equilibrium conditions. Alternatively, it can be dissolved from the solid phase of the rock matrix into the matrix pore water, diffuse into the flowing fluids of the fracture and is advected out of it. We show that such chemical interactions between the fluid and solid phases have significant impact on tracer transport in fractured rocks. We invoke the dual-porosity conceptualization to represent the fractured rocks and develop a semi-analytical solution to describe the transient transport of tracers in interacting fluid-rock systems. To test the accuracy and stability of the semi-analytical solution, we compare it with simulation results obtained with the TOUGHREACT simulator. We observe that, in a chemically interacting system, the tracer breakthrough curve exhibits a pseudo-steady state, where the tracer concentration remains more or less constant over a finite period of time. Such a pseudo-steady condition is not observed in a non-reactive fluid-rock system. We show that the duration of the pseudo-state depends on the physical and chemical parameters of the system, and can be exploited to extract information about the fractured rock system, such as the fracture spacing and fracture-matrix interface area.
Preparation of asymmetric porous materials
Coker, Eric N.
2012-08-07
A method for preparing an asymmetric porous material by depositing a porous material film on a flexible substrate, and applying an anisotropic stress to the porous media on the flexible substrate, where the anisotropic stress results from a stress such as an applied mechanical force, a thermal gradient, and an applied voltage, to form an asymmetric porous material.
Fabricating porous silicon carbide
NASA Technical Reports Server (NTRS)
Shor, Joseph S. (Inventor); Kurtz, Anthony D. (Inventor)
1994-01-01
The formation of porous SiC occurs under electrochemical anodization. A sample of SiC is contacted electrically with nickel and placed into an electrochemical cell which cell includes a counter electrode and a reference electrode. The sample is encapsulated so that only a bare semiconductor surface is exposed. The electrochemical cell is filled with an HF electrolyte which dissolves the SiC electrochemically. A potential is applied to the semiconductor and UV light illuminates the surface of the semiconductor. By controlling the light intensity, the potential and the doping level, a porous layer is formed in the semiconductor and thus one produces porous SiC.
Machado, A C; Teles, A P; Pepin, A; Bize-Forest, N; Lima, I; Lopes, R T
2016-04-01
Porous space characterization of carbonate rocks is an important aid in petroleum exploration from carbonate reservoir. In this study, X-ray microtomography technique was applied to evaluate total porosity of a coquina sample extracted from pre-salt reservoir, in Brazil, before and after acid injection. Two image processing program were used in order to assess performance. The results showed that microtomography has potential to compute porosity of coquina samples and provides information about rock porous network. PMID:26794261
Probing porous media with gas diffusion NMR.
Mair, R W; Wong, G P; Hoffmann, D; Hurlimann, M D; Patz, S; Schwartz, L M; Walsworth, R L
1999-10-18
We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks. PMID:11543587
Probing porous media with gas diffusion NMR
NASA Technical Reports Server (NTRS)
Mair, R. W.; Wong, G. P.; Hoffmann, D.; Hurlimann, M. D.; Patz, S.; Schwartz, L. M.; Walsworth, R. L.
1999-01-01
We show that gas diffusion nuclear magnetic resonance (GD-NMR) provides a powerful technique for probing the structure of porous media. In random packs of glass beads, using both laser-polarized and thermally polarized xenon gas, we find that GD-NMR can accurately measure the pore space surface-area-to-volume ratio, S/V rho, and the tortuosity, alpha (the latter quantity being directly related to the system's transport properties). We also show that GD-NMR provides a good measure of the tortuosity of sandstone and complex carbonate rocks.
Technology Transfer Automated Retrieval System (TEKTRAN)
Fractal and prefractal geometric models have substantial potential of contributing to the analysis of flow and transport in porous media such as soils and reservoir rocks. In this study, geometric and hydrodynamic parameters of saturated 3D mass and pore-solid prefractal porous media were characteri...
X-ray microtomography application in pore space reservoir rock.
Oliveira, M F S; Lima, I; Borghi, L; Lopes, R T
2012-07-01
Characterization of porosity in carbonate rocks is important in the oil and gas industry since a major hydrocarbons field is formed by this lithology and they have a complex media porous. In this context, this research presents a study of the pore space in limestones rocks by x-ray microtomography. Total porosity, type of porosity and pore size distribution were evaluated from 3D high resolution images. Results show that carbonate rocks has a complex pore space system with different pores types at the same facies. PMID:22264795
NASA Technical Reports Server (NTRS)
Peters, Gregory; Brown, Kyle; Fuerstenau, Stephen
2009-01-01
The rollerjaw rock crusher melds the concepts of jaw crushing and roll crushing long employed in the mining and rock-crushing industries. Rollerjaw rock crushers have been proposed for inclusion in geological exploration missions on Mars, where they would be used to pulverize rock samples into powders in the tens of micrometer particle size range required for analysis by scientific instruments.
Porous Organic Molecular Materials
Tian, Jian; Thallapally, Praveen K.; McGrail, B. Peter
2012-01-01
Most nanoporous materials with molecular-scale pores are extended frameworks composed of directional covalent or coordination bonding, such as porous metal-organic frameworks and organic network polymers. By contrast, nanoporous materials comprised of discrete organic molecules, between which there are only weak non-covalent interactions, are seldom encountered. Indeed, most organic molecules pack efficiently in the solid state to minimize the void volume, leading to non-porous materials. In recent years, a significant number of nanoporous organic molecular materials, which may be either crystalline or amorphous, have been confirmed by the studies of gas adsorption and they are surveyed in this Highlight. In addition, the possible advantages of porous organic molecular materials over porous networks are discussed.
Three-phase flow in porous media: A review of experimental studies on relative permeability
NASA Astrophysics Data System (ADS)
Alizadeh, A. H.; Piri, M.
2014-09-01
We present a detailed, synthesized review of experimental studies on three-phase relative permeability published since 1980. We provide comprehensive, yet highly focused, analysis of critical aspects of the field and their evolution over the last three decades. In particular, we review the effects of saturation history, wettability, spreading, and layer drainage on the measured flow properties. We also list all the processes, rock types, fluid systems, and measurement techniques in order to provide a clear map for future studies. Behavior of the measured three-phase relative permeabilities with respect to fluid saturations, saturation histories, wettability of rock samples, spreading characteristics, interfacial tensions, and other pertinent properties are carefully discussed. Studies that use a diverse set of experimental techniques and data analysis to deduce relative permeability are included. The experimental techniques that should be utilized to reduce uncertainty are also explored. We interpret the measured properties and outcomes of different studies and compare them to substantiate distinct trends at various saturation ranges and provide ideas for new studies. This is intended to distill a clear image of where the field stands and to allow composition of possible paths for future investigations. The areas of critical relevance that have not been investigated or require further studies are highlighted.