FORMATION EVALUATION PETE 663

Transcription

1 FORMATION EVALUATION PETE 663 SHALY SAND EVALUATION - B Summer 2010

2 SHALY FORMATION ISSUES Water Saturation

3 Obstacles with Shaly Sands Fresh formation water can cause conventional analysis to overestimate water saturation Salty formation water can cause low resistivity, meaning pay zones can be bypassed Thin beds may lead conventional log analysis to underestimate porosity and overestimate water saturation

4 SHALY FORMATION ISSUES LECTURE A Shales/clays have several origins and forms Shales/clays affect: Porosity Permeability Vshale Estimations Assumptions Log responses LECTURE B Shales conduct electricity Problems with Archie-based methods Rwa problem Sw errors

5 WELL LOG EFFECTS - 1 Well X Water leg 150 ft Shaly interval ft. Resist. increase Sonic Δ t increase Density ρ b increase

6 HC zone 150 ft Shaly interval ft. WELL LOG EFFECTS - 2 Resist. decrease Sonic Δt increase Density ρ b increase Neutron φ N increase

7 Without shale R o = FR w C o 1 or F C o = C w F C w

8 With shale C C = w + o F X C o X 1 F C w The factor X is the excess conductivity caused by the fact that the clays and shales are conductors of current

9 ELECTRICAL CURRENT FLOWING THROUGH BRINE CLEAN SANDSTONE Pore Current Pore Throat Grain Brine

10 SCHEMATIC - ELECTRICAL CURRENT FLOWING THROUGH DISPERSED CLAY COATING AND BRINE, SHALY SANDSTONE What is the effect of dispersed clay on resistivity when pores are filled with brine? Grain Brine Clay Coat Total Conductivity =

11 ELECTRICAL CURRENT FLOWING THROUGH DISPERSED CLAY COATING, SHALY SANDSTONE What is the effect of dispersed clay on resistivity when the porefilling fluid is oil rather than saline water? Electrical Current Grain Oil Clay Coat

12 SILICATE TETRAHEDRON SiO2 SILICATE MINERALS Silicates are the most abundant minerals Basic building block is the silicate tetrahedron (SiO4) Oblique View Si, + 4 O -2 Map View Modified from Grim, 1968

13 MONTMORILLONITE STRUCTURE n H2O & Mg, Na, Ca Modified from Halliburton, EL-1007 n H2O & Mg, Na, Ca

14 Montmorillonite

15 MUSCOVITE STRUCTURE (Similar to Illite) Electrical Current Aluminum replaces silicon Mobile cations; Includes Helmholtz Planes From Grim, 1968

16 Illite

17 WHAT IS SHALE? Clay + silt + other Clays Plate-like form Large surface area Contain Al +3 and Si +4 Substitution by Mg +2 Negative charge results Attraction by water and cations Clay Crystal x H Absorbed Water Sodium Ion Hydration Water Water Outer Helmholtz Plane Schematic Water Molecule

18 DIFFERENT MODELS OF DIFFUSE LAYER Clay Crystal Cl - Na + Ionic Concentration In NaCl Solution Saline Water 0 X d Distance From Clay Surface Gouy profile of diffuse layer, thickness X d = ( n) increases as salinity decreases. Clay Crystal x H Absorbed Water Sodium Ion Hydration Water Saline Water (X H / 6.18 A) O H + H + Model of exclusion layer (Helmholtz Plane)sodium ions excluded from surface layer by dielectric properties of water

19 SPECIFIC SURFACE AREAS OF SOME MINERALS Mineral Sand Kaolinite Illite Montmorollinite Ft^2/ft^ thousand 15.2 million 85.4 million 274 million Clays have extremely large surface areas Surface area varies greatly among clay minerals

20 SURFACE AREA vs CEC

21 WATER SATURATION SHALY SANDS

22 Shaly Formation Issues Water Saturation Shales/clays have several origins and forms Vshale Estimation Assumptions Log responses Shales/clays affect formation Permeability Porosity Shales conduct electricity Problems with Archie-based methods Rwa problem Sw errors

