The elastic properties and anisotropy of artificial compacted clay samples

Geophysics ◽  
2021 ◽  
Vol 86 (1) ◽  
pp. MR1-MR15
Author(s):  
Fei Gong ◽  
Bangrang Di ◽  
Lianbo Zeng ◽  
Jianxin Wei ◽  
Jiwei Cheng ◽  
...  

Clay minerals are a major component of hydrocarbon reservoir rocks, and they are known to play important roles in the physical and elastic properties of rocks. However, it is difficult to directly measure these properties of single-crystal clays due to their small particle size. Therefore, we have constructed three sets of artificial clay samples with different compaction stresses to investigate the effect of the compaction stress and clay mineralogy on their elastic properties and anisotropy. All of the dry samples are measured by the pulse-transmission method. The results indicate that the compaction stress and clay mineralogy have a significant influence on the physical and elastic properties of the clay samples. The microstructures of clay samples indicate that the clay platelets are aligned almost perpendicularly to the direction of compaction stress, and the ultrasonic velocity analysis validates the assumption of transverse isotropy of our clay samples. The velocities increase with the compaction stress, especially at low stress, which corresponds to the rapid porosity reduction at low stress levels. Velocity anisotropy parameters increase with increasing of compaction stress due to the increase of texture sharpness for clay minerals during the compaction process. The elastic moduli of the clay samples display a significant stress sensitivity and a strong directional dependence, with the Young’s moduli increasing and the Poisson’s ratios decreasing with the compaction stress. A simple theoretical template is used to quantify the orientation distribution functions (ODFs) of clay platelets, and the generalized Legendre coefficients of ODF increase with the increase of compaction stress, especially at low stress. Further, the compressional-wave (P-wave) and shear-wave anisotropy increase with the ODF coefficients [Formula: see text] and [Formula: see text], especially P-wave anisotropy.

Geophysics ◽  
1997 ◽  
Vol 62 (2) ◽  
pp. 521-532 ◽  
Author(s):  
Lev Vernik ◽  
Xingzhou Liu

Using ultrasonic velocity and anisotropy measurements on a variety of shales with different clay and kerogen content, clay mineralogy, and porosity at a wide range of effective pressure, we find that elastic anisotropy of shales increases substantially with compaction. The effect is attributed to both porosity reduction and smectite‐ to‐illite transformation with diagenesis. A means of kerogen content mapping using velocity versus porosity crossplot for shales is shown. Matrix anisotropy of shales dramatically increases with kerogen reaching the maximum values of about 0.4 at total organic carbon (TOC)=15–20%. A strong chemical softening effect was found in shales containing even minor amounts of swelling (smectite) clay when saturated with aqueous solution. This effect results in a significant P‐wave anisotropy reduction as compared to dry and oil‐saturated shales. Since mature black shales are normally oil wet, this effect can only have a local significance restricted to the wellbore wall. Accurate measurements of phase velocities, including velocities at a 45° direction to the bedding plane, allow us to immediately calculate elastic stiffnesses and anisotropic parameters. Intrinsic (high pressure) properties of shales display an ε > δ > 0 relation. Introduction of the bedding‐parallel microcracks in overpressured shales results in a δ decrease when fully fluid saturated and a δ increase when partially gas saturated, with a characteristic effect on the shape of the P‐wave velocity surface at small angles of incidence. Filtering the contribution of the intrinsic anisotropy of shales, it is possible to estimate the pore fluid phase, microcrack density, and aspect ratio parameters using seismic anisotropy measurements.


1962 ◽  
Vol 42 (2) ◽  
pp. 296-301 ◽  
Author(s):  
J. S. Clark ◽  
J. E. Brydon ◽  
H. J. Hortie

X-ray diffraction analysis was used to identify the clay minerals present in fourteen subsoil samples that were selected to represent some more important clay-bearing deposits in British Columbia. The clay mineralogy of the subsoils varied considerably but montmorillonitic clay minerals tended to predominate in the water-laid deposits of the south and illite in the soil parent materials of the Interior Plains region of the northeastern part of the Province.


Geosciences ◽  
2020 ◽  
Vol 11 (1) ◽  
pp. 8
Author(s):  
Mei Zhang ◽  
Hongfeng Lu ◽  
Qing Chen ◽  
Gayan Bandara ◽  
Hui Zhang ◽  
...  

In the northern South China Sea, pockmarks are widely distributed on the seabed offshore on the southwestern Xisha Uplift. The mineralogy and geochemistry of the clay minerals and surface sediments from the pockmark field were identified using X-ray diffraction (XRD) analysis and X-ray fluorescence (XRF) analysis to trace the provenance, weathering, and sediment transportation system in the area. The clay minerals are primarily comprised of illite, smectite, kaolinite, and chlorite, showing a distribution of average weight percentages of 35%, 35%, 18%, and 13%, respectively. Based on the surrounding fluvial drainage basins and various transport mechanisms (current or monsoon), illite and chlorite primarily originate from rivers in Taiwan and the Mekong and Red Rivers. Kaolinite primarily originates from the Pearl River, and smectite derived from the Luzon arc system is primarily transported by surface currents with significant influence from the Kuroshio intrusion.


Author(s):  
Suresh Dande ◽  
◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
◽  
...  

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.


Geophysics ◽  
2017 ◽  
Vol 82 (5) ◽  
pp. C175-C185 ◽  
Author(s):  
Ivan Pšenčík ◽  
Véronique Farra

