The paper describes a model to predict and analyze the time-dependent closure of a borehole drilled in a soft rock subjected to an initial isotropic stress, and the build-up of stress on a rigid casing after contact with the deforming rock. The rock behaves as a viscoplastic material characterized by a Mohr–Coulomb yield criterion and plastic potential, and by a time-dependent stress–strain deviatoric response akin to a Bingham rheology. The model formulation recognizes the particular structure of the solution, namely that the borehole is encircled by an evolving viscoplastic annulus, itself surrounded by an infinite domain, where the rock is either elastic or is unloading elastically. Noting that the solution outside the viscoplastic boundary is given explicitly by the Lamé solution, the evolution problem is formulated as an initial boundary value problem in the viscoplastic region only, but with a free boundary — the growing or shrinking interface between the viscoplastic and the elastic domains. The equations governing the mechanical fields and the evolution of viscoplastic boundary are spatially discretized on a moving mesh with a fixed number of nodes. The final system of equations is a set of first order ODEs that are efficiently solved using the MATLAB routine ODE45. The numerical simulations reveal (i) the time-dependent deformation and stress state of rock in the viscoplastic annulus before and after the contact, (ii) the time of first contact between the deforming rock and the casing, (iii) the effective duration of the stress build-up on the casing, and (iv) the large time contact pressure, which is well approximated by the elastoplastic solution.
Bibliographical noteFunding Information:
This work was sponsored by The Research Council of Norway, Aker BP, BP, ConocoPhillips, Equinor, Shell, and Total, through the PetroMaks 2 KPN-project: Shale as a Permanent Barrier after Well Abandonment, at SINTEF Industry (Grant no. 244420/E30 ).
- Borehole stability