ARTÍCULO DE INVESTIGACIÓN

ANÁLISIS TEÓRICO DEL EFECTO DE LA PRODUCCIÓN DE POZOS VECINOS EN PRUEBAS DE INTERFERENCIA

THEORETICAL ANALYSIS OF THE EFFECT OF PRODUCTION OF NEIGHBORING WELLS IN INTERFERENCE TESTS

Freddy Humberto Escobar Macualo* Felipe Hernando Repizo** Matilde Montealegre***

Resumen

Normalmente, un pozo productor de hidrocarburos emplea el 99% de su vida productiva operando bajo el régimen de estado pseudoestable y continúa su operación mientras se ejecutan pruebas de interferencia en otros pozos dentro del mismo campo. Esto ocasiona que la estimación de parámetros utilizando datos recogidos durante una prueba de interferencia, en la cual los pozos involucrados se asumen únicos dentro de un yacimiento infinito, puede no ser muy precisa, en virtud a que la producción de pozos vecinos dentro del mismo yacimiento puede ocasionar una caída de presión adicional, que afectará la presión inicial y la medida durante la prueba, haciendo que su interpretación no sea adecuada.

En este artículo se presenta un método para corregir o extrapolar la presión aplicando el principio de superposición, incluyendo la declinación de la presión causada por pozos vecinos en producción, al modelamiento tradicional para pruebas de interferencia. El análisis e interpretación de los datos obtenidos son realizados con la novedosa técnica Síntesis Directa de Tiab, que analiza el comportamiento de los datos de presión y su derivada, sin el uso de curvas tipo. Si la derivada de la presión se calcula de manera apropiada y ésta se define satisfactoriamente, la técnica permite estimar directamente los parámetros del yacimiento con bastante precisión mediante puntos característicos, ya que usa soluciones analíticas.

Al hallar la derivada a partir de los datos tomados directamente de la prueba, la curva resultante no se definirá apropiadamente y los parámetros estimados pueden presentar errores. El contraste entre los resultados determinados en la interpretación de los datos tomados de la prueba y los corregidos mediante el principio de superposición, se ilustra con dos ejemplos sintéticos.

Palabras clave: superposición, análisis de pruebas de presión, permeabilidad, técnica TDS.

Abstract

Normally, a hydrocarbon well spends 99% of its productive life under pseudosteady state flow regime and continues operating while conducting interference test on other wells inside the same field. This probably causes inaccuracy in the estimation of parameters using recorded raw data during the interference test, in which the involved wells are considered unique inside an infinite reservoir. The production of neighboring

Artículo recibido: 02/06/2009 Aprobado: 12/08/2009

Doctor en Ingeniería de Petróleos. Profesor Asociado, Facultad de Ingeniería Universidad Surcolombiana. E-mail: fescobar@usco.edu.co Ingeniero de Petróleos de la Universidad Surcolombiana. E-mail: fepizo@hotmail.com *** Magíster en Ingeniería de sistemas y computación. Profesora de la Facultad de Ingeniería, Universidad Surcolombiana. E-mail: matildemm@usco.edu.co

ENTORNOS, No. 22. Universidad Surcolombiana. Vicerrectoría de Investigación y Proyección Social, 2009, pp. 141-150

wells inside the same reservoir can cause additional pressure drops, which affect the initial and measure pressure during the test, leading to a possible misinterpretation of the pressure test.

In this article a method is presented to either correct or extrapolate the pressure by applying the superposition principle, including the pressure decline caused by production of neighboring wells on the traditional modeling for interference testing. The analysis and interpretation of the data recorded is carried out by means of Tiab Direct Synthesis Technique, TDS technique, which analyzes the behavior of the pressure and pressure derivative data, without using type-curve matching. If the pressure derivative is properly calculated and successfully defined, the technique allows us to estimate the reservoir parameters directly with enough precision through characteristic points, using analytic solutions.

If pressure derivative is calculated employing directly the raw recorded data, the resulting curve will not be properly defined and the estimation of parameters may be inaccurate. A comparison of the results from the interpretation of the raw data and the corrected/extrapolated by the application of the superposition principle is illustrated with two synthetic examples.

Key words: superposition, pressure transient analysis, permeability, TDS technique.

1. Introducción

Con la introducción de la derivada de presión por Tiab and Kumar en 1976, el desarrollo de modelos complejos de interpretación y la aparición de nuevas y mejores herramientas de medición de presiones de fondo, las pruebas de presión se han convertido en una herramienta poderosa utilizada para la caracterización de yacimientos de hidrocarburos. La determinación de las características de flujo y estimación de parámetros son de gran importancia en la toma de decisiones para el desarrollo y explotación de un campo como en la localización de nuevos pozos y la mejor implementación de procesos de recobro mejorado.

