Physical parameter estimation for sandstone reservoirs - crewes

7 downloads 0 Views 306KB Size Report
Canadian Occidental Petroleum Ltd.modified theFORTRAN program. The first author also would like to thank R. James Brown for showing the theory explained ...
291

Physical parameter estimation for sandstone reservoirs Taiwen Chen, Donald C. Lawton and Fred Peterson

ABSTRACT

Through theoretical derivation and numerical modeling, this report shows that BiotGassman's theory can be used for predicting some physical parameters (e.g., velocity, density, Poisson's ratio, bulk modulus) as functions of porosity or water saturation. Also, it is shown that oil and gas have different Poisson's ratio curve as a function of water saturation. Several plots in this report suggest the importance of obtaining good initial velocity, porosity and dry rock Poisson's ratio estimates when modeling gas-water or oil-water saturated reservoirs. INTRODUCTION

The Biot-Gassman theory has drawn many authors' attention (Geertsma and Smit, 1961: Domenico, 1974: Gregory, 1977: Hampson and Russell, 1990). Gregory (1977) suggests that the Biot-Gassman theory is a technique through which some physical parameters (e.g., velocity, density, Poisson's ratio, bulk modulus) as functions of porosity or water saturation can be obtained. Based on this theory, some modified equations were developed in order to obtain Poisson's ratio for AVO (amplitude versus offset) studies. Comparing the Biot-Gassman's equations (Gregory, 1977) with modified ones, one may feel that the modified equations are a more direct approach to get some parameters such as Poisson's ratio. Several numerical models show the application of this theory. THEORY

The main issue of the modified Biot-Gassman theory is to obtain Poisson's ratio ((r) as a function of porosity ((_)or water saturation (Sw). Following are the list of symbols: V Vp Vs

- velocity of the rock; - P-wave velocity; - S-wave velocity;

Vp0 M Md M0

-

initial P-wave velocity (at _0 and Sw0); space modulus (or P-wave modulus); dry rock space modulus; initial space modulus;

- rigidity modulus (or S-wave modulus);

292

lad 6 ¢_d P

- dry rock rigidity modulus; - Poisson's ratio; - dry rock Poisson's ratio; - bulk density;

Po

- initial bulk density (at _o and Swo);

Ps Pf

- bulk density of solid matter;, - bulk density of fluid;

Pro

- initial bulk density of fluid (at ¢0 and Sw0);

Pw Ph

- bulk density of water;, - bulk density of hydrocarbon (oil or gas);

¢

- porosity;

¢_0 Sw Sw0 K Kd Kd0 Ks

-

initial porosity; water saturation; initial water saturation; bulk modulus; dry rock bulk modulus; initial dry rock bulk modulus; solid matter bulk modulus;

Kp Kf Kw Kh C Cd Cs Cf Cw Ch

-

pore bulk modulus; fluid bulk modulus; water bulk modulus; hydrocarbon (oil and gas) bulk modulus; compressibility (reciprocal of the bulk modulus); dry rock compressibility; solid matter compressibility; fluid compressibility; water compressibility; hydrocarbon (oil or gas) compressibility;

It is well known that a can be obtained by:

(___p)2_ 2 6=_(_--_Ps)2-1]

(e.g., Sheriff, 1991),

(1)

VP M - C

(e.g., Gregory, 1977),

(2)

Vs = "_/g-¥ P

(e.g., Sheriff, 1991),

(3)

and Vp, Vs can be obtained:

Therefore, o can be obtained if M, I.tand p are a/l known.

293

Calculation

of p

It is known that: p = (l-¢)ps + Cpf, and pf = Swpw+ (1-Sw)Ph

(e.g., Gregory, 1977),

then we have: p = (1-¢)ps + ¢SwPw + ¢(1-Sw)Ph,

(4)

Ps, Pw and Ph are known. Equation (4) tells us that bulk density p is the function of porosity ¢ and water saturation Sw. Calculation

of M

Geertsma and Smit (1961) gave the equation to obtain space modulus (M).

KdZ (1-.v-) M=Md +

Ks

(5)

¢ +I-¢ Kf

Ks

Ks

Gregory (1977) gave the equation to get Kf: Cf = SwCw + (1-Sw)Ch,

and

1 =Sw.¢ Kf Kw Kf =

1-Sw Kh '

1 Sw + 1-Sw" Kw Kh

