Professional Documents
Culture Documents
Modeling Simultaneus Oil and Water Flow
Modeling Simultaneus Oil and Water Flow
SPE 39755
Modeling Simultaneous Oil and Water Flow in Reservoirs with Water Influx or Water
Injection Using Single-Phase Semi-Analytical Solutions
AK. Permadi, SPE, Institut Teknologi Bandung, D.D. Mamora, SPE, and W.J. Lee, SPE, Texas A&M U.
263
2 -... A. K. PERMADI, D. D. MAMORA, W, J. LEE SPE 39755
,.
264
SPE 39755 MODELING SIMULTANEOUS OIL AND WATER FLOW IN RESERVOIRS WITH WATER INFLUX OR 3
WATER INJECTION USING SINGLE-PHASE SEMI-ANALYTICAL SOLUTIONS
265
4 A. K. PERMADI, D. D. MAMORA, W. J. LEE SPE 39755
..................................... ...(A-l)
13. Perrine,R.L.: Analysisof PressureBuildupCurves: Drill. and
Prod Prac., API, Dallas (1956) 482-509.
14. Martin, J.C.: Simplified Equations of Flow in Gas Drive where constant al = 0.006327 in field units (time in days).
Reservoirs and the Theoreticrd Foundation of Multiphase For the composite reservoir system, the following initial and
Pressure Buildup Andysr.s~ JPT (Oct. 1959) 309-11; Trans., reservoir boundary conditions apply:
AIME (1959) 216.
15. Double4 L.E. and Blasingame,T.A.: Decline Curve Analysis initial condition,
Using Type Curves: Water influx/Waterflood Cases: paper
SPE 30774 presentedat tie I995 Annurd Technical Conference
md Exhibition, Dallas, TX, Oct. 22-25. p(r,t)=pj, rw<r<re, *=0. Y ...............................(A-2)
16. Levine, J.S. and Prats, M.: The Calculated Performance of
Solution-Gas-Drive Reservoirs, SP.W (Sept. 1961) 142-52. inner boundary condition,
17. Raghav~ R.: Well Test Analysis: Wells Producing by
SolutionGas Drive; SPU (Aug. 1976) 196-208.
kro
18. F- M.L. and Wattenbarger,R.A.: Decline Curve Analysis
for MuItiphaseFlow: paper SPE 18274 presented at the 1988
Annual Technicrd Conference and Exhibition, Houston, TX,
2tia4
(qB), [(o~+pw~ )]r: =], r=rw, t>o,
............................................... (A-3)
Oct. 2-5.
19. Vo, D.T.: WeII Test Analysis for Gas Condensate Reservoirs,
Ph.D. Dissertation, The University of Tuls~ TUISAOK (1989). where the constant total liquid rate (qB), is defined as
20. Fetkovich, M.J.: A Simplified Approach to Water Influx
Calculations - Finite Aquifer Systems, JPT (Jul. 1971)814-28.
21. Dake, L.P.: Fundamentals of Reservoir Engineering, Elsevier (q~)tq.B. +qwBw =constant; .....................(A-4)
Scientific Publishing Co,, New York, NY (1978).
22. C- B.C., Hawkins, M.: Applied Petroleum Reservoir outer boundary condition,
figineering, Revised by Terry, R.E,, Prentice Hall, Englewood
Cliff% NJ (1991).
23.
24.
tiIlins, RE.: Flow of Flui& through Porous Materials, The
PetroleumPublishing Company, TuIs4 OK (1976).
Perniadi, A,K.: Modeling Simultaneous Oil and Water Flow
2&a4
9ext [( kro
Po +pw
Po
km r~
Pw )] ~
=1
,=r
t>.
.......................... . . . . . ..(A-5)
titz SingIe-Phase Analytical Solutions, Ph.D. Dissertation,
Tmas A&M University, College Station, TX (1997).
25. SABRE Users Manual, S.A, Holditch & Associates, College where qti = specified flux of aquifer influx, and constant
~titiO~ TX (1989). a4 = 0.001127 in field units.
26. Standing, M.B.: Notes on Relative Permeability Relationships,
Div; of Pet. Eng. and AppI. Geoph., University of Trondheim,
Trondheim, Norway (1975), Method of Transformation. l13e mult~hase diffisivity
27. Muska~ M.: *The Production Histories of Oil Producing Gas equation @q. A-5) may be solved by fwst transforming it into
Drive Ke~-fioirs, ~ ofAppIied Physics (1945) 16, 147-59. a form of the single-phase diffisivity equation, Following
Perrine-Martin, the following dimensionless variables are used
to express Eq. A-5 in dimensionless form.
266
SPE 39755 MODELING SIMULTANEOUS OIL AND WATER FLOW IN RESERVOIRS WITH WATER INFLUX OR 5
WAIER INJEC710N USING SINGLE-PHASE SEMI-ANALYTICAL SOLUTIONS
9kt =
(qextB),
..........................................................(A-16)
(qB)*
,nd
Water Influx Model. In most practical cases, the Fetkovich
rD. L, .........................................................(A-8) water influx model is stilciently accurate compared to the
rw
van Everdingen-Hurst method.zzz The aquifer influx rate
equation based on the Fetkovich method is as follows.
where a2 = 141.2in field units, and total mobility 2t is
defined as
dwe
w = =dt J(p-~-)ap(-::at) (A-17)
At= Ao+Aw=&+k Y .$$.,,$.,,.,...,..... ...................... (A-9)
Po Pw
me incremental water influx A W& at time step n may be
and total compressibility c1 is defined as follows: dervied from Eq. A-17 and is as follows.
Ct=soco+swcw+cf .,.,.......................................(A-lO)
Wej
Aw:= -(F;-1-P;.a,me~-aP[-J;::tn)](A18)
Pj,a
Using dimensionless variab~es defined in Eqs. A-6 to A-8,
Eqs. A-1, A-2, A-3, and A-5 may transformed into the
following dimensionless forms. where pa
*1 is determined from material balance, and
p~_=,we is the average value between time steps n and (n-l).
8PD 2PD........................................(A-n)
Id
.
rD arD (-)
D
8?.D =
dtD
Cumulative water influx at time step n is obtained by
summing up AW: from time step n=l until n.
initial conditio~ For a finite aquifer model, the pseudo steady state
productivity index is defined as
PD(rD,tD) = O, I S rD < reD, tD = O ;................(A-12)
M
J= / _\ , ............................................(A-19)
inner boundary condition,
41n:-:J
where the drawdown is expressed as ~a - pr_a instead of
()
D=
dpD
=
_l
rD=l, tD > O;,........................(A-13)
Pip Pr-a .
ne dimensionless pressure ptit and dimensionless
outer boundary condition,
aquifer influx rate q~t are related to parameters in the
Fetkovich model as follows. Defining the dimensionless
PD
(D) drD
= q~~(tD) > rD = reD > tD ~ o ............(A-14)
water influx as
AW:D=
A w:
>.................... ......................................(A-20)
We
The dimensionless pressure at the reservoir-aquifer boundary
ptit is defined as Eq. A- 18 may be arranged to yield
n-l AW~DPi,a
P!-a,ove = Pa ............................(A-21)
c
267
6 A. K. PERMADI, D. D. MAMORA, W. J. LEE SPE 39755
[)
Jpi,a Atn
Cn=l.ap . ......................................(A-22) relative permeability data and the following equation.
Wei
.fw= 1 ...................................................T...(A.29)
Using p~.a~ from Eq. (A-21), we compute P;.= which is ~+
P. kro
then used in Eq. A-15 to calculate p~t. Furthermore, krw P.
and 2fJs.f,BT
Sw,BT=s~+ dfw
................................T.(A-30)
*ch_
- +ctrw2 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (A-25) ds w
al~t S~,BT
Knowing the average water saturation in the reservoir, water
Material Balance For a reservoir with no original
Model. saturation at the well after breakthrough may be determined in
gascap and where the reservoir pressure is always above the the following manner. Compute: (1) fw using k, vs. Sw
bubble point, the material balance equation may be expressed
df w
as follows. da@ (2) ,the slope with respect to &w= Sw - SWC
-idsw s%
[1
Swicw+cf
~ Boi ....................(A-26) the slope with respect to &w= SW,i_Sw,i_l
l-s.j 4=0 3) +dsw sw
where i is the index for ( Sw, kr ) data point.
FoIlowing Muska<27the instantmeous water-oil ratio WOR
3W is computed at each time step in the model. If
isgiven by:
~w < ~w,BT, the weIl is not allowed to produce water.
dwr 1
_ dso _ (l S.)
dBw Conversely, if ~w 2 ~W,BT,the well is allowed to produce
WOR = ~ = w p W2 p ,.............(A-2iI] water. For example, water saturation at the well at
dr?r 1 dso
---
so dBo
breakthrough would be given as. S .,Well= S.f = Swf,BT. Afier
B. dp B02 dp breakthrough, water saturation at the well may be calculated
USill~ Eq.A-3 1.
where the change in oil saturation with respect to pressure is
given by:
~-f~sw,wll
SW,WII = 5. df
....................(A-31)
S. kr. PQ dBQ
(~-so) dBw w
dso
= Bo kro Pw dp Bw dp = Sw,wll
..................(A-28)
dp ~~ krw P.
268
SPE 39755 MODELING SIMULTANEOUS OIL AND WATER FLOW IN RESERVOIRS WITH WATER INFLUX OR 7
. .. .. WA_~R INJECTION USING SINGLE-PHASE SEMI-ANALYTICAL SOLUTIONS
TABLE IRESERVOIR, ROCK, AND FLUID PROPERTIES
initial Conditions
Depth of gas-oil contaot = 4000 ft
Capillary pressure at gas-oil oontact = O psia
Depth of oil-water contaot = 5000 fi
Capillaty pressure at oil-water contact = O psia
Resetvoir Propeties
Net sand thickness = 30 R.
Porosity, fraction = 0.15
lm3d@le water saturation = 0.22
Nominal well spacing = 31.4 acres
Horizontal absolute permeability = 500 md
Vertical absolute permeability = 500 md
Original Oil-in-Place = 855 MSTB
Constant oil viscosity = 3.0 Cp
Constant water vismsity = I.ocp
Production Parameters
Location of Producing Well (1)
Constant total production rate = 100 STB/D
Initial reservoir pressure = 2500 psia
Initial saturation pressure = 1000 psia
f---- ~\.,
TQDViW. . FiniteAquifer
Stie View:
269
8 A. K. PERMADI, D. D. ~MORA, W. J. LEE SPE 39755
. -.
I reservok prassurs 14
I
+ +
Using Fetkovti model unnpule wrnulatiie
Using s~le singl~se sOlutiix
weter infius Fcf the cumuhtiie period ci time
end a wssure dq at h mservoir+fer -@e Well preseure aml pressure at the
reaerwir aqufer interfsce
Wary
1
I
&
T Y
Yes
Lul
0.4
0.2
0
0 0.2 U.4 0.6 0.8 1
tm
* S-8 lhkwork
0 5,Em5
o Shubtbn
80 0
0 4.E+05
m 0
3,eo5
0 0
0
40
2,E+05
20 1.E+05
t
, oo~o
n Q . . 0 O.E+OO
o
1.2M 0,[ too 2,0-03 4.oEto3 8.0E+03 a.ow3 t ,0EM4 1,ZBW
0.0500 20- 4.0E+02 6.om now 1.OW
Froducini th., days
~dUCh# tic, da~
Fig. 4-Oil production rate versus time Fig. 7-Cumulative water production versus time
120 120
100
0
m
~ 80 0
do
40
,
0
i 0
j 40
0
mkwti
It ~ whtbn
20 0 20
0
t
Fig. 5-Water production rate versus time Fig. S-Water influx rate versus time
271
2.E+05
o.~
t
1.506 t ,&M
Fig. 9--Cumulative
i ,s02 1,*
1
2700
O,a 2400
mbwti
o~
210U
0.6
4800
04
1500
02 Thkwoti
1200 o Shhtbn
I !
o~
<Ecu IEM i.Ero2 1.E+OO 1.EM2 lH 1.E+od
ml
1.EW l.&M i.koz 1.2+00 1 .E+02 ! .E+04 l.E+m
Produzing thug, day,
7koducing tim , days
Fig. 10-Average reservoir water saturation versus time Fig. 13-Average resewoir pressure versus time
2700 2700
I
2400 2400
Z*M 2100
1am 1800
\
1SW 15W
_miiwork
~ tiktbn
12M
i 200
\
2UOI J 200 J
1.s06 i .6M l&02 1,- 1 ,Wz 1,- 1 ,m 1.E-04 l&02 1 .W 1 .E+02 1- lE4m
~OdUCiIIS Ykll , dvs
Produchg tlmc, days
Fig. 1lAverage aquifer pressure versus time Fig. 14-Bottomhole pressure versus time
2?C4
51
\ x
10
Soo
0
l.~m 1.E04 l.Ea? 1.E4M I.m l,sun 1.-
O.om 2.0E403 4.osto3 6.0- hom 1.H 1=
Wddudngtlm., W
Producing tima, dap
Fw. I>Reservoir-aqui fer boundary pressure versus time Fig. 1S-Water-oil ratio versus time
272