Effect of Fluids on the Elastic Properties of 3D-Printed Anisotropic Rock Models

2021 ◽  
Vol 62 (5) ◽  
pp. 537-552
Author(s):  
Suresh Dande ◽  
◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
◽  
...  
Keyword(s):  
Wave Propagation ◽  
Elastic Properties ◽  
P Wave ◽  
Engine Oil ◽  
Rock Properties ◽  
Anisotropic Rock ◽  
Wave Velocities ◽  
Inclusion Model ◽  
3D Printed ◽  
Primary Porosity

Laboratory physical models play an important role in understanding rock properties and wave propagation, both theoretically and at the field scale. In some cases, 3D-printing technology can be adopted to construct complex rock models faster, more inexpensively, and with more specific features than previous model-building techniques. In this study, we use 3D-printed rock models to assist in understanding the effects of various fluids (air, water, engine oil, crude oil, and glycerol) on the models’ elastic properties. We first used a 3D-printed, 1-in. cube-shaped layered model. This model was created with a 6% primary porosity and a bulk density of 0.98 g/cc with VTI anisotropy. We next employed a similar cube but with horizontal inclusions embedded in the layered background, which contributed to its total 24% porosity (including primary porosity). For air to liquid saturation, P-velocities increased for all liquids in both models, with the highest increase being with glycerol (57%) and an approximately 45% increase for other fluids in the inclusion model. For the inclusion model (dry and saturated), we observed a greater difference between two orthogonally polarized S-wave velocities (Vs1 and Vs2) than between two P-wave velocities (VP0 and VP90). We attribute this to the S2-wave (polarized normal to both the layering and the plane of horizontal inclusions), which appears more sensitive to horizontal inclusions than the P-wave. For the inclusion model, Thomsen’s P-wave anisotropic parameter (ɛ) decreased from 26% for the air case to 4% for the water-saturated cube and to 1% for glycerol saturation. The small difference between the bulk modulus of the frame and the pore fluid significantly reduces the velocity anisotropy of the medium, making it almost isotropic. We compared our experimental results with theory and found that predictions using Schoenberg’s linear slip theory combined with Gassmann’s anisotropic equation were closer to actual measurements than Hudson’s isotropic calculations. This work provides insights into the usefulness of 3D-printed models to understand elastic rock properties and wave propagation under various fluid saturations.

3D-printed rock models: Elastic properties and the effects of penny-shaped inclusions with fluid substitution

Geophysics ◽  
2016 ◽  
Vol 81 (6) ◽  
pp. D669-D677 ◽  
Author(s):  
Long Huang ◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
Jose Baez-Franceschi
Keyword(s):  
3D Printing ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Bedding Plane ◽  
P Wave ◽  
Fluid Substitution ◽  
S Wave ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
3D Printed

3D printing techniques (additive manufacturing) using different materials and structures provide opportunities to understand porous or fractured materials and fluid effects on their elastic properties. We used a 3D printer (Stratasys Dimension SST 768) to print one “solid” cube model and another with penny-shaped inclusions. The 3D printing process builds materials, layer by layer, producing a slight “bedding” plane, somewhat similar to a sedimentary process. We used ultrasonic transducers (500 kHz) to measure the P- and S-wave velocities. The input printing material was thermoplastic with a density of [Formula: see text], P-wave velocity of [Formula: see text], and S-wave velocity of [Formula: see text]. The solid cube had a porosity of approximately 6% and a density of [Formula: see text]. Its P-wave velocity was [Formula: see text] in the bedding direction and [Formula: see text] normal to bedding. We observed S-wave splitting with fast and slow velocities of 879 and [Formula: see text], respectively. Quality factors for P- and S-waves were estimated using the spectral-ratio method with [Formula: see text] ranging from 15 to 17 and [Formula: see text] from 24 to 27. By introducing penny-shaped inclusions along the bedding direction in a 3D printed cube, we created a more porous volume with density of [Formula: see text] and porosity of 24%. The inclusions significantly decreased the P-wave velocity to 1706 and [Formula: see text] parallel and normal to the bedding plane. The fast and slow S-wave velocities also decreased to 812 and [Formula: see text]. A fluid substitution experiment, performed with water, increased (20%–46%) P-wave velocities and decreased (9%–10%) S-wave velocities. Theoretical predictions using Schoenberg’s linear-slip theory and Hudson’s penny-shaped theory were calculated, and we found that both theories matched the measurements closely (within 5%). The 3D printed material has interesting and definable properties and is an exciting new material for understanding wave propagation, rock properties, and fluid effects.

