maggiolo, r. - optimización de la producción mediante análisis nodal

Optimizacin de la Produccin mediante Anlisis Nodal

1

Programa de cursos 2008

ESP Oil_Copyright_01-P35/The document is property of ESP Oil Consultants, any reprodu

ction is strictly forbidden and will be prosecuted

Ing. Ricardo Maggiolo

OOppttiimmiizzaacciinn ddee llaaPPrroodduucccciinn mmeeddiiaannttee

AAnnlliissiiss NNooddaall Dictado por:

Msc. Ricardo Maggiolo

Del 07 al 11 de Julio de 2008 Instalaciones del Hotel El Condado

Lima - Per


ESP Oil_Copyright_01-P35/The document is property of ESP Oil Consultants, any reproduction is strictly forbidden and will be prosecuted

2



OPTIMIZACIN DE LA PRODUCCIN MEDIANTE

ANLISIS NODAL

INTRODUCCIN Las compaas productoras de petrleo y gas realizan continuamente grandes esfuerzos por agregar valor a sus corporaciones y mejorar as sus resultados financieros. Estos esfuerzos estn dirigidos a mediano y largo plazo a maximizar el factor de recobro de los yacimientos y a corto plazo a acelerar el recobro de las reservas recuperables, la primera es una meta de aos para el equipo multidisciplinario de personas que laboran en la Optimizacin Integrada del Yacimiento, la segunda es el da a da del equipo multidisciplinario de personas que laboran en la Optimizacin Total del Sistema de Produccin. Esta ltima, aunque es un subproceso de la primera, constituye el Ncleo del Negocio (Core Business) de la Corporacin ya que permite maximizar la produccin total diaria de hidrocarburos y/o el beneficio neto (M$$$) producto de la venta de los mismos. Una de las tcnicas mas utilizadas para optimizar sistemas de produccin, dada su comprobada efectividad y confiabilidad a nivel mundial, es el Anlisis Nodal; con la aplicacin de esta tcnica se adecua la infraestructura tanto de superficie como de subsuelo, para reflejar en el tanque el verdadero potencial de produccin de los pozos asociados a los yacimientos del sistema total de produccin. En otras palabras, se logra cerrar la brecha existente entre la produccin real de los pozos y la produccin que debera exhibir de acuerdo a su potencial real de produccin. El Anlisis Nodal bsicamente consiste en detectar restricciones al flujo y cuantificar su impacto sobre la capacidad de produccin total del sistema. Existen en el mercado varios simuladores comerciales que permiten aplicar dicha tcnica, entre los ms conocidos se tienen, por ejemplo: PERFORM-PIPESOFT2 de IHS, PIPESIM-PIPESIM GOAL y NET de Schlumberger (BJ), PROSPER-GAP de Petroleum Expert, WELLFLO-FIELDFLO-ReO de Weatherford (EPS), etc. El presente curso tiene como objetivo: Describir y aplicar una metodologa para optimizar Sistemas de Produccin de Hidrocarburos utilizando la tcnica del Anlisis Nodal. Para el cumplimiento de este objetivo se estructur un contenido programtico de cuatro captulos: en el captulo 1 se describe el sistema de produccin haciendo nfasis en el balance de energa requerido entre el yacimiento y la infraestructura instalada para establecer la capacidad de produccin del pozo. Adicionalmente se describe el principio fundamental de funcionamiento de los mtodos de levantamiento artificial. En el captulo 2 se detallan las ecuaciones y modelos matemticos simplificados para cuantificar la capacidad de aporte de fluidos de las formaciones productoras incluyendo el dao a la formacin y la forma de completacin del pozo (empaque con grava, caoneo convencional, etc.), se ilustra el uso de las mismas a travs de ejemplos numricos.



3



En el captulo 3 se describen algunas correlaciones de flujo multifsico en tuberas para cuantificar las prdidas y transformacin de la energia en las diferentes tuberas instaladas tanto en el pozo como en la superficie. Se ilustra el uso de las correlaciones a travs de ejemplos numricos. En el captulo 4 se determina la capacidad de produccin de pozos que producen por flujo natural, por levantamiento artificial por gas y por bombeo electro-sumergible. En el captulo 5 se describe la metodologa de optimizacin donde despus de cotejar el comportamiento actual del pozo se optiman los componentes del sistema tanto en superficie como en el subsuelo. Finalmente se describen diferentes escenarios de optimizacin del sistema incluyendo la distribucin ptima de gas en un sistema de gas-lift. A pesar de que solo se utilizar un simulador comercial como herramienta de optimizacin, no se sacrificar la generalidad de la aplicacin de la metodologa con otros simuladores disponibles en el mercado. La siguiente figura seala el marco de referencia donde se aplicar la metodologa de optimizacin.

Expansin lquida

Expansin gas en solucin

Expansin de una capa de gas

VOLUMTRICO

P

T

P

T

Expansin de un acufero

Cf = - 1/V . dV/dP

V = Cf . V . P

Marco de referencia



4



CONTENIDO

CAPTULO 1 EL SISTEMA DE PRODUCCIN

1.1 El Sistema de produccin y sus componentes 1.2 Proceso de produccin

Recorrido de los fluidos en el sistema 1.3 Capacidad de produccin del sistema.

Curvas de oferta y demanda de energa en el fondo del pozo. Balance de energa y capacidad de produccin Optimizacin del sistema Mtodos de produccin: Flujo natural y Levantamiento artificial

CAPTULO 2 COMPORTAMIENTO DE AFLUENCIA DE FORMACIONES PRODUCTORAS

2.1 Flujo de fluidos en el yacimiento. Estados de flujo Flujo de petrleo

