The Dual-Reciprocity Boundary Element Analysis for Hydraulically Fractured Shale Gas Reservoirs Considering Diffusion and Sorption Kinetics

This work presents a new application of boundary element method (BEM) to model fluid transport in unconventional shale gas reservoirs with discrete hydraulic fractures considering diffusion, sorption kinetics and sorbed-phase surface diffusion. The fluid transport model consists of two governing partial differential equations (PDEs) written in terms of effective diffusivities for free and sorbed gases, respectively. Boundary integral formulations are analytically derived using the fundamental solution of the Laplace equation for the governing PDEs and Green’s second identity. The domain integrals arising due to the time-dependent function and nonlinear terms are transformed into boundary integrals employing the dual-reciprocity method. This transformation retains the domain-integral-free, boundary-integral-only character of standard BEM approaches. In the proposed solution, the free- and sorbed-gas flow in the shale matrix is solved simultaneously after coupling the fracture flow equation of free gas. Well production performance under the effect of relaxation phenomenon due to delayed responses of sorbed gas under nonequilibrium sorption condition is rigorously captured by imposing the zero-flux condition at fracture–matrix interface for the sorbed-gas transport equation. The validity of proposed solution is verified using several case studies through comparison against a commercial finite-element numerical simulator.