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

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

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

2015-11-30

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

Laminar film condensation along a vertical plate embedded in an anisotropic porous medium with oblique principal axes

NASA Astrophysics Data System (ADS)

Degan, Gérard; Sanya, Arthur; Akowanou, Christian

2016-10-01

This work analytically investigates the problem of steady film condensation along a vertical surface embedded in an anisotropic porous medium filled with a dry saturated vapor. The porous medium is anisotropic in permeability whose principal axes are oriented in a direction which is oblique to the gravity vector. On the basis of the generalized Darcy's law and within the boundary layer approximations, similar solutions have been obtained for the temperature and flow patterns in the condensate. Moreover, closed form solutions for the boundary layer thickness and heat transfer rate have been obtained in terms of the governing parameters of the problem.

Theoretical and numerical aspects of fluid-saturated elasto-plastic soils

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

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,more » 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.« less

Frequency-dependent processing and interpretation (FDPI) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs

DOEpatents

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.

Frequency-dependent processing and interpretation (FDPI) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs

DOEpatents

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.

Measurements of acoustic surface waves on fluid-filled porous rocks

NASA Astrophysics Data System (ADS)

Adler, Laszlo; Nagy, Peter B.

1994-09-01

Novel experimental techniques to measure ultrasonic velocity and attenuation of surface waves on fluid-filled porous natural rocks are presented. Our experimental results are consistent with the theoretical predictions of Feng and Johnson (1983). Depending on the interface conditions, i.e., whether the surface pores are open or closed, pseudo-Rayleigh, pseudo-Stoneley, and/or Stoneley surface waves may exist on fluid-saturated rocks with closed 'slow' surface wave (true Stoneley mode) on fluid-filled porous rocks with closed surface pores. The velocity and attenuation of the 'slow' surface mode may be used to assess the dynamic permeabilty of porous formations.

Anisotropic dispersion and attenuation due to wave-induced fluid flow: Quasi-static finite element modeling in poroelastic solids

NASA Astrophysics Data System (ADS)

Wenzlau, F.; Altmann, J. B.; Müller, T. M.

2010-07-01

Heterogeneous porous media such as hydrocarbon reservoir rocks are effectively described as anisotropic viscoelastic solids. They show characteristic velocity dispersion and attenuation of seismic waves within a broad frequency band, and an explanation for this observation is the mechanism of wave-induced pore fluid flow. Various theoretical models quantify dispersion and attenuation of normal incident compressional waves in finely layered porous media. Similar models of shear wave attenuation are not known, nor do general theories exist to predict wave-induced fluid flow effects in media with a more complex distribution of medium heterogeneities. By using finite element simulations of poroelastic relaxation, the total frequency-dependent complex stiffness tensor can be computed for a porous medium with arbitrary internal heterogeneity. From the stiffness tensor, velocity dispersion and frequency-dependent attenuation are derived for compressional and shear waves as a function of the angle of incidence. We apply our approach to the case of layered media and to that of an ellipsoidal poroelastic inclusion. In the case of the ellipsoidal inclusion, compressional and shear wave modes show significant attenuation, and the characteristic frequency dependence of the effect is governed by the spatiotemporal scale of the pore fluid pressure relaxation. In our anisotropic examples, the angle dependence of the attenuation is stronger than that of the velocity dispersion. It becomes clear that the spatial attenuation patterns show specific characteristics of wave-induced fluid flow, implying that anisotropic attenuation measurements may contribute to the inversion of fluid transport properties in heterogeneous porous media.

Micro-poromechanics model of fluid-saturated chemically active fibrous media.

PubMed

Misra, Anil; Parthasarathy, Ranganathan; Singh, Viraj; Spencer, Paulette

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

Relationship between P-wave attenuation and water saturation in an homogeneous unconsolidated and partially saturated porous media : An experimental study

NASA Astrophysics Data System (ADS)

Barrière, J.; Sénéchal, P.; Bordes, C.; Perroud, H.

2010-12-01