Flujo No-Continuo o Transitorio (Unsteady State Flow Flujo Continuo o Estacionario (Steady State Flow): Ecuacin de Darcy para flujo continuo Flujo Semi-continuo (Pseudo-steady State Flow): ndice de productividad Eficiencia de flujo (EF) IPR (Inflow Performance Relationships). Ejercicios

Flujo de petrleo y gas en yacimientos saturados Ecuacin y Curva de Vogel para yacimientos saturados

Flujo de petrleo y gas en yacimientos sub-saturados Ecuacin de Vogel para yacimientos subsaturados

2.2 Flujo de fluidos en la completacin Tipos de completacin Hoyo desnudo Caoneo convencional Empaque con grava Cada de presin en la completacin Ecuaciones de Jones, Blount y Glaze Ejercicios Curva de oferta de energa o afluencia de fluidos que el yacimiento entrega en el fondo del pozo

CAPTULO 3 FLUJO MULTIFSICO EN TUBERAS

3.1 Flujo de fluidos en el pozo y en la lnea de flujo Algoritmo para calcular las prdidas de presin del fluido. Ecuacin general del gradiente de presin dinmica Clculo de la presin requerida en el cabezal Clculo de la presin requerida en el fondo del pozo



5



3.2 Consideraciones tericas del flujo multifsico en tuberas

Clculo del factor de friccin Definiciones bsicas: factor Hold-Up, densidad y viscosidad bifsica, etc. Patrones de flujo

3.3 Descripcin de correlaciones de flujo multifsico en tuberas

Correlacin de Hagedorn & Brown Correlacin de Duns & Ros Correlacin de Orkiszewski Correlacin de Beggs and Brill Ejemplos numricos Ejemplos con curvas de gradiente ya graficadas

3.4 Construccin de Curva de Demanda de energa Rangos caractersticos de la curva de demanda

CAPTULO 4 CAPACIDAD DE PRODUCCIN DEL SISTEMA

4.1 Capacidad de produccin del pozo en flujo natural Tasa de produccin posible o de equilibrio. Ejercicio Uso de reductores para controlar la produccin del pozo en FN Ecuaciones para estimar el comportamiento de estranguladores o reductores

4.2 Capacidad de produccin del pozo de Levantamiento Artificial por Gas Curva de rendimiento del pozo de LAG

4.3 Capacidad de produccin del pozo con bombeo electrocentrfugo sumergible (BES)

Curva de rendimiento del pozo en funcin de las RPM del motor CAPTULO 5 OPTIMIZACIN DEL SISTEMA DE PRODUCCIN

5.1 Cotejo del comportamiento actual del pozo Seleccin y Ajuste de las correlaciones empricas para calcular las propiedades del petrleo Seleccin y Ajuste de las correlaciones de Flujo Multifsico en Tuberas Cotejo del Comportamiento actual de Produccin

5.2 Optimizacin del sistema de produccin Anlisis Nodal del pozo: Oportunidades de aumentar la Oferta de energa y

fluidos del Yacimiento. Anlisis Nodal del pozo: Oportunidades de disminuir la Demanda de energa

para levantar fluidos del Yacimiento. Casos de estudio con utilizando un simulador de anlisis nodal.



6



CAPTULO I

El Sistema de Produccin

1.1 El Sistema de produccin y sus componentes El sistema de produccin est formado por el yacimiento, la completacin, el pozo y las lneas de flujo en la superficie. El yacimiento es una o varias unidades de flujo del subsuelo creadas e interconectadas por la naturaleza, mientras que la completacin (perforaciones caoneo), el pozo y las facilidades de superficie es infraestructura construida por el hombre para la extraccin, control, medicin, tratamiento y transporte de los fluidos hidrocarburos extrados de los yacimientos. 1.2 Proceso de produccin El proceso de produccin en un pozo de petrleo, comprende el recorrido de los fluidos desde el radio externo de drenaje en el yacimiento hasta el separador de produccin en la estacin de flujo. En la figura se muestra el sistema completo con cuatro componentes claramente identificados: Yacimiento, Completacin, Pozo, y Lnea de Flujo Superficial. Existe una presin de partida de los fluidos en dicho proceso que es la presin esttica del yacimiento, Pws, y una presin final o de entrega que es la presin del separador en la estacin de flujo, Psep.

YACIMIENTOYACIMIENTOCOMPLETACIN

Pesttica promedio (Pws)PRESIN DE ENTRADA: Pesttica promedio (Pws)PRESIN DE ENTRADA:

Pseparador (Psep)PRESIN DE SALIDA:

Pseparador (Psep)PRESIN DE SALIDA:

LINEA DE FLUJO

OP

OZ

LINEA DE FLUJO

OP

OZOP

OZOP

OZ

PwsPwfsPwf

Pwh Psep



7



Recorrido de los fluidos en el sistema Transporte en el yacimiento: El movimiento de los fluidos comienza en el

yacimiento a una distancia re del pozo donde la presin es Pws, viaja a travs del medio poroso hasta llegar a la cara de la arena o radio del hoyo, rw, donde la presin es Pwfs. En este mdulo el fluido pierde energa en la medida que el medio sea de baja capacidad de flujo (Ko.h), presente restricciones en la cercanas del hoyo (dao, S) y el fluido ofrezca resistencia al flujo (o). Mientras mas grande sea el hoyo mayor ser el rea de comunicacin entre el yacimiento y el pozo aumentando el ndice de productividad del pozo. La perforacin de pozos horizontales aumenta sustancialmente el ndice de productividad del pozo.

Transporte en las perforaciones: Los fluidos aportados por el yacimiento atraviesan

la completacin que puede ser un revestidor de produccin cementado y perforado, normalmente utilizado en formaciones consolidadas, o un empaque con grava, normalmente utilizado en formaciones poco consolidadas para el control de arena. En el primer caso la prdida de energa se debe a la sobrecompactacin o trituracin de la zona alrededor del tnel perforado y a la longitud de penetracin de la perforacin; en el segundo caso la perdida de energa se debe a la poca rea expuesta a flujo. Al atravesar la completacin los fluidos entran al fondo del pozo con una presin Pwf.

Transporte en el pozo: Ya dentro del pozo los fluidos ascienden a travs de la

tubera de produccin venciendo la fuerza de gravedad y la friccin con las paredes internas de la tubera. Llegan al cabezal del pozo con una presin Pwh.

Transporte en la lnea de flujo superficial: Al salir del pozo si existe un reductor de

flujo en el cabezal ocurre una cada brusca de presin que depender fuertemente del dimetro del orificio del reductor, a la descarga del reductor la presin es la presin de la lnea de flujo, Plf, luego atraviesa la lnea de flujo superficial llegando al separador en la estacin de flujo, con una presin igual a la presin del separador Psep, donde se separa la mayor parte del gas del petrleo.

En las siguientes figuras se presentan los componentes del sistema de una manera mas detallada as como el perfil de presin en cada uno de ellos.



8



Componentes del Sistema y Perfil de presiones

La perdida de energa en forma de presin a travs de cada componente, depende de las caractersticas de los fluidos producidos y, especialmente, del caudal de flujo transportado en el componente.



9



1.3 Capacidad de produccin del sistema. La perdida de energa en forma de presin a travs de cada componente, depende de las caractersticas de los fluidos producidos y, especialmente, del caudal de flujo transportado, de tal manera que la capacidad de produccin del sistema responde a un balance entre la capacidad de aporte de energa del yacimiento y la demanda de energa de la instalacin para transportar los fluidos hasta la superficie. La suma de las prdidas de energa en forma de presin de cada componente es igual a la prdida total, es decir, a la diferencia entre la presin de partida, Pws, y la presin final, Psep: Pws Psep = Py + Pc + Pp + Pl Donde: Py = Pws Pwfs = Cada de presin en el yacimiento, (IPR). Pc = Pwfs- Pwf = Cada de presin en la completacin, (Jones, Blount & Glaze). Pp = Pwf-Pwh = Cada de presin en el pozo. (FMT vertical). Pl = Pwh Psep = Cada de presin en la lnea de flujo. (FMT horizontal) Tradicionalmente el balance de energa se realiza en el fondo del pozo, pero la disponibilidad actual de simuladores del proceso de produccin permite establecer dicho balance en otros puntos (nodos) de la trayectoria del proceso de produccin: cabezal del pozo, separador, etc. Para realizar el balance de energa en el nodo se asumen convenientemente varias tasas de flujo y para cada una de ellas, se determina la presin con la cual el yacimiento entrega dicho caudal de flujo al nodo, y la presin requerida en la salida del nodo para transportar y entregar dicho caudal en el separador con una presin remanente igual a Psep. Por ejemplo, s el nodo esta en el fondo del pozo:



10



Presin de llegada al nodo: Pwf (oferta) = Pws - Py Pc Presin de salida del nodo: Pwf (demanda)= Psep + Pl + Pp

Pws

PsepPsepPsep PwsPwsNODO En cambio, si el nodo esta en el cabezal del pozo: Presin de llegada al nodo: Pwh (oferta) = Pws py pc - Pp Presin de salida del nodo: Pwh (demanda) = Psep + Pl

Pws

PsepPsepNODO

Pws



11



Curvas de oferta y demanda de energa en el fondo del pozo: Curvas VLP / IPR.

La representacin grfica de la presin de llegada de los fluidos al nodo en funcin del caudal o tasa de produccin se denomina Curva de Oferta de energa del yacimiento (Inflow Curve), y la representacin grfica de la presin requerida a la salida del nodo en funcin del caudal de produccin se denomina Curva de Demanda de energa de la instalacin (Outflow Curve). Si se elige el fondo del pozo como el nodo, la curva de oferta es la IPR (Inflow Performance Relationships) y la de demanda es la VLP (Vertical Lift Performance)

VLP

IPR

Pwf qliq. Como realizar el balance de energa? El balance de energa entre la oferta y la demanda puede obtenerse numrica o grficamente. Para realizarlo numricamente consiste en asumir varias tasas de produccin y calcular la presin de oferta y demanda en el respectivo nodo hasta que ambas presiones se igualen, el ensayo y error es necesario ya que no se puede resolver analticamente por la complejidad de las formulas involucradas en el calculo de las Ps en funcin del caudal de produccin.



12



Pyacimiento Pcompletacin

Ppozo Plnea

Donde: qo= Tasa de produccin, bbpd. o= Viscosidad, cps Factor volumtrico del petrleo, by/bn. re= Radio de drenaje, pies. rw= Radio del pozo, pies. S= Factor de dao, adim. Ko= Permeabilidad efectiva al petrleo, md. h= Espesor de arena neta petrolfera, pies. = Coeficiente de velocidad para flujo turbulento, 1/pie. o= Densidad del petrleo, lbm/pie3 rp= Radio de la perforacin, pulg. rc= Radio de la zona triturada alrededor del tnel perforado, pulg. Lp= Longitud del tnel perforado, pies. Kp= Permeabilidad de la zona triturada, md. TPP= Densidad de tiro, tiros/pie. hp= Longitud del intervalo caoneado, pies. g= Aceleracin de la gravedad, 32,2 pie/seg2 gc= Constante gravitacional, 32,2 pie/seg2. lbm/lbf. g/gc= Conversin de maas en fuerza, 1 lbf/lbm. At= Area seccional de la tubera, pie2. Z= Longitud del intervalo de tubera, pies.m= Densidad de la mezcla multifsica gas-petrleo, lbm/pie3 = Angulo que forma la direccin de flujo con la horizontal. fm= Factor de friccin de Moody de la mezcla multifsica gas-petrleo, adim. Vm= Velocidad de la mezcla multifsica gas-petrleo,pie/seg.

LgLLm HH 1Densidad:tA

oBoq,Vm

86400

6155 tA

gBsRRGPoq

86400

Velocidad:

h.Ko,

S,)rw/re(LnBo.o.qoPws007080

750 qo . Kp.Lp10 3- .,

)rprcLn(.Bo.o

qo 2 . Lp 2

)rc1-

rp1(. o.Bo 2..10 14-.30,2

007080

Psep)( Z. gc 2

Vm .m + d. gc 2Vm .m .fm +

gc

sen .m.g 144

Z n

22

1

) Z. gc 2


gc

sen .m.g 144

Z m

22

1

TPP2 . hP2 TPP . hP

(

Pyacimiento Pcompletacin

Ppozo Plnea

LgLLm HH 1Densidad: LgLLm HH 1Densidad:tA

oBoq,Vm

86400

6155 tA

gBsRRGPoq

86400

Velocidad:tA

oBoq,Vm

86400

6155 tA

gBsRRGPoq

86400tA

oBoq,Vm

86400

6155 tA

gBsRRGPoq

86400

Velocidad:

h.Ko,

S,)rw/re(LnBo.o.qoPws007080

750 qo . Kp.Lp10 3- .,

)rprcLn(.Bo.o

qo 2 . Lp 2

)rc1-

rp1(. o.Bo 2..10 14-.30,2

007080

Psep)( Z. gc 2


gc

sen .m.g 144

Z n

22

1

) Z. gc 2


gc

sen .m.g 144

Z m

22

1

TPP2 . hP2 TPP . hP

(

h.Ko,S,)rw/re(LnBo.o.qoPws

007080750 qo .

Kp.Lp10 3- .,

)rprcLn(.Bo.o

qo 2 . Lp 2

)rc1-

