An anisotropic viscoplasticity model for shale based on layered microstructure homogenization

TL;DR

This paper introduces a layered microstructure-based anisotropic viscoplasticity model for shale, capturing its time-dependent deformation behavior by homogenizing soft and hard micro-layers, validated with laboratory and simulation results.

Contribution

A novel two-scale homogenization approach deriving anisotropic viscoplastic behavior from layered microstructure, integrating a Modified Cam-Clay model for soft layers and elastic hard layers.

Findings

Model accurately predicts laboratory creep data.

Simulation demonstrates anisotropic borehole closure behavior.

Microstructure-based approach captures layered shale deformation.

Abstract

Viscoplastic deformation of shale is frequently observed in many subsurface applications. Many studies have suggested that this viscoplastic behavior is anisotropic---specifically, transversely isotropic---and closely linked to the layered composite structure at the microscale. In this work, we develop a two-scale constitutive model for shale in which anisotropic viscoplastic behavior naturally emerges from semi-analytical homogenization of a bi-layer microstructure. The microstructure is modeled as a composite of soft layers, representing a ductile matrix formed by clay and organics, and hard layers, corresponding to a brittle matrix composed of stiff minerals. This layered microstructure renders the macroscopic behavior anisotropic, even when the individual layers are modeled with isotropic constitutive laws. Using a common correlation between clay and organic content and magnitude of creep, we apply a viscoplastic Modified Cam-Clay plasticity model to the soft layers, while treating the hard layers as a linear elastic material to minimize the number of calibration parameters. We then describe the implementation of the proposed model in a standard material update subroutine. The model is validated with laboratory creep data on samples from three gas shale formations. We also demonstrate the computational behavior of the proposed model through simulation of time-dependent borehole closure in a shale formation with different bedding plane directions.