Comments
Transcript
Evidence for a variable Archie porosity exponent “ m
Evidence for a variable Archie porosity exponent “m” and impact on saturation calculations for Mesaverde tight gas sandstones: Piceance, Uinta, Green River, Wind River, and Powder River basins Robert M Cluff, The Discovery Group Inc. Alan P. Byrnes, Kansas Geological Survey1 1presently with Chesapeake Energy Stefani Whittaker, The Discovery Group Inc. Dan Krygowski, The Discovery Group Inc. AAPG Rocky Mountain Section meeting, Denver, Colorado 10 July 2008 US DOE Project Summary z DOE Contract # DE-FC26-05NT42660 z z completion date 30 June 2008 Organizations z z z University of Kansas Center for Research, Inc. Kansas Geological Survey, Lawrence, KS The Discovery Group Inc., Denver, CO z Principal Investigators: Alan P. Byrnes, KGS; Bob Cluff, Discovery Group z project website is http://www.kgs.ku.edu/mesaverde Objectives of this task z Characterize Mesaverde electrical properties as a function of porosity and salinity z Archie porosity (cementation) exponent “m” z Investigate behavior at low porosity end (<6%) not previously studied z Evaluate excess conductivity effects z Methods to compute accurate Sw from logs Sampling z z z z z systematic characterization of Kmv lithofacies over entire Rocky Mtn region 44 wells/6 basins Described ~7000 ~ ft core (digital) 2200 core samples 120-400 advanced properties samples Powder River Wind River Wyoming Green River N Washakie Utah Colorado Uinta Piceance Permeability vs Porosity z Samples collected over a wide range of porosity and permeability across 6 basins 0-24% porosity, spanning 1 nD to >100 mD 1000 Klinkenberg Permeability (4,000 psi, mD) z 100 10 1 0.1 Green River Piceance Powder River Uintah Washakie Wind River logK=0.3Phi-3.7 logK=0.3Phi-5.7 0.01 0.001 0.0001 0.00001 0.000001 0.0000001 0 2 4 6 8 10 12 14 In situ calc Porosity (%) 16 18 20 22 24 Archie’s equation Sw = (a / φ n z z * (Rw / Rt ) completely empirical – no theoretical basis “m” is the porosity or cementation exponent z z m) loosely related to tortuosity of the current flow path, better thought of as electrical efficiency of the path “n” is the saturation exponent z related to change in conductivity path with changing saturation Archie porosity exponent z for a simple bundle of capillary tubes oriented parallel to current flow direction: m→1 z insensitive to cross section shape, so fractures act like capillary tubes as porosity increases and less of it participates in the conductive path, m ↑ z for an “average” sandstone comprised of spherical grains, m → 2 z Resistivity of a simple rock model 0 ∞ For straight capillary tubes: Porosity (φ) Resistivity (Ro) F = Ro/Rw = 1/φ For rock with tortuous pores: F = Ro/Rw = 1/φ m (after Doveton, 2005) 1 Rw Capillary tube model for m m 1.0 >1 ~2 >2 m=1 after Herrick & Kennedy, 1993, SPWLA Paper HH A Φ Rw L high -P core holder electric insulator then we measure R0. F = R0/Rw. brine in Resistance Reference Cell Frequency Generator brine of known Rw & φ Core Plug Start with core plug saturated with high-P fluid Micropipette Core measurement of the formation factor, F Plated electrodes When F and φ are plotted log-log 1000 m= 2 m= 3 100 m= 1 F 10 1 0.01 0.1 φ 1 log F = -m log φ Observed porosity dependence of “m” R2 = 0.63 (RMA) 40Kppm dataset is largest and used for base case cap m at 1.95 40K ppm brine data 2.2 2.2 2.1 2.0 In situ Archie Cementaiton Exponent Empirical: m = 0.676 log φ + 1.22 each salinity is different In situ Archie m y = 0.5377x + 1.3313 R2 = 0.6331 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.0 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 log In situ Porosity (%) 1.4 2.0 1.8 1.6 1.4 1.2 1.0 0 4 8 12 16 In situ Porosity (%) 20 24 Dual porosity model z m = log[(φ-φ2)m1 + φ2m2]/log φ z φ expressed as V/V φ2 = 0.0035, m1=2, m2=1; SE both = 0.11 z rock behaves like a mixture of matrix porosity and cracks, with cracks dominating low porosity end cap at m = 1.95 (φ ~ 16%) 40K ppm brine data both models fit data 2.2 z z φ = bulk porosity φ2 = fracture porosity m1 = matrix cementation exponent m2 = fracture cementation exponent In situ Archie Cementaiton Exponent z 2.0 1.8 1.6 1.4 1.2 1.0 0 4 8 12 16 In situ Porosity (%) 20 24 And a third way to look at it.... z z z z z Why is the minimum m ~ 1.2, instead of 1? A – for a distribution of cracks of different crosssectional area, the largest (widest) cracks will dominate the conductivity The high tail of the distribution determines the bulk conductivity, while the rest of the cracks act like “excess” porosity that do not participate (significantly) in the conductivity. Therefore m ↑ And are the “cracks” all fractures? z z z Probably not.............. Slot-like pores oriented preferentially parallel to bedding also act like conductive cracks Thin parallel laminae of slightly coarser, more permeable sand will be crack-like Salinity dependence of “m” tested plugs with 20K, 40K, 80K, and 200K ppm brines Nearly all cores exhibit some salinity dependence 2.3 n=335 0.9 Core Conductivity (mho/m) In situ Archie Cementation Exponent, (m, A=1) 1.0 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 2.2 2.1 2.0 1.9 1.8 1.7 1.6 1.5 1.4 1.3 1.2 1.1 1.0 0 2 4 6 8 10 12 14 16 18 Brine Conductivity (mho/m) 20 22 0.01 0.1 Brine Resistivity (ohm-m) 1 All data, all salinities Archie Cementaiton Exponent (m, a=1) 2.40 2.20 2.00 1.80 1.60 1.40 1.20 200K 80K 1.00 40K 20K 0.80 0 2 4 6 8 10 12 14 16 18 20 22 In situ Porosity (% ) Salinity dependence of “m” z 20K ppm y = 0.2267Ln(x) + 2.2979 2 R = 0.6619 z 1.50 Series1 Log. (Series1) 1.00 40K ppm z 3.00 0.50 y = 0.2328Ln(x) + 2.409 2 R = 0.6547 2.50 0.050 0.100 0.150 0.200 0.250 2.00 insitu porosity (%) Series1 1.50 Log. (Series1) 80K ppm 1.00 0.50 3.00 0.00 0.000 2.50 y = 0.2149Ln(x) + 2.4354 2 0.050 0.100 0.150 0.200 R = 0.5132 0.250 insitu porosity (%) 2.00 Axis Title 0.00 0.000 Axis Title Axis Title 2.00 200K ppm Series1 1.50 Log. (Series1) 1.00 3.00 0.50 2.50 y = 0.1621Ln(x) + 2.3222 0.00 0.000 2 R = 0.3633 2.00 0.050 0.100 0.150 insitu porosity (%) 0.200 0.250 Axis Title 2.50 m = a log φ + b intercept b drops with decreasing salinity slope is ~ constant Series1 1.50 Log. (Series1) 1.00 0.50 0.00 0.000 0.050 0.100 0.150 insitu porosity (%) 0.200 0.250 Simple procedure to compute Sw z determine Rw @ Tf conventionally Pickett plots – focus on the lower porosity, wetter sandstones z produced waters z your best guess....... z z convert Rw to 75°F by chart lookup or Arps equation Pickett Plot example Rw = 0.306 pick m at low porosity end, where BVWirr ~ BVW Williams PA 424-34 Piceance basin Kmv above “top gas” Pickett plot Rw 0.306 ohmm @ 160°F = 0.7 @ 75°F (9K ppm) Our new procedure z compute m at 40K ppm from RMA regression: m40k = 0.676 log φ + 1.22 e.g. for 10% φ : m = 0.676 + 1.22 = 1.896 z correct m for salinity effect by m = m40k + ((0.0118 φ – 0.355) * (log Rw + 0.758)) z e.g. for 10% φ, Rw = 0.7 @ 75°F m = 1.896 + ((0.0118 * 10 – 0.355) * (log 0.7 + 0.758)) m = 1.896 + (-0.237 * 0.603) = 1.753 z cap m at 1.95 (~12% porosity) Practical impact Nominally, most of us use an m close to 2, but usually slightly less, for tight gas sand evaluations (e.g. 1.85, 1.90) z Variable m that DECREASES with decreasing porosity leads to lower Sw’s z Therefore, there is more gas in the tight rocks than we thought. z Above 10% porosity there is very little difference z Example: Low porosity, wet zone Moderate porosity, wet “High” porosity gas zone m is HIGHER than base case, so Sw is higher 20Kppm example, Natural Buttes 30K ppm example, Wamsutter Summary & Conclusions 335 Kmv samples run at multiple salinities z Archie porosity exponent m varies with z porosity z salinity z z m ↓ as porosity ↓ m ↓ as salinity ↓ behavior is consistent with increasing electrical efficiency with decreasing porosity, whatever the pore scale architecture variable m model can be implemented with a simple equation relating m to porosity and formation water salinity z m is constant above ~12% porosity at 1.95 z lowering m at 5-12% φ increases GIP z see no impact below ~5% porosity z BVWirr is typically 3-5% z no longer calculate Sw’s >> 1 z Sw = 1 at low φ validates Rw z Visit our project website http://www.kgs.ku.edu/mesaverde Questions?