rp1(. o.Bo 2..10 14-.30,2

007080

Psep)( Z. gc 2


gc

sen .m.g 144

Z n

22

1

) Z. gc 2


gc

sen .m.g 144

Z m

22

1

(

h.Ko,S,)rw/re(LnBo.o.qoPws

007080750 qo .

Kp.Lp10 3- .,

)rprcLn(.Bo.o

qo 2 . Lp 2

)rc1-

rp1(. o.Bo 2..10 14-.30,2

007080

Psep)( Z. gc 2


gc

sen .m.g 144

Z n

22

1

) Z. gc 2


gc

sen .m.g 144

Z m

22

1

TPP2 . hP2 TPP . hPTPP . hP

(



13



Para obtener grficamente la solucin, se dibujan ambas curvas en un papel cartesiano y se obtiene el caudal donde se interceptan. La figura muestra el procedimiento paso a paso: YACIMIENTOCOMPLETACIN

Pws

LINEA DE FLUJO

OP

OZOP

OZ

2.- Se repite el paso anterior para otros valores asumidos de ql, y se construye la curva de Oferta de energa del Sistema.

Como estimar la Capacidad de Produccin del Sistema ?

PwfsPwf

1.- Dado un valor de ql en superficie se determina Pwfs y Pwf a partir de la Pws, luego se tabula y grafica Pwf vs. ql.

ql Pwfs Pwf

ql

Pwf

Oferta

Para obtener la curva de oferta en el fondo del pozo es necesario disponer de un modelo matemtico que describa el comportamiento de afluencia de la arena productora, ello permitir computar Py y adicionalmente se requiere un modelo matemtico para estimar la cada de presin a travs del caoneo o perforaciones (Pc) y para obtener la curva de demanda en el fondo del pozo es necesario disponer de correlaciones de flujo multifsico en tuberas que permitan predecir aceptablemente Pl y Pp.

Demanda

3.- Similarmente para cada valor de ql en superficie se determina Pwh y Pwf a partir de la Psep y se construye la curva de Demanda.

PsepPwh

PwfPwfPwfPwf

Pwh Pwf

ql

Pwf

Capacidad de Produccin del Sistema.

ql = ?



14



- Optimizacin Global del Sistema Una de las principales aplicaciones de los simuladores del proceso de produccin es optimizar globalmente el sistema lo cual consiste en eliminar o minimizar las restricciones al flujo tanto en superficie como en el subsuelo, para ello es necesario la realizacin de mltiples balances con diferentes valores de las variables ms importantes que intervienen en el proceso, para luego, cuantificar el impacto que dicha variable tiene sobre la capacidad de produccin del sistema.

Qliq.

Pwf

AUMENTANDOOFERTA

La tcnica puede usarse para optimizar la completacin del pozo que aun no ha sido perforado, o en pozos que actualmente producen quizs en forma ineficiente. Para este anlisis de sensibilidad la seleccin de la posicin del nodo es importante ya que a pesar de que la misma no modifica la capacidad de produccin del sistema, si interviene en el tiempo de ejecucin del simulador. El nodo debe colocarse justamente antes (extremo aguas arriba) o despus (extremo aguas abajo) del componente donde se modifica la variable. Por ejemplo, si se desea estudiar el efecto que tiene el dimetro de la lnea de flujo sobre la produccin del pozo, es ms conveniente colocar el nodo en el cabezal o en el separador que en el fondo del pozo. La tcnica comercialmente recibe el nombre de Anlisis Nodal (Nodal Systems AnalysisTM) y puede aplicarse para optimar pozos que producen por flujo natural o por levantamiento artificial

Marca registrada por Dowell-Schlumberger

DEMANDADEMANDA

OFERTAOFERTA

DISMINUYENDOLA DEMANDA

q3q3q1q1q1 q2q2

PwsPws

PsepPsepPsep

Pwfcrit. Pwfcrit.

Ing. de YacimientoIng. de Produccin qL = J ( Pws - Pwf )sinergiaIng. de YacimientoIng. de Produccin qL = J ( Pws - Pwf )sinergia



15



- Mtodos de produccin: Flujo Natural y Levantamiento Artificial Cuando existe una tasa de produccin donde la energa con la cual el yacimiento oferta los fluidos, en el nodo, es igual a la energa demandada por la instalacin (separador y conjunto de tuberas: lnea y tubera de produccin) sin necesidad de utilizar fuentes externas de energa en el pozo, se dice entonces que el pozo es capaz de producir por FLUJO NATURAL. A travs del tiempo, en yacimientos con empuje hidrulico, los pozos comienzan a producir con altos cortes de agua la columna de fluido se har mas pesada y el pozo podra dejar de producir. Similarmente, en yacimientos volumtricos con empuje por gas en solucin, la energa del yacimiento declinar en la medida en que no se reemplacen los fluidos extrados trayendo como consecuencia el cese de la produccin por flujo natural.

Pws

20 %

30 %

50 %

NO FLUYE

Empuje Hidrulico

0 %

Qliq.

Pwf % AyS

q1 q2 q3

Pwf NO FLUYE

Pws1

Pws2Pws3Pws4

Empuje por gas en solucin

Qliq.

RGL (pcn/bn)

400 600 800 1000

q2q3 q1

Cuando cesa la produccin del pozo por flujo natural, se requiere el uso de una fuente externa de energa para lograr conciliar la oferta con la demanda; la utilizacin de esta fuente externa de energa en el pozo con fines de levantar los fluidos desde el fondo del pozo hasta el separador es lo que se denomina mtodo de LEVANTAMIENTO ARTIFICIAL.

LAG

qL

q

BOMBEO

NO FLUJO

DISMINUYENDO DEMANDA EN LA VLVULA

AUMENTANDO OFERTA EN LA DESCARGA DE

LA BOMBA

qL

Demanda (Outflow )

Oferta (Inflow)



16



Entre los mtodos de Levantamiento Artificial de mayor aplicacin en la Industria Petrolera se encuentran: el Levantamiento Artificial por Gas (L.A.G), Bombeo Mecnico (B.M.C) por cabillas de succin, Bombeo Electro-Centrifugo Sumergible (B.E.S), Bombeo de Cavidad Progresiva (B.C.P) y Bombeo Hidrulico tipo Jet ( B.H.J). El objetivo de los mtodos de Levantamiento Artificial es minimizar los requerimientos de energa en la cara de la arena productora con el objeto de maximizar el diferencial de presin a travs del yacimiento y provocar, de esta manera, la mayor afluencia de fluidos sin que generen problemas de produccin: migracin de finos, arenamiento, conificacin de agua gas, etc.

qliq

Pwf

IPR

qliq

Pwf

IPR

Pwf

IPRqliqPwf



17



- Principio de funcionamiento de cada mtodo de produccin: A travs de recursos audiovisuales presentados por la compaa Weatherford se explicaran el principio de levantamiento/funcionamiento de cada mtodo. El ingeniero de produccin debe participar en el desarrollo del plan de explotacin del yacimiento para realizar una adecuada seleccin del mtodo o mtodos de levantamiento en los pozos, acorde con la estrategia de explotacin establecida.



18



CAPTULO II

Comportamiento de afluencia de formaciones productoras 2.1 Flujo de fluidos en el yacimiento. Estados de flujo. La simulacin del flujo de fluidos en el yacimiento debe considerar la composicin de los fluidos presentes, y las condiciones de presin y temperatura para establecer si existe flujo simultneo de petrleo, agua y gas, las heterogeneidades del yacimiento, etc. Para describir el flujo de fluidos en el yacimiento a travs del tiempo, se debe utilizar el modelaje matemtico de yacimientos y las soluciones numricas de la ecuacin de difusividad obtenidas con los simuladores comerciales (Familia Eclipse, por ejemplo). La simulacin numrica de yacimientos es materia que no ser tratada en este curso. La capacidad de aporte del yacimiento hacia el pozo se cuantificar en este curso a travs de modelos matemticos simplificados como por ejemplo: la ecuacin de Vogel, Fetckovich, Jones Blount & Glace, etc. rea de drenaje

Con fines de simplificar la descripcin del flujo de fluidos en el yacimiento se considerar el flujo de petrleo negro en la regin del yacimiento drenada por el pozo, comnmente conocida como volumen de drenaje, y adicionalmente, se asumir homogneo y de espesor constante (h) por lo que en lo sucesivo se hablar de rea de drenaje del yacimiento.

Flujo de petrleo en el yacimiento

El movimiento del petrleo hacia el pozo se origina cuando se establece un gradiente de presin en el rea de drenaje y el caudal o tasa de flujo depender no solo de dicho gradiente, sino tambin de la capacidad de flujo de la formacin productora, representada por el producto de la permeabilidad efectiva al petrleo por el espesor de arena neta petrolfera (Ko.h) y de la resistencia a fluir del fluido representada a travs de su viscosidad (o). Dado que la distribucin de presin cambia a travs del tiempo es necesario establecer los distintos estados de flujo que pueden presentarse en el rea de drenaje al abrir a produccin un pozo, y en cada uno de ellos describir la ecuacin que regir la relacin entre la presin fluyente Pwfs y la tasa de produccin qo que ser capaz de aportar el yacimiento hacia el pozo.



19



Estados de flujo:

Existen tres estados de flujo dependiendo de cmo es la variacin de la presin con tiempo:

1. Flujo No Continuo: dP/dt 0 2. Flujo Continuo: dP/dt = 0 3. Flujo Semicontinuo: dP/dt = constante

1) Flujo No-Continuo o Transitorio (Unsteady State Flow):