Nowadays, it is well known that hydrogeological properties of the porous media (porosity, fluid saturation and permeability) can influence seismic properties. The major theory which links hydrogeological and seismic parameters is poroelasticity proposed by Biot (1956) for saturated porous media in a wetting phase fluid. However the Biot relaxation process can't explain the level of attenuation of seismic waves generally measured on field from seismic to sonic frequency range in the case of partially saturated media. Laboratory experiments are necessary to better understand the effects of fluids on the attenuation of waves but few ones are done in the low frequency range (1Hz to 10 kHz) where the wavelength is greater than heterogeneities size. We propose an experimental study to determine the attenuation of propagative P-wave in the sonic frequency range on unconsolidated and partially saturated porous media, typical of near surface hydrogeological media. 10 accelerometers (0.0001-17kHz) and 6 capacitance probes (soil moisture sensors) are placed in a container (107 cm x 34 cm x 35cm) full of homogeneous sand (99% silica). An acoustic source (0 - 20 kHz) generate seismic waves which are recorded by the accelerometers during three cycles of imbibition-drainage (corresponding to a water saturation range from 0% to 95%). Values of attenuation (quality factor Q) versus water saturation and frequency are calculated with the well-known spectral ratio method. The spectrum of each recorded P-wave is obtained by a continuous wavelet transform, more adapted than Fourier transform for a non-stationary signal, such as seismic signal, whose frequency content varies with time. The first analyses show a strong dependence of the quality factor with frequency and water saturation, notably at high water saturation (above 60 %) where the attenuation is maximum. Knowing some important parameters of the studied media such as porosity and permeability, we interpret physically our

Benchmarking variable-density flow in saturated and unsaturated porous media

NASA Astrophysics Data System (ADS)

Guevara Morel, Carlos Roberto; Cremer, Clemens; Graf, Thomas

2015-04-01

In natural environments, fluid density and viscosity can be affected by spatial and temporal variations of solute concentration and/or temperature. These variations can occur, for example, due to salt water intrusion in coastal aquifers, leachate infiltration from waste disposal sites and upconing of saline water from deep aquifers. As a consequence, potentially unstable situations may exist in which a dense fluid overlies a less dense fluid. This situation can produce instabilities that manifest as dense plume fingers that move vertically downwards counterbalanced by vertical upwards flow of the less dense fluid. Resulting free convection increases solute transport rates over large distances and times relative to constant-density flow. Therefore, the understanding of free convection is relevant for the protection of freshwater aquifer systems. The results from a laboratory experiment of saturated and unsaturated variable-density flow and solute transport (Simmons et al., Transp. Porous Medium, 2002) are used as the physical basis to define a mathematical benchmark. The HydroGeoSphere code coupled with PEST are used to estimate the optimal parameter set capable of reproducing the physical model. A grid convergency analysis (in space and time) is also undertaken in order to obtain the adequate spatial and temporal discretizations. The new mathematical benchmark is useful for model comparison and testing of variable-density variably saturated flow in porous media.

Wave propagation through elastic porous media containing two immiscible fluids

NASA Astrophysics Data System (ADS)

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

2005-02-01

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

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

NASA Astrophysics Data System (ADS)

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

2017-12-01

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

The impact of fluid topology on residual saturations - A pore-network model study

NASA Astrophysics Data System (ADS)

Doster, F.; Kallel, W.; van Dijke, R.

2014-12-01

In two-phase flow in porous media only fractions of the resident fluid are mobilised during a displacement process and, in general, a significant amount of the resident fluid remains permanently trapped. Depending on the application, entrapment is desirable (geological carbon storage), or it should be obviated (enhanced oil recovery, contaminant remediation). Despite its utmost importance for these applications, predictions of trapped fluid saturations for macroscopic systems, in particular under changing displacement conditions, remain challenging. The models that aim to represent trapping phenomena are typically empirical and require tracking of the history of the state variables. This exacerbates the experimental verification and the design of sophisticated displacement technologies that enhance or impede trapping. Recently, experiments [1] have suggested that a macroscopic normalized Euler number, quantifying the topology of fluid distributions, could serve as a parameter to predict residual saturations based on state variables. In these experiments the entrapment of fluids was visualised through 3D micro CT imaging. However, the experiments are notoriously time consuming and therefore only allow for a sparse sampling of the parameter space. Pore-network models represent porous media through an equivalent network structure of pores and throats. Under quasi-static capillary dominated conditions displacement processes can be modeled through simple invasion percolation rules. Hence, in contrast to experiments, pore-network models are fast and therefore allow full sampling of the parameter space. Here, we use pore-network modeling [2] to critically investigate the knowledge gained through observing and tracking the normalized Euler number. More specifically, we identify conditions under which (a) systems with the same saturations but different normalized Euler numbers lead to different residual saturations and (b) systems with the same saturations and the same

Numerical simulation of convective motion in an anisotropic porous medium and cosymmetry conservation

NASA Astrophysics Data System (ADS)

Abdelhafez, M. A.; Tsybulin, V. G.

2017-10-01

The onset of convection in a porous anisotropic rectangle occupied by a heat-conducting fluid heated from below is analyzed on the basis of the Darcy-Boussinesq model. It is shown that there are combinations of control parameters for which the system has a nontrivial cosymmetry and a one-parameter family of stationary convective regimes branches off from the mechanical equilibrium. For the two-dimensional convection equations in a porous medium, finite-difference approximations preserving the cosymmetry of the original system are developed. Numerical results are presented that demonstrate the formation of a family of convective regimes and its disappearance when the approximations do not inherit the cosymmetry property.