scholarly journals Upscaling of elastic properties in carbonates: A modeling approach based on a multiscale geophysical data set

2020 ◽  
Author(s):  
Jerome Fortin ◽  
Cedric Bailly ◽  
Mathilde Adelinet ◽  
Youri Hamon
Keyword(s):  
Elastic Properties ◽  
Wave Velocity ◽  
Representative Elementary Volume ◽  
P Wave ◽  
Elementary Volume ◽  
Data Set ◽  
Wave Velocities ◽  
P Wave Velocity ◽  
Ultrasonic Measurements ◽  
Seismic Scale

<p>Linking ultrasonic measurements made on samples, with sonic logs and seismic subsurface data, is a key challenge for the understanding of carbonate reservoirs. To deal with this problem, we investigate the elastic properties of dry lacustrine carbonates. At one study site, we perform a seismic refraction survey (100 Hz), as well as sonic (54 kHz) and ultrasonic (250 kHz) measurements directly on outcrop and ultrasonic measurements on samples (500 kHz). By comparing the median of each data set, we show that the P wave velocity decreases from laboratory to seismic scale. Nevertheless, the median of the sonic measurements acquired on outcrop surfaces seems to fit with the seismic data, meaning that sonic acquisition may be representative of seismic scale. To explain the variations due to upscaling, we relate the concept of representative elementary volume with the wavelength of each scale of study. Indeed, with upscaling, the wavelength varies from millimetric to pluri-metric. This change of scale allows us to conclude that the behavior of P wave velocity is due to different geological features (matrix porosity, cracks, and fractures) related to the different wavelengths used. Based on effective medium theory, we quantify the pore aspect ratio at sample scale and the crack/fracture density at outcrop and seismic scales using a multiscale representative elementary volume concept. Results show that the matrix porosity that controls the ultrasonic P wave velocities is progressively lost with upscaling, implying that crack and fracture porosity impacts sonic and seismic P wave velocities, a result of paramount importance for seismic interpretation based on deterministic approaches.</p><p>Bailly, C., Fortin, J., Adelinet, M., & Hamon, Y. (2019). Upscaling of elastic properties in carbonates: A modeling approach based on a multiscale geophysical data set. Journal of Geophysical Research: Solid Earth, 124. https://doi.org/10.1029/2019JB018391</p>

scholarly journals Simulation of the effect of stress-induced anisotropy on borehole compressional wave propagation

Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D205-D216 ◽  
Author(s):  
Xinding Fang ◽  
Michael C. Fehler ◽  
Arthur Cheng
Keyword(s):  
Wave Propagation ◽  
Stress Concentration ◽  
Elastic Properties ◽  
Wave Velocity ◽  
Wave Amplitude ◽  
P Wave ◽  
Azimuthal Variation ◽  
P Wave Velocity ◽  
Velocity Region

Formation elastic properties near a borehole may be altered from their original state due to the stress concentration around the borehole. This can lead to an incorrect estimation of formation elastic properties measured from sonic logs. Previous work has focused on estimating the elastic properties of the formation surrounding a borehole under anisotropic stress loading. We studied the effect of borehole stress concentration on sonic logging in a moderately consolidated Berea sandstone using a two-step approach. First, we used an iterative approach, which combines a rock-physics model and a finite-element method, to calculate the stress-dependent elastic properties of the rock around a borehole subjected to an anisotropic stress loading. Second, we used the anisotropic elastic model obtained from the first step and a finite-difference method to simulate the acoustic response of the borehole. Although we neglected the effects of rock failure and stress-induced crack opening, our modeling results provided important insights into the characteristics of borehole P-wave propagation when anisotropic in situ stresses are present. Our simulation results were consistent with the published laboratory measurements, which indicate that azimuthal variation of the P-wave velocity around a borehole subjected to uniaxial loading is not a simple cosine function. However, on field scale, the azimuthal variation in P-wave velocity might not be apparent at conventional logging frequencies. We found that the low-velocity region along the wellbore acts as an acoustic focusing zone that substantially enhances the P-wave amplitude, whereas the high-velocity region caused by the stress concentration near the borehole results in a significantly reduced P-wave amplitude. This results in strong azimuthal variation of P-wave amplitude, which may be used to infer the in situ stress state.

A predictive anisotropic rock-physics model for estimating elastic rock properties of unconventional shale reservoirs