(6)

In equation (5), Md=4_-_ + Kd

(e.g., Sheriff, 1991),

3(1-2Od)KA I_ = 2(]_a)_'_

(e.g., Sheriff, 1991),

and

then (5) becomes:

Kd 2 (I- _--_-s ) M =_d

+Kd t

¢+__1-¢ Kd Kf K_ K_

294

(1 M- 4x3 l'2_dKd + Ka I

Ks)2

-72 l+od

e 1-¢ Kd' Kf

(1 - K_?2

M=SKi•

Ks

K2

,

1-¢ Kd Kf

Ks

(7)

K_

here S = 3(1-t_d) l+_d " In equation (6), Kw, Kh are known. For example, 2.38 x 10a°dynes/cm 2 and 0.0208 x 1010

dynes/cm 2 are suggested values for bulk moduli of water(Kw), and gas, respectively

(Hilterman; Hampson and Russell, 1990). Hence, Kf is a function only of Sw. In equation (7), Ks, _d are known.

Values of 40.0 x 10l° dynes/cm 2 and 0.12 are suggested for sandstone

Ks and ad, respectively (Hilterman; Hampson and Russell, 1990). From equation (5) and (6), it is clear that space modulus M is a function of porosity _, water saturation Sw and bulk modulus of dry rock Kd. However,

Kd can be obtained by knowing Kd0. From Appendix A, we can obtain:

Cd = 0Cp + (l-O)Cs, then 1=¢ +1-¢ Ka Kp K's

Kd-

1 ¢ +--__¢ Kp

(8)

Ks

similarly, Kd0 =

I

¢o-I--1-Oo

(9)

Ks

Kp

From equation (9), we can obtain Kp.

_0 Kp=

.

Kao

Ks

(10)

295

Substituting (10) into (8): Kd=

1

1-¢"

_._..(1 ¢0 _

(11)

1-*0) ._ Ks Ks

Equation (7) can be used to solvc for Kdo. For the initial porosity _o, initial fluid hulk modulus Kf0, initial bulk modulus of the dry rock Kdo and initial space modulus Mo, equation (7) becomes:

(1 - ._--_°) 2 Mo = SKd0 t

Ks

,

(12)

¢o _1-¢o Kao Kro

Ks

Ks2

to separate the factor Kao in equation (12):

MO(K_ m + 1-¢o. __, ¢0 + 1-¢0 Ks Kd0) K_ =o°Kao___.fo Ks Kd0) K_ +(1 - Kd0)2, Ks further written as: S-_(Kdo)2

Mo S¢o S(1-0o)],_

+[MoOo+Mo(1-¢o) 1 =0,

then we have: A(I_o) 2+BI_o+C=0, where: A=S-1.

K_' B=

2 .Mo.S¢° S(1-¢o). Ks Ks2 Kf0 Ks

C = M0¢0 4 M0(1-_) -1; Kfo Ks 3 1-_d;

S = y-+b--rid Kf0 =

1

Swo + 1-Swo" Kw Kh

(13)

296

Mo = (Vpo)2po; and Po = (1-¢o)ps + ¢oSwoPw+ ¢o(1-Swo)Ph. By solving the equation (13), we can get Kd0: Kd0-

-B ___/-_- 4AC 2A

04)

Knowing Kd0 from equation (14), and by substituting (11) into (7), we can see clearly that M is only the function of porosity _ and water saturation Sw.

Calculation

of p.

Gregory (1977) gave the equation: 3(1-2(_d)v p.=gd=_X,.d.

(15)

¢_dis known, and from equation (11), Kd is the function of porosity _. Hence, p. is the function of porosity _ only. In equations (2) and (3), M, p and }.tcan be obtained by equations (7), (4)and (15), respectively.

From the above, it is clear that P-wave velocity (Vp) and S-wave velocity (Vs) are

the functions of porosity _ and water saturation Sw. By giving one value of ¢ or Sw, we can get the curves of Vp (Vs) versus Sw or _. Then from equation (1), Poisson's ratio (¢_)versus Sw or curves can be obtained. DISCUSSION

Figure 1 and Figure 2 are two examples for the application of the Biot-Gassman theory for water-gas and water-oil saturated sandstone. Parameters (most values were suggested by Hampson and Russell, 1990) are listed in Table 1. Table

Water-gas Water-oil

1:

Physical

