A fully implicit and consistent finite element framework for modeling reservoir compaction with large deformation and nonlinear flow model. Part II: verification and numerical example

A fully implicit, fully coupled, and fully consistent finite element framework has been formulated in part I of this work for modeling reservoir compaction through linearizing coupled solid and flow field equations and constructing a local material integrator. In part II of this work, we focus on verification and performance analysis of our numerical formulation and computer implementation using several numerical examples. First, we design a cube problem in triaxial compression to verify our numerical formulation and computer code implementation especially for rock formation in compaction using cap plasticity models. The finite element prediction on stresses is compared with the analytical solution. The second problem we select is a strip footing problem popular in the geotechnical area where the evolution of soil consolidation degrees following the diffusion of pore pressure is the main interest. In this example, we demonstrate a good performance of the proposed numerical formulation on solving different shear and compaction-dominated deformation behaviors by varying the footing length. Importantly, an extremely sharp cap model based on real experimental data for Leda clays, a challenging cap model, is successfully applied in this footing problem. Our focus in this work is to model field reservoirs undergoing serious compaction. A reservoir with complex payzone geometries, multiple horizontal wells, and cap plasticity models with sharp cap surfaces has been successfully solved using our fully implicit formulation. The last example is to model a horizontal wellbore damage problem. Finally, the sensitivity of predicted subsidence to nonlinear flow model, cap hardening parameters, and Lode angles have been systemically investigated and documented in detail, which can provide a constructive guidance on how to successfully model field reservoir compaction problems with cap plasticity models.