Abstract

Ubiquitous observations of channelised fluid flow in the form of pipes or chimney-like features in sedimentary sequences provide strong evidence for significant transient permeability-generation in the subsurface. Understanding the mechanisms and dynamics for spontaneous flow localisation into fluid conductive chimneys is vital for natural fluid migration and anthropogenic fluid and gas operations, and in waste sequestration. Yet no model exists that can predict how, when, or where these conduits form. Here we propose a physical mechanism and show that pipes and chimneys can form spontaneously through hydro-mechanical coupling between fluid flow and solid deformation. By resolving both fluid flow and shear deformation of the matrix in three dimensions, we predict fluid flux and matrix stress distribution over time. The pipes constitute efficient fluid pathways with permeability enhancement exceeding three orders of magnitude. We find that in essentially impermeable shale (10−19 m2), vertical fluid migration rates in the high-permeability pipes or chimneys approach rates expected in permeable sandstones (10−15 m2). This previously unidentified fluid focusing mechanism bridges the gap between observations and established conceptual models for overcoming and destroying assumed impermeable barriers. This mechanism therefore has a profound impact on assessing the evolution of leakage pathways in natural gas emissions, for reliable risk assessment for long-term subsurface waste storage, or CO2 sequestration.

Highlights

  • The buoyant pore-fluid triggers local decompaction of the porous medium and enables upward-migration within self-organised chimneys. This contrasts with Darcian flow models in non-deforming porous media that predict diffusive fluid flow and spreading of fluids

  • The second invariant of the deviatoric stress tensor quantifies the magnitude of shear deformation recorded by the porous matrix (Fig. 4d)

  • Our findings suggest that the dominant mechanisms responsible for the spontaneous formation of fluid escape pipes in the subsurface are viscous creep of the porous matrix, decompaction weakening, and hydro-mechanical coupling

Read more

Summary

Introduction

The hydro-mechanical simulation initial configuration consists of a rectangular box with a dimensionless extent of 30 × 30 × 60 in the x, y (horizontal) and z (depth) directions, respectively. Description Buk viscosity Permeability Fluid viscosity Chimney width Propagation speed Symbol ηφ κφ μf. Permeability values within the source region were nine times higher compared to normalised background values of 1 (Fig. 6). The computational domain was subjected to the downward-pointing gravity field and affected by a background horizontal strike-slip shear-deformation of similar magnitude than buoyancy forces. The mechanical problem was solved using free-slip (no shear stress) boundary conditions on all sides of the box. For the fluid flow problem, we applied no flux boundary conditions on all vertical sides of the box, and fixed flux value at the bottom and top boundaries to satisfy the condition pe = 0 (no compaction or decompaction). The dimensionless parameters used in the code (Table 1) allowed us to optimally converge the numerical simulation

Methods
Results
Conclusion
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call