A New Explicit Sequentially Coupled Technique for Chemo-Thermo-Poromechanical Modelling and Simulation in Shale Formations

Lateef T. Akanji (), Adamu Ibrahim, Hossein Hamidi, Stephan Matthai and Alfred Akisanya
Additional contact information
Lateef T. Akanji: School of Engineering, University of Aberdeen, Aberdeen AB24 3FX, UK
Adamu Ibrahim: School of Engineering, University of Aberdeen, Aberdeen AB24 3FX, UK
Hossein Hamidi: School of Engineering, University of Aberdeen, Aberdeen AB24 3FX, UK
Stephan Matthai: Peter Cook Centre for CCS, The University of Melbourne, Parkville, VIC 3010, Australia
Alfred Akisanya: School of Engineering, University of Aberdeen, Aberdeen AB24 3FX, UK

Energies, 2023, vol. 16, issue 3, 1-28

Abstract: A new explicit sequentially coupled technique for chemo-thermo-poromechanical problems in shale formations is developed. Simultaneously solving the flow and geomechanics equations in a single step is computationally expensive with consequent limitations on the computations involving well or reservoir-scale geometries. The newly developed solution sequence involves solving the temperature field within the porous system. This is followed by the computation of the chemical activity constrained by the previously computed temperature field. The pore pressure is then computed by coupling the pore thermal and chemical effects but without consideration of the volumetric strains. The geomechanical effect of the volumetric strain, stress tensors, and associated displacement vectors on the pore fluid is subsequently computed explicitly in a single-step post-processing operation. By increasing the borehole pressure to 20 MPa, it is observed that the rock displacement and velocities concurrently increase by 50%. However, increasing the wellbore temperature and chemical activities shows only a slight effect on the rock and pore fluid. In the chemo-thermo-poroelasticity steady-state simulation, the maximum displacements recorded in the H m i n and H m a x are 0.00633 m and 0.0035 m, respectively, for 2D and 0.21 for the 3D simulation. In the transient simulation, the displacement values are observed to increase gradually over time with a corresponding decrease in the maximum pore-fluid velocity. A comparison of this work and the partial two-way coupling scheme in a commercial simulator for the 2D test cases was carried out. The maximum differences between the computed temperatures, displacement values, and fluid velocities are 0.33 % , 0.7 % , and 0 % , respectively. The analysed results, therefore, indicate that this technique is comparatively accurate and more computationally efficient than running a full or partial two-way coupling scheme.

Keywords: coupled flow and geomechanics; wellbore stability; hydraulic fracturing; geomechanics; finite element; COMSOL multiphysics; CSMP++ (search for similar items in EconPapers)
JEL-codes: Q Q0 Q4 Q40 Q41 Q42 Q43 Q47 Q48 Q49 (search for similar items in EconPapers)
Date: 2023
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/1996-1073/16/3/1543/pdf (application/pdf)
https://www.mdpi.com/1996-1073/16/3/1543/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jeners:v:16:y:2023:i:3:p:1543-:d:1057405

Access Statistics for this article

Energies is currently edited by Ms. Agatha Cao

More articles in Energies from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().