Las técnicas tradicionales de interpretación de pruebas de interferencia asumen que la presión inicial del yacimiento se encuentra en equilibrio y consideran los pozos involucrados en la prueba como únicos dentro del yacimiento. No se tiene en cuenta la producción de pozos vecinos que afecten la presión de fondo inicial y la medida durante la prueba, debido al efecto de superposición. Además, los software comerciales están diseñados, en su mayoría, para modelos de yacimiento infinitos y podrían presentar ajustes de datos erróneos.

En este trabajo se plantean los errores que se pueden presentar en el análisis de los datos recogidos durante las pruebas de interferencia, al no tomar en cuenta el efecto de superposición de pozos vecinos en producción. Para mostrar las discrepancias entre los resultados obtenidos de la interpretación con datos tomados directamente de la prueba y los corregidos o extrapolados mediante el principio de superposición. Se utilizan datos simulados de una prueba de interferencia cuya interpretación se hace usando la técnica TDS, Tiab (1995), empleando los datos tomados directamente de la prueba y los corregidos mediante el principio de superposición, para efectos comparativos.

2. Efecto de superposición causada por la producción de pozos vecinos sobre las pruebas de interferencia

La caída de presión ocasionada por pozos vecinos en producción pueden afectar la presión inicial y la medida durante una prueba si existe interferencia; lo que busca este trabajo es plantear que al verse afectados los datos medidos en un pozo de observación al realizar una prueba, debido al efecto de superposición, la interpretación de los datos y la determinación de los parámetros del yacimiento presentan desviaciones. Los parámetros adimensionales están dados por (Earlougher, 1977):


Si consideramos un pozo el cual presenta una caída adicional de presión dada la interferencia de otros pozos operando dentro del mismo yacimiento (ver Figura 1), el comportamiento de la presión de fondo fluyendo, de acuerdo con Earlougher (1977), está dado por:

De la línea recta hallada al graficar en coordenadas cartesianas APwo (t) vs. t, que corresponde al régimen de flujo del estado pseudoestable, podemos determinar la pendiente m*, la cual es estimada antes de que la declinación de presión inicie:

dPws (Pws)2~(Pws)

(7)

dt

t2 ti


p = p

wf 1 i


'A'

K'-,


+ 1 ln

2


2,5458 C\


(8)


Asumiendo que todos los pozos producen a una rata constante a un tiempo t = 0, la caída de presión esperada, A Pwo (t), medido a un tiempo t = t1 + At puede ser estimada haciendo uso del principio de superposición en espacio.


141 2^ [Pfa r =li„.,)+s]-APow(t) (4)


kh


De acuerdo con Earlougher (1977), el comportamiento de la presión en un pozo bajo estado pseudoestable está dado por:


PD = 2ntDA + ^ ln


A


De las ecuaciones 4 y 7 obtenemos la siguiente expresión:

(5)

Normalmente todos los pozos de un yacimiento operan bajo condiciones de estado pseudoestable; la ecuación 4 quedaría (Escobar and Montealegre, 2007):

nr =

(9)

Es de aclarar que las pruebas de declinación no se efectúan para determinar el área del yacimiento, sin embargo, teóricamente es posible determinar dicho

Pnw (t) = b-m*t

(6)

ti

Tiempo de flujo Figura 1. Comportamiento de una prueba de declinación en un pozo repletado.


parámetro usando la ecuación 9 mediante el método convencional o la ecuación 10 mediante la técnica TDS, utilizando el punto de intersección entre las derivadas de la línea de comportamiento infinito y la línea de estado pseudoestable, trpil como lo presentó Tiab (1994):

kt .

rpi


A


301.77^^


(10) 4. Ejemplos sintéticos


3. Determinación de la permeabilidad del yacimiento

Si los datos de la prueba de presión están corregidos o no con la ecuación 6, la permeabilidad del yacimiento utilizando la técnica TDS puede estimarse mediante las ecuaciones 11 a 14, proporcionadas en el estudio de Escobar et al. (2007). Ver figura 2.

Pi = 4000 TVD = 4000 h = 30 ft Ct = 1x10-6

f = 20 % So = 1

V = 2 cp k = 30 md

Figura 2. Solución de la línea fuente, Escobar et al. (2007). --144


Ejemplo 1. Caso de tres pozos. Un yacimiento de forma circular con radio de 7000 ft (3533,9 Ac) tiene un pozo A localizado en el centro, un pozo B a 500 ft al este y otro pozo C situado a 600 ft exactamente al norte. El pozo C produce 200 BPD por 10000 hrs, mientras que el pozo B es abierto luego de 6000 hrs de producción del pozo C. El pozo A se mantiene cerrado y en él se registra la caída de presión ocasionada por el pozo B que produce 250 BPD. Datos adicionales:

