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***
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.
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.
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
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.
|
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 |
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. 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 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 / 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
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