Es un tipo de flujo donde la distribucin de presin a lo largo del rea de drenaje cambia con tiempo, (dP/dt 0). Este es el tipo de flujo que inicialmente se presenta cuando se abre a produccin un pozo que se encontraba cerrado viceversa. La medicin de la presin fluyente en el fondo del pozo (Pwf) durante este perodo es de particular importancia para las pruebas de declinacin y de restauracin de presin, cuya interpretacin a travs de soluciones de la ecuacin de difusividad, permite conocer parmetros bsicos del medio poroso, como por ejemplo: la capacidad efectiva de flujo (Ko.h), el factor de dao a la formacin (S), etc. La duracin de este perodo normalmente puede ser de horas das, dependiendo fundamentalmente de la permeabilidad de la formacin productora. Dado que el diferencial de presin no se estabiliza no se considerarn ecuaciones para estimar la tasa de produccin en este estado de flujo.

Transicin entre estados de flujo

Despus del flujo transitorio este perodo ocurre una transicin hasta alcanzarse una estabilizacin pseudo-estabilizacin de la distribucin de presin dependiendo de las condiciones existentes en el borde exterior del rea de drenaje.

2) Flujo Continuo o Estacionario (Steady State Flow):

Es un tipo de flujo donde la distribucin de presin a lo largo del rea

de drenaje no cambia con tiempo, (dP/dt = 0). Se presenta cuando se estabiliza la distribucin de presin en el rea de drenaje de un pozo perteneciente a un yacimiento lo suficientemente grande, asociado a un gran acufero, de tal forma que en el borde exterior de dicha rea existe flujo para mantener constante la presin (Pws). En este perodo de flujo el diferencial de presin a travs del rea de drenaje es constante y est representado por la diferencia entre la presin en el radio externo de drenaje, Pws a una distancia re del centro del pozo, y la presin fluyente en la cara de la arena, Pwfs a una distancia rw radio del pozo; ambas presiones deben ser referidas a la misma profundidad y por lo general se utiliza el punto medio de las perforaciones caoneo. Para cada valor de este diferencial (Pws-Pwfs), tradicionalmente conocido como Draw-down, se establecer un caudal de flujo del yacimiento hacia el pozo.



20



.- Estado de flujo continuo de un lquido (Pws constante en el lmite exterior)

.- Estado de flujo semicontinuo de un lquido (Pws constante en el lmite exterior)



21



Ecuaciones de flujo para estado continuo.

A continuacin se presenta la ecuacin de Darcy para flujo radial que permite estimar la tasa de produccin de petrleo que ser capaz de aportar un rea de drenaje de forma circular hacia el pozo productor bajo condiciones de flujo continuo.

Ecuacin 1.1

dpBooKroqoaSrwreLn hKqPws

Pwfs

o .')/( .00708,0 Donde: qo = Tasa de petrleo, bn/d K = Permeabilidad absoluta promedio horizontal del rea de drenaje, md h = Espesor de la arena neta petrolfera, pies Pws = Presin del yacimiento a nivel de las perforaciones, a r=re, lpcm Pwfs = Presin de fondo fluyente al nivel de las perforaciones, a r=rw lpcm re = Radio de drenaje, pies rw = Radio del pozo, pies S = Factor de dao fsico, S>0 pozo con dao, S



22



2.1 Flujo de fluidos en el yacimiento (continuacin) Simplificaciones de la ecuacin de Darcy:

La integral de la ecuacin 1.1 puede simplificarse para yacimientos sub-saturados con presiones fluyentes en el fondo del pozo, Pwfs, mayores que la presin de burbuja, Pb. Primeramente para presiones mayores a la presin de burbuja el producto o.Bo es aproximadamente constante y por lo tanto puede salir de la integral. En segundo lugar, dado que no existe gas libre en el rea de drenaje, toda la capacidad de flujo del medio poroso estar disponible para el flujo de petrleo en presencia del agua irreductible Swi, es decir, el valor de Kro debe ser tomado de la curva de permeabilidades relativas agua-petrleo a la Swi, este valor es constante y tambin puede salir de la integral. Normalmente el trmino de turbulencia aqo solo se considera en pozos de gas donde las velocidades de flujo en las cercanas de pozo son mucho mayores que las obtenidas en pozos de petrleo. Bajo estas consideraciones la ecuacin 1.1, despus de resolver la integral y evaluar el resultado entre los lmites de integracin, quedar simplificada de la siguiente manera:

Ecuacin 1.2 SrwreLnBoo PwfsPwshKoqo )/(. .00708,0 La misma ecuacin puede obtenerse con la solucin P(r,t) de la

ecuacin de difusividad bajo ciertas condiciones iniciales y de contorno, y evalundola para r=rw. En trminos de la presin promedia en el rea de drenaje Pws, la ecuacin quedara despus de utilizar el teorema del valor medio:

Ecuacin 1.3 SrwreLnBoo PwfsPwshKoqo 5,0)/(. .00708,0

Propiedades del petrleo

Las propiedades del petrleoo y Bo se deben calcular con base al anlisis PVT, en caso de no estar disponible, se deben utilizar correlaciones empricas apropiadas. En el CD anexo se presentan, en una hoja de Excel, algunas de las correlaciones ms importantes que se utilizaran en este curso para el clculo de la solubilidad del gas en el petrleo (Rs), factor volumtrico del petrleo (Bo), la viscosidad (o) y densidad del petrleo (o) para presiones tanto por encima como por debajo de la presin de burbuja. La Tabla 2.1 muestra las correlaciones mencionadas.




Tabla 1.1 Propiedades del petrleo

23 Ing. Ricardo Maggiolo

Bo, R

s, o

y o

, pa

ra p

etr

leo

satu

rado

(P<

=

Pb).

Sta

ndin

g

Sta

ndin

g

od :

sin

gas

en s

oluc

in

o : c

on g

as e

n so

luci

n

Con

:a =

10.

715

(Rs+

100)

- 0.5

15

b =

5.44

(Rs+

150)

- 0.

338

B

eggs

& R

obin

son

Co=

Com

Bo,

o y o

, par

a pe

trle

o su

bsat

urad

o (P

>Pb)

.

pres

ibili

dad

del p

etr

leo

(apr

ox. 1

5 x

10 -6

lpc

-1)

ob y

Bob

= o

y B

o @

P=P

b

o =

1.0

008

ob +

0.00

1127

(P-P

b) (0

.038

ob 1

.59 -

0.0

0651

7 ob

1.81

48)

ob= o

@ P

=Pb

Kar

toat

mod

jo y

Schm

idt

Vict

or P

opn

(Z)

Bg

(bls/

pcn)

= 0

.005

03*Z

.T(R

) / P

(lpca

) g

(lbs/

pc) =

2.7

g . P

(lpca

)/Z.T

(R)

Fact

or Z

, Bg

y g

par

a el

gas

.

20

48

14

12

181

000

091

00

125

0.

x.

.

)(

Pg

Rs

)F

(T

.A

PI

.

lp

ca

21

251

0001

20

9759

0.

og)

F(

T.

Rs

..

Bo

.1T

1010

163

.1)

API

0202

3.0

0324

.3(.

od

b

.a

odo

.1

825

.3

g7

85.1

)R

(T

.1

0).

lpca

(P

.3

44

40

0.1

Z

e)

PbP

.(C

o.

Bob

Bo

Pb

Rs

Pb

Bo

Pb

oe

)Pb

P.(

Co

0.

ob

Bo

./

Rs

..

.g

oo

615

507

640

462

Pb

oRs Bo

Bo, R

s, o

y o

, pa

ra p

etr

leo

satu

rado

(P<

=

Pb).

Sta

ndin

g

Sta

ndin

g

od :

sin

gas

en s

oluc

in

o : c

on g

as e

n so

luci

n

Con

:a =

10.

715

(Rs+

100)

- 0.5

15

b =

5.44

(Rs+

150)

- 0.

338

B

eggs

& R

obin

son

Co=

Com

Bo,

o y o

, par

a pe

trle

o su

bsat

urad

o (P

>Pb)

.

pres

ibili

dad

del p

etr

leo

(apr

ox. 1

5 x

10 -6

lpc

-1)

ob y

Bob

= o

y B

o @

P=P

b

o =

1.0

008

ob +

0.00

1127

(P-P

b) (0

.038

ob 1

.59 -

0.0

0651

7 ob

1.81

48)

ob= o

@ P

=Pb

Kar

toat

mod

jo y

Schm

idt

Vict

or P

opn

(Z)

Bg

(bls/

pcn)

= 0

.005

03*Z

.T(R

) / P

(lpca

) g

(lbs/

pc) =

2.7

g . P

(lpca

)/Z.T

(R)

Fact

or Z

, Bg

y g

par

a el

gas

.

20

48

14

12

181

000

091

00

125

0.

x.

.

)(

Pg

Rs

)F

(T

.A

PI

.

lp

ca

21

251

0001

20

9759

0.

og)

F(

T.

Rs

..

Bo

.1T

1010

163

.1)

API

0202

3.0

0324

.3(.

od

b

.a

odo

.1

825

.3

g7

85.1

)R

(T

.1

0).

lpca

(P

.3

44

40

0.1

Z

e)

PbP

.(C

o.

Bob

Bo

Pb

Rs

Pb

Bo

Pb

oe

)Pb

P.(

Co

0.

ob

Bo

./

Rs

..

.g

oo

615

507

640

462

Pb

oBo Rs



24



2.1 Flujo de fluidos en el yacimiento (continuacin) 3) Flujo Semi-continuo (Pseudo-steady State Flow):

