Sediment with porous grains: Rock-physics model and application to marine carbonate and opal

Geophysics ◽  
2009 ◽  
Vol 74 (1) ◽  
pp. E1-E15 ◽  
Author(s):  
Franklin Ruiz ◽  
Jack Dvorkin

We offer an effective-medium model for estimating the elastic properties of high-porosity marine calcareous sediment and diatomite. This model treats sediment as a pack of porous elastic grains. The effective elastic moduli of the porous grains are calculated using the differential effective-medium (DEM) model, whereby the intragranular ellipsoidal inclusions have a fixed aspect ratio and are filled with seawater. Then the elastic moduli of a pack of these spherical grains are calculated using a modified (scaled to the critical porosity) upper Hashin-Shtrikman bound above the critical porosity and modified lower (carbonates) and upper (opal) Hashin-Shtrikman bounds below the critical porosity. The best match between the model-predicted compressional- and shear-wave velocities and Ocean Drilling Program (ODP) data from three wells is achieved when the aspect ratio of intragranular pores is 0.5. This model assigns finite, nonzero values to the shear modulus of high-porosity marine sediment, unlike the suspension model commonly used in such depositional settings. The approach also allows one to obtain a satisfactory match with laboratory diatomite velocity data.

ICIPEG 2014 ◽  
2015 ◽  
pp. 313-321
Author(s):  
Ida Ayu Purnamasari ◽  
Wan Ismail Wan Yusoff ◽  
Chow Weng Sum

Geophysics ◽  
1996 ◽  
Vol 61 (2) ◽  
pp. 444-455 ◽  
Author(s):  
Jack Dvorkin ◽  
Amos Nur ◽  
Caren Chaika

Our observations made on dry‐sandstone ultrasonic velocity data relate to the variation in velocity (or modulus) with effective stress, and the ability to predict a velocity for a rock under one effective pressure when it is known only under a different effective pressure. We find that the sensitivity of elastic moduli, and velocities, to effective hydrostatic stress increases with decreasing porosity. Specifically, we calculate the difference between an elastic modulus, [Formula: see text], of a sample of porosity ϕ at effective pressure [Formula: see text] and the same modulus, [Formula: see text], at effective pressure [Formula: see text]. If this difference, [Formula: see text], is plotted versus porosity for a suite of samples, then the scatter of ΔM is close to zero as porosity approaches the critical porosity value, and reaches its maximum as porosity approaches zero. The dependence of this scatter on porosity is close to linear. Critical porosity here is the porosity above which rock can exist only as a suspension—between 36% and 40% for sandstones. This stress‐sensitivity pattern of grain‐supported sandstones (clay content below 0.35) practically does not depend on clay content. In practical terms, the uncertainty of determining elastic moduli at a higher effective stress from the measurements at a lower effective stress is small at high porosity and increases with decreasing porosity. We explain this effect by using a combination of two heuristic models—the critical porosity model and the modified solid model. The former is based on the observation that the elastic‐modulus‐versus‐porosity relation can be approximated by a straight line that connects two points in the modulus‐porosity plane: the modulus of the solid phase at zero porosity and zero at critical porosity. The second one reflects the fact that at constant effective stress, low‐porosity sandstones (even with small amounts of clay) exhibit large variability in elastic moduli. We attribute this variability to compliant cracks that hardly affect porosity but strongly affect the stiffness. The above qualitative observation helps to quantitatively constrain P‐ and S‐wave velocities at varying stresses from a single measurement at a fixed stress. We also show that there are distinctive linear relations between Poisson’s ratios (ν) of sandstone measured at two different stresses. For example, in consolidated medium‐porosity sandstones [Formula: see text], where the subscripts indicate hydrostatic stress in MPa. Linear functions can also be used to relate the changes (with hydrostatic stress) in shear moduli to those in compressional moduli. For example, [Formula: see text], where [Formula: see text] is shear modulus and [Formula: see text] is compressional modulus, both in GPa, and the subscripts indicate stress in MPa.


1995 ◽  
Vol 22 (5) ◽  
pp. 555-558 ◽  
Author(s):  
Tapan Mukerji ◽  
Jim Berryman ◽  
Gary Mavko ◽  
Patricia Berge

Geophysics ◽  
2017 ◽  
Vol 82 (2) ◽  
pp. D135-D156 ◽  
Author(s):  
Paul Sayar ◽  
Carlos Torres-Verdín