2019 ◽  
Vol 38 (5) ◽  
pp. 358-365 ◽  
Author(s):  
Colin M. Sayers ◽  
Sagnik Dasgupta
Keyword(s):  
Accurate Determination ◽  
Rock Physics ◽  
Seismic Inversion ◽  
P Wave ◽  
Rock Properties ◽  
Anisotropic Rock ◽  
Rock Classification ◽  
Shale Reservoirs ◽  
Physics Model ◽  
Rock Physics Model

This paper presents a predictive rock-physics model for unconventional shale reservoirs based on an extended Maxwell scheme. This model accounts for intrinsic anisotropy of rock matrix and heterogeneities and shape-induced anisotropy arising because the dimensions of kerogen inclusions and pores are larger parallel to the bedding plane than perpendicular to this plane. The model relates the results of seismic amplitude variation with offset inversion, such as P- and S-impedance, to the composition of the rock and enables identification of rock classes such as calcareous, argillaceous, siliceous, and mixed shales. This allows the choice of locations with the best potential for economic production of hydrocarbons. While this can be done using well data, prestack inversion of seismic P-wave data allows identification of the best locations before the wells are drilled. The results clearly show the ambiguity in rock classification obtained using poststack inversion of P-wave seismic data and demonstrate the need for prestack seismic inversion. The model provides estimates of formation anisotropy, as required for accurate determination of P- and S-impedance, and shows that anisotropy is a function not only of clay content but also other components of the rock as well as the aspect ratio of kerogen and pores. Estimates of minimum horizontal stress based on the model demonstrate the need to identify rock class and estimate anisotropy to determine the location of any stress barriers that may inhibit hydraulic fracture growth.

Elastic properties of rock salt: Laboratory measurements and Gulf of Mexico well-log analysis

Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. D303-D317 ◽  
Author(s):  
Jingjing Zong ◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
Michael T. Myers
Keyword(s):  
Gulf Of Mexico ◽  
Elastic Properties ◽  
Rock Salt ◽  
Model Building ◽  
Complex Geometry ◽  
Rock Physics ◽  
P Wave ◽  
Elastic Behavior ◽  
Laboratory Measurements ◽  
Wave Velocities

Rock salt (essentially halite) is a special type of sedimentary rock that has played a large role throughout tectonic and economic history. The unique physical properties of halite (ductility, low density, flowability, and impermeability) can be critical factors in hydrocarbon traps and underground storage. However, seismic imaging and interpretation can be challenging when salt structures are present due to their complex geometry and large impedance contrasts relative to surrounding rocks. To investigate the properties of rock salt in terms of elastic parameters, we use ultrasonic laboratory measurements and well logs. In the laboratory, we have analyzed the effects of composition, crystalline structure, pressure, and temperature on the elastic behavior of a variety of rock salt samples. The samples include pure halite (>95 wt%) from the Gulf of Mexico (GOM) area, argillaceous rock salt from the Zipaquirá Mine, Colombia, and crystalline salt from the Goderich Mine, Canada. Current measurements suggest that the GOM salt cores behave isotropically in general. The Zipaquirá salt samples show velocity and density variations on account of their heterogeneous composition. The Goderich halite crystals display distinct cubic anisotropy. Measurements on the GOM samples at varying confining pressures and temperatures indicate that increasing pressure elevates velocity whereas increasing temperature decreases velocity. From the analysis of 145 log suites from boreholes drilled through rock salt in the northern GOM, we found that, within the salt formations, P-wave velocities increased slightly with depth (approximately [Formula: see text] per km). The S-wave velocities from three wells range from 2280 to [Formula: see text]. Bulk densities from all the wells cluster at [Formula: see text]. These laboratory and log measurements provide new values for the elastic properties of rock salt, which can assist in velocity model building, synthetic seismogram generation, and the understanding of the rock physics of halite.

Petroacoustics and composition of presalt rocks from Santos Basin

2019 ◽  
Vol 38 (5) ◽  
pp. 342-348 ◽  
Author(s):  
Guilherme Fernandes Vasquez ◽  
Marcio José Morschbacher ◽  
Camila Wense Dias dos Anjos ◽  
Yaro Moisés Parisek Silva ◽  
Vanessa Madrucci ◽  
...  
Keyword(s):  
Clay Minerals ◽  
Chemical Elements ◽  
Pore Space ◽  
Rock Physics ◽  
P Wave ◽  
Rock Properties ◽  
Tectonic Activity ◽  
S Wave ◽  
Wave Velocities ◽  
Lacustrine Carbonate