We have developed approximate nonhyperbolic P-wave moveout formulas applicable to weakly or moderately anisotropic media of arbitrary anisotropy symmetry and orientation. Instead of the commonly used Taylor expansion of the square of the reflection traveltime in terms of the square of the offset, we expand the square of the reflection traveltime in terms of weak-anisotropy (WA) parameters. No acoustic approximation is used. We specify the formulas designed for anisotropy of arbitrary symmetry for the transversely isotropic (TI) media with the axis of symmetry oriented arbitrarily in the 3D space. Resulting formulas depend on three P-wave WA parameters specifying the TI symmetry and two angles specifying the orientation of the axis of symmetry. Tests of the accuracy of the more accurate of the approximate formulas indicate that maximum relative errors do not exceed 0.3% or 2.5% for weak or moderate P-wave anisotropy, respectively.


2004 ◽  
Vol 5 (1) ◽  
pp. 43
Author(s):  
K.G. PEHLIVANOGLOU ◽  
G. TRONTSIOS ◽  
A. TSIRAMBIDES

The Thermaikos Gulf constitutes the NW part of the North Aegean Sea and is limited eastward from the Chalkidiki Peninsula and westward from the Pieria Prefecture. Its plateau covers an area of 3,500 km2. The mechanisms responsible for the grain size distribution into the Gulf, the clay mineralogy and the chemistry of some bottom sediments from the outer Thermaikos Gulf, are examined. Source mixing during transportation, flocculation, differential settling processes and organic matter appear to be the main mechanisms for the distribution of clay minerals in shallow waters. All grain size fractions studied present a wide range of values confirming the extreme variations of the discharged load and the variability in marine processes. Plagioclases predominate over K-feldspars, while quartz is the most abundant mineral present. In addition, micas, chlorites, amphiboles and pyroxenes exist as primary and/or accessory minerals in all samples. Among clay minerals, illite predominates over smectite and smectite over chlorite (+ kaolinite). The ordered interstratified phase of I/S, with 30-35% S layers, is present in the 2-0.25µm fraction. The randomly interstratified phase of I/S, with 50% S layers, is present in the <0.25& micro; m fraction. On average the clay mineral content of the studied samples is: 48% I, 23% S, 17% Ch (+K) and 12% others for the 2-0.25µm fraction and 50% I, 30% S and 20% Ch (+K) for the <0.25 µm fraction. All these minerals are the weathering products of the rocks from the drainage basins of the rivers flowing into the Gulf, as well as of the Neogene and Quaternary unconsolidated sediments of the surrounding coasts. The terrigenous input, the water mass circulation and, to a lesser extent, the quality of the discharged material and the differential settling of grains, control the grain size distribution within the outer Thermaikos Gulf. The chemical composition of the analysed samples is generally in agreement with their mineral composition and signifies their terrigenous origin presenting discretely clastic character.


2020 ◽  
Author(s):  
Jerome Fortin ◽  
Cedric Bailly ◽  
Mathilde Adelinet ◽  
Youri Hamon

&lt;p&gt;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.&lt;/p&gt;&lt;p&gt;Bailly, C., Fortin, J., Adelinet, M., &amp; 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&lt;/p&gt;


2021 ◽  
Vol 18 (4) ◽  
pp. 567-582
Author(s):  
Shikai Jian ◽  
Li-Yun Fu ◽  
Chenghao Cao ◽  
Tongcheng Han ◽  
Qizhen Du

Abstract Natural fractures are usually multiscale in size, orientations and distribution, resulting in complex anisotropic characteristics. Analytical methods for quantifying the associated effective elastic properties are based on some assumptions, such as dilute fracture concentration and regular-shaped fractures, which do not occur in actual reservoirs. Here, we conduct anisotropic finite-element modeling of effective elastic properties of complex fractured rocks using the least-square fitting method. The algorithm is developed for a 3D case and validated by classical effective medium theories for models with different fracture densities. The results of the 3D numerical method agree well with the theoretical predictions at low fracture density. The model also considers the interactions among fractures to calculate equivalent stiffness tensors at high fracture density. Three 2D fracture models are simulated to demonstrate the basic behavior of stress interactions and their effect on the overall elasticity under different fracture densities. We applied the developed model to 3D natural fractures built from a real outcrop, and we found that the fracturing pattern significantly affects the effective anisotropy properties. The resultant P-wave phase velocities as functions of the incidence angle and frequency are anisotropic. This study provides a great potential to calculate equivalent stiffness tensors and anisotropic properties of 3D multiscale natural fractures.


Geophysics ◽  
2014 ◽  
Vol 79 (4) ◽  
pp. D205-D216 ◽  
Author(s):  
Xinding Fang ◽  
Michael C. Fehler ◽  
Arthur Cheng

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.


Clay Minerals ◽  
2006 ◽  
Vol 41 (4) ◽  
pp. 863-877 ◽  
Author(s):  
S. Guggenheim ◽  
J. M. Adams ◽  
D. C. Bain ◽  
F. Bergaya ◽  
M. F. Brigatti ◽  
...  

Brindley et al. (1951) reported the earliest efforts to obtain international collaboration on nomenclature and classification of clay minerals, initiated at the International Soil Congress in Amsterdam in 1950. Since then, national clay groups were formed, and they proposed various changes in nomenclature at group meetings of the International Clay Conferences. Most of the national clay groups have representation on the Nomenclature Committee of the Association Internationale pour l'Etude des Argiles (AIPEA, International Association for the Study of Clays), which was established in 1966. The precursor committee to the AIPEA Nomenclature Committee was the Nomenclature Subcommittee of the Comité International pour l'Etude des Argiles (OPEA, International Committee for the Study of Clays).


Sign in / Sign up

Export Citation Format

Share Document