Multiscale Investigation of Fluid Transport in Gas Shales Rob Heller and Mark Zoback
Multiscale Fluid Flow Process Control Production July 5 July 6 Valko and Lee, 1 Production Rate July 4
Hypotheses: 3 Production Rate 1 Hypothesis 1: desorption and diffusion may be responsible for the flat production tails characteristic of many gas shale reservoirs.
Hypotheses: 4 Production Rate 1 Hypothesis 1: desorption and diffusion may be responsible for the flat production tails characteristic of many gas shale reservoirs. Hypothesis : achieving a large, conductive percolation area is essential for production.
1 Adsorption and Gas Shales 5 Conceptual model for flow in gas shales Desorption From Internal Surfaces Flow Through Intact Matrix Flow Through Fracture Network Adsorbed Gas Pore Pressure Adsorption is the pressure dependent attraction of gas molecules to the surface of a solid, resulting in a dense phase of gas at the surface The release of adsorbed gas is pressure dependent Rocks rich in minerals with a high surface area (organic matter and clay) tend to adsorb more
Adsorption Measurement Methodology 6 Adsorption Measurement Details All tests done at 4 C Boyle s Law gas expansion method: 1. First measure porosity of sample using nonadsorbing gas (Helium). Then expand adsorbing gas into sample and calculate adsorption amount based on measured pressure deficit 3. Repeat procedure at range of pressures Sample Preparation Samples from Barnett, Montney, Eagle Ford and Marcellus shale reservoirs Crushed samples (1 5 μm) Dried under vacuum at 5 C
Methane Adsorption Results: All Samples 7 7 Absolute Adsorption (SCF/ton) 6 5 4 3 1 Barnett Montney Marcellus Eagle Ford 5.3% TOC, 37.4% Clay.% TOC, 3.7% Clay 1.% TOC, 51.4% Clay 1.8% TOC, 4.9% Clay 5 1 15 Pressure (psia) Samples with higher TOC and clay adsorb considerably more.
Very little production comes from adsorbed gas 8 35 Barnett 31 35 Marcellus Gas Produced (scf/ton) 3 5 15 1 5 Total Gas Produced Free Gas Produced Adsorbed Gas Produced Gas Produced (scf/ton) 3 5 15 1 5 Total Gas Produced Free Gas Produced Adsorbed Gas Produced Gas Produced (scf/ton) 35 3 5 15 1 1 3 4 5 Pressure (psia) 5 Eagle Ford 17 Total Gas Produced Free Gas Produced Adsorbed Gas Produced 1 3 4 5 Pressure (psia) Gas Produced (scf/ton) 35 3 5 15 1 1 3 4 5 Pressure (psia) 5 Montney Total Gas Produced Free Gas Produced Adsorbed Gas Produced 1 3 4 5 Pressure (psia)
Matrix Flow in Gas Shales 9 Conceptual model for flow in gas shales Permeability Desorption From Internal Surfaces Flow Through Intact Matrix Flow Through Fracture Network Pore Pressure. How does matrix permeability evolve during production? Stress effects
Matrix Flow in Gas Shales 1 Conceptual model for flow in gas shales Permeability Desorption From Internal Surfaces Flow Through Intact Matrix Flow Through Fracture Network Pore Pressure. How does matrix permeability evolve during production? Stress effects Flow regime effects
Matrix Permeability Influenced by: 11 Stress Effects Flow Regime Effects Zoback and Byerlee, 1975 Klinkenberg, 1941 Permeability (md) Permeability (md) Effective Stress, C p P p (bars) Reciprocal Pore Pressure, 1/P p (bars 1 ) At a given pore pressure, permeability decreases with confining stress eff k ( eff C p ) xp p Apparent increase in permeability at low pore pressure due to gas slippage b k ka k 1 p
Under what conditions is slip flow important? 1 As pore pressure decreases, distance between molecular collisions (mean free path) increases Diffusion (slip-flow) becomes increasingly more important
Permeability System Setup 13 Gas Cylinder QX-6 Pump Hydrostatic Pressure Vessel Pressure Generator
Sequence of Confining Pressure and Pore Pressure Steps 5 4.5 Ln(P up -P down ) 14 4 Pressure (psia) P Upstream Natural Log P up -P down 3.5 3 Slope = α P Downstream 195 4 6 8 1 Time (hr).5 4 6 8 1 Time (hr) Pressure on one side of sample increased Downstream pressure monitored as pulse travels through Natural log of P linear in time Permeability calculated from slope t P( t) P e ka V L down Brace, 1968
Samples Tested 15 Sample Photo Perm Range (nd) Orientation Eagle Ford 17 15-35 Horizontal Eagle Ford 174 5-9 Horizontal Marcellus -18 Vertical All measured made with Helium data points/day ~6 points per experiment ~1 month/sample Measurement repeatability was monitored
Permeability vs. Cp Pp 16 6 Marcellus 5 Eagle Ford 174 13 Eagle Ford 17 Pp=1 Pp= Pp=3 Pp=4 Permeability (nd) 5 4 3 Permeability (nd) 4 3 1 Permeability (ud) 1 11 1 9 8 7 4 6 8 Effective Stress: Cp-Pp (psi) 4 6 8 Effective Stress: Cp-Pp (psi) 6 4 6 8 Effective Stress: Cp-Pp (psi) Permeability decreases with increasing simple effective stress. Can we describe this behavior with an effective stress law?
Permeability vs. Cp χpp 17 6 Marcellus 5 Eagle Ford 174 13 Eagle Ford 17 Pp=1 Pp= Pp=3 Pp=4 Permeability (nd) 5 4 3 χ =.15 Permeability (nd) 4 3 1 χ =.4 Permeability (ud) 1 11 1 9 8 7 χ = =.6 4 6 8 Effective Stress: Cp-.15*Pp (psi) 4 6 8 Effective Stress: Cp-.4*Pp (psi) 6 4 6 8 Effective Stress: Cp-.6*Pp (psi) For each rock, successfully fit all measurements to a single trend So far, we have been able to explain all permeability variation with stress effects What about lower pore pressures?
Permeability vs. Cp χpp 18 Marcellus Vertical 1 Eagle Ford 174 Eagle Ford 17 Pp=1 Pp= Pp=3 Pp=4 Permeability (nd) 15 1 5 Permeability (nd) 8 6 4 Permeability (ud) 18 16 14 1 1 8 4 6 8 Effective Stress: Cp-.15*Pp (psi) 4 6 8 Effective Stress: Cp-.4*Pp (psi) 6 4 6 8 Effective Stress: Cp-.6*Pp (psi)
Permeability vs. Cp χpp 19 Pp=5 Pp=5 Pp=75 Pp=1 Pp= Pp=3 Pp=4 Permeability (nd) 15 1 5 Marcellus Vertical Permeability (nd) 1 8 6 4 Eagle Ford 174 Permeability (ud) 18 16 14 1 1 8 Eagle Ford 17 4 6 8 Effective Stress: Cp-.15*Pp (psi) 4 6 Effective Stress: Cp-.4*Pp (psi) 8 6 4 6 8 Effective Stress: Cp-.6*Pp (psi)
Permeability vs. Cp χpp Pp=5 Pp=5 Pp=75 Pp=1 Pp= Pp=3 Pp=4 Permeability (nd) 15 1 5 Marcellus Vertical Permeability (nd) 1 8 6 4 Eagle Ford 174 Permeability (ud) 18 16 14 1 1 8 Eagle Ford 17 4 6 8 Effective Stress: Cp-.15*Pp (psi) 4 6 Effective Stress: Cp-.4*Pp (psi) 8 6 4 6 8 Effective Stress: Cp-.6*Pp (psi) eff = psi 16 Marcellus Vertical 7 Eagle Ford 174 18 Eagle Ford 17 eff =3 psi eff =4 psi Permeability (nd) 14 1 1 8 K =531.4 psi b K =854.7 psi b Permeability (nd) 6 5 4 3 K =688. psi b K =916.7 psi b Permeability (ud) 16 14 1 K =195.9 psi b K =.8 psi b 6 K b =1811.3 psi K b =1176.3 psi 1 K b =7. psi 4 1 3 4 1 1 3 4 8 1 3 4 1/Pore Pressure (psi -1 ) x 1-3 1/Pore Pressure (psi -1 ) x 1-3 1/Pore Pressure (psi -1 ) x 1-3
1 Assumptions: Total flow is sum of viscous (Darcy) flow and Knudsen/slip flow Slit shaped pore geometry Model viscous flow using Poiseuille equation P L w RT M P c P L w Q 3 1 4 3 4 1 P L ka Q 1 4 16 1 1 M RT wp c A w k 1 16 M RT w c k b 1 16 M RT k c w b P b k k k a 1 Effective Pore Size from Klinkenberg Slope
Effective Pore Size vs. Effective Stress 14 1 Eagle Ford Pore size 1- nm Effective Pore Width (nm) 1 8 6 4 Eagle Ford 17 Marcellus Eagle Ford 174 5 3 35 4 Effective Stress (psi) Pore width decreases with increasing effective stress Pore widths range from -4nm in Marcellus samples, ~13nm in Eagle Ford Klinkenberg pore sizes consistent with SEM images Pore size 1 s of nm Kohli and Zoback, 1 Image from Sondergeld, 1
To what extent does diffusion contribute to total flow? Eagle Ford 174 Marcellus Vertical 3 3 Diffusive Flux/Darcy Flux (-).5 1.5 1.5 eff = psi eff =3 psi eff =4 psi Diffusive Flux/Darcy Flux (-) 1.5 1.5 eff = psi eff =3 psi eff =4 psi Diffusive Flux/Darcy Flux (-) 3 4 5 6 7 8 Pore Pressure (psi) Eagle Ford 17.7.6.5.4.3..1 eff = psi eff =3 psi eff =4 psi 3 4 5 6 7 8 Pore Pressure (psi) 3 4 5 6 7 8 Pore Pressure (psi) Diffusive flow contributes appreciably to total flow at pore pressures < 8 psi Diffusive flow is sometimes more important than Darcy flow at pore pressure < 5 psi As we increase effective stress for a given pore pressure, we narrow the pore aperture and the relative contribution of diffusion decreases
Matrix Flow Conclusions 4 Gas slippage seems to enhance permeability at low pore pressure Effective pore widths are estimated to be 1-15nm, consistent with SEM images At low pore pressures, Knudsen diffusion (or slip flow ) becomes increasingly more important, in some cases surpassing Darcy flow Heller, R.J., Vermylen, J.P., & Zoback, M.D. (13). Experimental Investigation of Matrix Permeability of Gas Shales. AAPG Bulletin, In Press. Heller, R.J, and Zoback, M.D. (14). Adsorption of Methane and Carbon Dioxide on Gas Shale and Pure Mineral Samples. In Review. Heller, R. J., & Zoback, M. D. (11, June). Adsorption, Swelling And Viscous Creep of Synthetic Clay Samples. In 45th US Rock Mechanics/Geomechanics Symposium. Heller, R., & Zoback, M. (13, August). Laboratory Measurements of Matrix Permeability and Slippage Enhanced Permeability in Gas Shales. In Unconventional Resources Technology Conference.