Long-wave equivalent viscoelastic solids for porous rocks saturated by two-phase fluids

NASA Astrophysics Data System (ADS)

Santos, J. E.; Savioli, G. B.

2018-07-01

Seismic waves travelling across fluid-saturated poroelastic materials with mesoscopic-scale heterogeneities induce fluid flow and Biot's slow waves generating energy loss and velocity dispersion. Using Biot's equations of motion to model these type of heterogeneities would require extremely fine meshes. We propose a numerical upscaling procedure to determine the complex and frequency-dependent Pwave and shear moduli of an effective viscoelastic medium long-wave equivalent to a poroelastic solid saturated by a two-phase fluid. The two-phase fluid is defined in terms of capillary pressure and relative permeability flow functions. The Pwave and shear effective moduli are determined using harmonic compressibility and shear experiments applied on representative samples of the bulk material. Each experiment is associated with a boundary value problem that is solved using the finite element method. Since a poroelastic solid saturated by a two-phase fluid supports the existence of two slow waves, this upscaling procedure allows to analyse their effect on the mesoscopic loss mechanism in hydrocarbon reservoir formations. Numerical results show that a two-phase Biot medium model predicts higher attenuation than classic Biot models.

Mechanical behavior of concrete and related porous materials under partial saturation: The effective stress and the viscous softening due to movement of nanometer-scale pore fluid

NASA Astrophysics Data System (ADS)

Vlahinic, Ivan

It has been said that porous materials are like music: the gaps are as important as the filled-in bits. In other words, in addition to the solid structure, pore characteristics such as size and morphology play a crucial role in defining the overall physical properties of the porous materials. This work goes a step further and examines the behaviors of some porous media that arise when the pore network is occupied by two fluids, principally air and water, as a result of drying or wetting. Such a state gives rise to fluid capillarity which can generate significant negative fluid pressures. In the first part, a constitutive model for drying of an elastic porous medium is proposed and then extended to derive a novel expression for effective stress in partially saturated media. The model is motivated by the fact that in a system that is saturated by two different fluids, two different pressure inherently act on the surfaces of the pore network. This causes a non-uniform strain field in the solid structure, something that is not explicitly accounted for in the classic formulations of this problem. We use some standard micromechanical homogenization techniques to estimate the extent of the 'non-uniformity' and on this basis, evaluate the validity of the classic Bishop effective stress expression for partially saturated materials. In the second part, we examine a diverse class of porous materials which behave in an unexpected (and even counterintuitive) way under the internal moisture fluctuations. In particular, during wetting and drying alike, the solid viscosity of these materials appears to soften, sometimes by an order of magnitude or more. Under load, this can lead to significantly increased rates of deformations. On account of the recent experimental and theoretical findings on the nature of water flow in nanometer-size hydrophillic spaces, we provide a physical explanation for the viscous softening and propose a constitutive law on this basis. To this end, it also

Infiltration behaviour of elemental mercury DNAPL in fully and partially water saturated porous media

NASA Astrophysics Data System (ADS)

D'Aniello, Andrea; Hartog, Niels; Sweijen, Thomas; Pianese, Domenico

2018-02-01

Mercury is a contaminant of global concern due to its harmful effects on human health and for the detrimental consequences of its release in the environment. Sources of liquid elemental mercury are usually anthropogenic, such as chlor-alkali plants. To date insight into the infiltration behaviour of liquid elemental mercury in the subsurface is lacking, although this is critical for assessing both characterization and remediation approaches for mercury DNAPL contaminated sites. Therefore, in this study the infiltration behaviour of elemental mercury in fully and partially water saturated systems was investigated using column experiments. The properties affecting the constitutive relations governing the infiltration behaviour of liquid Hg0, and PCE for comparison, were determined using Pc(S) experiments with different granular porous media (glass beads and sands) for different two- and three-phase configurations. Results showed that, in water saturated porous media, elemental mercury, as PCE, acted as a non-wetting fluid. The required entry head for elemental mercury was higher (from about 5 to 7 times). However, due to the almost tenfold higher density of mercury, the required NAPL entry heads of 6.19 cm and 12.51 cm for mercury to infiltrate were 37.5% to 20.7% lower than for PCE for the same porous media. Although Leverett scaling was able to reproduce the natural tendency of Hg0 to be more prone than PCE to infiltrate in water saturated porous media, it considerably underestimated Hg0 infiltration capacity in comparison with the experimental results. In the partially water saturated system, in contrast with PCE, elemental mercury also acted as a nonwetting fluid, therefore having to overcome an entry head to infiltrate. The required Hg0 entry heads (10.45 and 15.74 cm) were considerably higher (68.9% and 25.8%) than for the water saturated porous systems. Furthermore, in the partially water saturated systems, experiments showed that elemental mercury displaced