Es un tipo de flujo donde la distribucin de presin a lo largo del rea de drenaje cambia con tiempo pero a una tasa constante, (dP/dt = cte). Se presenta cuando se seudo-estabiliza la distribucin de presin en el rea de drenaje de un pozo perteneciente a un yacimiento finito de tal forma que en el borde exterior de dicha rea no existe flujo, bien sea porque los lmites del yacimiento constituyen los bordes del rea de drenaje o por que existen varios pozos drenando reas adyacentes entre s. Las ecuaciones homlogas a las anteriores pero bajo condiciones de flujo semicontinuo son las siguientes:

Ecuacin 1.4 SrwreLnBoo PwfsPwshKooq 5,0)/(. .00708,0 En trminos de la presin promedia en el rea de drenaje Pws, la

ecuacin quedara: Ecuacin 1.5 SrwreLnBoo PwfsPwshKoqo 75,0)/(. .00708,0 Este es el estado de flujo mas utilizado para estimar la tasa de

produccin de un pozo que produce en condiciones estables.

Uso importante de las ecuaciones

Para estimar el verdadero potencial del pozo sin dao, se podran utilizar las ecuaciones 1.2 y 1.5 asumiendo S=0 y compararlo con la produccin actual segn las pruebas, la diferencia indicara la magnitud del dao seudodao existente.

Modificacin de las ecuaciones para los casos donde la forma del rea de drenaje no sea circular:

Los pozos difcilmente drenan reas de formas geomtricas definidas, pero con ayuda del espaciamiento de pozos sobre el tope estructural, la posicin de los planos de fallas, la proporcin de las tasas de produccin de pozos vecinos, etc. se puede asignar formas de reas de drenaje de los pozos y hasta, en algunos casos, la posicin relativa del pozo en dicha rea. Para considerar la forma del rea de drenaje se sustituye en la ecuacin 1.5 el trmino Ln (re/rw)" por Ln (X) donde X se lee de la tabla 2.2 publicada por Mathews & Russel, el valor de X incluye el factor de forma desarrollado por Dietz en 1965.



25



Tabla 2.2 Factores X de Mathews & Russel



26



2.1 Flujo de fluidos en el yacimiento (continuacin) A continuacin se definen algunas relaciones importantes muy utilizadas en Ingeniera de Produccin, para representar la capacidad de aporte de fluidos del yacimiento: Indice de productividad

Se define ndice de productividad (J) a la relacin existente entre la tasa de produccin, qo, y el diferencial entre la presin del yacimiento y la presin fluyente en el fondo del pozo, (Pws- Pwf). Para el caso de completaciones a hoyo desnudo, la Pwf es igual a Pwfs, luego (Pws- Pwf)= (Pws- Pwfs) De las ecuaciones 1.2 y 1.5 se puede obtener el ndice de productividad, despejando la relacin que define al J, es decir: Para flujo continuo: Ecuacin 1.6 SrwreLnBoo hKoPwfsPws qolpcbpdJ )/(. .00708,0)/( Para flujo semi-continuo: Ecuacin 1.7 SrwreLnBoo hKoPwfsPws qolpcbpdJ 75,0)/(.. ..00708,0)/( En las relaciones anteriores la tasa es de petrleo, qo, ya que se haba asumido flujo solo de petrleo, pero en general, la tasa que se debe utilizar es la de lquido, ql, conocida tambin como tasa bruta ya que incluye el agua producida. Escala tpica de valores del ndice de productividad en bpd/lpc: Baja productividad: J < 0,5 Productividad media: 0,5 < J < 1,0 Alta Productividad : 1,0 < J < 2,0 Excelente productividad: 2,0 < J

Eficiencia de flujo (EF)

Cuando no existe dao (S=0) el ndice J reflejar la verdadera productividad del pozo y recibe el nombre de Jideal y en lo sucesivo se denotara J para diferenciarlo del ndice real J. Se define eficiencia de flujo a la relacin existente entre el ndice de productividad real y el ideal, matemticamente: EF= J/ J



27



2.1 Flujo de fluidos en el yacimiento (continuacin) IPR (Inflow Performance Relationships)

La curva IPR es la representacin grfica de las presiones fluyentes, Pwfs, y las tasas de produccin de lquido que el yacimiento puede aportar al pozo para cada una de dichas presiones. Es decir para cada Pwfs existe una tasa de produccin de lquido ql, que se puede obtener de la definicin del ndice de productividad: ql= J.(Pws- Pwfs) o tambin Pwfs = Pws - ql/ J Obsrvese que la representacin grfica de Pwfs en funcin de ql es una lnea recta en papel cartesiano. La IPR representa una foto instantnea de la capacidad de aporte del yacimiento hacia el pozo en un momento dado de su vida productiva y es normal que dicha capacidad disminuya a travs del tiempo por reduccin de la permeabilidad en la cercanas del pozo y por el aumento de la viscosidad del crudo en la medida en que se vaporizan sus fracciones livianas.

Ejercicio para ilustrar el clculo de J, EF, qo y Pwfs.

Un pozo de dimetro 12 y bajo condiciones de flujo semicontinuo drena un rea cuadrada de 160 acres de un yacimiento que tiene una presin esttica promedio de 3000 lpcm y una temperatura de 200 F, el espesor promedio del yacimiento es de 40 pies y su permeabilidad efectiva al petrleo es de 30 md. La gravedad API del petrleo es de 30 y la gravedad especifica del gas 0,7. La presin de burbuja es de 1800 lpcm y de una prueba de restauracin de presin se determin que el factor de dao es 10. Se pregunta: 1) Cul seria la tasa de produccin para una presin fluyente de 2400

lpcm? 2) El pozo es de alta, media o baja productividad? 3) Si se elimina el dao, a cuanto aumentara el ndice de

productividad? 4) Cunto es el valor de la EF de este pozo? 5) Cunto producira con la misma presin fluyente actual si se elimina

el dao? 6) Cul seria Pwfs para producir la misma tasa actual si se elimina el

dao? Nota: Utilice para las propiedades de los fluidos las correlaciones indicadas en la hoja de Correl_PVT y para el Bo con P>Pb use una compresibilidad del petrleo de 15x 10-6 lpc-1.



28



2.1 Flujo de fluidos en el yacimiento (continuacin)

Solucin : De la tabla 1.2 para un rea de drenaje cuadrada con el pozo en el centro se tiene el siguiente factor de forma: ( re/rw)= X = 0,571 A1/2/rw es decir, que el re equivalente si el rea fuese circular seria: re equiv. = 0,571 A1/2 = 0,571x (43560x160) 1/2 = 1507 pies (rea circular = 164 acres) Con el valor de la Pb se obtiene la solubilidad de gas en el petrleo Rs,utilizando la correlacin de Standing que aparece en la Tabla1.1, luego se evalan el factor volumtricoBo y la viscosidad o tanto a Pws como a Pb para luego promediarlos. Los resultados obtenidos son los siguientes: Rs = 311 pcn/bn Bo = 1,187 by/bn o = 0,959 cps Despus de obtener los valores de las propiedades se aplican la ecuacin para determinar qo, J, EF,y Pwfs. 1) 1075,0))24/25,12/(1507(187,1.959,0 1800300040.30.00708,0 Lnqo = 260 bpd 2) J = 0,433 bpd/1pc, luego es de baja productividad 3) J = 1,03 bpd/1pc 4) EF = 0,42 5) q1 = 618 bpd 6) Pwfs = 2790 1pcm



29



2.1 Flujo de fluidos en el yacimiento (continuacin) Flujo de petrleo y gas en yacimientos saturados

En yacimientos petrolferos donde la presin esttica, Pws, es menor que la presin de burbuja, Pb existe flujo de dos fases: una liquida (petrleo) y otra gaseosa (gas libre que se vaporizo del petrleo). El flujo de gas invade parte de los canales de flujo del petrleo disminuyendo la permeabilidad efectiva Ko, a continuacin se describen las ecuaciones utilizadas para obtener la IPR en caso de tener flujo bifsico en el yacimiento. La ecuacin general de Darcy establece que:

Pws

Pwfsooo dpBKrSrwreLn

