Robust and accurate formulation for modeling of acid stimulation

Abstract

Accurate representation of processes associated with energy extraction from subsurface formations often requires models which account for chemical interactions between different species in the presence of multiphase flow. In this study, we focus on modeling of acid stimulation in the near-well region. For the chemical processes which include a dissolution of rock material, an issue arises with the predictive representation of flow. Taking into account the spatial scale of discretization, some of simulation control volumes can have values of porosity close to 1, which makes an application of Darcy's law inconsistent and requires employing a true momentum equation such as the Darcy-Brinkman-Stokes (DBS) equation. The DBS equation automatically switches the description between Darcy equation in control volumes with low porosity and Stokes equation in grid blocks with high porosity. For chemical reactions, we propose a local nonlinear solution technique that allows solving the balance of solid species separately yet retaining the full coupling with rest of the equations. Finally, we study the impact of multiphase flow. The DBS approach is not well established for multiphase flow description. Therefore we employ a hybrid approach, where we assume that the single-phase DBS flow and the multiphase Darcy flow occur in separate regions. We test the accuracy and performance of both approaches on realistic models of practical interest.