Infiltration behaviour of elemental mercury DNAPL in fully and partially water saturated porous media.

PubMed

D'Aniello, Andrea; Hartog, Niels; Sweijen, Thomas; Pianese, Domenico

2018-02-01

Mercury is a contaminant of global concern due to its harmful effects on human health and for the detrimental consequences of its release in the environment. Sources of liquid elemental mercury are usually anthropogenic, such as chlor-alkali plants. To date insight into the infiltration behaviour of liquid elemental mercury in the subsurface is lacking, although this is critical for assessing both characterization and remediation approaches for mercury DNAPL contaminated sites. Therefore, in this study the infiltration behaviour of elemental mercury in fully and partially water saturated systems was investigated using column experiments. The properties affecting the constitutive relations governing the infiltration behaviour of liquid Hg 0 , and PCE for comparison, were determined using P c (S) experiments with different granular porous media (glass beads and sands) for different two- and three-phase configurations. Results showed that, in water saturated porous media, elemental mercury, as PCE, acted as a non-wetting fluid. The required entry head for elemental mercury was higher (from about 5 to 7 times). However, due to the almost tenfold higher density of mercury, the required NAPL entry heads of 6.19cm and 12.51cm for mercury to infiltrate were 37.5% to 20.7% lower than for PCE for the same porous media. Although Leverett scaling was able to reproduce the natural tendency of Hg 0 to be more prone than PCE to infiltrate in water saturated porous media, it considerably underestimated Hg 0 infiltration capacity in comparison with the experimental results. In the partially water saturated system, in contrast with PCE, elemental mercury also acted as a nonwetting fluid, therefore having to overcome an entry head to infiltrate. The required Hg 0 entry heads (10.45 and 15.74cm) were considerably higher (68.9% and 25.8%) than for the water saturated porous systems. Furthermore, in the partially water saturated systems, experiments showed that elemental mercury

Long-wave equivalent viscoelastic solids for porous rocks saturated by two-phase fluids

NASA Astrophysics Data System (ADS)

Santos, J. E.; Savioli, G. B.

2018-04-01

Seismic waves traveling across fluid-saturated poroelastic materials with mesoscopic-scale heterogeneities induce fluid flow and Biot's slow waves generating energy loss and velocity dispersion. Using Biot's equations of motion to model these type of heterogeneities would require extremely fine meshes. We propose a numerical upscaling procedure to determine the complex and frequency dependent P-wave and shear moduli of an effective viscoelastic medium long-wave equivalent to a poroelastic solid saturated by a two-phase fluid. The two-phase fluid is defined in terms of capillary pressure and relative permeability flow functions. The P-wave and shear effective moduli are determined using harmonic compressibility and shear experiments applied on representative samples of the bulk material. Each experiment is associated with a boundary value problem that is solved using the finite element method. Since a poroelastic solid saturated by a two-phase fluid supports the existence of two slow waves, this upscaling procedure allows to analyze their effect on the mesoscopic-loss mechanism in hydrocarbon reservoir formations. Numerical results show that a two-phase Biot medium model predicts higher attenuation than classic Biot models.

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

Effect of Carreau-Yasuda rheological parameters on subcritical Lapwood convection in horizontal porous cavity saturated by shear-thinning fluid

NASA Astrophysics Data System (ADS)

Khechiba, Khaled; Mamou, Mahmoud; Hachemi, Madjid; Delenda, Nassim; Rebhi, Redha

2017-06-01

The present study is focused on Lapwood convection in isotropic porous media saturated with non-Newtonian shear thinning fluid. The non-Newtonian rheological behavior of the fluid is modeled using the general viscosity model of Carreau-Yasuda. The convection configuration consists of a shallow porous cavity with a finite aspect ratio and subject to a vertical constant heat flux, whereas the vertical walls are maintained impermeable and adiabatic. An approximate analytical solution is developed on the basis of the parallel flow assumption, and numerical solutions are obtained by solving the full governing equations. The Darcy model with the Boussinesq approximation and energy transport equations are solved numerically using a finite difference method. The results are obtained in terms of the Nusselt number and the flow fields as functions of the governing parameters. A good agreement is obtained between the analytical approximation and the numerical solution of the full governing equations. The effects of the rheological parameters of the Carreau-Yasuda fluid and Rayleigh number on the onset of subcritical convection thresholds are demonstrated. Regardless of the aspect ratio of the enclosure and thermal boundary condition type, the subcritical convective flows are seen to occur below the onset of stationary convection. Correlations are proposed to estimate the subcritical Rayleigh number for the onset of finite amplitude convection as a function of the fluid rheological parameters. Linear stability of the convective motion, predicted by the parallel flow approximation, is studied, and the onset of Hopf bifurcation, from steady convective flow to oscillatory behavior, is found to depend strongly on the rheological parameters. In general, Hopf bifurcation is triggered earlier as the fluid becomes more and more shear-thinning.