The deposition of the presalt section from Santos Basin began when Gondwana started to break up and South America and Africa were separating. Initial synrift carbonate deposits affected by relatively severe tectonic activity evolved to a lacustrine carbonate environment during the later stages of basin formation. Although the reservoirs are composed of carbonate rocks, the occurrence of faults and the intense colocation of igneous rocks served as a source of chemical elements uncommon in typical carbonate environments. Consequently, beyond the presence of different facies with complex textures and pore geometries, the presalt reservoir rocks present marked compositional and microstructural variability. Therefore, rock-physics modeling is used to understand and interpret the extensive laboratory measurements of P-wave velocities, S-wave velocities, and density that we have undertaken on the presalt carbonate cores from Santos Basin. We show that quartz and exotic clay minerals (such as stevensite and other magnesium-rich clay minerals), which have different values of elastic moduli and Poisson's ratio as compared to calcite and dolomite, may introduce noticeable “Poisson's reflectivity anomalies” on prestack seismic data. Moreover, although the authors concentrate their attention on composition, it will become clear that pore-space geometry also may influence seismic rock properties of presalt carbonate reservoirs.

Estimation of sea‐floor wave velocities and density from pressure and particle velocity by AVO analysis

Geophysics ◽  
1995 ◽  
Vol 60 (5) ◽  
pp. 1575-1578 ◽  
Author(s):  
Lasse Amundsen ◽  
Arne Reitan
Keyword(s):  
Numerical Study ◽  
Offshore Structures ◽  
P Wave ◽  
Rock Properties ◽  
Sea Floor ◽  
S Waves ◽  
Interface Waves ◽  
Wave Velocities ◽  
Sea Bottom ◽  
Sea Bed

Sea‐bottom properties play an important role in fields as diverse as underwater acoustics, earthquake and geotechnical engineering, and marine geophysics. Water‐column acousticians study shear and interface waves in the nearbottom sediments with the aim of inferring sea‐bed geoacoustic parameters for predicting reflection and absorption of waves at the sea floor. On the other hand, geotechnical engineers working on design and siting of offshore structures focus on these waves to characterize soil and rock properties. In the field of geophysics, sea‐bottom parameters are of interest for several reasons. In conventional marine acquisition, these parameters determine the partitioning of the incident P‐wave energy from the source into transmitted P‐waves and mode‐converted S‐waves (Tatham and Goolsbee, 1984; Kim and Seriff, 1992). The sea‐floor P‐ and S‐wave velocities and density are also necessary inputs for decomposing multicomponent sea‐floor data into P‐ and S‐waves (Amundsen and Reitan, 1995a and b), as well as in the numerical study of wave propagation phenomena.

Laboratory ultrasonic measurements: Shear transducers for compressional waves

2019 ◽  
Vol 38 (5) ◽  
pp. 392-399 ◽  
Author(s):  
Alexey Yurikov ◽  
Nazanin Nourifard ◽  
Marina Pervukhina ◽  
Maxim Lebedev
Keyword(s):  
Elastic Properties ◽  
Elastic Waves ◽  
High Energy ◽  
P Wave ◽  
P Waves ◽  
S Waves ◽  
Wave Velocities ◽  
Dipole Structure ◽  
Ultrasonic Measurements ◽  
P Wave Velocities

The ultrasonic measurements technique is well established to measure the elastic properties of rocks in the laboratory for seismic and well-log data interpretation. The key components of every laboratory ultrasonic setup are piezoelectric transducers, which generate and register elastic waves in rock samples. The elastic properties of rocks are determined through the velocities of elastic waves, which are measured by the times of the waves' travel from the source to the receiver transducer. Transducers can be specifically designed to generate P-waves (P-transducers) or S-waves (S-transducers). In limited studies, the measurement of P-wave velocities with S-transducers is mentioned. Such measurement is possible due to specific aspects of the operation of S-transducers. Namely, S-transducers are known to emit parasitic low-energy P-waves, which travel faster than high-energy S-waves and hence can be registered. However, no justification or elaboration of this method of measuring P-wave velocities was reported. To fill this gap, we first compare P-wave velocities measured with S-transducers against P-wave velocities measured with P-transducers in different rocks and materials. We show that the discrepancy between velocities measured with the two methods in homogeneous materials is less than 1% and can be up to 4% for natural rocks. Second, we numerically simulate the operation of S-transducers, show that parasitic P-waves have a dipole structure, and explain how the receiver transducer can register this compressional dipole. Finally, we use laser doppler interferometry to measure the displacement of the free surface of a sample caused by elastic waves emitted by the source S-transducer. We observed the dipole structure of the sample's surface displacement upon P-wave arrival on the surface.