23 WATER SATURATION EQUATIONS Many different water saturation equations have been developed Archie s model for a clean formation is: n w S = FR R t All other models are for shaly formations where the rock is not a perfect insulator w

24 With shale C C = w + o F X C o X 1 F C w The factor X is the excess conductivity caused by the fact that the clays and shales are conductors of current

25 Commonly used formulas to account for shale: Vsh Models Simandeaux (better with saline fm water) Indonesia (developed for fresher fm water) Double-layer model (Attempt to avoid using Vsh) Waxman-Smits Dual water

26 CO vs CW in Shaly and Clean Sands Non-Linear Region (Indonesian Equation) Linear Region (Simandoux Equation) CO Slope = 1 / F CW Modified from Halliburton, EL-1007

27 Cw/Co Variation with Cw and Vcl Clean Sand, F = aφ -m CW CO CW Modified from Halliburton, EL-1007

28 Archie n w m e w t S A C C φ = ( ) sh sh n w sh n w m e w t C V S V A S C C = φ Simandeaux (better with saline water) All V sh models are similar: total C = clean C + shale C Indonesia (better with fresh water) 2 / 2) / ( 1 2 / n w sh V sh n w w t S C V S F C C sh + =

29 Waxman-Smits model C t = C ' w φ m t A S n wt Where: C ' = C + w w BQ S v wt Note independent conduction paths by free water and bound water

30 New terms BQ V : conductivity of bound water Q V : cation exchange capacity (meq/gm dry clay) 1 meq = 6E20 atoms measures how many cations are present different clays have different CEC s kaolinite 0.03 to 0.06 chlorite 0 to 0.1 illite 0.1 to 0.4 montmorillonite 0.8 to 1.5 B is: specific counterion conductivity (mho/m per meq/cc) Counterions are the charge-balancing Na cations B is a per unit measure Measures how effectively cations conduct electricity

31 Waxman-Smits Swt obtained by iteration S wti+ 1 = 1 R w + ( BQ S ) wti where Swt 0 is the initial guess, Swt 1 is the next guess, etc., and F B Note: R w in B equation is at 75 F. F R t A = φ m t [ ( 0. 5 R )] = e w B max v 1 n

32 Maximum Equivalent Conductance of Sodium-Exchange Ions, λ NA or B max vs Temperature λ NA or B max, mho - cm -2 mca Temperature, C

33 Graphing values of B max vs. temperature ( R and/or F) on various types of graph paper, one finds that B max vs. log (T R) is more or less linear: B max = (51.31)ln T R ( )

34 Approximate values for Q v are: Very shaly Q v =1.5 Moderate shale Q v =1.0 Medium shale Q v =0.5 Low shale Q v =0.25 No shale Q v =0 CEC or Q v should, however, be lab measured Q v may correlate with logs (e.g., GR)

35 Dual water model C t = C w m t φ A S n wt where C w = S S b wt C wb + 1 S S The dual-water model is a more general form of the Waxman-Smits model. The free water salinity can be different than the salinity of the claybound water. To determine Sw, use iterative method, like W-S b wt C wf

36 Case study 1. Pennsylvanian Cottage Grove Sandstone The Cotton Grove is identified as the interval from 6542 to Sample shows calcareous shaly sandstone. Calculate Sw at 6566 and 6680 ft. Conventional Analysis Rw = Tf φ e = 2 D φ + φ 2 2 N Shaly sand F = 0.81/ф2 Rmf = 0.47@Tf φ D = ρ ρ ma ma ρ ρ b fl ρf = 1.0 gm/cc ρma = 2.68 gm/cc

37 At depths 6566 and

38 Calculate Sw at 6566 and 6680 ft. Conventional Analysis (Archie s method) gives the results below Sw Results are high (60%), which could make one cautious about developing the well Depth φ d φ n φ xp ρ b φ d Sw

39 sh sh app corr V φ φ φ = Ncorr Dcorr corr φ φ φ + = ρ ρ ρ ρ φ = ma fl ma b D MIN MAX MIN SH GR GR GR GR I = 1) 0.33(2 2* = sh I V sh Shaly Sand Analysis - Steps = SH SH SH SH t w e e w we R V R V R R R S * 5*.4* 0 φ φ 5. Simandoux Equation

40 Using Simandoux Equation for shaly sand analysis we have: Depth φ d ρ b GR I GR V cl Depth φ nc φ dc φ e Swe Comparison of Sw values: with shaly sand analysis Sw = without shaly sand analysis Sw =

41 Case Study 2. Permian Basin, Spraberry Sandstone, Midland Basin Deep Spraberry sandstone was encountered at a depth of 7720 to The formation is not very clean, as is evident from the log sandstone

42 Log Analysis of Depths 7724,7732,

43 Conventional log analysis produced the following results: ρ f = 1.0 gm/cc ρ ma = 2.68 gm/cc Depth φ d φ n φ xp ρ b φ d Rt Sw The saturations obtained are 0.50 to 0.66, which is fairly high

44 Shaly sand analysis produced the following: Depth ρ b φ d GR LOG I GR V SH Depth φ nc φ dc φ e Swe Comparison of saturation values: Without shaly sand analysis - range With shaly sand analysis - range

45 Detection of Secondary Porosity

46 Vertical, Mineralized Fracture: 1U Payzone Shackelford 1-38A

47 Mineralized Fracture: 1U Payzone Shackelford 1-38A

48 Shackelford 1-38A (1-U in the Upper Spraberry) water saturation with different m & n compared with measured water saturation from whole core analysis. Sharp contrast between pay and non-pay is observed, by fluorescence, at a depth of 7092 ft. 1 Water saturation (Sw) Pay Sw (a=0.81, m=2, n=2) Sw (a=1, m=1.66, n=1.46) Sw (core) Non-pay Depth, ft.

49 Fractured Zone Identification 0.3 Fractured Zone (pay) Non-fractured zone (Non-pay) 0.2 Porosity (neutron) Porosity 0.1 Porosity (sonic) Porosity (core) Depth, ft.

50 Fracture Detection

51 Sponge Core - 5U Zone O Daniel #37

52 FMI O Daniel #37 5U Zone

53 Detection of Secondary Porosity

54 SHALY SAND ANALYSIS - INDONESIA EQUATION WELL X EXAMPLE

55 WELL X EXAMPLE Shale layer : 0-12 feet Clean layer: Approx from feet Required : Water saturation at 225 feet and 47 feet using the Indonesia equation The matrix is sandstone

56 INDONESIA EQUATION C t = C F w S n / 2 1 ( V / 2) n / 2 w + V sh sh C sh S w

57 Porosity Estimation Using Vsh 1 φ corr = φ app V φ sh sh Effective porosity = φ corr Apparent porosity, matrix adjusted = φ app Apparent porosity in shale = φ sh Example...

58 SHALY SAND ANALYSIS - INDONESIA EQUATION WELL X EXAMPLE Depth Rhob Ø D Ø NLS Ø NSS GR Vsh Ø D (corr) Ø N (corr) Rt Sw

59 SHALY SAND ANALYSIS - INDONESIA EQUATION WELL X EXAMPLE Depth ρ b φ D φ NLS φ NSS CGR V SH φ D corr φ N corr R t S w % % %

60 New terms W-S and Dual Water models depend on CEC Without CEC, have to use nearby shale S b - bound water saturation S b = f(cec, C wf ) S b =V sh φ sh /φ t C wb - bound water conductivity C wb = g(cec, S b )

61 SUMMARY Clays are conductive complex resistivity responses, Sw determination Two types of Sw models for shaly sand Vsh models (e.g., Simandeaux) Double-layer models Vsh models Empirical Assume shale properties same as nearby shale All shales have same effect Double-layer models Do not use Vsh Use electrical properties of clays (CEC)