Ultrasonic elastic‐wave velocities, as well as the stress‐strain behavior, were measured on dry dune sand from Saudi Arabia under hydrostatic loading from 0 to 50 MPa. The samples came from the dune's crest, limb, and base. The experiments included several loading and unloading cycles. All samples exhibited irreversible deformation and porosity reduction during the tests due to grain repacking in low‐stress regimes and grain breakage at higher stress. The loading behavior was plastic, while the unloading/reloading behavior was elastic. In spite of noticeable porosity variations during the experiments, the elastic‐wave velocities mainly depended on the stress rather than on porosity. Our velocity versus stress data is very close to earlier published results, although the provenance of the sand is very different, pointing at a universality of such behavior in loose sand. Our results were matched by the contact Hertz‐Mindlin theory with varying coordination number and shear‐stiffness adjustment. We also offer a new theoretical model for the static bulk modulus measured in these experiments. This model is based on the concept of stress heterogeneity (stress chains or patches) in a granular pack and uses elastic bound averaging of the respective moduli of stressed and unstressed parts of the pack. The same model is used to describe the velocity versus stress behavior. These results and theory constitute a new installment into the published sparse data of the mechanical properties of unconsolidated sand as well as theoretical analyses thereof and are the first ever such data generated in Saudi Arabia.

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 V p , V s , 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.

A carbonate sample extracted from the depth of about 10 kft was subjected to uniaxial loading while the confining stress remained constant. Post-experiment inspection of the sample showed an inclined crack at an angle less than 20° to the horizontal. This subhorizontal crack orientation was contrary to the expected 45° inclination, the plane of the maximum shear stress. Coincidentally, as shown by CT-scan prior to loading, there was a boundary between two layers of different density inside the sample located almost exactly where the crack appeared. This density difference has arguably translated into the contrast in the elastic properties at the boundary. The hypothesis is that because of this elastic heterogeneity, an incipient crack developed at the boundary due to the unavoidable tensile stressing of the sample as it was brought to the benchtop from its original state of high confining stress at depth. Controlled uniaxial compression made the sample slip along this crack, which then developed into a prominent feature. This assumption was corroborated by a numerical experiment showing a strong von Mises stress concentration at the elastic contrast boundary during hydrostatic tensile loading. Another sample, from the same formation, but without strong density heterogeneity, exhibited a classic 45° crack after uniaxial loading. These results provide a novel and important insight into the mechanics, breakage, and strength of natural rock.

No abstract

Sonic and dipole wireline tools measure Vp and Vs along the vertical direction. The state of stress in the subsurface is predominantly anisotropic, while most laboratory experiments measuring the dynamic elastic properties are conducted under hydrostatic stress. The question we ask is whether such laboratory experiments provide the velocities that are close to those measured in the vertical wellbore where the stresses are anisotropic. To address this question, we conducted ultrasonic pulse transmission experiments on several room-dry rock samples. The comparison was made between the P- and S-wave velocities obtained at pure hydrostatic loading conditions and those at a smaller hydrostatic stress with added axial stress, so that the total stress along the axis of a cylindrical plug was the same as under pure hydrostatic loading. These differences were significant in the extreme case of only 1 MPa hydrostatic confining stress with the axial stress increasing up to 40 MPa. However, as the hydrostatic (confining) stress increased, the differences between the velocities along the axis of the sample became smaller and smaller. For example, at 1 MPa confining and 30 MPa axial stress, the relative difference in Vp was about 10%, while that in Vs was about 20%. However, at 10 MPa confining stress, these differences became about 3% and 6%, respectively, and further decreased as the confining stress increased. This means that even at strong in-situ contrasts between the vertical and horizontal stresses, the results of laboratory hydrostatic experiments can be used for in-situ velocity estimates. These results also appear to be consistent with a theoretical model that predicts the directional velocities at any triaxial stress conditions from those measured versus hydrostatic stress.