Khqo ./)/(

00708,0

Asumiendo que se conoce Pws, S=0, el limite exterior es cerrado y Pws



30



2.1 Flujo de fluidos en el yacimiento (continuacin) Trabajo de Vogel

Dado un yacimiento con K, h, re, rw, curvas de permeabilidades relativas y anlisis PVT conocidos, se podran calcular para cada valor Pwfs el rea bajo la curva de Kro/o.Bo desde Pwfs hasta Pws y estimar la tasa de produccin qo con la ecuacin anterior. De esta forma en un momento de la vida productiva del yacimiento se puede calcular la IPR para yacimientos saturados. Inclusive a travs del tiempo se podra estimar como vara la forma de la curva IPR a consecuencia de la disminucin de la permeabilidad efectiva al petrleo por el aumento progresivo de la saturacin gas, en el rea de drenaje, en la medida que se agota la energa del yacimiento. Para obtener la relacin entre la presin del yacimiento y el cambio de saturacin de los fluidos es necesario utilizar las ecuaciones de balance de materiales. Este trabajo de estimar curvas IPR a distintos estados de agotamiento del yacimiento fue realizado por Vogel en 1967 basndose en las ecuaciones presentadas por Weller para yacimientos que producen por gas en solucin, lo ms importante de su trabajo fue que obtuvo una curva adimensional vlida para cualquier estado de agotamiento despus que el yacimiento se encontraba saturado sin usar informacin de la saturacin de gas y Krg. La siguiente ilustracin indica esquemticamente el trabajo de Vogel

qmax1

Pws1

(q , Pwf)

2

max PwsPwfs8.0

PwsPwfs2.0.1

qq

1.

1.q/qmax

PwfPws

qmax1qmax1

Pws1Pws1

(q , Pwf)

2

max PwsPwfs8.0

PwsPwfs2.0.1

qq

1.

1.q/qmax

PwfPws

1.

1.q/qmax

PwfPws

1.

1.q/qmax

PwfPws

1.q/qmax

PwfPws

q/qmax

PwfPwsPwfPws



31



2.1 Flujo de fluidos en el yacimiento (continuacin) Ecuacin y Curva de Vogel para yacimientos saturados

Como resultado de su trabajo Vogel public la siguiente ecuacin para considerar flujo bifsico en el yacimiento:

2

8.02.01max/

PwsPwfs

PwsPwfsqqo

La representacin grfica de la ecuacin anterior es la curva IPR adimensional presentada por Vogel, y que se muestra a continuacin:

Validez de la ecuacin de Vogel

La solucin encontrada ha sido ampliamente usada en la prediccin de curvas IPR cuando existen dos fases (lquido y gas) y trabaja razonablemente segn Vogel para pozos con porcentajes de agua hasta 30%.

Ejercicio para ilustrar el uso de la ecuacin de Vogel

Dada la siguiente informacin de un pozo que produce de un yacimiento saturado: Pws= 2400 lpc qo= 100 b/d Pwf= 1800 lpc Pb = 2400 lpc.

Calcular la tasa esperada para Pwf = 800 lpc



32



Solucin : Primero se debe resolver la ecuacin de Vogel para obtener el qomax

28.02.01

max

PwsPwf

PwsPwf

qoqo

Sustituyendo: bpdqo 250

240018008.0

240018002.01

100max2

Luego para hallar qo para Pwf = 800 lpc se sustituye Pwf en la misma ecuacin de Vogel:

bpdqo 21124008008.0

24008002.01250

2

Construccin de la IPR para Yacimientos Saturados

Para construir la IPR para yacimientos saturados se deben calcular con la ecuacin de Vogel varias qo asumiendo distintas Pwfs y luego graficar Pwfs v.s. qo. Si se desea asumir valores de qo y obtener las correspondientes Pwfs se debe utilizar el despeje de Pwfs de la ecuacin de Vogel, el cual quedara:

max/80811125.0 qoqoPwsPwfs Esta curva representa la capacidad de aporte de fluidos del yacimiento hacia el pozo en un momento dado. Como ejercicio propuesto construya la IPR correspondiente al ejercicio anterior.

La siguiente figura muestra la IPR resultante.

C U R V A S D E O F E R T A VALORES Jreal= 0,188 Jideal= 0,188 Jfutura= 0,188

ASUMIDOS EF= 1,00 EF= 1,00 EF= 1,00

Pwf / Pws ql IPR Real ql IPR Ideal ql IPR Futura0 2400 0 2400 0 2400

1,00 0 2400 0 2400 0 24000,90 43 2160 43 2160 43 21600,80 82 1920 82 1920 82 19200,70 117 1680 117 1680 117 1680 00,60 148 1440 148 1440 148 1440 0,20,50 175 1200 175 1200 175 1200 0,40,40 198 960 198 960 198 960 0,60,33 211 800 211 800 211 800 0,80,20 232 480 232 480 232 480 10,10 243 240 243 240 243 2400,00 250 0 250 0 250 0

qmax-qb= 250 qmax-qb= 250 qmax-qb= 250qmax= 250 qmax= 250 qmax= 250

CURVAS DE OFERTA EN EL FONDO DEL POZO

0

500

1000

1500

2000

2500

3000

0 50 100 150 200 250 300

ql (bpd)

Pw

f (lp

c)

IPR Real

IPR Ideal

IPR Futura

Pwf_prueba



33



2.1 Flujo de fluidos en el yacimiento (continuacin) Flujo de gas y petrleo en yacimientos sub-saturados

En yacimientos subsaturados existir flujo de una fase liquida (petrleo) para Pwfs> Pb y flujo bifsico para Pwfs < Pb. En estos casos la IPR tendr un comportamiento lineal para Pwfs mayores o iguales a Pb y un comportamiento tipo Vogel para Pwfs menores a Pb tal como se muestra en la siguiente figura.

Ntese que la tasa a Pwfs= Pb se denomina qb Ecuacin de Vogel para yacimientos subsaturados

Dado que la IPR consta de dos secciones, para cada una de ellas existen ecuaciones particulares:

En la parte recta de la IPR, q qb Pwfs Pb, se cumple: )(. PwfsPwsJq de donde, J se puede determinar de dos maneras: 1) Si se conoce una prueba de flujo (Pwfs, ql) donde la Pwfs > Pb.

)()(

pruebaPwfsPwspruebaqJ

2) Si se dispone de suficiente informacin se puede utilizar la ecuacin

de Darcy:

SrwreLnoBohKoJ 75.0/

.00708,0

Pwfs Pb

Pwfs Pb

qmax

Pws

Pb qb, Pb

qb



34



2.1 Flujo de fluidos en el yacimiento (continuacin)

En la seccin curva de la IPR, q < qb Pwfs > Pb, se cumple:

2

8,02,01maxPb

PwfsPb

Pwfsqbqqbq

)(. PbPwsJqb

8,1.max PbJqbq

La primera de las ecuaciones es la de Vogel trasladada en el eje X una distancia qb, la segunda es la ecuacin de la recta evaluada en el ltimo punto de la misma, y la tercera se obtiene igualando el ndice de productividad al valor absoluto del inverso de la derivada de la ecuacin de Vogel, en el punto (qb, Pb). Las tres ecuaciones anteriores constituyen el sistema de ecuaciones a resolver para obtener las incgnitas J, qb y qmax. Introduciendo las dos ltimas ecuaciones en la primera y despejando J se obtiene:

28,02,01

8,1 PbPwfs

PbPwfsPbPbPws

qJ

El valor de J, se obtiene con una prueba de flujo donde la Pwfs est por debajo de la presin de burbuja, una vez conocido J, se puede determinar qb y qmax quedando completamente definida la ecuacin de q la cual permitir construir la curva IPR completa. Otra manera de calcular el ndice de productividad es con la ecuacin de Darcy cuando se dispone de suficiente informacin del rea de drenaje del yacimiento. A continuacin se presentan dos ejercicios para ilustrar el uso de la ecuacin de Vogel para yacimientos subsaturados.



35



Ejercicio usando la ecuacin de Darcy

Dada la informacin de un yacimiento subsaturado: Pws = 3000 lpc h = 60 pies Pb = 2000 lpc re = 2000 pies o = 0,68 cps rw = 0,4 pies Bo = 1,2 md. Ko = 30 md. Calcular: 1.- La tasa de flujo (qb) a una Pwfs= Pb. 2.- La qmax total. 3.- La q para una Pwf = a) 2500 lpc y b) 1000 lpc Solucin: 1) Inicialmente se aplica la ecuacin de Darcy:

075.04.0/200068.02.1 200030006010)30(08.74/3/1008.7

33

LnSrwreLnBouoPwfsPwsKhqb

evaluando se obtiene dbqb /2011

Luego ...... lpcbpdPbPwsqbJ /011.2

200030002011

2) Aplicando la ecuacin de qmax en funcin de J se tiene: bpdJPbqbq 4245

8.12000011.2

20118.1

max 3.a) bdpPwfsPwsJqo 100525003000011.2

3.b) dosustituyenPb

PwfsPb

Pwfsqbqqbqo

2

8.02.01max

dbqo /3575200010008.0

200010002.01)20114245(2011

2

Si se desea obtener la curva IPR se asumen otros valores de Pwfs y se calculan sus correspondientes qo para luego graficar Pwfs vs. qo.



36



Ejercicio usando los resultados de una prueba de flujo.

Dada la informacin de un yacimiento subsaturado: Pws = 4000 lpc Pb = 3000 lpc y qo = 600 b/d para una Pwfs = 2000 lpc. Calcular: 1.- La qmax. 2.- La qo para Pwfs= 3500 lpc. 3.- La qo para Pwfs= 1000 lpc. Procedimiento: Para resolver este problema, primero se determina el ndice de productividad utilizando la solucin obtenida para J al resolver el sistema de ecuaciones para la parte curva de la IPR ya que Pws>Pb y Pwfs



37



En resumen

Para cada presin fluyente en el fondo del pozo (en la cara de la arena) el rea de drenaje del yacimiento quedar sometida a un diferencial de presin que depender de la energa del yacimiento (Pws-Pwfs), este diferencial provocar el flujo de fluidos del yacimiento hacia el pozo y la mayor o menor tasa de produccin aportada depender fundamentalmente del ndice de productividad del pozo. La IPR se considerar en lo sucesivo como una curva de oferta de energa o afluencia de fluidos que el yacimiento entrega al pozo (Pwfs v.s. q).



38



2.2 Flujo de fluidos en la completacin Descripcin La completacin representa la interfase entre el yacimiento y el pozo, y a

travs de ella el fluido sufre una prdida de presin la cual depender del tipo de completacin existente:

Tipo de completacin Ilustracin

1) Hoyo desnudo: son completaciones donde existe una comunicacin directa entre el pozo y el yacimiento, normalmente se utilizan en formaciones altamente consolidadas y naturalmente fracturadas.