K(x 1010 dynes/cm 2) Ks Kw Kh 40.0 2.38 0.0208 40.0 2.38 1.0

Parameters

Ps 2.65 2.65

used for Sandstone

P(g/cm3) Pw Ph 1.089 0.103 1.089 0.75

(_d

Vp0

Sw0

1_0

0.12 0.12

(m/s) 2500 3000

1.0 1.0

0.33 0.33

3000

"

7000 I_rmlty

P-wave

,

= :_1%

-

6000 _,

2000

Wat_

Saturation.

100%

_._ 4000

1500 •

1000 2500

oo

1000

012

014 Water

016

018

Lo

"

'

0.0

0.2

0.4

0.6

0.8

i.o

Saturation Porosity

04"

0.4 •

I_rmlly.

_

0

33%

0.2'



0.2-

0.3"_

i

03-

Oi

o.o

o12

_,_ MD

oi,_

olo

o18

_.o

OI

o.o

o12

Water Saturation

Figure 1: Water-gas Saturated Sandstone

o'.4

oi_

Porosity

o'.8

,.o

6000,

Por_lty

3500

5000 "

= 33%



_

P-wave

_ 2500

3000 -

2000

1500 0.0

0:2

, 0.4 Water

016

0:8

1000 0.0

1.0

_

0.22 '

0.18

'

0.0

0:4

0:6

0]8

1.0

Porosity

0.3

0.24 •

0.20

0:2

saturation

0,28 •

"_

-

012

: ,

0.4

Water

0:6

_

ha

0.2

.

o'.8

,.o

O,

oo

0:2

,

o._

Saturation

Porosity

Figure 2: Water-oil Saturated Sandstone

0:6

0:8

i.o

299

Shown in the left side diagrams of Figure 1, it is clear that a small percentage of gas in a sand has a strong effect on the Vp and a, which is consistent with Ostrander's (1984) and Hilterman's results. Oil, however, does not have the same effect. In the left side diagrams of Figure 2, Vp and ff change smoothly with the change of Sw; P increases with the increasing of Sw, thus causes the Vs to decrease with increasing Sw, which can be seen in the plots on the upper left hand side both in Figure 1 and Figure 2. Generally Vp and Vs decrease with the increasing _, increases with the increasing of _. These can be seen in the fight-hand side diagrams of both Figure 1 and Figure 2. Also, it can be noticed that o changes nonlinearly with the change of _ for small values of 0.

Effect of Initial

Estimation

of Velocity

From the theoretical derivation, it is clear that initial velocity Vp0 is one of the important input parameters for the calculation of Poisson's ratio. Therefore, it is needed to know how would the initial estimation of velocity affect the final results. Figure 3 shows the plots of three models with Vpo changed (arbitrary values) and other parameters kept constant. were suggested by Hampson and Russell, 1990) are listed below:

