Partial Saturation Fluid Substitution with Partial Saturation

Transcription

1 Fluid Substitution with 261

2 5 4.5 ρ fluid S w ρ w + S o ρ o + S g ρ g Vp (km/s) K fluid S w K w + S o K o + S g K g Patchy Saturation Drainage Fine-scale mixing 1 = S w + S o + S g K fluid K w K o K g Sw (fraction) Imbibition Knight and Nolen-Hoeksema (GRL, 1990) found saturation hysteresis at ultrasonic frequencies. We know now that velocities depend, not just on saturation, but also on the scales at which the phases are mixed. The curve labeled imbibition is typical when phases are mixed at a fine scale. The curve labeled drainage is typical when the phases are mixed at a coarse scale -- which we call patchy. K.1 262

3 Vp (km/s) sandstone porosity = 30% patchy homogeneous Oil Saturation 263

4 Increasing water saturation a. b. c. d. Endres and Knight (The Log Analyst, 1989) modeled different microdistributions of pore fluids and gas in the stiff and soft portions of the pore space. They concluded that the scale and distribution of fluids influence velocities. K.4 264

5 Patchy Saturation High frequency response Isolated patches some stiff, some softer Overall a higher effective velocity Very low frequency response Gassmann behavior with a single "effective fluid" 1 K eff.fl = Σ i S i K i Overall the softest, lowest velocity Critical frequency: f visc = kk f l 2 η 265

6 Patchy Diffusion Scales Characteristic diffusion time for a pressure disturbance with length scale L to relax: 1 f = τ L2 4D Inverting this, we can find the characteristic diffusion length over which pressure differences can relax at seismic frequency f L 4κK fl ηf D is the hydraulic diffusivity, K is the fluid bulk modulus, κ is the permeability, and η is the viscosity. f 10 Hz 1000 md 100 md 10 md 1 md 1000 md 100 md 100 Hz 10 md 1 md 10 5 Hz κ 1000 md 100 md 10 md 1 md.1 md L 1 m.3 m.1 m.03 m.3 m.1 m.03 m.01 m.01 m.003 m.001 m.0003 m.0001 m 266

7 K.5 Thierry Cadoret studied velocity vs. saturation using the resonant bar and found the coarse-scale and fine-scale behavior. 267

8 Estaillades Limestone Cadoret s velocity and attenuation vs. saturation. The fine scale distribution gives relaxed viscoelastic behavior, and the coarse scale gives unrelaxed. Therefore, we expect the largest attenuation when the velocity dispersion is largest. Hence, we get the important result that P-wave attenuation in a partially saturated rock can be much larger than in the dry or fully saturated case. 268

9 The problem that we address is the nonunique response of seismic velocity to fluid saturation. What are the physical conditions that cause patchy behavior? When do we use the patchy model and when do we use the homogeneous model? Our approach is to use flow simulation to study the parameters that control fluid distributions at a fine scale. Vp (km/s) sandstone porosity = 30% patchy homogeneous Oil Saturation study by Madhumita Sengupta and G. Mavko 269

10 Typical Eclipse Cell 10m x 10m Our cells: 1m x 1m ~ L c Our approach is simply to run flow simulations at a fine scale to help discover the reservoir and fluid properties that control the saturation scale. The fine scales are chosen to be approximately the critical diffusion length, so that any mix of fluids within the cell can be represented as a Reuss average fluid. 270

11 Porosity and permeability models for flow simulation 271

12 We will consider two important cases: water flood into oil, and gas flood into oil. The parameters that we consider are: relative permeability wettability density contrast permeability heterogeneity capillary pressure 272

13 Water Injection in Oil Relative Permeability Curves for Oil and Water 1 Rel Perm oil oil water water Sw 273

14 Saturations obtained from flow simulations using the dashed (top) and solid (bottom) relative permeability curves. The irreducible saturations are critical controls on the saturation extremes. Sw Sw Sw Sw 274

15 The patchy and uniform saturation curves are upper and lower bounds. They describe the range of velocity signatures that we can achieve by mixing the end members. Finite irreducible saturations drastically narrow the range of uncertainty. 275

16 Wettability The saturation distribution depends on the wettability of the rock. Most sandstone reservoirs are water wet and most carbonate reservoirs are oil wet. Water saturation map and histogram of saturation in an oil wet rock 276

17 Wettability Vp Oil Wet (Drainage) * * * **** Water Wet (Imbibition) Sw Wettability tends to determine whether the velocities fall high or low in the allowable range. 277

18 Mobility Ratios MR = k max wr η w / 2500 k or max η o Low Mobility Ratio High Mobility Ratio

19 V P V p 3150 Low MR Low MR High MR Water Saturation S w High MR 279

20 The uniform saturation model is good enough for waterflood oil-water cases. Exceptions: when the irreducible oil is very low in an oil wet rock. The main control is the finite irreducible saturations. 280

21 Gas Injection Into Oil Sg Sg Sg Sg Sg Sg 281

22 Gas Injection Effect of Mobility (a) Low Mobility Ratio Sg (b) High Mobility Ratio Sg Sg Sg 282

23 Vp High Mobility Ratio Low Mobility Ratio So 283

24 Heterogeneity of Perm Perm Models Saturation Satn Histograms

25 Vp 2950 Large Scale Heterogeneities Small Scale Heterogeneities So 285

26 Summary of Mixing Rules K Voigt = S w K w + S o K o + S g K g K Brie = ( K liquid K g )( 1 S g ) e + K g 1/K Reuss = S w / K w + S o / K o + S g / K g Brie, et al.: SPE

27 Conclusions Reservoirs with gas are very likely to show patchy behavior. The uniform saturation model may be good enough for reservoirs with only oil and water The main mechanism that causes patchy behavior at the field scale is gravity. 287