Subcritical crack propagation due to chemical rock weakening: macroscale chemo-plasticity and chemo-elasticity modeling

Abstract:

Crack propagation in a subcritically stressed rock subject to chemically aggressive environment is analyzed and numerically simulated. Chemically induced weakening is often encountered in hydraulic fracturing of low-permeability oil/gas reservoirs and heat reservoirs, during storage of CO2 and nuclear waste corroding canisters, and other circumstances when rock matrix acidizing is involved. Upon acidizing, mineral mass dissolution is substantially enhanced weakening the rock and causing crack propagation and eventually permeability changes in the medium.

The crack process zone is modeled mathematically via a chemo-plastic coupling and chemo-elastic coupling model. In plasticity a two-way coupling is postulated between mineral dissolution and a yield limit of rock matrix. The rate of dissolution is described by a rate law, but the mineral mass removal per unit volume is also a function of a variable internal specific surface area, which is in turn affected by the micro-cracking (treated as a plastic strain). The behavior of the rock matrix is modeled as rigid-plastic adding a chemical softening capacity to Cam-Clay model. Adopting the Extended Johnson’s approximation of processes around the crack tip, the evolution of the stress field and deformation as a function of the chemically enhanced rock damage is modeled in a simplified way. In addition, chemical reactive transport is made dependent on plastic strain representing micro-cracking. Depending on mechanical and chemical boundary conditions, the area of enhanced chemical softening is near or somewhat away from the crack tip.

In elasticity, chemo-mechanical effect is postulated via a chemical volumetric shrinkage strain proportional to mass removal variable, conceived analogously to thermal expansion. Two versions are considered: of constant coefficient of shrinkage and a variable one, coupled to deviatoric strain. Airy Potential approach used for linear elasticity is extended considering an extra term, which is uncoupled or coupled to strain. The later case requires iterations with solution of reactive transport equation. A decrease of stress intensity factor with time of reaction is well reproduced.