Parameters (most values

p(g/cm3): Pw = 1.089, Ph = 0.75, Ps = 2.65; K(x 101° dynes/cm2): Kw = 2.38, Kh = 1.0, Ks = 40.0; (_d= 0.12, Sw0 = 0.30, _0 = 0.15;

Vp0(m/s): Vpm = 4000, Vp02 = 3600, Vp03= 3200. The plot on the upper left hand side of Figure 3 represents the Vp change with the change of Sw for different Vpo, and the plot on the upper right shows the Vs change. From these two, it is clear that for a given Sw, decreasing Vp0 will decrease both Vp and Vs- t_ increases with the decreasing Vpo for a given Sw, which can be seen on the lower left hand side plot. The lower right hand plot shows that with the decreasing of Vp0, Kd decreases. compressible, in other words, less compatible. Effect of Initial Estimation

This means that the rock is more

of Porosity

Initial estimation of porosity (¢0) is another important factor for the calculation.

Figure 4

shows the plots of three models with O0 changed (arbitrary values) and other parameters

5000 3200 -

----o---VpO =4000 VpO =3600 _, 4500

Pore_ty=15%

3000-

-_

_

------o---VpO = 4000 _ VpO=3600

Por_s_ty=15%

2800 -

VpO=3200

"--'_--

VpO= 3200

2400 _

__.c_.c._e...o.__e--.e-.-c_c.r" 3000

0.0

012

'

04

'

0.6

018

20001800

1.0

oo

WaterSaturation

,

012

04

0'.6

018

10

WaterSaturation

0.34.

2OO £.o

0.30"

o---- VpO=4oO0 * VpO= 3600 , ,[7---VpO= 3200

Poro_ty = 15%

._.$

cr"

_" _

_ 180

_

160

o

o 0

VpO=4oo0 VpO= 3600 VpO 320O

0

i

0.26" 120

o.18.

"g

0,22' ,.._.,._,__

_,_

0.14

oo

012

014

016

Water Saturation

&

1.0

_

ioo 140 80

0.0

o12

o14

o16

Water Saturation

Figure 3: Effect of Initial Estimation of Velocity

o18

,0

2320 ----o---

3800 •

--

{

_

I.l_-asity=2O%



/

• LPormlty=15%

_-

_

jCj

_ 3600 3700 y

0'.6

018

_

2180 2240'! 2160" "-°" _+_d"_'_ 2140

014

.

LPereslty=15% L pe.re_ty = 20% LPorostty=10%

2200

_

012

*

2220 - _ 2260-

_

].l:'oresity:lO%

_,

_00 0.0

-_

2300 2280

_.0

0.0

_'_"a"'>"t_'t_'_

0'.2

014

0'.6

0'.8

,.0

Water Saturation

Water Saturation 0.28 0.26'

_ ----o---

L porosl_ pere*l_ == 20% _% L

,_.f,_.,,1

_

0"22I

_

I 17

_

L Porosity = 20%

--o--e --

1. por_ty = 10% L poro61ty= 15%

o

I 15

0.20' 0A8 2 0.24- _ 0.16 0,0

"

_

0',2

" 0.4

I

0.6

I

1.0

0"0

0.8

0_2

014

Water Saturation

Figure 4: Effect of Initial Estimation

0"6

Water Saturation

of Porosity

. 0"8

I'O

302

unchanged. below:

Parameters (most values were suggested by Hampson and Russell, 1990) are listed

p(g/cm3): Pw = 1.089, Ph = 0.75,Ps= 2.65; K(× 10I°dynes/cm2): Kw = 2.38,Kh = 1.0, Ks = 40.0; crd = 0.12,Swo = 0.30,Vpo = 3600m/s;

(_ot = 0.20, 002 = 0.15, (_03= 0.10. From the two upper plots, we can see that there is a crossing of the Vp curve with Sw, and for a given Sw, Vs decreases with the decreasing value of t_o. The change in Vs is greater than that of Vp, therefore for a given Sw, _ increases with the decreasing value of _0. Effect of initial estimation

of dry rock Poisson's

ratio

Figure 5 are plots of three models with (_dvarying (arbitrary values) while all other parameters remain constant. Parameters (most values were suggested by Hampson and Russell, 1990) are listed below: p(g/cm3): Pw = 1.089, Ph = 0.75, Ps = 2.65; K(x 101° dynes/cm2): Kw = 2.38, Kh = 1.0, Ks = 40.0; Swo = 0.30, Vpo = 3600 m/s, _0 = 0.15; t_d: Odl

=

0.12, t_d2 = 0.14, t_d3 = 0.16.

From the two upper plots in Figure 5, it is obvious that Vp remains almost constant with the change of ad, while Vs decreases with the increasing _d. Hence, for a given Sw, _ increases with the increasing value of od. The lower right hand side plot shows the bulk modulus increasing with the increasing t_a. This infers that the rock is more compacted (or less compressible). Assumptions

for Biot-Gassman

Theory

The Biot-Gassman theory is only used for water-gas and water-oil saturated reservoirs, but it is not valid for water-oil-gas saturatedreservoirs. Also, as Hilterman cited in his notes, the Biot-

2280.

3750

i_"

----o---•

P.R. (dry rock)= 0.12 P.R.(dryro_k)=0.14

_

_

2200.

-----o----

P.R. (dry rock) = 0.12

2240.

_



P.R. (dry rock) = 0.14 P.R. (dry rock) = 0.16

os 2200. _

3650

- -

2100.

3550

0.0

012

014

0;6

018

2140

_.0

0.0

_-_.__

012

Water Saturation

,

0._

01_

,

0.6

_.0

Water Saturation

Iz;0 0.26.

o---•

P.R. (dry rock) = 0.12 • ----c_---

135,

P.R. (dry rock) =0.14 P.l_ (dry rock) = 0.16

0.24"

_

_ R

J

to o t.o

P.R. (dry rock) = 0.12 P.R. (dry rock) = 0.14

----_--"P.R. (dry rock) = 0.16 130'

.¢_125, i_

0.22, --_

:::: :

: : : : : : : : : : : : : : ::

]15'

_

0.20. _ 0.16

0.0

"_ _ =

0:2

, 0.4

0;6

Water Saturation

0:6

_.o

_

IlO

o.o

0:2

, 0.4

0;6

Water Saturation

Figure 5: Effect of Initial Estimation of Dry Rock Poisson's Ratio

o;8

,.o

304

Gassman's theory assumes a homogeneous rock in which the pore fluid is uniformly distributed in the pores, the shear modulus is not affected by the pore fluid, and the pore shapes are spheroidal. CONCLUSION

Four conclusions can be obtained from this report: (1): The Biot-Gassman theory can be used for predicting physical parameters (e.g., veloci!y, density, Poisson's ratio, bulk modulus, etc.) as a function of porosity or water saturation. (2): Oil and gas have different Poisson's ratio change curve as a function of water saturation for a given porosity. (3): When modeling gas-water and oil-water saturated reservoirs, it is important to obtain good estimation of initial velocity, initial porosity and dry rock Poisson's ratio. (4): Generally, decreasing the initial P-wave velocity or initial porosity, or increasing dry rock Poisson's ratio, and keeping the rest of the parameters constant, will cause Poisson's ratio to increase for a given water saturation. ACKNOWLEDGEMENT

This study was initially suggested by Brian Russell of Hampson-Russell Software Services Ltd. We wish to give our special thanks to him for his suggestions and kind help. Marco Leal of Canadian Occidental Petroleum Ltd.modified the FORTRAN program. The first author also would like to thank R. James Brown for showing the theory explained in Appendix A. REFERENCES

Domenico,S.N.,1974,Effectof watersaturationonseismicreflectivityof sandreservoirsencasedin shale: Geophysics, 39, 759-769. Gardner,G.H.F.,Gardner,LW., and Gregory,A.R., 1974,Formationvelocityand density- the diagnosticbasics for stratigraphic traps: Geophysics, 39, 770 - 780. Geerstma, J., and Smit, D.C., 1961, Some aspects of elastic wave propagation in fluid - saturated porous solids:

Geophysics, 26, 169-181. Gregory,A.R., 1977,Aspectsof rockphysicsfromlaboratoryand logdatathat are importantto seismic interpretations: AAPG Memoir, 26, 15-46. Hampson, D., and Russell, B., 1990, AVO: Amplitude versus offset analysis: handbook for AVO version 4.1. Hilterman,F., Gassman-Biot-Geertma Equation: unpublishedcoursenotes. Russell, B., 1990, The theory of the amplitude versus offset analysis program: unpublished notes. Sheriff,R.E., 1991,Encyclopedicdictionaryof explorationgeophysics,third edition:SEG.

305

APPENDIX

A

According to the concept of the compressibility: C = AT?Z = AT., AP ZAP

(A-l)

i.e., the compressibility (C) is the relative change in volume (ATJZ) with pressure (AP). AZ is the absolute change in volume, and Z is the total volume of the rock. The total change in volume AZ is equal to the change in volume of the pore space (AZp) plus the change in volume of the solid matter (AZs), i.e.: AZ = AZp + AZs,

therefore:

Am_Amp AZs ZZ 4 Z •

(A-2)

According to the concept of the porosity: Zp = 0Z, Zs= (I-,)Z, Z= Zp, Z=

¢

Zs

1-¢'

where Zp is the total volume of pore space and Zs is the total volume for solid matter. equation (A-2) becomes:

Then

AZ = CAZp _ (1-t_)AZs Z Zp Z_ ' divided

by Ap:

AZ _ OAZp + (1-_b)AZs ZAP ZpAP ZsAP comparing the above equation with equation (A-I), we can get: C d = _bCp +

(1-¢)Cs).

(A-3)

Equation (A-3) was used to obtain Kd in equation (8). Equation (A-3) can also be written as:

Ca= ¢Cp+ Cs- ¢cs. Cs is quite small compared to Cp, and t)Cs is very small compared with both @Cpand Cs. Therefore the above equation can be simplified as: Cd = ¢Cp + Cs. Equation (A-4) is the same equation given by Domenico (1974).

(A-4)