Hydrodynamic cavitation in Stokes flow of anisotropic fluids.

PubMed

Stieger, Tillmann; Agha, Hakam; Schoen, Martin; Mazza, Marco G; Sengupta, Anupam

2017-05-30

Cavitation, the nucleation of vapour in liquids, is ubiquitous in fluid dynamics, and is often implicated in a myriad of industrial and biomedical applications. Although extensively studied in isotropic liquids, corresponding investigations in anisotropic liquids are largely lacking. Here, by combining liquid crystal microfluidic experiments, nonequilibrium molecular dynamics simulations and theoretical arguments, we report flow-induced cavitation in an anisotropic fluid. The cavitation domain nucleates due to sudden pressure drop upon flow past a cylindrical obstacle within a microchannel. For an anisotropic fluid, the inception and growth of the cavitation domain ensued in the Stokes regime, while no cavitation was observed in isotropic liquids flowing under similar hydrodynamic parameters. Using simulations we identify a critical value of the Reynolds number for cavitation inception that scales inversely with the order parameter of the fluid. Strikingly, the critical Reynolds number for anisotropic fluids can be 50% lower than that of isotropic fluids.

Hydrodynamic cavitation in Stokes flow of anisotropic fluids

NASA Astrophysics Data System (ADS)

Stieger, Tillmann; Agha, Hakam; Schoen, Martin; Mazza, Marco G.; Sengupta, Anupam

2017-05-01

Cavitation, the nucleation of vapour in liquids, is ubiquitous in fluid dynamics, and is often implicated in a myriad of industrial and biomedical applications. Although extensively studied in isotropic liquids, corresponding investigations in anisotropic liquids are largely lacking. Here, by combining liquid crystal microfluidic experiments, nonequilibrium molecular dynamics simulations and theoretical arguments, we report flow-induced cavitation in an anisotropic fluid. The cavitation domain nucleates due to sudden pressure drop upon flow past a cylindrical obstacle within a microchannel. For an anisotropic fluid, the inception and growth of the cavitation domain ensued in the Stokes regime, while no cavitation was observed in isotropic liquids flowing under similar hydrodynamic parameters. Using simulations we identify a critical value of the Reynolds number for cavitation inception that scales inversely with the order parameter of the fluid. Strikingly, the critical Reynolds number for anisotropic fluids can be 50% lower than that of isotropic fluids.

Stability of anisotropic self-gravitating fluids

NASA Astrophysics Data System (ADS)

Ahmad, S.; Jami, A. Rehman; Mughal, M. Z.

2018-06-01

The aim of this paper is to study the stability as well as the existence of self-gravitating anisotropic fluids in Λ-dominated era. Taking a cylindrically symmetric and static spacetime, we computed the corresponding equations of motion in the background of anisotropic fluid distributions. The realistic formulation of energy momentum tensor as well as theoretical model of the scale factors are considered in order to describe some physical properties of the anisotropic fluids. To find the stability of the compact star, we have used Herrera’s technique which is based on finding the radial and the transverse components of the speed of sound. Moreover, the behaviors of other physical quantities are also discussed like anisotropy, matching conditions of interior metric and exterior metric and compactness of the compact structures are also discussed.

A new numerical benchmark for variably saturated variable-density flow and transport in porous media

NASA Astrophysics Data System (ADS)

Guevara, Carlos; Graf, Thomas

2016-04-01

In subsurface hydrological systems, spatial and temporal variations in solute concentration and/or temperature may affect fluid density and viscosity. These variations could lead to potentially unstable situations, in which a dense fluid overlies a less dense fluid. These situations could produce instabilities that appear as dense plume fingers migrating downwards counteracted by vertical upwards flow of freshwater (Simmons et al., Transp. Porous Medium, 2002). As a result of unstable variable-density flow, solute transport rates are increased over large distances and times as compared to constant-density flow. The numerical simulation of variable-density flow in saturated and unsaturated media requires corresponding benchmark problems against which a computer model is validated (Diersch and Kolditz, Adv. Water Resour, 2002). Recorded data from a laboratory-scale experiment of variable-density flow and solute transport in saturated and unsaturated porous media (Simmons et al., Transp. Porous Medium, 2002) is used to define a new numerical benchmark. The HydroGeoSphere code (Therrien et al., 2004) coupled with PEST (www.pesthomepage.org) are used to obtain an optimized parameter set capable of adequately representing the data set by Simmons et al., (2002). Fingering in the numerical model is triggered using random hydraulic conductivity fields. Due to the inherent randomness, a large number of simulations were conducted in this study. The optimized benchmark model adequately predicts the plume behavior and the fate of solutes. This benchmark is useful for model verification of variable-density flow problems in saturated and/or unsaturated media.