Calculation of the velocity of elastic waves that propagate in saturated rocks is complicated by many dispersion and attenuation mechanisms. Wave-induced fluid flow at pore and mesoscopic scales is responsible for substantial dispersion and attenuation effects. In addition, dispersion and attenuation phenomena in rocks are due to spatial heterogeneities in density and elastic properties. Whenever the wavelength of a propagating wave is of the same scale as the size of the pores, dispersion takes place in the form of multiple scattering. We have developed an effective medium model that reproduced the combined effects on acoustic velocities of four wave-attenuation mechanisms by replicating results predicted by existing models in saturated isotropic rocks, Biot’s and squirt flow, and acoustic scattering from spherical pores and randomly oriented penny-shaped cracks. First, the strengthening of rocks and wave velocity increase with frequency associated with fluid-related dispersion was modeled by introducing the concept of complex-valued and frequency-dependent equivalent moduli for the solid background. Effective properties were calculated by invoking the newly defined dynamic solid moduli and the classic self-consistent approximation theory. Second, we replicated the geometric effect of idealized spherical and spheroidal porous inclusions on propagating waves with two existing dynamic self-consistent embedding schemes. In doing so, we explicitly derived the equations of the self-consistent model previously developed for the analysis of waves that propagated in a rock containing randomly oriented cracks. Finally, we established an original procedure to combine all dispersion effects into one simple workflow. We examined the performance of the model for a variety of cases, and consideration was given to the possibility of quantifying petrophysical and rock-physics parameters, such as the presence and shape of porous inclusions, permeability, and fluid viscosity, from seismic, sonic, and ultrasonic measurements.


2021 ◽  
Author(s):  
Phillip Cilli ◽  
Mark Chapman

<p>The combination of electrical and elastic measurements can be useful in lowering subsurface characterisation uncertainty. The majority of rock physics relations which relate a rock’s electrical and elastic properties, however, rely on the estimation of porosity as an intermediate step. By combining differential effective medium schemes which relate a rock’s electrical and elastic properties to porosity and pore shape, we obtain cross-property expressions which are independent of porosity, depending only on pore aspect ratio. Analysing published joint electrical-elastic measurements shows the cross-property model works well for clean sandstones, and models V<sub>p</sub>/V<sub>s</sub> ratios as a function of resistivity without porosity. Although clay-bearing sandstones are more complex, our model can still identify the correct trends. On theoretical grounds, it seems our approach has the potential to produce additional cross-property relations; a topic upon which we speculate.</p>


Geophysics ◽  
2006 ◽  
Vol 71 (1) ◽  
pp. N11-N19 ◽  
Author(s):  
Ayako Kameda ◽  
Jack Dvorkin ◽  
Youngseuk Keehm ◽  
Amos Nur ◽  
William Bosl

Numerical simulation of laboratory experiments on rocks, or digital rock physics, is an emerging field that may eventually benefit the petroleum industry. For numerical experimentation to find its way into the mainstream, it must be practical and easily repeatable — i.e., implemented on standard hardware and in real time. This condition reduces the size of a digital sample to just a few grains across. Also, small physical fragments of rock, such as cuttings, may be the only material available to produce digital images. Will the results be meaningful for a larger rock volume? To address this question, we use a number of natural and artificial medium- to high-porosity, well-sorted sandstones. The 3D microtomography volumes are obtained from each physical sample. Then, analogous to making thin sections of drill cuttings, we select a large number of small 2D slices from a 3D scan. As a result, a single physical sample produces hundreds of 2D virtual-drill-cuttings images. Corresponding 3D pore-space realizations are generated statistically from these 2D images; fluid flow is simulated in three dimensions, and the absolute permeability is computed. The results show that small fragments of medium– to high-porosity sandstones that are statistically subrepresentative of a larger sample will not yield the exact porosity and permeability of the sample. However, a significant number of small fragments will yield a site-specific permeability-porosity trend that can then be used to estimate the absolute permeability from independent porosity data obtained in the well or inferred from seismic techniques.


2008 ◽  
Vol 75 (14) ◽  
pp. 4104-4116 ◽  
Author(s):  
Dan Su ◽  
Michael H. Santare ◽  
George A. Gazonas

Sign in / Sign up

Export Citation Format

Share Document