An Adaptive Stratified Joint PP and PS AVA Inversion Using Accurate Jacobian Matrix

Geophysics ◽  
2021 ◽  
pp. 1-145
Author(s):  
Xiaobo Liu ◽  
Jingyi Chen ◽  
Jing Zeng ◽  
Fuping Liu ◽  
Handong Huang ◽  
...  
Keyword(s):  
Seismic Data ◽  
Jacobian Matrix ◽  
Incidence Angle ◽  
P Wave ◽  
Rock Properties ◽  
Inversion Method ◽  
S Wave ◽  
Zoeppritz Equations ◽  
Wave Velocities ◽  
Strong Contrast

Amplitude variation with incidence angle (AVA) analysis is an essential tool for discriminating lithology in the hydrocarbon reservoirs. Compared with the traditional AVA inversion using only P-wave information, joint AVA inversion using PP and PS seismic data provides better estimation of rock properties (e.g., density, P- and S-wave velocities). At present, the most used AVA inversions depend on the approximations of Zoeppritz equations (e.g., Shuey and Aki-Richards approximations), which are not suitable for formations with strong contrast interfaces and seismic data with large incidence angles. Based on the previous derivation of accurate Jacobian matrix, we find that the sign of each partial derivative of reflection coefficient with respect to P-, S-wave velocities and density changes across the interface, represents good indicator for the reflection interfaces. Accordingly, we propose an adaptive stratified joint PP and PS AVA inversion using the accurate Jacobian matrix that can automatically obtain the layer information and can be further used as a constraint in the inversion of in-layer rock properties (density, P- and S-wave velocities). Due to the use of the exact Zoeppritz equations and accurate Jacobian matrix, this proposed inversion method is more accurate than traditional AVA inversion methods, has higher computational efficiency and can be applied to seismic wide-angle reflection data or seismic data acquired for formations with strong contrast interfaces. The model study shows that this proposed inversion method works better than the classical Shuey and Aki-Richards approximations at estimating reflection interfaces and in-layer rock properties. It also works well in handling a part of the complex Marmousi 2 model and real seismic data.

scholarly journals On The Low-Frequency Elastic Response of Pierre Shale During Temperature Cycles

2021 ◽  
Author(s):  
Stian Rørheim ◽  
Andreas Bauer ◽  
Rune M Holt
Keyword(s):  
Wave Velocity ◽  
Low Frequency ◽  
P Wave ◽  
Rock Properties ◽  
Fluid Viscosity ◽  
S Wave ◽  
Temperature Cycles ◽  
Wave Velocities ◽  
Young’S Moduli ◽  
P Wave Velocity

Summary The impact of temperature on elastic rock properties is less-studied and thus less-understood than that of pressure and stress. Thermal effects on dispersion are experimentally observed herein from seismic to ultrasonic frequencies: Young’s moduli and Poisson’s ratios plus P- and S-wave velocities are determined by forced-oscillation (FO) from 1 to 144 Hz and by pulse-transmission (PT) at 500 kHz. Despite being the dominant sedimentary rock type, shales receive less experimental attention than sandstones and carbonates. To our knowledge, no other FO studies on shale at above ambient temperatures exist. Temperature fluctuations are enforced by two temperature cycles from 20 via 40 to 60○C and vice versa. Measured rock properties are initially irreversible but become reversible with increasing number of heating and cooling segments. Rock property-sensitivity to temperature is likewise reduced. It is revealed that dispersion shifts towards higher frequencies with increasing temperature (reversible if decreased), Young’s moduli and P-wave velocity moduli and P-wave velocity maxima occur at 40○C for frequencies below 56 Hz, and S-wave velocities remain unchanged with temperature (if the first heating segment is neglected) at seismic frequencies. In comparison, ultrasonic P- and S-wave velocities are found to decrease with increasing temperatures. Behavioural differences between seismic and ultrasonic properties are attributed to decreasing fluid viscosity with temperature. We hypothesize that our ultrasonic recordings coincide with the transition-phase separating the low- and high-frequency regimes while our seismic recordings are within the low-frequency regime.

Sign in / Sign up

Export Citation Format

Share Document