2) Caoneo convencional: son completaciones donde se perfora caonea la tubera de revestimiento, el cemento y la formacin productora para crear tneles que comuniquen el pozo con el yacimiento, normalmente se utilizan en formaciones consolidadas.



39



2.2 Flujo de fluidos en la completacin (continuacin.)

3) Empaque con grava: son completaciones donde se coloca un filtro de arena de granos seleccionados (grava) por medio de una tubera ranurada para controlar la entrada de arena al pozo, normalmente se utilizan en formaciones poco consolidadas. El empaque puede realizarse con la tubera de revestimiento perforada con el hoyo desnudo.

Cada de presin en la completacin

A continuacin se presenta la manera de calcular la prdida de presin en cada tipo de completacin:

1) Cada de presin en completaciones a hoyo desnudo

En este tipo de completaciones la cada de presin es cero ya que la comunicacin entre el yacimiento y el pozo es directa, luego: Pc= Pwfs Pwf = 0 Pwf= Pwfs

2) Cada de presin en completaciones con caoneo convencional

La ecuacin presentada por Jones, Blount y Glaze puede ser utilizada para evaluar la prdida de presin a travs de la completacin con caoneo convencional.

bq + qa= Pwf-PwfsPc 2 La completacin se dice, con base a la experiencia, que no es restrictiva cuando la cada de presin a travs del caoneo est entre 200 a 300 lpc. Antes de definir los coeficientes a y b se deben describir algunas premisas establecidas por los autores.



40



2.2 Flujo de fluidos en la completacin (continuacin) Premisas para las ecuaciones de Jones, Blount y Glaze

Se ha demostrado que alrededor del tnel caoneado, durante una perforacin normal, existir siempre una zona triturada o compactada que exhibe una permeabilidad sustancialmente menor que la del yacimiento. A fin de analizar los efectos de este caoneo y su efecto restrictivo sobre la capacidad de flujo se han realizado varias suposiciones basndose en el trabajo de numerosos autores. La siguiente figura muestra que mediante un giro de perforacin de 90 el tnel caoneado puede ser tratado como un pozo miniatura sin dao.

Otras suposiciones

1. La permeabilidad de la zona triturada o compactada es: a) El 10% de la permeabilidad de la formacin, si es perforada en

condicin de sobre-balance. b) El 40% de la permeabilidad de la formacin si es perforada en

condicin de bajo-balance. Mcleod especific un rango de valores pero se trabajara con estos promedios.

2. El espesor de la zona triturada es de aproximadamente 1/2 pulgada. 3. El pequeo pozo puede ser tratado como un yacimiento infinito: es

decir, Pwfs permanece constante el lmite de la zona compacta, de este modo se eliminan el -3/4 de la ecuacin de Darcy para la condicin de flujo radial semicontinuo.



41



2.2 Flujo de fluidos en la completacin (continuacin) Ecuacin de Jones, Blount & Glaze para caoneo convencional

La ecuacin de Jones, Blount & Glaze establece que

bq + qa= Pwf-PwfsPc 2 Donde:

Lp

)rc1-

rp1( oBo10302

= a 2

214- .., y

KpLp 08

)rprc(ooB

= b..007,0

Ln con Kp

10 332= 201110

,,

(Firoozabadi y Katz, presentaron una correlacin de en funcin de K,

ver grfico en la prxima pgina)

q = tasa de flujo/perforacin, b/d/perf = factor de turbulencia, pie-1 Bo= factor volumtrico del petrleo, by/bn o = densidad del petrleo, lb/pie3 Lp = longitud del tnel caoneado, pie o = viscosidad del petrleo, cp. Kp = permeabilidad de la zona triturada, md. (Kp= 0.1 K para caoneo con sobrebalance y Kp= 0.4 K para caoneo con bajobalance) rp = radio del tnel caoneado, pie rc = radio de la zona triturada, pie



42



2.2 Flujo de fluidos en la completacin (continuacin) Ecuacin de Jones, Blount & Glaze para caoneo convencional (continuac)

Sustituyendo a y b la ecuacin de Jones, Blount & Glaze quedara:

KpLp0,00708

)rprc(oo

+ q Lp

)rc1-

rp1( oBo10302

=Pc 22

214-

..

Ln...

.....,

La informacin acerca de los caones de perforacin debe ser solicitada a la contratista de servicio quienes podran suministrar la longitud estimada de la perforacin Lp ya corregida y adaptada a las condiciones del caoneo. La grfica presentada por Firoozabadi y Katz de vs. K, es la siguiente:



43



2.2 Flujo de fluidos en la completacin (continuacin)

Ejercicio propuesto para calcular Pc en una completacin con caoneo convencional

Dada la siguiente informacin de un pozo caoneado convencionalmente: K = 5 md Pws = 3500 1pc Ty = 190F Pb = 2830 1pc re = 1500 pies h = 25 pies g = 0,65 rw = 0,36 pies Densidad de tiro = 2 tpp hoyo = 8,75 RGP = 600 pcn/bl Bo = 1,33 by/bn hp = 15 pie casing = 5-1/2" Pwh = 200 1pc o = 0,54 cp API = 35 tubera = 2-3/8" OD Perforado con sobrebalance utilizando can de casing de 4" (dimetro de la perforacin= 0,51", longitud de la perforacin = 10,6 pulg.) Determine la prdida de presin a travs de la completacin para una tasa de produccin de 100 bpd.

3) Cada de presin en completaciones con empaque con grava

La ecuacin presentada por Jones, Blount y Glaze puede ser utilizada para evaluar la prdida de presin a travs del empaque:

bq + qa= Pwf-PwfsPc 2 Al igual que en el caso anterior la completacin, con base a la experiencia, es ptima cuando la cada de presin a travs del caoneo est entre 200 a 300 lpc. Antes de definir los coeficientes a y b se deben describir algunas premisas establecidas por los autores.

Premisas para las ecuaciones de Jones, Blount y Glaze

Los fluidos viajan a travs de la formacin a la regin cercana que rodea el pozo, entran por las perforaciones de la tubera de revestimiento hacia el empaque de grava y luego pasar el interior del "liner" perforado o ranurado. Las siguientes premisas se consideran para utilizar las ecuaciones de Jones, Blount & Glaze:



44



2.2 Flujo de fluidos en la completacin (continuacin) 1) Tipo de flujo a travs del empaque: Se asume que el flujo a travs del empaque es lineal y no radial, de all que se utiliza la ecuacin de Darcy para flujo lineal.

2) Longitud lineal de flujo L: es la distancia entre la pared del liner ranurado y la pared del hoyo del pozo. En las siguientes figuras se indica la longitud L lineal del flujo a travs del empaque.

3) Permeabilidad de la grava: La grava posee una permeabilidad sustancialmente mayor que la del yacimiento, el tamao de las ranuras de la tubera liner ranurado depende de la grava utilizada y el tamao de los granos de grava debe ser seleccionado segn el tamao promedio de los granos de arena de la roca de yacimiento. Para cada tamao de grava existe un estimado de su permeabilidad suministrado por el proveedor , por ejemplo: Tamao Permeabilidad 20-40 Mesh 100.000,0 md 40-60 Mesh 45.000,0 md



45




Ecuacin de Jones, Blount & Glaze para completaciones con empaque con grava

La ecuacin de Jones, Blount & Glaze establece que bq + qa= Pwf-PwfsPc 2 Donde:

A

Lo.Bo..10089=a 2-13 ..,

y

AKg101271LBoo=b 3- ..,

.. con

K10 471= 550

g

7

,., (segn Firoozabadi y Katz)

Ntese que aqu se utiliza la ecuacin de para formaciones no consolidadas

q = Tasa de flujo, b/d Pwf = Presin fluyente en el fondo del pozo, 1pc Pwfs= Presin de fondo fluyente del pozo a nivel de la cara

de la arena, lpc = Coeficiente de turbulencia para grava, pie-1. Bo = Factor volumtrico de formacin, by/bn o = Densidad del petrleo, lbs/pie 3 L = Longitud de la trayectoria lineal de flujo, pie A = rea total abierta para flujo, pie2 (A = rea de una perforacin x densidad de tiro x

longitud del intervalo perforado). Kg = Permeabilidad de la grava, md. (Para 20-40 mesh 100 Darcies y para 40-60 mesh 45

Darcies) Sustituyendo a y b la ecuacin de Jones, Blount & Glaze quedara:

q AK10 1271

LBoo + q A

LoBo10 089= Pc g

3-2

2

2-13

...,......,



46




Ejercicio propuesto para calcular Pc

Dada la siguiente informacin de un pozo con empaque con grava: Pwh = 280 1pc Pws = 3500 1pc Dw = 8000 pies Ko = 170 md h = 25' pies re = 1500 pies hoyo = 12-1/4" revestidor = 9-5/8" "liner" = 5-1/2" OD rw = 0,51 pies tubera = 4" Tamao de grava 40-60 (45000 md) g = 0,65 API=35 T = 190F RGP = 600 pcn/bl Bo = 1,33 b/bn Densidad de tiro=4 tpp (perf 0,51") hp = 15 pies Pb = 2380 1pc o = 0,54 cps AyS= 0 % Determine: 1) La cada de presin a travs del empaque de grava para una

tasa de 500 bpd 2) Cual ser la tasa de produccin para generar una cada de

presin a travs del empaque de 200 1pc.

Nota importante

Debe recalcarse que las completaciones con empaques con grava se utilizan en formaciones no consolidadas y de all el inters en mantener suficiente rea abierta al flujo. En formaciones compactadas el inters no est solamente en el rea abierta a flujo, sino tambin en la longitud del tnel caoneado, ambas tienen sus efectos sobre la cada de presin a travs de la completacin.