Hydrodynamic cavitation in Stokes flow of anisotropic fluids

PubMed Central

Stieger, Tillmann; Agha, Hakam; Schoen, Martin; Mazza, Marco G.; Sengupta, Anupam

2017-01-01

Cavitation, the nucleation of vapour in liquids, is ubiquitous in fluid dynamics, and is often implicated in a myriad of industrial and biomedical applications. Although extensively studied in isotropic liquids, corresponding investigations in anisotropic liquids are largely lacking. Here, by combining liquid crystal microfluidic experiments, nonequilibrium molecular dynamics simulations and theoretical arguments, we report flow-induced cavitation in an anisotropic fluid. The cavitation domain nucleates due to sudden pressure drop upon flow past a cylindrical obstacle within a microchannel. For an anisotropic fluid, the inception and growth of the cavitation domain ensued in the Stokes regime, while no cavitation was observed in isotropic liquids flowing under similar hydrodynamic parameters. Using simulations we identify a critical value of the Reynolds number for cavitation inception that scales inversely with the order parameter of the fluid. Strikingly, the critical Reynolds number for anisotropic fluids can be 50% lower than that of isotropic fluids. PMID:28555615

Capillary Imbibition of Hydraulic Fracturing Fluids into Partially Saturated Shale

NASA Astrophysics Data System (ADS)

Birdsell, D.; Rajaram, H.; Lackey, G.

2015-12-01

Understanding the migration of hydraulic fracturing fluids injected into unconventional reservoirs is important to assess the risk of aquifer contamination and to optimize oil and gas production. Capillary imbibition causes fracturing fluids to flow from fractures into the rock matrix where the fluids are sequestered for geologically long periods of time. Imbibition could explain the low amount of flowback water observed in the field (5-50% of the injected volume) and reduce the chance of fracturing fluid migrating out of formation towards overlying aquifers. We present calculations of spontaneous capillary imbibition in the form of an "imbibition rate parameter" (A) based on the only known exact analytical solution for spontaneous capillary imbibition. A depends on the hydraulic and capillary properties of the reservoir rock, the initial water saturation, and the viscosities of the wetting and nonwetting fluids. Imbibed volumes can be large for a high permeability shale gas reservoir (up to 95% of the injected volume) or quite small for a low permeability shale oil reservoir (as low as 3% of the injected volume). We also present a nondimensionalization of the imbibition rate parameter, which facilitates the calculation of A and clarifies the relation of A to initial saturation, porous medium properties, and fluid properties. Over the range of initial water saturations reported for the Marcellus shale (0.05-0.6), A varies by less than factors of ~1.8 and ~3.4 for gas and oil nonwetting phases respectively. However, A decreases significantly for larger initial water saturations. A is most sensitive to the intrinsic permeability of the reservoir rock and the viscosity of the fluids.

A note on flow reversal in a wavy channel filled with anisotropic porous material

NASA Astrophysics Data System (ADS)

Karmakar, Timir; Raja Sekhar, G. P.

2017-07-01

Viscous flow through a symmetric wavy channel filled with anisotropic porous material is investigated analytically. Flow inside the porous bed is assumed to be governed by the anisotropic Brinkman equation. It is assumed that the ratio of the channel width to the wavelength is small (i.e. δ2≪1). The problem is solved up to O(δ2) assuming that δ2λ2≪1, where λ is the anisotropic ratio. The key purpose of this paper is to study the effect of anisotropic permeability on flow near the crests of the wavy channel which causes flow reversal. We present a detailed analysis of the flow reversal at the crests. The ratio of the permeabilities (anisotropic ratio) is responsible for the flow separation near the crests of the wall where viscous forces are effective. For a flow configuration (say, low amplitude parameter) in which there is no separation if the porous media is isotropic, introducing anisotropy causes flow separation. On the other hand, interestingly, flow separation occurs even in the case of isotropic porous medium if the amplitude parameter a is large.

Numerical study at moderate Reynolds number of peristaltic flow of micropolar fluid through a porous-saturated channel in magnetic field

NASA Astrophysics Data System (ADS)

Ahmed, Bilal; Javed, Tariq; Ali, N.

2018-01-01

This paper analyzes the MHD flow of micropolar fluid induced by peristaltic waves passing through the porous saturated channel at large Reynolds number. The flow model is formulated in the absence of assumptions of lubrication theory which yields the governing equations into a non-linear set of coupled partial differential equations which allows studying the peristaltic mechanism at non-zero Reynolds and wave numbers. The influence of other involved parameters on velocity, stream function and microrotation are discussed through graphs plotted by using Galerkin's finite element method. Besides that, the phenomena of pumping and trapping are also analyzed in the later part of the paper. To ensure the accuracy of the developed code, obtained results are compared with the results available in the literature and found in excellent agreement. It is found that the peristalsis mixing can be enhanced by increasing Hartmann number while it reduces by increasing permeability of the porous medium.