Bo = 1.2 bb/STB rw = 0.5 ft

s = 0

Los datos de presión obtenidos de la prueba y los corregidos mediante el principio de superposición son ilustrados en la Figura 3 que muestra la relación de la presión con respecto al tiempo de producción de los pozos.


Solución. La pendiente m* obtenida de los últimos datos de la línea correspondiente al estado pseudoestable, antes de abrir el pozo B, es m* = -0.029272 psi/ciclo.

^ 3800

<1

*

Ph

Figura 3. Prueba de interferencia para el caso con tres pozos, ejemplo 1.


La Figura 4 presenta las presiones corregidas aplicando el principio de superposición y sus derivadas. Allí se observa el intercepto entre las dos curvas, (APf int = (íAP'J nt = 29 psi, con el cual se puede determinar una permeabilidad de 31,52 md por medio de la ecuación 12. En esta figura difícilmente se observa el régimen de flujo radial puesto que se desarrolla el flujo pseudoestable de pendiente unitaria rápidamente. Adicionalmente, del valor de la derivada durante flujo radial, (fAP)r = 51,1761 psi, con la cual se obtiene una permeabilidad de 27,59 md por medio de la ecuación 14. Nótese en la Figura 4 que es difícil determinar con exactitud la línea de flujo radial.


Figura 4. Gráfico log-log de presión y derivada de presión contra el tiempo a partir de los datos tomados directamente de la prueba para el ejemplo 1.

f = 30 % h = 100 ft s = 0


P, = 4000 w = 0,5 ft = 3 x10-6


1 w ct


El área del yacimiento, aclarando que estas pruebas no buscan obtener este parámetro, puede determinarse usando el punto de intercepto entre el flujo radial y el estado pseudoestable, trpi = 753 hrs, el cual da un valor de 1978,7 Ac usando la ecuación 10.

La figura 5 presenta la presión y la derivada en función del tiempo para los datos corregidos y extrapolados para el presente ejemplo; en forma similar al caso de los datos sin corregir se observa el intercepto entre las dos curvas, (APJ int = (t*AP'w) int = 30 psi, con el cual se halla una permeabilidad de 30,47 md por medio de la ecuación 12. El flujo radial está perfectamente definido con un valor de derivada de (t*A P)r =45,1 psi, con la cual se determina una permeabilidad de 31,30 md con la Ecuación 14 y un área de yacimiento de 3528,3 Ac usando un valor de trpi =1183 hrs en la ecuación 10. Para este caso, los parámetros obtenidos se ajustan muy bien a los empleados en la simulación.

Ejemplo 2. Caso con cuatro pozos. Considere cuatro pozos verticales dentro de un yacimiento circular de un área de 3535,6 Ac. El pozo 1 de observación estará cerrado durante la prueba de interferencia y se encuentra en el centro de dicho yacimiento. Un pozo 3 ubicado


Para este ejemplo se realizó la extrapolación al estado pseudoestable en dos ocasiones, con anterioridad a la apertura del pozo 2 y el pozo 4. Se obtuvo, m1* = -0,024884 psi/ ciclo. La pendiente m1* fue determinada con los últimos puntos antes de ser abierto el pozo 2. Luego, se obtuvo m2* = -0,02342 psi/ciclo. La pendiente m2* fue determinada con los últimos puntos antes de ser abierto el pozo 4.

Los datos obtenidos corresponden al análisis del transiente de presión desde las 10000 hrs a las 15000 hrs en las cuales produjo el pozo 4.


a 600 ft al norte del pozo de observación y produce a una rata constante de 850 BPD por 15000 hrs. El pozo 2 se encuentra a 500 ft exactamente al este del pozo 1 y es abierto luego de que el pozo 3 lleva 5000 hrs de producción, a una rata de 800 BPD. El pozo activo de la prueba de interferencia es el pozo 4 y se encuentra ubicado a 45° al noreste del pozo de observación a 781 ft y es abierto cuando el pozo 3 lleva 10000 hrs fluyendo, para producir 900 BPD (ver figuras 6 y 7).


Bo = 1.2 bb/STB |j = 2 cp k= 100 md


At, hrs

Figura 5. Gráfico log-log de presión y derivada de presión contra el tiempo a partir de los datos corregidos para el ejemplo 1.

Figura 6. Esquema ilustrativo del ejemplo 2.

isd ‘iv*ui - JMd ^JMd

At, hr

Figura 7. Prueba de interferencia para el caso con cuatro pozos, ejemplo 2.

Del intercepto entre las dos curvas, (APJ int = (t*AP'J ¡nt = 10.1 psi, de la figura 8, se determina una permeabilidad de 97,74 md usando la ecuación 12 y de 77,7 md con la ecuación 14 empleándose un valor de (fÁP)r= 19,6 psi. Con el valor de trpj = 260 hrs, el área estimada del yacimiento, Ec. 10, fue de 1281,9 Ac. Nuevamente se observa una desviación de los parámetros obtenidos y los reales.

La Figura 9 presenta la presión y la derivada en función del tiempo para los datos corregidos y extrapolados para el ejemplo 2. El intercepto entre las dos curvas, (APJ


Figura 8. Gráfico log-log de presión y derivada de presión contra el tiempo a partir de los datos tomados directamente de la prueba caso con cuatro pozos, Ejemplo 2.

At, hrs

Figura 9. Gráfico log-log de presión y derivada de presión contra el tiempo a partir de los datos corregidos, caso con cuatro pozos, ejemplo 2.

int = (tAP'f ¡nt = 9,8 psi, valor que permite hallar una permeabilidad de 110,7 md usando la ecuación 12 y de 112,3 md usando la ecuación 14 con (fAP)r = 13,5 psi. El área de yacimiento estimada fue de 3561,5 Ac usando un valor

de trpi=500 hrs en la ecuación 10.

5. Análisis de resultados

En los ejemplos podemos apreciar que cuando los datos de la prueba no se corrigen es muy difícil diferenciar el régimen de flujo radial, línea de pendiente cero en la derivada, ocasionando una determinación precisa de la permeabilidad de la formación. Aunque las pruebas de interferencia, por efectos económicos y prácticos, no se corren con el fin de determinar el área del yacimiento, los resultados difieren significativamente cuando no se corrigen los datos. No se establece el valor de desviación alguna, en virtud a que existe subjetividad a la hora de determinar la línea recta de flujo radial. Sin embargo, los resultados de los parámetros resultan muy satisfactorios cuando se efectúa la corrección.

ti

■int

6. Conclusiones

1.    Los resultados del análisis de los datos obtenidos directamente de la prueba presentan diferencias con respecto a los que se obtuvieron, con los datos depurados mediante el principio de superposición. Esto en virtud a que en el gráfico de la derivada de presión construido a partir de los datos tomados directamente de la prueba, se puede apreciar indefinición del régimen de flujo radial caracterizado por la pendiente cero, desarrollándose el flujo pseudoestable tempranamente.

K Permeabilidad, md q Rata de producción o inyección, STB/D Distancia radial entre pozos, ft D Distancia radial adimensional Radio del pozo, ft Tiempo, horas Tiempo adimensional P Presión, psi

P' Derivada de presión con respecto al tiempo, psi/hr

PD Presión adimensional

tD*P D Derivada de presión adimensional

AP Cambio de presión

Tiempo de intercepción, horas

tD    Tiempo adimensional

m    Viscosidad, cp

f    Porosidad, % o fracción.

Pwo(t)    Caída de presión esperada psi

Pwf    Presión de fondo fluyendo psi

A    Cambio, caída

Pws    Presión de cierre de pozo psi

APa    Caída de presión esperada psi

P*    Diferencia de presión entre la presión observada y la presión extrapolada, psi

Pi    Presión inicial psi

m    Pendiente psi/*ciclo

m*    Pendiente en el plano cartesiano durante el estado pseudoestable psi/hr

Subíndices

/    Nombre del pozo

int    Intercepción

r    Flujo radial

ext    Corregida o extrapolada

rpi    Intersección entre el estado radial y pseudoestable


A

B

C

h

Referentes bibliográficos

Escobar, F.H. and Montealegre, M. "Application of TDS Technique to Developed Reservoirs". Journal of Petroleum Science and Engineering. February 2007; 55 (Issues 3-4): 252-258.

Escobar, F.H., Cubillos, J. and Montealegre-M, M. "Estimation of Horizontal Reservoir Anisotropy without Type-Curve Matching". Journal of Petroleum Science and Engineering. 2008; 60: 31-38.

Earlougher, R.C., Jr. "Advances in Well Test Analysis". SPE Monograph 1977; 5.

Tiab, D. and Kumar, A. "Application of PD' Function to Interference Analysis" (paper SPE presented at the 51st Annual Fall Technical Conference and Exhibition of the Society of Petroleum Engineers of AIME, held in New Orleans, Oct. 3-6, 1976). Also, JPT Aug. 1980; 1465-1470.

Tiab, D. "Analysis of Pressure and Pressure Derivative without Type-Curve Matching: 1- Skin and Wellbore Storage". Journal of Petroleum Science and Engineering 1995; 12: 171-181.

Tiab, D. "Analysis of Pressure Derivative without Type-Curve Matching: Vertically Fractured Wells in Closed Systems". Journal of Petroleum Science and Engineering 1994; 11: 323-333.


150