The Elastic Stiffness C13 in Transversely Isotropic Shales
Abstract
Shales are often anisotropic due to fine-scale layering and preferred alignment of clay particles with bedding. They may be approximated as transversely isotropic (TI) media with symmetry axis (x3) perpendicular to bedding. TI media have five independent elastic stiffnesses C11, C33, C55, C66 and C13, with C12 = C11−2C66. While C11, C33, C55 and C66 may be estimated from velocities measured parallel and perpendicular to x3, C13 requires measurements at oblique angles, resulting in greater uncertainty. Despite this, C13 is important in seismic velocity analysis, imaging, amplitude variation with offset (AVO), and geomechanical applications including wellbore stability, stress estimation, hydraulic fracture design, and reservoir compaction. In isotropic media C13 = C33 − 2C55 and C13 = C12 , but shale microstructure may break these equalities due to preferred orientation of clay particles, fine-scale layering, and bedding-parallel microcracks. For a TI orientation distribution (ODF) of clay particles, the elastic properties may be characterized by two coefficients W200 and W400 in an expansion of the ODF in generalized Legendre functions. C13 may be greater than C12 if W200 is significantly larger than W400, whereas C13 may be lower than C33 − 2C55 if W200 is significantly smaller than W400. Fine-scale layering also influences C13. Depending on layer properties, sequences with C12 less than C13 and C13 less than C33 − 2C55 may occur. Bedding-parallel microcracks typically increase C12 − C13 and C13 −C33 +2C55 , but C13 −C33 +2C55 may decrease if the ratio of normal-to-shear crack compliance is small, as expected for isolated fluid-filled cracks. These observations provide insight into the physical controls on C13 and improve understanding of its variability in shale formations.
Related
- Vanadium Occurrence and Enrichment Mechanisms in Lower-Middle Carboniferous Black Shales of the Central Taurus Belt, Türkiye: Insights into Critical Element Potential
- Separation of Source, Site, and Path Effects Based on Elastic Radiative Transfer Theory. Part 1: Mapping of Scattering and Absorption Properties and Application to Metropolitan France
- Stiffness Formula and Conditioning Analysis for General Brown–Korringa Fluid Substitution
- Separation of Source, Site, and Path Effects Based on Elastic Radiative Transfer Theory, Part 2: Generalized Inversion of Source and Site Parameters and Application to Metropolitan France
- Upper Limits on the Isotropic Gravitational-Wave Background from the first part of LIGO, Virgo, and KAGRA's fourth Observing Run
BibTeX
@misc{sayers2026elastic,
title = {{The Elastic Stiffness C13 in Transversely Isotropic Shales}},
author = {Colin M Sayers},
howpublished = {EarthArXiv},
year = {2026},
month = sep,
doi = {10.31223/x5tv3m},
url = {https://eartharxiv.org/repository/view/14768/},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - The Elastic Stiffness C13 in Transversely Isotropic Shales
AU - Colin M Sayers
PY - 2026
DA - 2026/09/03
PB - EarthArXiv
DO - 10.31223/x5tv3m
UR - https://eartharxiv.org/repository/view/14768/
AB - Shales are often anisotropic due to fine-scale layering and preferred alignment of clay particles with bedding. They may be approximated as transversely isotropic (TI) media with symmetry axis (x3) perpendicular to bedding. TI media have five independent elastic stiffnesses C11, C33, C55, C66 and C13, with C12 = C11−2C66. While C11, C33, C55 and C66 may be estimated from velocities measured parallel and perpendicular to x3, C13 requires measurements at oblique angles, resulting in greater uncertainty. Despite this, C13 is important in seismic velocity analysis, imaging, amplitude variation with offset (AVO), and geomechanical applications including wellbore stability, stress estimation, hydraulic fracture design, and reservoir compaction. In isotropic media C13 = C33 − 2C55 and C13 = C12 , but shale microstructure may break these equalities due to preferred orientation of clay particles, fine-scale layering, and bedding-parallel microcracks. For a TI orientation distribution (ODF) of clay particles, the elastic properties may be characterized by two coefficients W200 and W400 in an expansion of the ODF in generalized Legendre functions. C13 may be greater than C12 if W200 is significantly larger than W400, whereas C13 may be lower than C33 − 2C55 if W200 is significantly smaller than W400. Fine-scale layering also influences C13. Depending on layer properties, sequences with C12 less than C13 and C13 less than C33 − 2C55 may occur. Bedding-parallel microcracks typically increase C12 − C13 and C13 −C33 +2C55 , but C13 −C33 +2C55 may decrease if the ratio of normal-to-shear crack compliance is small, as expected for isolated fluid-filled cracks. These observations provide insight into the physical controls on C13 and improve understanding of its variability in shale formations.
ER -