Scattering rings in optically anisotropic porous silicon

NASA Astrophysics Data System (ADS)

Oton, C. J.; Gaburro, Z.; Ghulinyan, M.; Pancheri, L.; Bettotti, P.; Negro, L. Dal; Pavesi, L.

2002-12-01

We report the observation of strongly anisotropic scattering of laser light at oblique incidence on a (100)-oriented porous silicon layer. The scattered light forms cones tangent to the incident and reflected beams. The conical pattern is caused by scattering on the vertical walls of pores, which are straight along the layer thickness. The light cone defines structured light rings onto a screen normal to the cone axis. We explain the various structures by optical anisotropy of porous silicon. For the sample under analysis, we directly measure from the ring patterns a value of Δn/nord=8% of positive birefringence.

On Thermal Instability of Kuvshiniski Fluid with Suspended Particles Saturated in a Porous Medium in the Presence of a Magnetic Field June 13, 2017

NASA Astrophysics Data System (ADS)

Singh, M.

2017-12-01

The thermal instability of a Kuvshiniski viscoelastic fluid is considered to include the effects of a uniform horizontal magnetic field, suspended particles saturated in a porous medium. The analysis is carried out within the framework of the linear stability theory and normal mode technique. For the case of stationary convection, the Kuvshiniski viscoelastic fluid behaves like a Newtonian fluid and the magnetic field has a stabilizing effect, whereas medium permeability and suspended particles are found to have a destabilizing effect on the system, oscillatory modes are introduced in the system, in the absence of these the principle of exchange of stabilities is valid. Graphs in each case have been plotted by giving numerical values to the parameters, depicting the stability characteristics. Sufficient conditions for the avoidance of overstability are also obtained.

Direct laboratory observation of fluid distribution and its influence on acoustic properties of patchy saturated rocks

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

Probing the water distribution in porous model sands with two immiscible fluids: A nuclear magnetic resonance micro-imaging study

NASA Astrophysics Data System (ADS)

Lee, Bum Han; Lee, Sung Keun

2017-10-01

The effect of the structural heterogeneity of porous networks on the water distribution in porous media, initially saturated with immiscible fluid followed by increasing durations of water injection, remains one of the important problems in hydrology. The relationship among convergence rates (i.e., the rate of fluid saturation with varying injection time) and the macroscopic properties and structural parameters of porous media have been anticipated. Here, we used nuclear magnetic resonance (NMR) micro-imaging to obtain images (down to ∼50 μm resolution) of the distribution of water injected for varying durations into porous networks that were initially saturated with silicone oil. We then established the relationships among the convergence rates, structural parameters, and transport properties of porous networks. The volume fraction of the water phase increases as the water injection duration increases. The 3D images of the water distributions for silica gel samples are similar to those of the glass bead samples. The changes in water saturation (and the accompanying removal of silicone oil) and the variations in the volume fraction, specific surface area, and cube-counting fractal dimension of the water phase fit well with the single-exponential recovery function { f (t) = a [ 1 -exp (- λt) ] } . The asymptotic values (a, i.e., saturated value) of the properties of the volume fraction, specific surface area, and cube-counting fractal dimension of the glass bead samples were greater than those for the silica gel samples primarily because of the intrinsic differences in the porous networks and local distribution of the pore size and connectivity. The convergence rates of all of the properties are inversely proportional to the entropy length and permeability. Despite limitations of the current study, such as insufficient resolution and uncertainty for the estimated parameters due to sparsely selected short injection times, the observed trends highlight the first

A lattice Boltzmann investigation of steady-state fluid distribution, capillary pressure and relative permeability of a porous medium: Effects of fluid and geometrical properties

NASA Astrophysics Data System (ADS)

Li, Zi; Galindo-Torres, Sergio; Yan, Guanxi; Scheuermann, Alexander; Li, Ling

2018-06-01