47



Curva de oferta de energa o afluencia de fluidos que el yacimiento entrega en el fondo del pozo (Pwf v.s. q)

Para obtener la curva de oferta de energa en el fondo del pozo, Pwf vs ql, se le debe sustraer a la IPR para cada tasa de produccin, la cada de presin que existe a travs de la completacin, es decir: Pwf (oferta) = Pwfs - Pc donde Pc se estima por las ecuaciones sugeridas por Jones, Blount & Glaze bien sea para caoneo convencional o para empaque con grava, y Pwfs es la presin fluyente obtenidas en los clculos de la IPR. La siguiente figura muestra la grafica de Pwf y Pwfs en funcin de la tasa de produccin q.

Ilustracin

Pwf vs q, Oferta en el fondo del pozo

Pc

Pwfs vs q, Oferta en la cara de la arena

P, lpc

q, bpd



48



CAPTULO 3

Flujo Multifsico en Tuberas

El estudio del flujo multifsico en tuberas permite estimar la presin requerida en el fondo del pozo para transportar un determinado caudal de produccin hasta la estacin de flujo en la superficie. El objetivo del presente del capitulo es determinar, mediante correlaciones de flujo multifsico en tuberas (FMT), la habilidad que tiene un pozo para extraer fluidos del yacimiento.

Flujo de fluidos en el pozo y en la lnea de flujo

Durante el transporte de los fluidos desde el fondo del pozo hasta el separador en la estacin de flujo existen prdidas de energa tanto en el pozo como en la lnea de flujo en la superficie. Las fuentes de prdidas de energa provienen de los efectos gravitacionales, friccin y cambios de energa cintica. Algoritmo para calcular las prdidas de presin del fluido.

1. Determinar un perfil de temperaturas dinmicas tanto en la lnea como en el pozo.

(Ecuacin de Ramey en el pozo, por ejemplo) 2. Dividir tanto la lnea de flujo como la tubera de produccin en secciones de 200

a 500 pies de longitud. 3. Considerar el primer tramo y asignar P1= Psep y asumir un valor de P2a 4. Calcular P y T promedio para el tramo y determinar las propiedades de los

fluidos: petrleo, agua y gas. 5. Calcular el gradiente de presin dinmica (P/Z) utilizando la correlacin de

FMT mas apropiada. 6. Calcular: P = Z.[P/Z] y P2c = P1 + P; luego compararlo con P2a, si

satisface una tolerancia pre-establecida se repite el procedimiento para el resto de los intervalos hasta el fondo, de lo contrario se repiten los clculos en el mismo intervalo tomando como asumido el ltimo valor de P2 calculado

P en la lnea de flujo= Pl =

n

i iZPZ

1

.

P en el pozo = Pp =

m

i iZPZ

1

.

Donde n representa el nmero de secciones de la lnea de flujo y m representa el nmero de secciones de la tubera en el pozo.

3 12

Psep



49



Clculo de la presin requerida en el cabezal

Una vez conocida para una determinada tasa de produccin las prdidas de energa en la lnea de flujo, Pl, se puede obtener la presin requerida en el cabezal, Pwh, de la siguiente manera: Pwh = Psep + Pl

Clculo de la presin requerida en el fondo del pozo

Similarmente, una vez conocida para una determinada tasa de produccin las prdidas de energa en el pozo, Pp, se puede obtener la presin requerida en el fondo, Pwf, de la siguiente manera: Pwf = Pwh + Pp

Ecuacin general del gradiente de presin dinmica

El punto de partida de las diferentes correlaciones de FMT es la ecuacin general del gradiente de presin la cual puede escribirse de la siguiente manera (ver deduccin en el anexo A):

Grad.total (lpc/pie) = )( Z gc 2V +

d gc 2V fm +

gc

sen g 144

1= ZP

..

.

.... 22

Siendo:

gc

sen g =)elevZP(

144..

= gradiente de presin por gravedad (80-90%).

d g 2

V fm=)ZP(

cfricc. )(144

2 = gradiente de presin por friccin (5-20%).

Z g 2

V =)ZP(

cacel. )(144

2

= gradiente de presin por cambio de energa

cintica aceleracin. La componente de aceleracin es muy pequea a menos que exista una fase altamente compresible a bajas presiones (menores de 150 lpcm).

En las ecuaciones anteriores: = ngulo que forma la direccin de flujo con la horizontal, ( =0 para flujo horizontal e =90 en flujo vertical) = densidad de la mezcla multifsica, lbm/pie3 V = velocidad de la mezcla multifsica, pie/seg. g = aceleracin de la gravedad, 32,2 pie/seg2 g/g = constante para convertir lbm a lbf fm = factor de friccin de Moody, adimensional. d = dimetro interno de la tubera, pie. Es indispensable el uso de un simulador de flujo multifsico en tuberas en el computador ya que el clculo es iterativo en presin y en algunos casos ms rigurosos iterativos en temperatura y presin.



50



Correlaciones de flujo multifsico mas utilizadas en tuberas

Entre las correlaciones para flujo multifsico que cubren amplio rango de tasa de produccin y todos los tamaos tpicos de tuberas se encuentran, para flujo horizontal: Beegs & Brill, Duckler y colaboradores, Eaton y colaboradores, etc. y para flujo vertical: Hagedorn & Brown, Duns & Ros, Orkiszewski, Beggs & Brill, Ansari, etc.

3.2 Consideraciones tericas del flujo monofsico y multifsico en tuberas A continuacin se presentan algunas consideraciones tericas requeridas para comprender el clculo del flujo monofsico y multifsico en tuberas, para luego describir las correlaciones de Hagedorn & Brown y la de Beggs & Brill. Clculo del Factor de Friccin

El clculo del gradiente de presin por friccin requiere determinar el valor del factor de friccin, fm. El procedimiento requiere evaluar si el flujo es laminar o turbulento. Para ello es necesario calcular el nmero de Reynolds.

No. de Reynolds

Est definido como:

..

Re

VdN En unidades prcticas .

..0,1488

Re

VdN Donde: d = dimetro interno de la tubera, pie. V = velocidad de la mezcla multifsica, pie/seg. = densidad de la mezcla multifsica, lbm/pie3 = viscosidad del fluido Existe flujo laminar si el nmero de Reynolds es menor de 2100 en caso contrario es turbulento.



51



Consideraciones tericas del flujo multifsico en tuberas (continuacin)

n en

Laminar. trmino del gradiente de friccin (ecuacin de

rcy Weisbach).

Ecuacin de Poiseuille

Factor de fricciFlujo

Para determinar el factor de friccin en flujo laminar, se utiliza una expresin analtica derivada igualando el gradiente de presin de Poiseuille con el Da

dL32(Obtenida integrando el perfil

dPdV2

de velocidad para este tipo de flujo en bos capilares horizontales)

Combinando esta ecuacin con la componente de friccin, se tiene:

tu

Rem N

64dv

64ff

En adelante se considerara el factor de friccin de Moody con la letra f nicamente.

n en

o. as

lisas.

ara predecir el

las ecuaciones ms utilizadas en sus ngos de aplicabilidad son:

rew, Koo y McAdams5

Factor de fricciFlujo TurbulentTuber

Numerosas ecuaciones empricas han sido propuestas pfactor de friccin bajo condiciones de flujo turbulento. En el caso de tuberas lisasra D :

3000 < NRe < 3x106

lasius6

32.0ReN5.00056.0

f B

NRe < 105

,

la edad de la tubera y del medio ambiente a la cual esta xpuesta.

onsideraciones tericas del flujo multifsico en tuberas (continuacin)

25.0ReN316.0

f Como las paredes internas de una tubera no son normalmente lisases necesario utilizar ecuaciones que consideren la rugosidad de la pared interna de la tubera. En flujo turbulento, la rugosidad puede tener un efecto significativo sobre el factor de friccin. La rugosidad de la pared es una funcin del material de la tubera, del mtodo del fabricante,e

C



52



n en

o.

rugosas.

to, sino a su valor relativo al dimetro terno de la tubera, /d.

relacin para tubera de pared completamente rugosa es la siguiente:

Factor de fricciFlujo TurbulentTuberas

El anlisis dimensional sugiere que el efecto de la rugosidad no es debido a su valor absoluin El experimento de Nikuradse genera las bases para los datos del factor de friccin a partir de tuberas rugosas. Su cor

d

2Log274.1f

1

La regin donde el factor de friccin vara con el nmero de Reynolds y la rugosidad relativa es llamada la regin de transicin o

ared parcialmente rugosa.

una ecuacin emprica para describir la variacin e f en esta regin:

p Colebrook propuso d

fN7.18

d2Log274.1

f1

Re Note que para nmeros de Reynolds grandes correspondientes a flujo completamente turbulento esta ecuacin puede reducirse a la cuacin de Nikuradse.

e de un proceso e ensayo y error por lo que puede expresarse como:

e La ecuacin propuesta por Colebrook, para f requierd

2

sRec

fN7.18

d2Log274.1f

Valores de f son supuestos (fs) y luego calculado (fc), hasta que ellos se aproximen dentro de una tolerancia aceptable. El valor inicial para fs, puede ser obtenido a partir de una de las ecuaciones explicitas para

bera lisa.

onsideraciones tericas del flujo multifsico en tuberas (continuacin)

tu

C


ESP Oil_Copyright_01-P35/The document is property of ESP Oil Consultants, any reproduction is strictly f

maggiolo, r. - optimización de la producción mediante análisis nodal

Documents