Search

Published laboratory elastic-wave velocity versus porosity data in carbonate rocks exhibit significant scatter even at a fixed mineralogy. This scatter is usually attributed to the strong variability in the rock-frame or pore-space geometry, which, in turn, is driven by the richness and complexity of diagenetic alteration in these very reactive sediments. Yet, by examining wireline data from oil-bearing high-to-medium porosity chalk deposits, we find surprisingly tight velocity–porosity trends. Moreover, these trends are continued into the low-porosity domain by data from a location thousands of miles away from the chalk field. This congruence implies a universality of diagenetic trends, at least in the massive deposits under examination. We also find that the elastic bulk and shear moduli of the pure-calcite end member are somewhat smaller than such values reported in the literature. Using the end-member elastic constants relevant to the data under examination, we establish a theoretical rock physics model to match and generalize these data. [ABSTRACT FROM AUTHOR]

ABSTRACT: We explore a package of parallel porous layers, each filled with a different fluid. Assume that this package is sampled by an elastic wave with the wavelength much larger than the thickness of an individual layer. Also assume that the layers are hydraulically isolated from each other, meaning that the diffusion length is smaller than that of the individual layer. This assumption is relevant to a patchy saturation scenario. Suppose that we wish to conduct the fluid substitution operation on this package treated as a single porous elastic body. What is the effective bulk modulus of the pore fluid to be used in this operation that will result in the same elastic modulus as computed by Backus averaging the individual moduli of the layers? We address this question analytically by assuming that the porosity, dry frame, and the mineral matrix properties of the individual layers are the same for all layers. The only difference between the layers is the pore fluid. We find that the resulting effective bulk modulus of the fluid thus derived falls between the arithmetic and harmonic averages of the fluid bulk moduli in the layers. It can be approximated by a linear combination of these two bounds where the weights are 0.50 and 0.50 or 0.75 for the arithmetic average and 0.25 for the harmonic average, depending on the elastic moduli of the dry frame, the mineral, and the pore fluids. This solution also provides a relation between the effective bulk modulus of the pore fluid in the system under examination and water saturation to be used in the fluid substitution operation at a coarse spatial scale. [ABSTRACT FROM AUTHOR]

A patchy saturated two‐layered porous rock, with each layer filled with a different fluid, is examined. We assume that this system is probed by an elastic wave having a wavelength much larger than the layers' thicknesses. We also assume that the diffusion length is smaller than the thickness of an individual layer, which implies hydraulic disconnection between layers. In the context of Gassmann's fluid substitution theory, we analytically derived an exact expression for the effective fluid bulk modulus assuming that both layers have the same porosity, dry frame, and mineral matrix properties. In addition we derived an approximate solution that works well at relatively high porosities. Both solutions are expressed as a weighted average of the arithmetic and harmonic averages of individual bulk moduli of the pore fluid. These weights are explicitly given as functions of the porosity, the fractional thicknesses of both layers, and the elastic moduli of the constituents. For the approximate solution, one does not require explicit knowledge of the shear modulus of the rock. The comparison with laboratory data showed that, in the case where a porous, isotropic rock is filled with water and gas, the approximate solution can be used to model the measured data for high values of water saturation. Key Points: We derive exact and approximate analytic solutions to model the effective fluid bulk modulus of patchy saturated porous rocksWe compare and validate the solutions with stochastic simulation and laboratory data [ABSTRACT FROM AUTHOR]

Sand and glass bead samples were monotonically loaded from 0 to 50 MPa and then monotonically unloaded, all under hydrostatic stress conditions. Sand deformed irreversibly and its porosity loss was permanent. In contrast, the glass bead deformation was almost fully reversible with the initial porosity recovered. During loading, the static bulk modulus in all sand samples increased from 0 to about 0.5 GPa between 0‐ and 15‐MPa stress and remained constant for the remainder of the loading. Such behavior, somewhat unexpected in a granular system, means that sand deforms as a linear elastic body during loading from 15 to 50 MPa. The unloading behavior was very different—the static bulk modulus almost linearly decreased from about 4 GPa at 50 MPa to zero as the stress was reduced. Once again, this behavior was very different from that observed in glass beads, where the static bulk modulus during loading was essentially the same as during unloading and appeared to be a unique function of stress, independent of the direction of its variation. Moreover, the static and dynamic bulk moduli in glass beads, the latter computed from Vp, Vs, and density measured during the loading/unloading cycle, were very close to each other. In contrast, the dynamic bulk modulus in dry sand exceeded the static modulus during unloading by a factor of 1.5 to 2.0. We used a quantitative theory based on contact mechanics and the concept of stress heterogeneity in particulates to analytically match the observed results. Key Points: Measurements of static moduli of granular materialsDeformation of sand and glass beadsTheoretical explanation of results [ABSTRACT FROM AUTHOR]