Simulations of simultaneous steady-state two-phase flow in the capillary force-dominated regime were conducted using the state-of-the-art Shan-Chen multi-component lattice Boltzmann model (SCMC-LBM) based on two-dimensional porous media. We focused on analyzing the fluid distribution (i.e., WP fluid-solid, NP fluid-solid and fluid-fluid interfacial areas) as well as the capillary pressure versus saturation curve which was affected by fluid and geometrical properties (i.e., wettability, adhesive strength, pore size distribution and specific surface area). How these properties influenced the relative permeability versus saturation relation through apparent effective permeability and threshold pressure gradient was also explored. The SCMC-LBM simulations showed that, a thin WP fluid film formed around the solid surface due to the adhesive fluid-solid interaction, resulting in discrete WP fluid distributions and reduction of the WP fluid mobility. Also, the adhesive interaction provided another source of capillary pressure in addition to capillary force, which, however, did not affect the mobility of the NP fluid. The film fluid effect could be enhanced by large adhesive strength and fine pores in heterogeneous porous media. In the steady-state infiltration, not only the NP fluid but also the WP fluid were subjected to the capillary resistance. The capillary pressure effect could be alleviated by decreased wettability, large average pore radius and improved fluid connectivity in heterogeneous porous media. The present work based on the SCMC-LBM investigations elucidated the role of film fluid as well as capillary pressure in the two-phase flow system. The findings have implications for ways to improve the macroscopic flow equation based on balance of force for the steady-state infiltration.

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

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.

Stability Analysis and Internal Heating Effect on Oscillatory Convection in a Viscoelastic Fluid Saturated Porous Medium Under Gravity Modulation

NASA Astrophysics Data System (ADS)

Bhadauria, B. S.; Singh, M. K.; Singh, A.; Singh, B. K.; Kiran, P.

2016-12-01

In this paper, we investigate the combined effect of internal heating and time periodic gravity modulation in a viscoelastic fluid saturated porous medium by reducing the problem into a complex non-autonomous Ginzgburg-Landau equation. Weak nonlinear stability analysis has been performed by using power series expansion in terms of the amplitude of gravity modulation, which is assumed to be small. The Nusselt number is obtained in terms of the amplitude for oscillatory mode of convection. The influence of viscoelastic parameters on heat transfer has been discussed. Gravity modulation is found to have a destabilizing effect at low frequencies and a stabilizing effect at high frequencies. Finally, it is found that overstability advances the onset of convection, more with internal heating. The conditions for which the complex Ginzgburg-Landau equation undergoes Hopf bifurcation and the amplitude equation undergoes supercritical pitchfork bifurcation are studied.

Klinkenberg effect in hydrodynamics of gas flow through anisotropic porous materials

NASA Astrophysics Data System (ADS)

Wałowski, Grzegorz; Filipczak, Gabriel

2017-10-01

This study discusses results of experiments on hydrodynamic assessment of gas flow through backbone (skeletal) porous materials with an anisotropic structure. The research was conducted upon materials of diversified petrographic characteristics, both natural origin (rocky, pumice) and process materials (char and coke). The study was conducted for a variety of hydrodynamic conditions, using air, as well as for nitrogen and carbon dioxide. The basis for assessing hydrodynamics of gas flow through porous material was a gas stream that results from the pressure forcing such flow. The results of measurements indicate a clear impact of the type of material on the gas permeability, and additionally - as a result of their anisotropic internal structure - to a significant effect of the flow direction on the value of gas stream.

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.

Spin echo SPI methods for quantitative analysis of fluids in porous media.

PubMed

Li, Linqing; Han, Hui; Balcom, Bruce J

2009-06-01

Fluid density imaging is highly desirable in a wide variety of porous media measurements. The SPRITE class of MRI methods has proven to be robust and general in their ability to generate density images in porous media, however the short encoding times required, with correspondingly high magnetic field gradient strengths and filter widths, and low flip angle RF pulses, yield sub-optimal S/N images, especially at low static field strength. This paper explores two implementations of pure phase encode spin echo 1D imaging, with application to a proposed new petroleum reservoir core analysis measurement. In the first implementation of the pulse sequence, we modify the spin echo single point imaging (SE-SPI) technique to acquire the k-space origin data point, with a near zero evolution time, from the free induction decay (FID) following a 90 degrees excitation pulse. Subsequent k-space data points are acquired by separately phase encoding individual echoes in a multi-echo acquisition. T(2) attenuation of the echo train yields an image convolution which causes blurring. The T(2) blur effect is moderate for porous media with T(2) lifetime distributions longer than 5 ms. As a robust, high S/N, and fast 1D imaging method, this method will be highly complementary to SPRITE techniques for the quantitative analysis of fluid content in porous media. In the second implementation of the SE-SPI pulse sequence, modification of the basic measurement permits fast determination of spatially resolved T(2) distributions in porous media through separately phase encoding each echo in a multi-echo CPMG pulse train. An individual T(2) weighted image may be acquired from each echo. The echo time (TE) of each T(2) weighted image may be reduced to 500 micros or less. These profiles can be fit to extract a T(2) distribution from each pixel employing a variety of standard inverse Laplace transform methods. Fluid content 1D images are produced as an essential by product of determining the

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

NASA Astrophysics Data System (ADS)

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

2017-08-01

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