Comportamiento de afluencia IPR - Parte 1
Genera la curva de IPR con el método de Darcy.
La curva IPR es una representación gráfica de la relación entre la presión en el fondo del pozo y el gasto de producción de líquido. La magnitud de la pendiente de la curva IPR se denomina índice de productividad (\(PI\) or \(J\)), expresado como:
\[J=\frac{q}{P_e-P_{wf}} \ \ (1)\]
Para yacimientos de petróleo bajosaturados o regiones de yacimientos donde la presión es mayor que la presión del punto de burbuja podemos utilizar las siguientes ecuaciones
““Flujo radial transitorio en un pozo vertical”” \[J^*=\frac{q}{(P_e-P_{wf})}=\frac{kh}{162.6B_o\mu_o(logt+log \frac{k}{\phi\mu_oc_tr^{2}_w}-3.23+0.87S)}\] ““Flujo radial estacionario en un pozo vertical”” \[J^*=\frac{q}{(P_e-P_{wf})}=\frac{kh}{141.2B_o\mu_o(In\frac{r_e}{r_w}+S)}\] ““Flujo radial pseudo-estacionario en un pozo vertical”” \[J^*=\frac{kh}{141.2B_o\mu_o(In\frac{r_e}{r_w}-\frac{3}{4}+S)}\\J^*=\frac{q}{(P-P_{wf})}=\frac{kh}{141.2B_o\mu_o(\frac{1}{2}In\frac{4A}{\gamma{C_A}r^{2}_w}+S)}\] ““Flujo pseudo-estacionario en un pozo fracturado”” \[J^*=\frac{q}{(P_e-P_{wf})}=\frac{kh}{141.2B_o\mu_o(In\frac{r_e}{r_w}+S_f)}\] ““Flujo pseudo-estacionario en un pozo horizontal”” \[j^{*}=\frac{q}{(P_e-P_{wf})}=\frac {k_Hh} {141.2B_o\mu_o\left[In\left[\frac{a+\sqrt{a^{2}-(L/2)^{2}}}{L/2}\right]+\frac{1_{ani}h}{L}In\left[\frac{1_{ani}h}{r_w(1_{ani}+1)}\right]\right]}\]
Example in R
Usando los siguientes datos podemos calcular y graficar el IPR para un pozo vertical de un yacimiento de aceite considerando flujo pseudo-estacionario.
\(Porosity = 0.19\)
\(Permeability: k = 8.2 \ md\)
\(Pay \ zone \ thickness: h = 53 \ ft\)
\(Reservoir \ pressure: P_e = 5651 \ psi\)
\(Bubble-point \ pressure: p_b=5651 \ psia\)
\(Oil \ formation \ volume \ factor: B_o=1.1\)
\(Oil \ viscosity: \mu _o = 1.7 \ cp\)
\(Total compressibility: c_t = 0.0000129 \ psi^{-1}\)
\(Drainage \ area: A = 640 \ acres \ (re=2980 \ ft)\)
\(Wellbore \ radius: r_w=0.328 \ ft\)
\(Skin \ factor: S=0\)
Ver código
py <- 5651
k <- 8.2
h <- 53
Bo <- 1.1
viso <- 1.7
re <- 2980
rw <- 0.328
S <- 0
#Indice de productividad
J = (k*h)/(141.2*Bo*viso*(log(re/rw)-3/4+S))
print(J)[1] 0.1967785Los puntos de la curva IPR puede calcularse de la siguiente manera
Ver código
#Se crea un vector de presión desde la presion de yacimiento a 0
pwf <- seq(0,py, length.out = 10)
#Usando la ecuacion 1 y el valor de indice de productividad calculamos gasto de aceite (QO)
qo <- J*(py-pwf)
#Usando los vector de Pwf y Qo creamos un dataframe
IPR <- data.frame(Pwf = pwf, Qo = qo)
print(IPR) Pwf Qo
1 0.0000 1111.9951
2 627.8889 988.4401
3 1255.7778 864.8851
4 1883.6667 741.3301
5 2511.5556 617.7751
6 3139.4444 494.2201
7 3767.3333 370.6650
8 4395.2222 247.1100
9 5023.1111 123.5550
10 5651.0000 0.0000Usando ggplot2 graficamos la curva IPR
Ver código
library(ggplot2)
ggplot(IPR, aes(x = Qo, y = Pwf)) +
geom_point() + geom_line()
Usando el código anterior podemos definir una función para agilizar los caculos.
Ver código
Darcy_PSS <- function(py,k,h,Bo,viso,re,rw,S){
J = (k*h)/(141.2*Bo*viso*(log(re/rw)-3/4+S))
pwf <- seq(0,py, length.out = 10)
qo <- J*(py-pwf)
IPR <- data.frame(Pwf = pwf, Qo = qo)
return(IPR)
}
IPR <- Darcy_PSS(4500, 8.2,53,1.1,1.7,2980,0.328,0)
print(IPR) Pwf Qo
1 0 885.50310
2 500 787.11387
3 1000 688.72464
4 1500 590.33540
5 2000 491.94617
6 2500 393.55693
7 3000 295.16770
8 3500 196.77847
9 4000 98.38923
10 4500 0.00000Ahora podemos comparar facilmente el IPR con y sin factor de daño
Ver código
IPR_nS <- Darcy_PSS(5000, 8.2,53,1.1,1.7,2980,0.328,0)
IPR_S <- Darcy_PSS(5000, 8.2,53,1.1,1.7,2980,0.328,2)
ggplot() +
geom_point(data = IPR_nS, aes(x = Qo, y = Pwf), color = "red") +
geom_line(data = IPR_nS, aes(x = Qo, y = Pwf), color = "red")+
geom_point(data = IPR_S, aes(x = Qo, y = Pwf), color = "green") +
geom_line(data = IPR_S, aes(x = Qo, y = Pwf), color = "green")
Referencias:
Guo, Boyun (2008) Well productivity handbook