Pore Bore 2.70 > ax'abx C70 Be con OF O ' ANALISIS NODAL DEL SISTEMA DE PRODUCCION Ca4 ye et NIVEL BASICO-NTERMEDIO lar 2 (a-b) = - + Ing. Richard Marquez, Ph.D ISBN 978-980-12.2581-2 Deposito Legal Ne LF06120075002073 Febrero 2012Indice General INTRODUCCION. INDICE GENERAL Capitulo! ACTIVIDADES DEL INGENIERO DE PRODUCCION.... 14 20 Prod ele peloleo i de guy 12 Actividades del Ingeniero de Produccién. 4.4.1 Seguimiento y Control. 1.1.2 Disefio de Facilidades, 1.1.3. Manejo de Informacién y Optimizacién del Sistema... 114 Estudios. 1.1.5 Elaboracién de Informes. Desarrollo de Nuevas Tecnologias. Capitulo Il COMPORTAMIENTO DE AFLUENCIA. 24 Comportamiento de Afluencia en Pozos Productores de Petréleo, 244 LeydeDary..... 24.2 Métodos Empiricos para Construir ta Curve de fc vas do Gacleale de presiom —Aftuencia PR. 22 23 Comportamiento de Afluencia en Pozos Productores de Gas. Casos Particulares. 2.3.1 JPR Compuesta. a 23.2 Prediccion de Presién de Yacimiento a partir de la JPR Compuesta.. | Capitulo ll FLUJO MULTIFASICO EN TUBERIAS. 34 3.2 Ecuacién General de Gradiente de Presién.. Curvas de Gradiente 3.2.1 Curvas de Gradiente Estatico. 3.2.1 Curvas de Gradiente Dinamico. Construccion de Curvas de Gradiente. INDIGO EIT Dep. Legal NL F061 2ORTSDIRATSINDICE GENERAL 33.1 Curvas de Gradiente en Tuberia Vertical...... 85 3.3.2 Curvas de Gradiente en Tuberia Horizontal... 129 3.33 Curvas de Gradiente en Tuberia Inclinada.... 143 Capitulo lV ANALISIS NODAL PARA POZOS DE PETROLEO 44 Nodo Solucin... ' 4.4.4. Procedimiento Solucion Peeks Selveeo% 4.2 Nodo Funcional no. nesnsee 421 Estrangulador de Superficie. retlo Rolo furciore| 4.2.2 Empaque con Grave Abierto/Entubado)...... 423 Pozos Cafioneados... CapituloV POZO DE INYECCION............. 5.1 Pozos de Inyeccién de Agua... 5.2 Pozos de Inyeccién de Gas...... Bibliografia. Nomenclatura... Anexo... fl ao ANEXO A: Permeabilidad... ANEXO B: Propiedades Fisicas de los Fluidos. ANEXO C: Problemario. Bro Aagiclor cle agua ISON SIRT ISTT Dep. 154 154 188 179 179 (Hueco 186 195 202 202 212 24 217 223 224 231 239 FLATTERINTRODUCCION Un sistema pozo-yacimiento se encuentra conformado por una serie de elementos, entre los cuales tenemos: tuberla de produccién, linea de flujo, valvulas de seguridad, perforaciones, empaques con grava, estranguladores de flujo, estaciones de separacién, equipos de levantamiento artificial, entre otros. Cada uno de estos elementos afecta la produccién de un pozo de petréleo y gas, y la determinacion del grado de influencia de estos elementos puede ser medida a partir de la técnica de analisis nodal. Esta técnica permite no solo definir los componentes 0 secciones principales que conforman estos sistemas, sino también la interrelacion que éstos poseen entre si, La aplicacién de este método a sistemas productores de hidrocarburos fue inicialmente propuesta por Gilbert en 1954, y el procedimiento consiste en elegir un nodo solucién para dividir el sistema en este punto y poder analizar la interrelacién entre ambos componentes. Mediante un andlisis nodal se puede definir el didmetro 6ptimo de estranguladores de fiujo a utilizar en pozos con problemas de arenamiento y conificacién por agua 0 gas, también se puede identificar restricciones en lineas de flujo, definir presiones de separacién, optimizar instalaciones de levantamiento artificial por inyeccién continua de gas, realizar disefos de ‘equipos de levantamiento artificial por bombeo, efectuar andlisis de sensibilidad de densidad de tiros/pie, evaluar la eficiencia de empaques con grava, definir ef tamafio de la tuberia de produccién, cuantificar el impacto en la produccién por el uso de las técnicas de estimulacién y fracturamiento hidraulico, entre otros, TENE 990-1225 Dep egal NLR GTAOOTSDODOTSLa técnica de andlisis nodal no solo se encuentra dirigido a la evaluacion mecanica de la infraestructura, También se puede evaluar el impacto que sobre la produccién de un pozo de petr6leo, por ejemplo, ocurre cuando inorementa el corte de agua y sedimentos, la relacién gas-liquido o la presencia de particulas solidas como arena. Elementos indeseables que afectan considerablemente la eficiencia de produccién de un pozo. Adicionalmente, la técnica de andlisis nodal presenta facil adaptacién a sistemas informaticos. En la actualidad, existe una moderada variedad de simuladores los cuales de una manera rapida, eficiente y segura permiten no solo evaluar el sistema pozo-yacimiento, sino también tomar decisiones inmediatas que contribuyan a corregir cualquier eventualidad que afecte la eficiencia del sistema de produccién come tal. Este libro se encuentra estructurado en cinco capitulos, de la siguiente manera: el Capitulo | trata sobre las diversas actividades y responsabilidades que el ingeniero de produccién tiene en su rea y las cuales se encuentran referidas, entre otras cosas, al seguimiento, supervisin y control de la produccién y a la busqueda del desarrollo e implementacion de nuevas tecnologias; La construccion de curvas IPR para pozos verticales de petroleo 0 gas, mediante el uso de modelos y correlaciones existentes en fa literatura, es propuesto en el Capitulo Il; El Capitulo Ill presenta las diferentes técnicas existentes para la determinacién del gradiente de presion en tuberias verticales y horizontales, de una mezcla multifasica de petroleo y gas. Ya que el movimiento de fluidos a través de tuberias se encuentra directamente relacionado con el gradiente de presién, en la industria petrolera es de sumo interés determinar el mismo para el disefio de tuberias 0 facilidades de superficie que permita transportar los fluidos producidos por un pozo hasta los tanques de almacenamiento, de una manera eficiente y répida; En el Capitulo IV se define el término de nodo solucion y se ofrece toda ‘una metodologia para el andlisis de ciertos elementos, propios de un sistema pozo-yacimiento; La técnica de anglisis nodal puede también ser utilizada en pozos inyectores de agua o gas. El RTT - ap Lea LOST DOOTSOCSOTSIneTRODUCCION Capitulo V presenta toda una metodologia de aplicacién y ofrece las técnicas basicas necesarias para el andlisis de ciertos elementos de vital importancia que contribuyen al eficiente proceso de inyeccién en pozos ya perforados o suministraré la informacién basica requerida en el disefio de nuevos pozos.CAPITULO I ACTIVIDADES DEL INGENIERO DE PRODUCCION La misién de un Ingeniero de Produccién de Hidrocarburos es bastante amplia y compleja, ya que es responsable de realizar diversas actividades, tales como: Disefio y optimizacion de equipos de levantamiento artificial, Seguimiento y control del potencial de produccién; Elaboracién de estudios, informes y programas disponibles; Entre otros. También, todo ingeniero debe conocer criterios y conceptos claves sobre el potencial de produccién, declinacion y contribuciones, definiciones éstas muy utilizadas en la estimacién del maximo potencial de un pozo 0 grupo de pozos. No conforme con esto, debe conocer el desarrollo y resultado del uso de de nuevas tecnologias en el area, las cuales contribuirian al mantenimiento o incremento del potencial de un pozo © campo. El siguiente capitulo presenta un resumen de las diversas actividades llevadas a cabo por el ingeniero de produccién de hidrocarburos en la industria. Constituye una simple guia, cuyo objetivo no es mas que contribuir al fortalecimiento técnico-supervisorio de los futuros ingenieros. 4.1 Actividades del Ingeniero de Produccién 4.1.1 Seguimiento y Contro! La principal responsabilidad de un ingeniero de produccién es realizar un continuo seguimiento del comportamiento de produccién de un pozo o grupo de pozos, con el objeto de verificar si el pozo se comporta como Io esperado o si él sistema pozo-yacimiento requiere ser optimizado, Para lograr estos dos ONT 7 ‘ae Legal 10607 S00R0TS[ACTIVIDADES DEL INGENTERO DE FRODUCCION CAPTULOT objetivos, el ingeniero de produccién debe entender no solo e! comportamiento mismo de los pozos, sino también la interaccién que el mismo posee con otras importantes disciplinas envueltas en la produccién de petréleo, entre las cuales se puede destacar: Perforacién; Yacimientos; Completacion; Pruebas de presiones; Métodos de estimulacién; Flujo multifasico en tuberias; Facilidades de superficies; Entre otros. Ademas, resulta necesario realizar un continu seguimiento y control de produccién de sélidos o agua de formacién, elaborar cierres de produccién y reportarlos a los organismos correspondientes, elaborar y supervisar trabajos de mantenimiento en plantas, reportar la produccién del gas de formacién, entre otros. El seguimiento al potencial de un pozo 0 grupo de pozos conlleva entonces al conocimiento y uso de ciertos conceptos y definiciones propias de esta actividad, las cuales se detallan a continuacién: A-. Maximo Potencial de Produecién MPT Se refiere al nivel maximo de produccién que puede ser sustentado consistentemente en pozos conectados a facilidades de produccién, bajo condiciones normales de operacién, cumpliendo con las normas ambientales y de ufllizacion de gas vigente. Este es un ejercicio que se realiza mensualmente, a efecto de disponer de una cifra de maximo potencial de cada una de las segregaciones de crudo de una empresa, y la cual sirve como nivel de referencia para el cdlculo de la produccién disponible de los proximos meses. El valor de! MPT es discutido y acordado por un equipo de trabajo, integrado, ademas del ingeniero de produccién, por los ingenieros de yacimiento y facilidades de superficie. En el maximo potencial de produccién MPT debe incluirse todos los pozos capaces de producir por un periodo de tiempo determinado, es decir: Pozos activos; Pozos cerrados por problemas operacionales menores de superficie 0 subsuelo, con disponibilidad inmediata; Pozos cerrados por problemas operacionales mayores 0 menores de superficie 0 subsuelo, con disponibilidad no inmediata (siempre y cuando se tenga planificado repararlos o activarios en los tres meses siguientes). Ban OTE DO EDRT Ba Legal LRT SOORTTT_ ACTIVIDADES DEL INGENIERO DEMHODUCCION * CAPTULDT La revision mensual del MPT es un ejercicio post-mortem, En general, se trata de un balance mensual de pérdidas y ganancias de potencial, que se ° puode presenar mediante la siguiente ecuAcion, 60 ge Ao ON ae RO dF CESS re i apr l MPI Pendidas + Ganancias, (ty donde el MPY,,, representa e! maximo potencial al final del periodo en estudio (a fin de mes). MPT,., representa el maximo potencial a comienzos de afio. La variable “Pérdidas" considera todas aquellas acumuladas hasta la fecha. La variable “Ganancias’ es la sumatoria de las ganancias de potencial acumuladas hasta la fecha. Este ejercicio contable, a pesar de ser de frecuencia mensual, se realiza entonces con valores acumulados de pérdidas y ganancias hasta la fecha. Por esta raz6n, se parte de un valor inicial del MPT a comienzos de ano, ‘como base para obtener el MPT, La estimacin mensual de! MPT es un ejercicio que se realiza a nivel de pozo, a final de cada mes, y sobre la base de las ultimas medidas representativas que se dispongan. Durante esta revisién podria ocurrir cualquiera de los siguientes eventos. a. Las ultimas medidas muestran que el pozo se encuentra produciendo en linea con el potencial asignado: En este caso, simplemente se ratifica su MPT para ese periodo b. Las ultimas medidas muestran que el pozo est produciendo por debajo del potencial: En este caso, se debe revisar todos los parametros asociados al sistema de produccion (THP, CHP, %dyS, Reductor) y se asume que el potencial es recuperable, por lo que ia produccién faltante se mantiene como diferida y la duracion de este proceso dependera de las acciones acordadas y la capacidad de respuesta de la organizacién. Si agotadas todas las acciones econémicamente posibles y revisados los parametros de yacimiento asociados al pozo no permiten recuperar los niveles de produccién fijados, entonces se BONES Dep Legal LPUGTROTSOTSACTIVIDADES DEL INGENIERO DE PRODUCCION. ~ CARTULOT debe tomar la decision de ajustar el valor del MPT del pozo. La reduccién del maximo potencial de produccién del pozo podra asociarse a factores intrinsicos, ‘como: Yacimiento (debido a agotamiento de presién, incremento del corte de ‘agua 0 gas, deposicion de asféllenos 0 migracién de finos que contribuyan a la reduccién de la permeabilidad del yacimiento, entre otros); Completacion mecénica (lo que ocasionaria una pérdida permanente de la capacidad de produccion y cuya correccién requeriria un trabajo de reacondicionamiento mayor con taladro); Restricciones, impuestas por razones de control de yacimiento (pozos con alta RGP, etc.) y uso de estranguladores de flujo (control de arena, agua 0 gas); Ajustes administrativos, asociadas a la sinceracion de potenciales mal asignados por algiin tipo de error de interpretacion 0 problemas de medici6n; Entre otros. c-. Las tltimas medidas muestran que el pozo se encuentran produciendo por encima del potencial asignado: Entonces, se desprenden acciones dirigidas ala verificacién de los datos y si se estima procedente, se ajustaré el potencial, elevando el mismo, Estos casos ocurren generaimente cuando al pozo se le ha efectuado algun ajuste u optimizacién del gas de inyecci6n, algun cambio de reductor, cambio de zona, estimulaci6n de pozos, cafioneo adicional, entre otros. B+. _ Pérdidas Operacionales Normales o Diferida Ordinaria DO El valor del MPT es dificil de alcanzar en la practica puesto que la operacion normal de operacién es un tanto diferente 2 las condiciones de prueba. En términos reales, se tienen todos los pozos fluyendo a la estacién en forma conjunta, bajo un ambiente de interrupciones ordinarias por razones inherentes al proceso de produccion y en presencia de las ineficiencias naturales de la infraestructura de produccién, como filtraciones, mermas, etc. Debido a sto, la produccién obtenida en los puntos de medicién y fiscalizacion es menor que la sumatoria de los MPT de cada pozo. Esta diferencia entre produccién operada y MPT es lo que se conoce come “Produccién Diferida” 0 simplemente “Diferida” DO. TREN TRL Dep Lega CALI SOOTSORIOTSACTIVIDADES DEL INGENIFRO DE PRODLUCCION, " Las causas que ocasionan la produccién diferida son variadas y pueden ser clasificadas, por lo general, en dos renglones bésicos: Diferida Extraordinaria, cuando parte de la produccién faltante se encuentra asociada a causas extraordinarias, como paros de mantenimiento programados o siniestros; Diferida Ordinaria, originadas tanto por causas conocidas y cuantificables (como, interrupciones por trabajos de guaya, pruebas especiales, cierres de valvulas de seguridad, paras repetitives no programados, etc.), como por causas intangibles y no cuantificables (como, mermas en estaciones y oleoductos, efectos climatolégicos, variaciones en los volimenes de gas para uso en instalaciones de levantamiento artificial por gas LAG, cambios en las presiones de separacién, entre otros). C-. Potencial de Produccién PT Ya que el valor del MPT puede que no sea materializable en la practica, PT, el cual se establece entonces el concepto de “Potencial de Producciér representa el volumen de petréleo 0 potencial producible de un campo, o segregacion de crudo, en un momento dado. Sobre este valor de potencial producible es que se establecen los compromisos comerciales de produccion Su valor se obtiene en base a los valores estadisticos de diferida ordinaria observados en cada campo en particular, el cual a su vez es estimado sobre la base a las diferencias obtenidas entre la produccién operada en los patios de tanques y el MPT de los pozos que se encuentren produciendo en ese momento, restandole luego la diferida que por razones extraordinarias haya cocurrido en ese periodo. En ocasiones se extrapolan los valores de PT de un campo o segregacién, a nivel de pozo, mediante el uso de la siguiente ecuacién matematica: a pokacul Ob Produces joy eral Be. MPT xy (12) MPT J ape| | ONTETRTIST Fomox golercad fe peodueese™ACTIVIDADES DEL INGENIERO DE PRODUCCION D-. _ Produccién Disponible PD La produccién disponible PD se encuentra formada por la “Produccién Operada’ PO (referida a la produccién manejada en los patios de tanque y asociada a los pozos activos en un periodo dado) y la preduccién cerrada por restricciones de mercado CM (referida al volumen de petroleo que deja de producirse por razones de mercado 0 por estrategias de producsién, y forma parte de la produccién inmediatamente disponible). El valor de PD se diferencia de PT solo cuando haya ocurrido alguna diferida extraordinaria durante el mismo periodo. 1.1.2 Disefio de Facilidades El ingeniero de produccién debe poseer la destreza para seleccionar el diémetro 6ptimo del estrangulador, a fin de controlar la produccion de fluidos 0 sélidos indeseables e incrementar el recobro. Ademés, resulta imprescindible que este profesional conozca a detalle el disefio de cualquier equipo de levantamiento artificial, de manera que pueda interpretar los resultados obtenidos mediante el uso de algun simulador comercial disponible, asi como también resulta recomendable, siempre y cuando el tiempo se lo permita, en que participe y supervise la instalacién del equipo de fondo de algun método de levan‘amiento artificial. También, debe participar en el disefio y tendido de las lineas de flujo en superficie. Es responsabilidad del ingeniero de produccién visualizar las facilidades disponibles y hacer algunas recomendaciones a tiempo, sobre la base de algiin plan de perforacion/reparacion que la empresa plantea en el mediano y largo plazo. La Fig. 1.1 muestra una facilidad de superficie tipica encontrada en el sector petrolero. 4.1.3 Manejo de Informacién y Optimizacién del Sistema El ingeniero de produccin debe hacer seguimiento a las carpetas de produccién, en donde se recopila toda la informacién oficial referente a las actividades operacionales realizadas en el pozo. Esta informacién se encuentra plasmada en los reportes de operaciones diarias de perforacién, completacion original, servicios, trabajos de reacondicionamiento, estimulaciones, entre otros Dep- Legal WTACTIVIDADES DEL INGENIERO DE PRODUCCION 8 En la actualidad, se cuenta con sistemas en linea que permiten, por ejemplo consultar la produccién de un pozo 0 grupo de pozos, hacer seguimiento de las operaciones de subsuelo (cambios de valvulas/zonas, HUD, etc.), relacionar, validar y analizar datos (dependientes del tiempo) de yacimiento y produccion, entre otros, por lo que el ingeniero debe hacer uso de estas de manera permanente. Figura 1.1. Faclidades de Superficie Tipicas en la Industria Petrolera 1.1.4 Estudios Mediante técnicas como por ejemplo el andlisis nodal, el ingeniero de produccién puede analizar en detalle, mediante el uso de simuladores, un sistema pozo-yacimiento. De esta manera, el ingeniero puede: Definir el sistema de levantamiento artificial adecuado para el pozo; Determinar posibles cuellos de botellas presentes en el sistema en estudio; Definir el indice de productividad del ozo; Optimizar el gas de levantamiento; Disefiar equipos de levantamiento artificial mediante bombas de subsuelo, tales como: bombeo electrosumergible_ BES, bombeo por cavidad progresiva BCP, Bomibeo Mecdnico BM, Bombeo Hidrdulico BH, Entre otros. La Fig. 1.2 muestra el equipo de superficie de un sistema de levantamiento artificial por bombeo electrosumergible TaN OTE ETE ‘Dep Lean LFOGIDO0TSOORTTSACTIVIDADES DEL INGENIERO DE PRODUCCION CAPITULOT Figura 1.2. Equipo de Superficie de un Sistema de Levantamiento Arti 4.1.5 Elaboracién de Informes La elaboracion de informes es otra de las responsabilidades del ingeniero de produccién. Algunas veces, resulta necesario reportar las actividades desarrolladas en clertos pozos considerados como especiales, 0 simplemente elaborar un informe sobre el nivel de actividades ejecutadas durante cierta actividad de campo. 1.2. Desarrollo de Nuevas Tecnologias Un profesional responsable debe estar atento al desarrollo de nuevas y avanzadas tecnologias, determinando espeoificamente las ventajas y limitaciones de su uso, asi como también como estas tecnologias son o pueden ser afectadas por el comportamiento futuro del pozo o campo, y como estas tecnologias afectarian los costos y ganancias de la empresa durante el periodo de vida util del sistema. Ademds, un ingeniero de produccién debe entender que la busqueda de nuevos yacimientos de hidrocarburos _conlleva TaN TRIE Dep Legal LOG RATTIENTSACTIVIDADES DEL INGENIEKO DE PRODUCCION iB CARTULOT irremediablemente al establecimiento de nuevas fronteras, especialmente en el rea de levantamiento artificial. En los ullimos afios, la producsién en pozos costa afuera (aguas profundas) ha creado una serie de nuevos métodos de levantamiento artificial. Métodos que consideran el uso de equipo de levantamiento aguas arriba de! cabezal del pozo, el cual se encuentra ubicado (algunas veces) en el fondo del lecho marino. La versatilidad de estos métodos de levantamiento se hace elocuente en el uso de sistemas informaticos en linea, que facilitan su operatividad sin necesidad de estar presente en el area y permiten regular fa produccién, controlar sus datos, cuidar y supervisar el comportamiento de la infraestructura subsuelo/superficie de manera automatica y sin la supervision de un gran ndmero de profesionales en el area. La Fig. 1.3 muestra graficamente el uso de nuevas tecnologias para la produccién de hidrocarburos en campos ubicados en aguas profundas aN OTE: Dap Lage LFEPACAPITULOI COMPORTAMIENTO DE AFLUENCIA EI desarrollo de herramientas que permitian medir la presion en el fondo del pozo, a finales de la década de 1920, permitieron al ingeniero de campo medir la produccién del pozo en superficie y correlacionarla con la presion medida en el fondo de! mismo. Desde entonces, los esfuerzos se concentraron en obtener una expresién matematica que permitiese estimar la tasa de produccién, como una funcién de la presién de fondo fluyente, y que la misma fuese valida para un amplio rango de condiciones operacionales. Cominmente, la expresion “Inflow Performance Relationship IPR ( en espafiol Curva de Comportamiento de Afluencia) es utilizada para definir la relacién que existe entre la tasa de produccién de petréleo en superficie y su correspondiente presion de fondo fluyente. En el caso de pozos productores de gas, generalmente los ingenieros utilizan la expresién “Backpressure Curve" (Curva de Contrapresién) para referirse a la misma relacion pero para una tasa de gas en superficie. Indistintamente para pozos de petréleo 0 gas, resulta claro que la produccién se encuentra directamente asociada con una presién de fondo fluyente P,,, la cual se encuentra referida en la mayoria de los casos a la profundidad promedio de las perforaciones 0 cara de la arena productora. La construccién de curvas IPR (para pozos de petréleo o gas) son de relevante importancia en el andlisis de sistemas de produccién pozo-yacimiento. Este capitulo presenta solo algunas de las diferentes metodologias disponibles en la TINDER DIET 7 ‘Dep. Legal N° LFOGTRDISOOROTS_COMPORTAMIENTO DE AFLUENCI : TTULO literatura para fa elaboracion de las curvas de comportamiento de afluencia en pozos verticales solamente, completados en una arena o grupos de arenas y con produccion simultanea de petréleo, gas y agua 2.4 Comportamiento de Afluencia en Pozos Productores de Petréleo 2.1.1 Ley de Darcy (1856) Para estimar la caida de presién que oourre a través del medio poroso fue necesario desarrollar una expresion matematica que expresase las pérdidas de energia o presién, debido a fuerzas de corte viscose o friccién, como una funcion de la velocidad o tasa de flujo. Aunque la forma de la ecuacién puede ser bastante diferente dependiendo del tipo de fluido, la ecuacion base sobre la cual la mayoria de los modelos se sustenta es universalmente conocida como la Ley de Darcy. En 1856, mientras realizaba experimentos para el disefio de un filtro de arena para la purificacion del agua, Darcy establecié que para una tasa de flujo dada, la velocidad de flujo v era directamente proporcional a la caida de presién a través del filtro mismo, De acuerdo a la ley de Darcy, v se encuentra definida como: (2.1) La Ec. 2.1 puede también ser expresada en términos de tasa de flujo volumétrico g, como: {luxe yack al 4 en donde, v es la velocidad aparente del fluido. & y u representan la permeabilidad del medio poroso y la viscosidad del fiuido, respectivamente, El TaN FREI Dep Legal LPVGIDOOTSONDOTSCCOMPORTAMIENTO DE AFLUENCIA i término dp/dl. no es mas que el gradiente de presi6n en la direccién de flujo y es considerado negative ‘Aunque los experimentos fueron tlevados a cabo solo en la direccién vertical descendente, la expresion también resulté valida para flujo horizontal: condicién de particular interés en la industria petrolera mundial. Es importante resaltar que los experimentos considerados por Darcy involucraban el uso de un nico fluido (agua) y que el filtro de arena se encontraba saturado completamente por este fluido. Por lo tanto, los efectos de la densidad y viscosidad del fluide sobre el comportamiento de flujo no fueron considerados. La nica variacién considerada en los experimentos de Darcy fue el cambio del tipo de empaque de arena, el cual tenia el efecto de alterar el valor de la permeabilidad k Sobre la base de la forma geométrica presentada en Fig. 2.1 y asumiendo flujo lineal y que la permeabilidad, viscosidad y tasa de flujo volumétrica son independiente de presién, la Ec. 2.2 puede ser integrada para obtener el siguiente resultado: _CkA arr (2.3) BR PR Figura 2.1. Geometria Considerada para Flujo Lineal SEN FH ORIE ITT Ta Legal CF UGITOOSOOROTS‘COMPORTAMIENTO DE AFLUENCIA CARTULO La constante C de la Ec. 2.3 representa el factor de conversion: 1.0 para unidades Darcy y 1.127107 para unidades de campo. La Tabla 2.1 presenta las unidades coménmente utilizadas por la ley de Darcy. ‘Tabla 2.1. Unidades Comtinmente Utilizadas por la Ley de Darcy. Unidades Vari Simbolo Darey campo Tasa de Flujo om? Seq 8D Permeabilidad Dareys md Area om? Presiin atm Lpe Viscosided cos Longitud om pies Por razones practicas, la Ec. 2.2 ha sido adaptada a condiciones de fiujo radial, similares a las mostradas en Fig. 2.2, en donde el fluido presente en el yacimiento converge en el fondo de! pozo. En este caso, el area abierta al flujo no podré ser considerada constante y por lo tanto debe ser incluida en la integral de la Ee. 2.2 Figura 2.2. Geometria Considerada para Flujo Radial EAN T78 980 TIT ‘Lega LOI 2OOTSOOROTSCCOMPORTAMIENTO DE AI ~~ CARLTULOT JENCIA 2» Bajo condiciones de flujo monofasico, generalmente se asume que el petrdleo es ligeramente compresible. En consecuencia, pequefios cambios de q con presién podria ser manejado sencillamente por el factor de volumétrico de formacion B. Por otra parte, también podria asumirse que el término k/(j: B) es independiente de la presién. En consecuencia, la integracion de la Ec. 2.2, bajo las condiciones anteriormente mencionadas, permitiria definir la siguiente expresién matemética (valida para una fase cualquiera & y en unidades de Pz Res. Pramedso yore nes ford Pluye fre taco ale be 107 kh P= Pa w zyadeo A 7.08 10" , h ( 2 wy 3 (24) mB WGI Xen sadco de dienal en donde, g representa la tasa de flujo, en BD. k y h representan la permeabilidad efectiva y el espesor del yacimiento, en mD y pies, respectivamente. uy B son la viscosidad y el factor volumétrico de la formacién, en cps y Bbl/ Bn, respectivamente. r, y 7, son el radio de drenaje del yacimiento y del pozo, respectivamente, expresado en pies. P, y P., son la presién promedio del yacimiento y Ic presién de fondo fluyente frente a la cara de la arena, respectivamente, expresado en ipca. La Ec. 2.4 es valida para condiciones de estado estable y flujo laminar, en un pozo ubicado en el centro de un area de drenaje circular como el mostrado en Fig. 2.2. La permeabilidad & es una propiedad del medio poroso y representa una idad de medida de la capacidad que el medio posee para permitir la transmisi fluidos a través de éste. Por simple analogia con conductores eléctricos, la permeabilidad representa el reciproco de la resistencia que el medio poroso ofrece a que fluyan los fluidos. La unidad de medicién es el Darcy, la cual es la conductividad de la roca a un fluido que la satura en un 100%, que tiene una viscosidad de Icp y que fluye en un régimen laminar con una velocidad de TaN DED ERTT Dep Lesa WEFT SOOTSOOBDISCAPiTULOT lem/seg, bajo un gradiente de presion de lam/cm. Los valores de permeabilidad & de una formacion pueden obtenerse mediante pruebas de laboratorio a nuicleos convencionales 0 muestras de pared, aunque los valores ‘obtenidos de estos ultimos pueden resultar inexactos y por lo general no corresponden a los valores obtenidos de niicleos convencionales. Algunos intentos se han realizado para obtener valores de permeabilidad a partir de registros eléctricos, pero hasta el presente ninguno- de estos dispositivos de medicién ha ofrecido resultados aceptables. £1 Anexo A muestra en detalle las diferentes formas existentes para estimar la permeabilidad de una formacién, bajo condiciones de fiujo lineal y radial, a través de fracturas, entre otros El espesor de la zona productora puede ser estimado a partir de registros eléctricos 0, en algunos casos, de registros de perforacién y nucleos convencionales. La aplicacién de la ley de Darey considera el valor de / igual al espesor vertical de la arena productora. En aquellas arenas donde claramente se pueda diferenciar zonas contentivas de agua, el valor de h estara referido solamente al espesor de la arena contentiva de hidrocarburo. Por otra parte, la presién promedio del yacimiento P, generalmente es conocida. En caso contrario, su valor puede ser estimado mediante pruebas de restauracién de presién. Algunas veces, su valor también puede inferirse a partir de la medicién del nivel estatico del fluid almacenado en el espacio anular del pozo. En cuanto a los valores de viscosidad 1 y factor volumétrico B de los fluidos producides, éstos pueden ser obtenidos a partir de un anélisis de presion, Volumen y temperatura PVT realizada en un laboratorio a los fluidos producidos por el pozo. Sin embargo, algunas veces no es posible obtener muestras representativas de estos fluidos para su estudio respectivo, por lo que se hace necesario estimar algunas propiedades fisicas de los fluidos por analogias o mediante el uso de correlaciones empiricas. En el Anexo B, se presentan algunas correlaciones tipicas, disponibles en la literatura, para estimar ciertas TSN 978 990-122 ep. Legal LOSTS0GTSOOEOTS2 CAPTULOT propiedades de los fluidos, como la viscosidad y el factor volumétrico, entre algunas otras, Determinar exactamente el radio de drenaje r, del yacimiento es definitivamente una tarea dificil de realizar. Generalmente, se posee un estimado de su valor real, el cual puede depender del espaciamiento entre pozos y de la forma geométrica considerada en el yacimiento. Sin embargo y sobre la base de la Ec. 2.4, cualquier error que se cometa en la determinacién del mismo es amortiguado por el logaritmo natural de r, /r,, tal como puede apreciarse en Tabla 2.2. Tabla 2.2. Efecto del Logaritmo Natural Sobre la Relacién , /7, . Logaritmo Natural (r,/tw) ry = 0.5 Pie 1, (Pies) Ln (rfrw) 500 694 1000 7.60 2000 8.29 5000 ei 0000 2.80) , representa el radio del pozo y se encuentra referido Unicamente al diametro del hoyo perforado. Su valor puede ser obtenido directamente del registro del caliper. En caso de que no se disponga del registro, el tamafio de la barrena de perforacién puede asumirse perfectamente como el valor de Sobre la base de la Ec. 2.4, puede notarse claramente que la tasa volumétrica de flujo q es directamente proporcional a la caida de presién en el yacimiento (P, - P,,). Por lo tanto, la representacién grafica de Ec. 2.4 define una linea recta, tal como puede apreciarse en Fig. 2.3. TaN ORS Depa LF TOOOTSOOBOTSCCOMPORTAMIENTO DE AFLUENCIA CARTULOT Curva de Oferta IPR Presién Fluyente = Tasa de Flujo Figura 23. Curva de Oferta IPR, La constante de proporcionalidad, conocida comimmente como Indice de Productividad J, define la relacién entre la tasa de flujo a la caida de presién en el yacimiento. En otras palabras, J se encuentra dado por: yates Soovesaburado 7.0810 ky hh ar Hy By I(r, 1r) (P, — Pag) ey +B Sobre la base de la Fig. 2.3 0 la Ec. 2.5, se pueden establecer algunas importantes observaciones: El indice de productividad J representa el inverso de la pendiente de la recta Cuando la presién de fondo fluyente P,, es igual a la presi6n promedio del yacimiento P,, no existird afluencia de fluidos desde el yacimiento OTTO ‘Dep aga TRGIDOOTSOGBOTCOMPORTAMIENTO DE AFLUENCIA, u hasta el fondo del pozo, es decir q, =0, ya que no existiré algin diferencial de presién a través de la arena productora, ¥ La maxima tasa de flujo g,_ comesponderd al valor de P,, =0. En la practica, esta condicién jamas podrd ser generada ya que P,, debe ser igual a cero, justo frente a la cara de la arena. Sin embargo, su valor es una referencia obligada en la industria petrolera, particularmente al comparar el potencial de diferentes pozos en una misma area Y La experiencia de campo ha comprobado que bajo condiciones de yacimiento sobresaturado, la Ec. 2.5 aplica perfectamente. Y El Indice de productividad es un concepto muy util para describir el potencial de! pozo, ya que combina simulténeamente las propiedades del fluido y roca, asi como las formas geométricas del yacimiento, dentro de una simple constante, Generalmente, el indice de productividad J viene expresado en BD/ipca Por otra parte, resulta importante resaltar que el concepto de flujo estable involucra una constante de productividad en particular, lo que asume la i-estado estable. Es decir, se asume que el volumen entero de condicién de ser drenaje de un pozo contribuye totalmente a la produccién de! mismo. A-. Factores que Afectan el indice de Productividad ‘Se ha comprobado que el indice de productividad J cambia con tiempo, a medida que se desplaza 0 produce los fluidos confinados en el medio poroso de un pozo. Si J cambia, la pendiente de una curva JPR como la mostrada en Fig 2.3 cambiard y por lo tanto no existira jams una relacién lineal entre la Py y En secciones siguientes, se presentard un analisis de los principales factores que afectan a J y sus consecuencias sobre la prediccién de la curva de oferta de un pozo, particularmente bajo condiciones por debajo de la presion de burbuja (yacimiento saturado). Sobresavreda = es Zh sqiuteda amen TaN TRIE IT Dep Legal TFCOMPORTAMIENTO DE AFLUENCTA CAPTULOT Mecanismos de Empuje Representan la fuente de energia natural, responsables del movimiento de fluidos dentro del medio poroso, y por ende, a través del sistema de produccién, Entre los mecanismos de empuje mas comunes, se encuentran: Empuje por Gas en Solucisn, Capa de Gas y Empuje Hidrdulico EI empuje por gas en solucién asume que no existen cambios en el tamafo inicial del yacimiento y que no ocurre intrusion de agua. Para presiones por encima de la presin de burbuja, no existird la presencia de gas libre y la fuente de material que reemplaza los fluidos producidos es la expansién de los fluidos remanentes en el yacimiento. Bajo la accién de este mecanismo, la presién del yacimiento declinara rapidamente hasta alcanzar la presion de burbuja, debido a que el petréleo tinicamente se expande para reemplazar los fluidos producidos. Por debajo de la presién de burbuja, el gas libre también se expanderd y la presién del yacimiento disminuiré pero a menor velocidad. Bajo esta condicién, se produciré una mezcla multifasica donde el gas es producido junto al petroleo y no pasaré en algtin momento a formar parte de alguna capa de gas en el yacimiento mismo. En consecuencia, la produccién de petréleo es el resultado de una expansién volumétrica del gas en solucién y una expulsion volumétrica del petroleo, Como puede apreciarse en Fig. 2.4, tan pronto como la saturacién del gas exceda la saturacion de gas critica, la RGP se incrementaré y la presion disminuira aceleradamente. En el empuie por capa de gas, el petréleo se encuentra saturado por gas bajo condiciones iniciales de presién y por lo tanto existira la presencia de gas libre en el yacimiento. A medida que el petréleo sea producido, la capa de gas se expendera y ayudaré 2 mantener la presién del yacimiento. La presion del yacimiento disminuira en menor proporcién que la presién de un yacimiento bajo empuje por gas en solucién, pero a medida que la capa de gas libre se expanda, algunos de los pozos que se encuentren én la parte més alta de la estructura producirén con una alta RGP. La Fig. 2.5 muestra los cambios de presin, RGP, y %AyS que tradicionalmente ocurre en yacimientos de este tipo. TEEN SRDS Dap Legal FORT RIOTOODOTSCCOMPORTAMIPNTO DE_AFLUENCIA “CAPITULO Figura 2.4. Empuje por Gas en Solucién. PRESION Figura 2.5. Empuje por Capa de Gas. y, 26 BN TORE IST Dep Lag LRIGTIOOTSOONOTCOMPORTAMIENTO DE AFLUENCIA ‘CAPITULO En un empuje hidraulico, el petréleo y gas es desplazado por la expansion del agua solamente. De hecho, si el yacimiento posee un empuje hidréulico muy activo, la declinacion de presién puede resultar muy pequefia, incluso despreciable, por lo que J permanece constante durante la vida productiva del pozo. La Fig, 2.6, muestra el comportamiento de presién, RGL y %4yS de un yacimiento sometido a un empulje hidraulico. Figura 2.6. Empuje Hidraulico. Existen algunos yacimientos que pueden tener la combinacién simulténea de varios mecanismos de empuje durante su vida productiva. Predecir el recobro esperado y el comportamiento a futuro del yacimiento por la acci6n combinada de estos mecanismos puede resultar dificil, en especial, por la amplia variacion de tamafios de capa de gas y aculferos presentes en los mismos. Definitivamente, no hay sustituto que pueda reemplazar la toma de informacion al inicio de la explotacion de un yacimiento, si se desea realizar una prediccion aN TITAS Bap La LF ROOTSOOROTS__ COMPORTAMIENTO DE AFLUENCTA 2s CAPTTULO t acertada a futuro de su comportamiento. Prediccién que resulta recomendable ya que la variacién de presién durante la vida productiva del yacimiento afectara las propiedades de roca y fluidos, lo que se traduce en cambios substanciales del indice de productividad /. Como puede apreciarse en Fig. 2.7, J se vera afectado en menor o mayor proporcién, dependiendo del mecanismo de empuje presente en el yacimiento. Presién de Burbule al ‘Empuje Hidraulico Empuje Capa do Gas Indice de Productividad J Empule Gas en Soluclén Produccién Acumulada N, Figura 2.7. Empuje Hidrdulico b. Comportamiento de Fases Un hidrocarburo almacenado en un yacimiento puede existir como una simple fase 0 como una mezcla multifasica dependiente de la presion, temperatura y composici6n. La Fig. 2.8 muestra un diagrama de fase tipico de un sistema multicomponente de hidrocarburos. La forma de la envolvente de fases y los rangos de presion y temperatura para los cuales son generados, varian ampliamente con composicién. La linea superior de la envolvente de THON SR9RT ETT Dap Legal LF DGI ROOT SORES‘COMPORTAMIENTO DE AFLUENCIA » fases que se ubica a la izquierda del punto critico y se conecta con éste es denominada linea de puntos de burbuja. Asi mismo, la linea que parte de! punto critico de la envolvente y se desplaza a la derecha de éste se denomina linea de puntos de rocié. Las lineas internas de la envolvente denotan la region de dos fases y su valor porcentual representa la cantidad de liquid presente en la mezcla, a una determinada presion y temperatura Panto Bubuja Figura 2.8, Diagrama de Fases Tipico. El punto C de la Fig. 2.8 representa el sistema totalmente liquide a una presién relativamente alta. A medida que disminuya isotérmicamente la presion, se alcanzara el punto C, denominado presién de burbujeo o saturacin, en donde se liberara la primera burbuja de gas proveniente del gas disuelto en el petréleo. A menor presién, mayor presencia de gas libre en el sistema. El ico de yacimientos comportamiento anteriormente descrito resulta ti convencionales de petréleo. De igual manera y sobre la base de la Fig. 2.6, si TaN rs ISH Dep LaCOMPORTAMIENTO DE AFLUENCIA 30 ‘CAPITULO las condiciones iniciales del yacimiento se encuentran entre la temperatura critica y la cricondentérmica, y la presion sea tal que exista una sola fase en el yacimiento, se puede hablar entonces de yacimientos de condensado, destilado © de condensacion retrograda (un fenémeno en el cual se produce condensacién ‘0 vaporizacién en sentido inverso a lo que normalmente ocurre). Yacimientos de gas que presenten temperaturas a la derecha del punto cricondentérmico. Puede ocurrir que sean yacimientos de gas seco o himedo, siempre y cuando las condiciones de separacion en superficie generen la condicién 4, 0 4, sobre la envolvente, respectivamente. En el caso de yacimientos de gas seco, los hidrocarburos liquidos no se condesan en el yacimiento, asi como tampoco en la superficie. En consecuencia, el gas siempre permaneceré en su estado gaseoso. En el caso de yacimientos de gas hiimedo, las condiciones de separaci6n podrian generar un destilado o condensado. El concepto de presi6n de burbuja rocio resulta de vital importancia en el andlisis del comportamiento de produccién de un yacimiento. Por ejemplo, en un yacimiento donde la presién se encuentre por encima de la presién de burbuja P,, nunca existird la presencia de gas libre y por ende no se observara cambios en la permeabilidad del petréleo. Sin embargo, si la presion disminuye por debajo de P, , se permitird la presencia de gas libre en el medio poroso cuyo valor de saturasién incrementaré a medida disminuya dicha presion. El incremento en la saturaci6n de gas libre alrededor del pozo incrementara 1a permeabilidad al gas y disminuira la permeabilidad al petréleo, En otras palabras, el gas se moverd preferencialmente a través del medio poroso, afectando directamente el indice de productividad J c. Pe neabilidad Relativa ky La permeabilidad relativa k, de una fase cualauiera k ha sido el concepto frecuentemente utilizado como referencia para describir el flujo de dos o mas fases de fluidos inmiscibles a través de! medio poroso. SEN OE SIE Dep Lagat WL FOGLIDUTSDOROTSCCOMPORTAMIENTO DE AFLUENCIA BY CAPITULO T La permeabilidad efectiva es una medida de la conductividad del medio poroso a un determinado fiuido cuando existe mas de un fiuido en dicho medio que son méviles 0 no. Generalmente, se les designa como k, y donde el subindice & denota las posibles fases existentes en el yacimiento (petrdleo, gas © agua). &, es funcién del medio poroso y de Ia saturacién y distribucion del fluido en referencia. La suma de las permeabilidades de cada una de las fases es menor a la permeabilidad absoluta y ello se debe a la interferencia mutua entre los fluidos ‘que fluyen simulténeamente en el medio poroso. Debido al amplio campo de valores de permeabilidades que pueden tener las rocas de un yacimiento, resulta conveniente normalizar las permeabilidades efectivas. Esto puede ser posible, si se relacionan con la permeabilidad absoluta. En consecuencia, se define como permeabilidad relativa a! cociente entre la permeabilidad efectiva de un fluido, a un valor de saturacién dado, y la permeabilidad absoluta de! medio poroso. Por lo tanto, la permeabilidad relativa al petréleo, gas y agua, es decir kegs kay ¥ Rye Fespectivamente, se encuentra dado por: kg = (26) (27) (2.8) Las permeabilidades relativas son estimadas a partir de pruebas en el laboratorio y sus valores generalmente se presentan en formas de curvas para sistemas bifasicos agua-petréleo y gas-petréleo, como una funcién de la TUNER TEAHTT ‘Dep LT LFUOTHOTSOONOTSCCOMPORTAMIENTO DE AFLUENCIA Mm ‘CAPITULO T saturacién de la fase mojante. La Fig. 2.9 muestra una curva tipica de permeabilidad relativa Soturaci Residual de Potolno y Agus Connats SSmuracon Residue Gas 10 s Lo ‘Permesbilad Relais a Pteleo 100 ‘Saturacin Laude % Figura 2.9. Curva Tipica de Permeabilidad Relativa. Sistome Gas-Petrdleo. Sin embargo, algunas veces sus valores pueden ser estimados mediante correlaciones disponibles en la literatura o mediante ecuaciones basicas, desarrolladas a partir de simples suposiciones. El Anexo A muestra, en detalle, algunas de estas ecuaciones 0 correlaciones. Como se explicd en la seccién anterior, debido a cambios de presién o temperatura, se puede generar la presencia de gas libre, cuyos valores de saturacién puede incrementarse 2 medida que la presién disminuya en el yacimiento, permitiendo de esta manera la reduccién del indice de productividad .J e incrementando o permitiendo que la fase gaseosa se mueva con mayor facilidad a través del medio poroso, SON TTR SATIS = Dap Lega LRG DOOTEONEOT yo‘COMPORTAMIENTO DE AFLUENCIA 3 TARTULOT d-. Viscosidad del Petroleo 14, Representa una medida de la fricci6n interna o resistencia que ejercen las moléculas del petréleo a fluir (moverse) y es un importante parémetro en el célculo del indice de productividad ./. La viscosidad se encuentra principalmente afectada por presién y temperatura, Salvo algunas excepciones, la mayoria de los yacimientos exhiben un bajo grado de buzamiento, por lo que resultaria valido asumir que la viscosidad se encuentra principalmente afectada por presiOn, mas que por temperatura. Bajo esta condicién particular, la viscosidad del petréleo saturado con gas disminuiré, a medida que la presién disminuya desde la presién inicial del yacimiento hasta la presion de burbuja P, . Por debajo de la P, la viscosidad incrementaré nuevamente a medida que mas gas en solucion deje las moléculas mas pesadas en la fase liquida. La Fig. 2.10 muestra el comportamiento tipico de la viscosidad a medida que cambia la presién, a temperatura constante. En consecuencia, cambios de viscosidad con presion afectaran el indice de productividad J T = Constante Figura 2.10. Comportamiento de la Viscosidad del Petrdleo con Presin TL ROSTaTSOROTSCCOMPORTAMIENTO DE AFLUENCIA a CAMTULOT Siempre que sea posible, la viscosidad debe ser determinada en el laboratorio, bajo condiciones de presién y temperatura. Sin embargo, existen srentes correlaciones disponibles en la literatura para determinar su valor, en caso de que no se disponga informacién de taboratorio. El Anexo B presenta algunas de éstas correlaciones necesarias para estimar el valor de la viscosidad, a cualquier condicién de presién y temperatura. e. Factor Volumétrico det Petréleo 2, El factor volumétrico del petroleo B, es oro de los factores que afectan el indice de productividad de un pozo. A medida que la presién disminuye, el liquido se expande lo que se traduce en un aumento del factor volumétrico. Esta condicién se mantiene hasta alcanzar la presién de burbuja P,, a partir del cual el B, volverd a disminuir. Estos cambios de B, con presién puede detallarse en Fig. 2.11 T = Constante Figura 2.11, Comportamiento del Factor Voluméitico del Petréleo con Presi EN THORNE TATA Dep Legal N° LFGT2OOTSOORDTS‘COMPORTAMIENTO DE AFLUENCIA 3% . CAPITULO EI B, puede determinarse mediante el andlisis de una muestra de fluido en el laboratorio, Sin embargo, algunas veces puede no contarse con este tipo de anélisis y en consecuencia se debe hacer uso de ciertas correlaciones disponibles en la literatura para estimar B,. El Anexo B presenta algunas de éstas correlaciones necesarias para estimar el valor del factor volumétrico, a cualquier condicion de presién y temperatura f.. Efectos Colaterales La ley de Darcy asume que la permeabilidad del medio poroso es constante en todo el medio poroso y que ademds el fluido se desplaza bajo flujo laminar. Ambas suposiciones resultan no ser totalmente ciertas a nivel de campo por varias razones. Se ha demostrado que la permeabilidad cambia en todas las direcciones del yacimiento y esto se debe principalmente a la manera en que los sedimentos fueron depositados inicialmente. Pero no solo depende de la roca, también depende de la direccién del flujo, por lo que es una propiedad direccional. En consecuencia, los medics porosos son de naturaleza anisotropica, en cuanto a permeabilidad absoluta se refiere. En las cercanias del pozo, esta permeabilidad absoluta puede ser incrementada mediante algtin trabajo de estimulacién o fracturamiento hidrdulico, o puede disminuir debido a trabajos de completacién del pozo (invasion del filtrado del lodo, dispersion de arcillas, presencia de revoque de lodo 0 cemento, presencia de una alta saturacién de gas alrededor del pozo, penetracion parcial, limitada o una baja densidad de tiro, son algunos de los factores responsables en la reduccién de la permeabilidad). En ambos casos, se genera un cambio en la pendiente de la curva de presién, entre la cara de la arena y el radio para el cual la permeabilidad fue afectada r,. La Fig. 2.12 muestra los efectos que genera el cambio de permeabilidad sobre el perfil de presion en el yacimiento. Posiblemente, resulte dificil determinar 7, 0 en todo caso estimar el valor de la permeabilidad de la zona de dafio k,, denominada asi por van Everdingen y Hurst en 1983. Sin embargo, se ha sumido que el cambio de presion debido a la TENS ORTTEST A Dep TaN EFOETROOTROODOTSCCOMPORTAMIENTO DE AFLUENCIA 2% TAPTULO I alteracién de la permeabilidad puede ser perfectamente representada por el factor de dafio S q= Constente I Figura 2.12. Pfetos que Ocasiona la Alteracin de la Permeablidad sobre el Perfil de Presion en el Yacimiento El factor de dafio S, se encuentra formado a su vez por dos tipos de dafo: for S, yturbulencia S, S=5,+5, (2.9) El dafio de formacion S, define la resistencia al flujo de fiuidos sobre el jam? estrato, justamente en frente de la cara de la arena en el fondo del pozo, y Su y valor puede ser determinado mediante una prueba de restauracién de presién en | mediante una prueba de restauracion de ProSiOn Sr __ el pozo (b lup test). Si el valor de S, resulta positive, se dice entonces que existe dario en la formacién, justo alrededor del fondo del pozo. Este dafio puede reducirse mediante tratamiento con Acido. El tipo de acido a utilizar dependera del tipo de roca y del material de taponamiento que debe ser removido. Si la EN VeSR PITT ‘Dep Legal LFOGLROOTSTORNTS~cAPTULOT formacién es una caliza, el tratamiento con dcido clorhidrico removera el dario debido a la solubilidad de la roca misma. En el caso de areniscas, donde la matriz de la roca no es soluble, se hace uso de ciertos quimicos denominados lodos dcidos. Como resultado de un trabajo exitoso de acidificacion, el factor de dafo podria reducirse a cero o menor a este. Si S, es igual a cero, entonces no existe dario en la formacién. Finalmente, el valor de S, sera negativo cuando la permeabilidad de la zona de dafio sea mayor a la permeabilidad del yacimiento, como consecuencia de una acidificacion o fracturamiento. Pozos exitosamente 5. Por otra fracturados hidraulicamente exhiben valores de S, entre -3 es de efecto considerable en parte, el dafio generado por turbulencia de flujo pozos de gas o petréleo que presenten allas tasas de flujo. Se encuentra relacionado especificamente en la manera como fue completado el pozo y su valor puede ser estimado mediante técnicas 0 metodologias que serén explicadas a detalle en secciones posteriores. Estas técnicas requieren informacion adicional como algunas caracteristicas especificas del yacimiento y dos o mas pruebas de flujo estabilizado 0 isocronales. Como pudo notarse, el indice de productividad J se encuentra afectado por el factor de dao S. Sin ‘embargo, la metodologia propuesta por Darcy no ‘considera este efecto, y en consecuencia este fendmeno podria causar ciertas dificultades al momento de construir la JPR de un pozo. La siguiente seccion presenta una solucién a la limitacién por este concepto encontrada en la ecuacién de Darcy. B.. _Ecuacion de Difusividad Representa la eouaci6n diferencial basica de gran uso en la ingenieria de yacimiento para simular el fiujo de fluides a través de un medio poroso. Mediante un balance de masa y con la combinacién de la ecuacién de Darcy se puede obtener la ecuacién de difusividad, cuya solucién analitica depende en gran medida de las condiciones iniciales y de borde que sean consideradas. BONER ATT Dep. Lee LAGOSCOMPORTAMIENTO DE: AFLUENCIA a CAPITULO Sobre la base de la Fig, 2.13, la ecuaci6n de difusividad puede obtenerse asumiendo el flujo radial uniforme de un fluido a través del espesor total de una roca, considerada homogénea (en cuanto a propiedades de la roca) e isotropica con respecto a la permeabilidad. Se asume ademas que la roca se encuentra completamente saturada con un fluido monofésico. aa, Peau > SSO SEAN) $ Figura 2.13. Diagrama Tipico del Flujo de Fluidos en las Cercanias del fondo del Pozo. Considerando el flujo a través de un elemento diferencial de volumen, de espesor dr, situado a una distancia r del fondo del pozo, el principio de conservacién de la masa permite obtener la siguiente ecuacion diferencial: ap a 2.10) or a. % Volover POlaSe en donde qg y ¢ representan la tasa de flujo y la densidad del fluido que circula a través del medio poroso, respectivamente, h y ¢ definen el espesor y la porosidad de la roca en estudio, respectivamente. Combinando la Ec. 2.10 con la ecuacién de Darcy para flujo radial horizontal (Ec. 2.2) y haciendo uso de la definicion termo-dinamica de ta compresibilidad isotérmica, se puede obtener la siguiente ecuacién diferencia! TEN TET ITT ; ‘De Lagal * LFOGTROOTERRTSCOMPORTAMIENTO DE AFLUENCIA » ‘CAPITULOT uy representan la viscosidad y compresibilidad del fluido, respectivamente. La Ec. 2.11 representa la ecuacién diferencial basica para simular el flujo radial de un fluido monofasico a través de un medio poroso. Como bien puede notarse, la Ec. 2.11 es una ecuacién no lineal debido a la dependencia implicita de presién sobre algunos de sus parametros, tales como: k. p. #y €. No siempre sera posible encontrar una simple solucién analttica de la Ec. 2.11, en especial, debido a su no linealidad. Sin embargo y en funcién de las condiciones iniciales y de borde que se consideren, existe disponible en la literatura un variado numero de soluciones de multiple aplicacién en diversas reas de a ingenieria de yacimiento. Una de las mas comunes soluciones es denominada solucién de tasa terminal constante, en donde la condicién inicial asume que por un periodo de tiempo dado la presién del yacimiento se encuentra en equilibrio y por lo tanto el pozo se encuentra produciendo una tasa constante q, justo en el fondo del pozo (r=r,). Existen tres posibles condiciones a la solucién de la tasa terminal constante y estas son: Estado transitorio; Semi-estado estable; Estado estable. a. Condicién de Estado Transitovio pana Esta condicién resulta valida solamente por un periodo de tiempo relativamente corto, después de haberse creado alguna perturbacién de presion en el yacimiento. En términos del modelo de flujo radial, esta perturbacion seria causada por la alteracién de la tasa de produccién del pozo frente a la cara de la arena (r=r,). Durante este periodo de tiempo, se asume que la presién de respuesta del yacimiento no se encuentra afectada por la presencia de algun limite exterior o barrera. En consecuencia, el comportamiento de presion durante este periodo de tiempo resulta ser muy similar al comportamiento de presion esperado en un yacimiento infinitamente grande. TENET Dep Legal LFDSI2OOTSOOROTSCCOMPORTAMIENTO DE AFLUENCIA w Este tipo de solucién es comunmente utilizado en el andlisis de pruebas de pozos (well fests) en donde la tasa de produccién de un pozo es deliberadamente cambiada, y la respuesta de presién resultante en las cercanias del pozo es medida y analizada durante un breve periodo de tiempo, es decir, unas pocas horas después de haber ocurrido el cambio de tasa. La solucién de la ecuacién de difusividad bajo la condicién de estado transitorio ha permitido encontrar una solucién matematica para estimar el factor de dafio de la formacién. La Fig. 2.14 muestra el comportamiento tipico de la presién bajo condiciones de flujo radial en estado transitorio. La presién, asi como su derivada con tiempo, son funciones de tiempo y espacio. Ploy radial ech, \ anerdocid arto coninve p , 9=const SSCS SAAS +— afluencia de fluidos > Figura 2.14, Flujo Radial bajo Condiciones de Estado Transitorio, A con|inve Sen? eslecioracs® b-. Condicién de Semi-Estado Estable * seri! conitave Se ? eslaccora Esta condicion es aplicable a yacimientos que han estado aportando fluidos por un largo periodo de tiempo, de manera que la presion de respuesta del yacimiento se encuentra afectada por la presencia de algiin limite exterior 0 ~ Dap: Lega TRC z0OFSOOROTSOMPORTAMIENTO DE APLUENCIA CAPITULO barrera. En términos del modelo de flujo radial, la condicién de semi-estado estable puede ser descrita mediante la Fig. 2.15. Pp, 9 =const tam moratinat ae 8 y g 2 consi 4 p, a Goes yg : Zao ie — afluencia de fluidos Gg g g Z r a % Figura 2.15. Flujo Radial bajo Condiciones de Semi-Estado Estable. Se asume que el pozo se encuentra rodeado en su limite exterior por una pared sélida, impermeable, la cual impide el flujo de fluidos desde y hacia el yacimiento en cuestién. Por lo tanto, mas alld de este limite exterior cualquier cambio de presién con tiempo sera considerado igual a cero ((éP/a/)=0). En consecuencia, si el pozo se encuentra produciendo a una tasa constante de flujo, la presién en el yacimiento disminuira de una manera tal que (4P/ér) sera aproximadamente igual a una constante, para cualquier condicion de tiempo y espacio. Esta condicién ha permitido desarrollar ecuaciones de flujo para pozos de petroleo y gas. Sobre la base de la Fig. 2.15, la ecuacion de difusividad para flujo radial puede ser resuelta para flujo bajo condiciones de semi-estado estable. Esta solucién coméinmente utilizada, se encuentra dada por: egal LROTDOOTSONROTS TENE ONISHICCOMPORTAMIENTO DE AFLUENCIA CATO Pasa flo Semi efodo estable = 7.0810" kh las, = “te aaa iy aaa |e) oan La Ec. 2.12 describe el flujo de fluidos en el medio poroso, bajo condiciones de semi-estado estable. A diferencia de la Ec. 2.4, la solucién definida por la Ec. 2.12 considera el efecto por dafio a la formacién S q=const p Figura 2.16, Geometria y Distribucién de Presién Considerada para la Solucién de la Eeuacién de Difusividad para Flujo bajo Condiciones de Semi-Estado Estable, Note en Ec. 2.12 que el término encerrado entre corchetes representa el indice de productividad J, el cual es constante, siempre y cuando no cambien las propiedades de los fluidos y roca (1, B y S). Para una presi6n inferior a la presién de burbuja (P,, < P,), las propiedades fisicas de los fluidos cambiarian y por ende el indice de productividad 7, exhibiendo la curva de afluencia un comportamiento muy similar al mostrado en Fig. 2.17, Dap Lage LUG TDOTSOORTTS «de pede de BIG JE hein of CCOMPORTAMIENTO DE AFLUENCIA o CAPITULO Elconcepto de J constante no aplica, asus % 108107 Ap ka yp Ree 1 Ung i.) 31449) HB Figura 2.17, Comportamiento de la Curva de Oferta para Flujo bajo Condiciones de Semi-Estado Estable Matthews, Brons y Hazebroek (1954) sefialaron e! hecho que una vez que un yacimiento se encuentra produciendo bajo condiciones de semi-estado estable, cada pozo drenaré desde su propio area de drenaje e independientemente del resto de los pozos. La solucién de la ecuacién de difusividad, dada por la Ec. 2.12, asume que la forma de estos yacimien‘os puede ser del tipo radial o circular. Sin embargo, esta suposicién puede estar bastante lejana de la realidad. Para considerar esta pérdida de simetria, la Ec. 2.12 ha sido modificada por Odeh (1978), a fin de considerar las diferentes formas y posiciones del pozo en un area de drenaje dada. Para cualquier forma del drea de drenaje, la Ec. 2.12 ha sido reescrita como: C, Seri Ciode CEG ec - 7.08107 ky h “| my, By (In X =3/4 +S) Pag) CO g C 2 la Ew 7,\t en donde, la variable X puede ser definida de Fig. 2.18. SEND DROITS zai Dap Lega LFOGTROOTSOOIO7SCCOMPORTAMIENTO DE AFLUENCIA “CAPITULOT Tae aa 4 Jot AS a Mathews y Russel (1967) Figura 2.18. Factores para Considerar las Diferentes Formas Geométricas y Posiciones de un Pozo en una Determinada Area de Drenaje. . Condicién de Estado Estable ; fluo cortiava 5 @lac ona Esta condicién de flujo asume que el yacimiento es infinitamente grande y én todo caso, si se considerase algiin limite externo éste se encuentra abierto al fiujo, permitiendo que a una tasa constante cualquier volumen de fiuido producido sea exactamente balanceado por otro fluido que entra a través del TNT RITEISTT ep. Legal W LFOGT20OVSOMBOASCOMPORTAMIENTO DE AFLUENCIA 4s ace ‘CAPITULO limite extemo del yacimiento. Bajo esta condicién, la presion puede ser considerada independiente del tiempo. Estrictamente hablando, la condicion de estado estable ocurre solo sino existe una disminucién de masa en algiin punto del sistema de flujo. En la practica, esta condicién puede ser observada en yacimientos que se encuentran bajo la accién de un fuerte empuje hidraulico 0 que se encuentre sometide a un proyecto de inyeccion. La Fig. 2.19 muestra el comportamiento de la presién bajo esta condicién de flujo q=const pane | ee +— afluencia de fluidos Figura 2.19, Fiajo Radial bajo Condiciones de Estado Estable. Bajo flujo radial en estado estable, la solucién de la ec difusividad permite definir la siguiente expresion matemética para estimar la tasa de flujo a través, del oe poroso, en términos de la presion promedio del yacimiento fl 4o eae gobs nar hedicedt co 7.08107 k, ht Je-e Hy B, (ine, /7,)-1/2 +S). mh a é TENSE ISTT De Legal LPSSTROOTSOGEOTS‘COMFORTAMIENTO DE AFLUENCIA 4 CAPITULO 2.4.2 Métodos Empiricos para Construir la Curva de Afluencia JPR La determinacién de la JPR es posible si el valor de todas y cada una de las variables envueltas en las ecuaciones matematicas descritas en la seccién anterior es conocido. Sin embargo, rara vez suficiente y exacta informacion de nes de flujo multifasico, la solucién yacimiento existe. De hecho, bajo condi puede no ser sencilla y la forma mas exacta seria resolver las ecuaciones que gobiernan el flujo a través de! medio poroso mediante el uso de alguin simutader. En campo, algunos métodos empiricos pueden ser usados para estimar la PR de manera sencilla. La mayoria de estos métodos tan solo requieren de al menos una prueba estabilizada del pozo, como la tasa de produccién y su respectiva presién fluyente. ‘A-. Método de Vogel (1968) Mediante el uso de un computador, Vogel calculé curvas de IPR para pozos productores de petréleo pertenecientes a yacimientos por empuje de gas. en solucién. Estas curvas cubrieron un amplio rango de propiedades PVT’, permeabilidades relativas, espaciamiento de pozos y factores de dafio. Entre las suposiciones del método, se destacan: Eficiencia de flujo EF igual al 100% ; flujo radial uniforme; medio poroso uniforme e isotrépico; los efectos debido a segregacién gravitacional y compresibilidad de roca y agua son despreciados; P. P,), resultaria claro distinguir de la curva JPR una seccién recta y otra curva, tal como puede apreciarse en Fig, 2.22. Para la parte curva, se podria definir la siguiente ecuacién matematica, a partir de la ecuacién de Vogel y en funci¢n de las condiciones mostradas: ~ g=4% oun ~ M0 (2.18) Al derivar la Ec. 2.18 con respecto a la presién de fondo fluyente (d/dP,,), se obtendria la siguiente ecuacién valida para estimar el indice de y productividad 7 uJ Fics a Paro TPR os ~ 9b 4+(Qrox -F| Pore. tor =o bogy: J.(Rus- Po) ‘ BGnvx- gh = Se. SRE Huyeada Jo Ccuacton © r (eb J = Da TTS dS Q) edremos sovonde (4 x debaye di © CG & } >“ | eb bqCCOMPORTAMICNTO DE AFLUENCTA 50 CAPITULO 02 Saf) Sore ees (Frans a * So (2.19) Evaluando el indice de productividad J en el punto en comtn entre ambas curvas (P,, = P,), la E¢. 2.19 finalmente podria ser escrita como: a (2.20) Figura 2.22. Comportamiento Tipico de Ia Curva /PR, para un Yacimiento Sub-Saturado y Saturado. En consecuencia, el sistema de ecuaciones a resolver para fa construccién de la curva de comportamiento de afiuencia en un yacimiento sub-_ >P,) estarla definido por las Eos. 2.5, 2.18 y 2.20. La Fig. 2.23 (FB) turado muestra la forma tipica de una curva IPR de un yacimiento sub-saturado. SAN wR.CCOMPORTAMIENTO DE AFLUENCIA s CAPTULOT Figura 2.23. Combinacin dl indie de Productvidd y el Método de Vogel Una de las ventajas del método de Vogel es que la curva de afluencia IPR puede ser obtenida a partir de una simple prueba de produccién del pozo. Aunque el método fue desarrollado para yacimientos por empuje de gas en_ _solucién, la ecuacion ha mostrado también aceptables resultados en yacimientos con otros tipos de mecanismos de empuje. Por otra parte, el método puede también ser aplicado a pozos que producen simultaneamente agua, petréleo y gas, ya que el incremento en la saturacién del gas reduciria la permeabilidad relativa al agua. En este caso, ia tasa de produccién debe ser reemplazada por una tasa liquida (g, = 4, + 9,). La efectividad de la ecuacion propuesta por Vogel ha sido comprobada en pozos productores con un corte de agua de hasta 97%. por Vogel no considera los efectos del factor, El método original propues de dafio de la formacién. Sin embargo, una modificacién propuesta por Standing (1970) permitid extender la aplicacion del método a pozos dafiados o estimulados. Segtin Standing, el grado de alteracion de la permeabilidad puede TSR S90 SOIT(COMPORTAMIENTO DE AFLUENCTA _ ss CAPITULOT ser expresado en términos de una relacién de permeabilidad 0 eficiencia de flujo. denominado EF eee ah (00% (2.21) La condicion Pi’ se refiere al hecho de que un pozo no presente danio en la formacion (EF 00% ). El efecto que la eficiencia de flujo EF tiene sobre la produccién de un pozo puede ser analizado en Fig. 2.24. oe PO) 79 eshimulad es estimulada Leeadura kyhraalica Figura 2.24. Efectos del Datioy la Estimulacién Sobre la Presién de Fondo Fluyente P, Una eficiencia de flujo menor a 1.0 significa que el pozo presenta dafo y en consecuencia se tendra que efectuar algun trabajo en el fondo del pozo, como por ejemplo una acidificacién, para reestablecer el potencial del mismo. TON FRI TEISHT ‘Dp Taal N* LFOGTIOOTSOODOTSCOMPORTAMIENTO DE AFLUENCIA. 53 CAPITULO Una eficiencia mayor a 1.0 es el resultado de algtin cambio de la permeabilidad en el fondo del pozo y esto es posible mediante un fracturamiento hidraulico. Elefecto que la EF tiene sobre la prediccién de la tasa de flujo mediante ieacion el método de Vogel es posible, siempre y cuando se considere la mox de Standing. La £F puede ser estimada a partir del factor de dafio de la formacién S, (determinada mediante una prueba de restauracion de presién en el pozo (buildup test)), como ieee c dare =o para S In(y, /r,)—3/4 7 | in (r, /r,) 3/448] |7+S| 22) La Ec. 2.17 puede ser entonces utiizada para construir diferentes curvas de afluencia, valida para EF mayores 0 menores a 1.0, @ partir del uso de la Ec. 2.21, la cual puede también ser re-escrita como: Piel = P= (P, = Py) BF (2.23) Standing presenté una serie de curvas adimensionales de IPR para diferentes valores de EF entre 0.5-1.5, a partir del uso del método de Vogel y la Ec, 2.23. Estas curvas se presentan en Fig. 2.25. Pero el uso de la Ec. 2.23 presenta la particularidad que para bajos valores de P,, y altos valores de EF la modificacién propuesta por Standing genera valores negativos, tal como puede apreciarse en Fig, 2.26. Con el fin de corregir esta limitacidn, se ha propuesto el uso de dos posibles soluciones: La correlacién propuesta por Harrison y el uso de una ecuacién generalizada DONTECOMPORTAMIENTO DE AFLUENCIA st ‘CAPITULO oe EF =1.0 -» SinDaio EF>10 > Pozo Estimulado | oe EF<10 > Paxo Dodo | 00 04 02 03 04 08 08 OF 08 0.9 10 (al toners Figura 2.25. Comrelacién de Standing para diferentes EF . BF >0 y bajar valores de Pi 1P, Fs = Frroren ta Prediecién (a Goucdrs de Standing para diferentes EF . TSN OTESRO EIST ‘Dap Legal N*LFOGI2D07S0000754 Jauching- Waxswwso’ doch d iudnsurnieres pe Brcuiess ses ace Sees esenaasace ‘CAPITULO TI = ede ser utilizada en lugar de la La correlacién propuesta_por Harrison pt (4) 12-02 Ep), (2.24) ae) cast ‘Se ha comprobado que la correlacién propuesta por Harrison subestima las tasas de flujo al compararse con la ecuacién de Vogel. La Ec. 2.24 también puede ser re-escrita como una funcién de la EF . (2.25) 7 (2.25) (£] =12-02 Expl tra posible solucion a la limitacién de la correlacién de Standing puede ser obtenida mediante el uso de la siguiente ecuaci6n generalizada qe J, (PP Pay (2.26) Los valores positives obtenidos a partir de la ecuacién de Vogel representaran una linea recta si éstos se grafican en papel log-log, tal como puede observarse en Fig. 2.27. El coeficiente ./, y ef exponente n pueden féciImente ser obtenido del mencionado grafico, como el intercepto de la recta con la ordenada y el inverso de la pendiente, respectivamente. Una vez definidos los valores de J, y n, la Ec. 2.26 puede ser utiizada para completar los valores oorrespondientes de la IPR, tal como puede apreciarse en Fig. 2.28 aR TRISTE Dep Lag LAT ROTSEOTS(COMPORTAMIENTO DE AFLUENCIA in aati eee Rhy -pty Grifico de Bevacién Generalizada para Corregir la Limitacion de la Cortelacién de Standing. Figura 2.28. Extension de la Curva de Afluencia mediante Ecuacién Generalizada THON REIT ‘Dep Legal LPI DOOTSONETTCCOMPORFAMIENTO DE AFLUENCIA st CAPITULO Método de Fetkovich (1973) Fetkovich propuso un método para calcular la JPR en pozos de petréleo, ‘ando el mismo tipo de ecuacién que ha sido utilizada por afios para analizar pozos productores de gas. Esta técnica ha sido corroborada a través de pruebas isocronales (isochronal) y estabilizadas (flow-after-flow). Las pruebas de flujo isooronales son pruebas cortas de flujo, de igual tiempo de duracién, con cada prueba de flujo separada mediante un perfodo de cierre que permite la restauracién momentanea de la presién en las cercanias del pozo. Las pruebas estabilizadas presion fluyente estabilzada, la cual debera ser medida, Una vez logrado esto, ja tasa de produccién es nuevamente cambiada y el procedimiento se repite nsisten en colocar un pozo en produccién hasta aleanzar_una eee para diferentes tasas. Es igualmente valida para yacimientos saturados y sub- saturados, con permeabilidades entre 6-100 mD. El método de Fetkovich utiliza la misma ecuacién generalizada utilizada para pozos de gas y la cual se encuentra dada por: (2.27) donde C representa el coeficiente de flujo. El exponente » es el inverso de la pendiente y depende de las caracteristicas de! pozo. Los valores de C y » deben ser estimados a partir de pruebas de campo para de esta manera poder generar la curva JPR De acuerdo a la Ec. 2.27, se requiere al menos de dos pruebas para estimar C yn, asumiendo que la presién del yacimiento /, sea conocida, Sin embargo, ha sido costumbre el usar al menos cuatro pruebas de flujo para determinar ambas- constantes y reducir de esta manera la posibilidad de cometer errores en los calculos: consideracién igualmente recomendada para pozos de petréleo. La data obtenida de las pruebas isoc isocronales 0 estabilizadas representaran una linea recta si se grafican en papel log-log, tal como puede apreciarse en Fig. 2.31 Seal WEF DSTOTSOGROTSCOMPORTAMIENTO DE AFLUENCIA. ss ‘CAPTULOT Figura 2.31, Prueba de Cuatro Puntos. Ecuacién de Fetkovich, En la practica, el tipo de prueba de cuatro puntos a elegir para aplicar el método de Fetkovich dependera en gran medida del tiempo de estabilizacién det ozo, el cual es funcion directa de la permeabilidad del yacimiento. Si un pozo establliza_rapidamente, se recomienda considerar una prueba estabiizada_ convencional de cuatro puntos (flow-atter-flow). En aquellos pozos que se requiera largos periodos de tiempos para su estabilizacién, se recomienda “Considerar as pruebas isocron tipo convencional o modificada. €- Método de Jones, Blount y Glaze (1876) En 1976, Jones, Blount y Glaze publicaron un método que considera los efectos de turbulencia o flujo no-Darcy que ocurren en las cercanias del fondo del pozo y evaluaron su incidencia sobre la IPR. Estos investigadores suyirieron que el fiujo radial para gas y petréleo podia ser representado de otra forma, de manera que si alguna restriccién existiese esta pudiera ser considerada. El procedimiento propuesto utiliza pruebas de produccién (en sustitucién de las pruebas de restauracién de presién) para determinar si la turbulencia es el TaN STOTT ‘De -Legat *LPOTDOOVSUOROTS(COMPORTAMIENTO DE AFLUENCTA » CAPITULO principal factor que reduce la capacidad de flujo del pozo. Si la caida de presion por turbulencia es mayor a lo esperado, el pozo probablemente no habra sido completado eficientemente. En pozos carioneados, esto podria significar pocas perforaciones abiertas al flujo. Para un pozo fracturado, podria indicar que las fracturas son muy angostas. Aunque el método tiene mayor aplicacién en pozos productores de gas, resultados han mostrado igual efectividad en pozos con altas tasas de flujo de petréleo. Este método requiere de informacién adicional y de una prueba de flujo de al menos cuatro puntos. Sobre la base de la ecuacién de Forchheimer (1901), la siguiente ecuacién resulta valida para relacionar la caida de presién como una funcién de pe grade lurbulercia la tasa de flujo de petréleo. - : ion dedi ( “yecondee | ag+b. (2.28) =F) El coefciente « indica el grado de turbulencia en las cercantas del fondo del pozo. El coeficiente b indica las condiciones de dafio o no de la formaci6n. Si se dispone de al menos dos pruebas de fiujo estabilizadas, un grafico de una linea recta, a partir (2, ~P,,)/q Vs. q en coordenadas cartesianas defini del cual los coeficientes a y b podrian ser determinados. a representaria la pendiente de la recta y b el intercepto de ésta cuando q tienda a 0 Usualmente, se recomienda disponer de més pruebas de fiujo a fin de reducir posibles errores en las mediciones. La Fig, 2.34 muestra la representacion grafica de (P. - P,,)/q VS. 9 Un grafico de ia ecuacin de Jones, Blount y Glaze utilizando tres 0 cuatro pruebas de produccién permitiria no solo definir los coeficientes a y b, sino también distinguir las pérdidas de presién, ocasionadas por flujo no-Darcy de las pérdidas de presion ocasionadas por dafio (8). Este hecho seria un factor Den aN LTTElg = (COMPORTAMIENTO DE AFLUENCIA CAPITULO importante para decidir en realizar algtin trabajo de estimulaci6n o de otro tipo, para mejorar la productividad de un pozo, Otro parametro denominado relacién de b /& puede ser utilizado también como un buen indicador para determinar las pérdidas de presién causadas por flujo no-Darcy. EI valor de 6° puede ser calculado mediante la siguiente ecuacion: b+ 0 Gan (2.29) El valor de b’ también puede ser determinada gréficamente de Fig, 2.34, para la condicién de (P, -0)/g. Figura 2.34. Curva Tipica dela Ecuacién de Jones, Blount y Glaze La Fig. 2.36 ilustra algunas posibles conclusiones que pueden ser obtenidas de graficar la ecuacién de Jones, Blount y Glaze. De acuerdo a Fig. Dep. Lapa FORT EOOTSOORN SBN OTERO. IEST_COMPORTAMIENTO DE AFLUENCIA CAPITULO 2.35, las siguientes consideraciones resultan validas del método de Jones, Blount y Glaze. 4-. Sfel valor de b es bajo (< 0.05), se considera que la formacién esta libre de dano, Si la relacién b'/b es baja (<2.0), se puede asumir que no existe restricci6n de flujo por turbulencia. Silos valores de b y b'/b son bajos, el pozo esta bien completado. jel valor de b es bajo y el de 6'/b es allo, no se recomienda realizar una estimulacién al pozo. La baja productividad es causada por una insuficiente érea cafioneada. En este caso, se recomienda aumentar la densidad de tiro. Si el valor de b es alto y el de 5'/b es bajo, se recomienda realizar una estimulacién al pozo. Le wots + te wols ‘coro Turbulence Figura 2.35, Interpretacin de Varias Pruebas de Pozos, mediante el Método de Jones, Blount y Glaze TEN ITER ISTD Dep Leal LOSTCCOMPORTAMIENTO DE AFLUENCIA @ ‘CAPITULOT! La versatilidad del método de Jones, Blount y Glaze radica en el hecho que puede ser aplicado para el cdlculo de la caida de presién por efecto de: turbulencia de flujo; densidad de cafioneo; empaque por grava; entre otros. En esta secci6n, las ecuaciones discutidas seran validas para el flujo de petroleo, solamente. Las modificaciones respectivas para flujo de gas se discutiran en secciones posteriores en este libro. La caida de presién por "Efecto de Turbulencia” en pozos de petréleo generalmente se desprecia. Sin embargo, si se considera que su efecto es importante debido a las altas tasas, esta caida de presion puede ser estimada por el método de Jones, Blount y Glaze. Si se dispone de suficiente informacién de yacimiento, los coeficientes a y b pueden determinarse a partir de las siguientes Tne e. “que Pirie dubse coe g \ 310 f B p, aca (2.30) B, [in (r, /r,) - 3/4 +5] 231) 7.081 4, YB, representan la viscosidad y el factor volumétrico del petréleo, expresado en cps y Bbi/ Bn, respectivamente. S define el dafio de la formacién, solamente (S=S',). & representa la permeabilidad absoluta del yacimiento, en mD. El espesor de a zona productora es representado por hi, en pie. La densidad del petrleo a condiciones de yacimiento p, puede ser estimada $e, mediante: jesod eee . . 15.615 pOAn Corey, BUSES. (2.32) Po 7CCOMPORTAMIENTO DE AFLUENCIA “APTTULO | donde y, y 7, representan la densidad relativa de las fases petroleo y gas, respectivamente. R, define la solubilidad del gas en el petréleo, en pon/ Bn. El coeficiente de velocidad para flujo turbulento # puede obtenerse analitica 0 gréficamente. Para una formacién consolidada y no-consolidada, el coeficiente i puede estimarse analiticamente respectivamente como: 3.10" a (2.33) 4710" b= (2:34) Graficamente, los valores del coeficiente # pueden también obtenerse mediante Fig. 2.36. 5 5 3 rotor 9 ae war, xf ast ae a PemMeAILT, MIUDARCIES emt Brown (1854) Figura.36. Coeficiente de Velocidad para Flujo Turbulento, como una Funciin de Permeabilidad SORT IEa— i aaa LOST SOITCOMPORTAMIENTO DE AFLUENCIA os ‘CAPITULOT 2.2 Comportamiento de Afiuencia en Pozos Productores de Gas A. Ecuacién de Difusividad La ecuacién de difusividad para flujo radial (Ec. 2.11) puede ser linearizada para flujo real de gas. Esta linearizacion es posible gracias al uso de una funcién de seudo presién m (P), propuesta por Al- Hussainny et al. (1966) y la cual se encuentra definida como: A m(P)=2 @ ; @ ! agree: (2.38) donde P y P, representan la presién del sistema y la presién base, respectivamente. 1, y Z definen la viscosidad y el factor de compresibilidad del ‘gas. Por definicion, la densidad del gas_p, puede definirse como: = (2.36) donde, M, R y T definen el peso molecular, la constante universal de los gases y la temperatura, respectivamente. Combinando Ecs. 2.11, 2.35 y 2.38, se puede definir ia ecuacién de difusividad para una fase gaseosa como: out) BE He Om(P) (237) k a Como puede apreciarse, la Ec. 2.37 requiere de una condicién inicial y dos condiciones de borde. Bajo una condicién de semi-estado estable, la condicién inicial se encuentra dada por: BON RITE AST ep Legal N° LFOGIROO7 DODDSCOMPORTAMIENTO DE AFLUENCIA « ‘CAPTTULOT am(P)_ oo. 2Pa dee (FT a ar My hoe\Ty J (2.38) Por lo tanto, la Ec. 2.38 permite re-escribir la Ec. 2.37 como: Pa Ge. (1) ee ( TJ (2.39) En este caso, la condicién de borde estaria dada por: om(P) ar (2.40) para r Utiiizando una metodologia similar a la empleada para resolver la ecuacién de difusividad para una fase liquida, Ecs. 2.39 y 2.40 permiten definir la siguiente ecuacion de flujo de gas, valida bajo condiciones de semi-estado estable: (2.41) dq, fepresenta la tasa de gas a condiciones normals, en MPCND. ky hn representan la permeabilidad y el espesor, en md y pie, respectivamente. T define la temperatura del yacimiento, ’X. La funcién de seudo presion m(P) $2 encuentra dada en Lpe? /eps. La funcién de seudo presion m(P) seré definida en este caso como: TENGE STITT Dap Legal LF OGIROO7SDOOOTS‘COMPORTAMIENTO DE AFLUENCIA ‘CAMTULOT m(P)= j2e a, (2.42) 2 nba 2 donde, P y P, comesponderé a la presién del sistema y la presién base, respectivamente. Un grafico tipico de la relacion (4, Z) vs. P se muestra en Fig. 2.39 y 2.40. De acuerdo a estos graficos, puede concluirse que la relacion (ui, Z) sera constante para presiones comprendidas entre 0-100 Zc. En consecuencia, la relacién (1, Z) puede salir de la integral y la funcion de seudo presion m (P) podria ser definida como: m(P)-—- [Par (2.43) Zi ae (cas cRaviry = 065 L 1 L L L L L hee was S008 ea PRESSURE, psc omit Bren (1988) Figura 2.39. Comportamiento Tipico de la Relacién (44, Z) vs. Presién. SENTERO ISTE Dap Legal N LROSIROOTSUOROTSfit CCOMPORTAMIENTO DE AFLUENCIA a ‘CAPTULOT Z| TEMPERATURE © 200°F “cay PRESSURE, psi Koma rou 1886) Figura240,Comportanieno Tipco de a Relacin (ji, Z) vs. Presi, «Temperatura Constante, La solucion de la Ec. 2.43 permitira entonces definir la funcion de seudo presion m (P) como’ (2.44) donde jz, y Z representaran la viscosidad y el factor de compresibilidad del gas, respectivamente, estimadas a una presién promedio P, la cual sera determinada como: (2.45) NTs — Dap Lega LAST IOOTSOOIOTSCCOMPORTAMIENTO DE AFLUTENCIA CAMTULOT La Fig. 2.41 muestra una curva de comportamiento de afluencia tipica para pozos de gas, obtenida mediante Ec. 2.41 10" kh (om (2) = m (Py) Presion do Fondo Fluyente Figura 2.41. Curvas Tipica de Comportamiento de Afluencia para Pozos de Gas. Por otra parte, esta misma Ec 2.41 puede también ser resuelta en combinacién con la Ec. 2.44, lo que permitiria determinar la siguiente expresién matemética para estimar el flujo de gas que circula a través de un medio poroso, como una funcién de la presién de yacimiento y la presién de fondo fluyente, P, y P.;, respectivamente. pt 70310 kh (Fe (2.46) NEO a. Lagat CFT ROOT SOCAOTSCOMPORTAMIENTO DE AFLUENCIA o ~~ CAPTULOT De acuerdo a Figs. 2.39 y 2.40, la Ec. 2.46 es valida para presiones menores @ 1000 Zpe. Sin embargo, algunos autores sugieren utilizar la misma hasta presiones < 2500 /pc. Esta sugerencia aparenta ser valida al analizar la Fig. 2.42, en donde se presenta el comportamiento tipico de la relacién m (P) vs. Py (2 P*)Au, Z) vs. P, a temperatura constante zs00f PLOT OF MIP] va Yq 1 O70,T + 200°F & % ssooh joo at a “FP rose ico Sf e soo} 400 DSS TOGO 30004000 50506003 TOOT PRESSURE, PSIA| Kent Brown (1084) Figura 2.42. Comportamiento Tipico de la Relscion m (P) vs. P y (2 P*)/(t4, Z) ve. Py ‘Temperatura Constante, De acuerdo a la Fig. 2.42, la relacién (2 P?)/(u, Z) vs. P se desvia de la relacién m(P) vs. P, a presiones superiores a 2500 Lpe. Por otra parte y de acuerdo a Fig. 2.42, resultaria valido asumir que para presiones mayores a TaN TE PHOT ISI Dep. Leal LPOSTROOTSOOROTSCOMPORTAMIENTO DE AFLUENCIA 0 — enTULOT 2500 Lpe, la pendiente de la curva es una constante. Esta suposicion, permite establecer la siguiente consideracion: =c (2.47) Combinando Ecs. 2.42 y 2.47, se tiene que la funcién de seudo presion m(P) se encuentra dada por: m(P)=C | dP=C(P-F), (2.48) donde, 2.49) Hs, 2) oe Finalmente, al combinar Ecs. 2.41 y 2.49 se fiene que: 4. ween’ (2.50) La Ec. 2.50, 0 la Ec. 2.41 en conjunto con Ia Fig. 2.36 para estimar la tasa de gas a condiciones nomales 9,,, pata_presionés __ mayores a 2500 Lpe. cr latices B.. Prueba de Cuatro Puntos La técnica propuesta por Jones, Blount y Glaze puede también ser utlizada en pozos de gas. Esta técnica permitiria re-esoribir la Ec. 2.45 como: TSUN TORO TE ITT ‘Dep Legal NERDS TROOT=OOROTS‘COMPORTAMIENTO DE AFLUENCIA a ‘CARTULOT (251) La solucién de la Ec. 2.51 requeriré de al menos una prueba de produccién de cuatro puntos, de manera que un grafico de (7? ~P3.)/q, VS. dy» en coordenadas cartesianas, definiria una linea recta, tal como puede apreciarse en Fig. 2.43. El coeficiente @ representa la pendiente de la recta y 6 el intercepto de ésta cuando q tienda a 0 Figura 2.43. Curva de (P? — P2.)/q vs. q- Método de Jones, Blount y Glaze para Pozos de Gas. Similar que en pozos de petréleo, los coeficientes @ y b indican el grado de turbulencia y el grado de dafio de la formacién, respectivamente. Jones, Blount y Glaze sugirieron el uso de este grafico para determinar la presencia de alguna restriccién en las cercanias del pozo. Para ello, inicialmente deberia $i a aCOMPORTAMIENTO DE AFLUENCIA n CAPITULO obtenerse el valor maximo de (??—P2)/g, cuando P,, tienda a cero. En consecuencia, este valor deberé compararse con el valor obtenido de = P2)/q en el intercepto b. Si el valor maximo de P?/q, supera dos o tres veces el valor de 4, entonces es posible que ocurran algunas restricciones en las cercanias del pozo, En este caso, seria recomendable aumentar la densidad _ de cafioneo. De cualquier forma, en pozos de gas se utilizaran las mismas consideraciones formuladas, analizadas y discutidas para pozos de petrdleo. 2.3 Casos Particulares Todos los métodos presentados hasta ahora en este libro (necesarios para construir la curva de comportamiento de afluencia) han considerado: 1) Pozo verticalmente perforado; 2) El pozo produce desde una Unica zona; 3) Produccién bajo condiciones de semi-estado estable o fiujo estabilizado. Sin embargo, en la vida real un pozo produce bajo condiciones completamente diferentes a las ya estudiadas. Es practica comin completar un pozo con caracteristicas muy particulares, tales como: Pozo vertical y cuya produccién provenga desde miiltiples arenas; Pozos horizontales; Pozos productores de yacimientos sometidos a proyectos de inyeccién de agua o gas; Pozos de alto potencial lo que restringe algunas veces medir la presién de yacimiento; Entre otros. También, algunas veces resulta interesante y necesario predecir a futuro el comportamiento de la curva de afluencia en pozos de petréleo o gas. Por lo tanto, los fundamentos tedricos estudiados y considerados hasta ahora no permiten manejar las condiciones anteriormente mencionadas. En consecuencia, estos escenarios deben ser tratados como casos particulares. Como parte de un nivel basico, este libro solo trataré dos casos de especial importancia, y los cuales serén presentados en las secciones siguientes. El resto de las condiciones mencionadas deberén ser tratadas bajo un nivel de entrenamiento avanzado, debido a lo extenso y complejo del tema ISON STR SOTA Dag Legal LFUGTDOOTSOOSOTSCCOMPORTAMIENTO DE AFLUENCIA 2 ‘CAPITULOT 2.3.1 IPR Compuesta Algunos pozos son completados en dos o més arenas o lentes estratificados, por lo que el aporte de todas estas zonas es producida en Conjunto “commingled” en el fondo del pozo. Aun cuando estas arenas 0 estratos pertenezcan a un mismo yacimiento, las caracteristicas en cuanto al tipo de fiuido, mecanismo de empuje, permeabilidad y hasta presién pueden variar entre las mismas. Trayendo como consecuencia que ante cualquier cambio de la presion de fondo fluyente podria alterar el corte de agua y la relacién gas-liquida, %AyS y RGL, respectivamente. Bajo este escenario, el ingeniero de produccién estaré en la disyuntiva en como construir una curva de afluencia comin para las arenas que producen en conjunto. Con el objeto de establecer una metodologia para construir una curva de afluencia representativa de ambas arenas (JPR Compuesta), considérese la Fig. 2.44 en donde se ilustra la completacién de dos zonas 0 lentes estratificados, con caracteristicas de flujo ligeramente diferentes (2. JY da) ¥ Cuya produccién es aportada a un pozo en comin. Figura 2.44, Deseripeién Esquematica de dos Zonas Completadas en un, TEN TORTS Dep- Legal W LPOSIROOTSOODOTS(COMPORTAMIENTO DE AFLITENCIA 4 Lom Ante cualquier cambio de la P,,, se deberd estimar valores promedios de la fraccién total de agua y la relacién gas-liquida, f, y RGL. respectivamente, con el objeto de caloular las propiedades de la mezcla. Valores promedios de /, y RGL para una mezcla pueden estimarse haciendo uso de las siguientes ecuaciones matematicas: fa, (2.52) (2.53) En este caso, se asume que g, es la tasa de la mezcla liquida (4, =4 +9») por estrato o lente de arena. Experiencia de campo ha mostrado que el considerar la produccién de agua y petréleo en conjunto, como una mezcla, ha resultado valida en pozos con una produccién de agua de hasta 97%. Esto se debe a que el cambio de J es causado principalmente por la reduccién de la permeabilidad relativa al petréleo, a medida que la saturacién del gas aumenta en el yacimiento. EI gas reduce la permeabilidad al agua en la misma forma que reduce la permeabilidad al petréleo, Por otra parte, si se dispone de algunas pruebas de produccién individual por arena (haciendo uso de algtin registro de produccién) 0 si se conoce las caracteristicas de cada zona referente a la roca y los fluidos producidos (mediante algiin nicleo o registro eléctrico, y pruebas de laboratorio), ademas de la presion de yacimiento, se puede determinar la curva IPR individual por pozo, simplemente utilizando la ecuacién de semi-estado estable o cualquier otro metodo emprrico. SEN SANTINI ‘Da Leal LOT DOOTSORBTTSCCOMPORTAMIENTO DE AFLUENCTA CAPTULO La Fig. 2.45 presenta a manera de ejemplo, la curva de afluencia de cada arena presente en la Fig. 2.44, asumiendo el caso en que P, > 7,,. También se ha representado la curva de afluencia 0 IPR compuesta (representativa de ambas arenas consideradas). y E,+P, Out da) Tn) Yacimierto'8" aS > IPR Compuesta 4 Figura 245. IPR Individual por Pozo y 1a IPR Compuesta de dos Arenas Completadas en un mismo Pozo. De acuerdo a la Fig. 2.45, cualquier valor de P,, superior a P,, , permitira que los fluidos producides por el yacimiento “B" fluyan hacia dentro del yacimiento “A’. No existiré un aporte neto de produccién hasta tanto la presién en el fondo del pozo sea suficientemente baja pala que el aporte de fluidos de la zona de mayor presién (yacimiento A) sea superior a la tasa de fiuidos que fluyen hacia la zona de menor presién (yacimiento B). Este valor de presion a partir del cual comienza ha existir una produccién neta Pz, deberd ser estimada para poder construir la JPR compuesta. a2 a Dep Lea N° FSR(COMPOWTAMIENTO DE AFLUENCIA 6 CAPITULOTT El valor P;, puede ser determinado asumiendo q, =, (ya que el indice de productividad J es lineal en pequefios intervalos de disminucién de presién) En consecuencia, dn = Iu (Pry — Peg) = 44 =I 4 (Pep ~P,) (2.54) De manera que Ec. 2.54 permite establecer la siguiente relacion matematica: oF +P, Sal Ts) t= Tala) aca) Cualquier valor de P,, menor o inferior a P’, permitiré que ocurra un aporte neto de las arenas, en conjunto. Una vez definido el valor de P',, se procedera a construir la JPR compuesta, para lo cual se requiere asumir diferentes valores de P,, (<,) y determinar las tasas de flujo correspondientes a cada arena. La tasa total sera la suma de éstas (q = 4, +9,), correspondiente ala P, asumida. La Fig. 2.39 presenta la curva tipica de la JPR compuesta, valida para el ejemplo anteriormente considerado. 2.3.2 Prediccién de Presion de Yacimiento a partir de la JPR Compuesta Comunmente, la presién de yacimiento Pp, es un parametro frecuentemente medido en el campo. Aunque lo ideal seria medir la presién en aquellos pozos que hayan estado cerrados por un largo periodo de tiempo (para evitar pérdidas de produccién), en la préctica resulta exactamente lo contrario los pozos se encuentran actives en cuanto a produccién se refiere. Al cerrar la produccién de un pozo, los costos por medir la presién aumentan considerablemente ya que este tipo de pruebas algunas veces demandan Dep Lagat LOGIT SOIDOTSCOMPORTAMIENTO DE AFLUENCIA n ~ CAPITULO mantener el pozo cerrado por largos periodos de tiempo. Adicionalmente, se debe siempre tener presente el riesgo asociado que significa bajar alguna herramienta bajo condiciones operativas del pozo (debido a las altas probabilidades de generar alguna obstruccién al flujo de fluidos dentro de la tuberia de produccién). En consecuencia, la medicién y monitoreo de la presion .da por cada uno de los factores de yacimiento puede algunas veces verse lit anteriormente mencionados. Sobre la base de la informacién de campo disponible y haciendo uso de la técnica de anélisis nodal, podria determinarse la jones estaticas. Un curva JPR de un pozo asumiendo un amplio rango de pre posible valor de P, seria estimado una vez que las ecuaciones reproduzcan las condiciones reales de produccién del pozo en el campo. Este proceso es aleatorio y puede requerir largas jomadas de tiempo y trabajo, algunas veces sin éxito. Sin embargo, si se dispone de varias pruebas de fiujo estabilizadas, el método empitico propuesto por Vogel 0 Fetkovich podria ser utilizado para estimar P. En el caso de Vogel, se pueden obtener algunas relaciones matematicas para estimar P., Para ello se debe al menos disponer de dos pruebas de produccién del pozo. Considérese el caso de un yacimiento sub-saturado en donde se dispone de dos pruebas de produccién, las cuales se ubican por encima de la presion de burbuja del yacimiento. Este caso se encuentra ilustrado en la Fig. 2.48. Bajo esta condicién, ambas pruebas poseen igual indice de productividad Jy por lo tanto a siguiente condicién resultaria valida: —f____ (2.56) @-Pp) @-Py) por lo que, ?, podria ser estimada como: SDN TE 90-12 IS a Legal LFS ROOTRCOROTSCCOMPORTAMIENTO DE AFLUENCIA CAPTULOT lz (2.87) Este primer caso permite facilmente estimar a partir de la Ec. 2.57 la presién de yacimiento P., como una funcién de ambas pruebas de flujo, En aquellos casos en donde no se disponga de alguna prueba PVT, se puede hacer uso de alguna correlacién para estimar P,. Experiencia de campo ha mostrado resultados confiables cuando se le compara con los resultados de pruebas de presin medida en el pozo. 2 Prusbas do Fiyo Establizadas Presion de Fondo Fuyente "rasa de Filo Figura 2.48, Pruebas de Flujo Estabilizadas por encima de la Presin de Burbuja Un segundo caso seria el disponer de dos pruebas de fiujo estabilizadas: tuna por encima de la presion de burbuja 7, y la otra en la region de dos fases, tal como puede detallarse en Fig. 2.49. UNRATE ITT Dep. Legal WL ROOTROOTEOAITSCCOMPORTAMIENTO DE AFLUENCIA » CAPITULO 2 Pruebas de Fu Estabiteasoe Prosion de Fondo Fluyente esa de Fluo Figura 2.49, Pruebas de Flujo Estabilizadas: Encima de la Presign de Burbuja y en la RegiGn Bifasica Al combinar las Ecs. 2.5, 2.18 y 2.20, la siguiente ecuacién para estimar P, resultaria valida: itd Pa Sue (2.58) J donde, el indice de productividad J estaria definido como ie 18@, (2.59) Pa) 18 Py, 0.2 Py, ~ 0.8] P+ La determinacién del indice de productividad J depende de las dos pruebas estabilizadas de flujo y la presion de burbuja A, BONES ISTT ‘Bip Lagat NL DOORNCCOMPORTAMIENTO DE AFLUENCIA a0 Un tercer caso posible seria aquel en donde las dos pruebas estabilizadas disponibles se encuentren por debajo de la presion de burbuja P,, tal como se muestra en Fig. 2.50 2 Pras ce Flyo Estabizacat Presion de Fondo Fluyente Tes de Fisjo Figura 2.50, Pruebas d» Flujo Estabilizadas: Regién Bifisica Asumiendo que la presién de yacimiento se encuentra por encima de la presion de burbuja (?, > P,), la siguiente relacién matemética resultaria valida para estimar P,: nanrh (2.60) donde, el indice de productividad J’ puede ser estimado a partir de la siguiente ‘ecuacién: TaN oreo IST B ap-Lagal W CFDSTDOOTSUODOTSCCOMPORTAMIENTO DE AFLUENCIA a ‘CAPITULO 1.8 (G2 = 90) a{r-oa(% eal Como puede notarse en Ec. 2.61, se hace necesario determinar una (2.61) ‘ecuacién adicional para estimar la tasa de flujo a condiciones de presién de burbuja ?, . Esta ecuacién adicional se encuentra dada por: ifn or Un caso particular seria aquel en donde la presion de yacimiento se is ~ (2.62) S3[g? | oo|s> encuentre por debajo de la presién de burbuja (yacimiento saturado). En este caso, se puede hacer un andlisis similar pero considerando solamente las ecuaciones desarrolladas por Vogel. La presién de yacimiento P, podria ser estimada a partir de la solucién de la siguiente ecuacién cuadratica. eS)» ee vloGe)eCeie-* En el caso en que se considere el método propuesto por Fetkovich, se tendra, ademés de la presién de yacimiento P,, dos variables adicionales desconocidas (C y m). En consecuencia, se requerira de al menos tres pruebas de flujo estabilizadas para resolver el sistema de ecuaciones planteado SEN 7950-10258 ‘Dep Legal LFOGIZ007SO02075CCOMPORTAMIENTO DE AFLUENCIA 2 Tal (mediante el uso de algtin método iterativo 0 de ensayo y error). A partir de la Ec, 2.27 propuesta por Fetkovich, se puede demostrar que: a a 4s e e : 2.64 @-Pay Pay P= PY et La presién de yacimiento puede determinarse a partir de las pruebas de flujo | y 2, y Ly 3, mediante: gsor® (2.65) (2.66) En ambas soluciones (Ecs. 2.65 y 2.68) existira un nico valor para la constante 1 donde se cumple que (P.),»=(P,),s. Por lo tanto, debera establecerse un procedimiento iterativo o de ensayo y error que permita estimar el valor de Finalmente, puede darse el caso en donde solo se disponga de dos pruebas de flujo estabilizadas. Bajo esta condicion, puede asumir que n=1, lo que permitiria determinar de manera sencilla la presién de yacimiento P, y el coeficiente C presentes en el método propuesto por Fetkovich NRT Dep aga LFDTSCAPITULO III FLUJO MULTIFASICO EN TUBERIAS Se entiende por flujo multifasico, e! movimiento simultaneo de una fase libre de gas y liquido a través de una tuberia o conducto que le transporta. El gas y el liquido pueden existir como una mezcla multifasica 0 como dos fases perfectamente definidas. La distribuci6n fisica de estas fases en la tuberia es definida como patrén de flujo y las mismas se encuentran plenamente definidas para tuberias verticales, horizontales 0 inclinadas. Ya que el movimiento de _ s se encuentra directamente_relacionado_con fluidos a través de tube gradiente de presion, en la industria petrolera es de sumo interés determinar el mismo para el disefio de tuberias 0 facilidades de superficie que permita transportar los fluidos producidos por un pozo hasta los tanques de almacenamiento, de una manera eficiente y répida. Sin embargo, la determinacién de este gradiente de presién puede resultar compleja. En especial, por la naturaleza del flujo y la variacion de las propiedades de los fluidos producidos debido a cambios en su composicién por efecto de presion y temperatura. La determinacién del gradiente de presién es de suma utilidad para el disefio de instalaciones de levantamiento artificial, lineas de flujo y tuberias de produccién en pozos verticales y desviados, disefio de intercambiadores de calor, lineas de gas, entre otfos. Actualmente, se dispone de numerosas correlaciones y ecuaciones que han sido propuestas para el caloulo del gradiente de presién y las cuales seran presentadas en capitulo. SRN 978 980A ® ‘ae Lega NL 20073002073FLUJO MULTIFASICO EN TURERIAS a “APITULO UI 3.1 Ecuacién General de Gradiente de Presion La base tedrice para la mayoria de las ecuaciones de flujo de fiuidos esta basada sobre la ecuacién general de energia. Bajo condiciones de estado estable, Ia ley de conservacién no es mas que un balance de energia entre dos Puntos cualesquiera de un sistema, Considérese un balance de energia alrededor de un volumen de control, como el mostrado en Fig. 3.1. +2 Intercambiador 0 Fuente de Calor i Volumen de Control ' @ Figura 3.1. Volumen de Control en un Sistema de Flujo, Entre los puntos 1 y 2, la primera ley de la termodindmica sobre el sistema permitiria establecer lo siguiente: energia del trabajo ‘pérdidas] [energia del Sluido |+| adicional sobre|-| de |=) fluido que entra el flujo energia | | saliendo TaN TIT ‘Dap Legal PLFDSTDOOTSOODDISFLUO MULTIFASICO EN TUBERIAS, as CAPITULO La ecuacién de momento lineal y/o la ecuacién general de energia, que resultaria de aplicar la primera ley de termodinamica podria resultar dificil de resolver, en especial, por la limitacién en estimar el término de energia interna. Por tal motivo, la mayoria de los autores consideran aplicar simplemente la ley de Bemoulli entre dos puntos cualesquiera en la tuberia La ley de Bernoulli establece que la energia de presién, la energia potencial y de aceleracién es constante en un punto, y se mantendrd constante a lo largo de una misma linea de corriente. Sin embargo, el flujo de fluidos a través de un conducto experimentara pérdidas de energia que deberan ser consideradas. Aplicando este concepto entre los puntos | y 2, denotados en pro @.) donde 5 Perd|,, representa las pérdidas de energia entre los puntos 1 y 2, debido al sistema. El calor transferido @ es la energia en forma de calor que puede entrar o salir del sistema y W representa el trabajo. P, y, Y = representan la presién, el peso especifico del fluido y la posicién con respecto al sistema de referencia, respectivamente. V es la velocidad del fluido. Aplicando: la Ec. 3.1 a un sistema como el mostrado en Fig.3.2, se podria obtener la siguiente ecuacién de gradiente de presion dl 2g, dL fa haves ap _pdV? , pa Send, f pV 2) TN STS ITT Dep Legal NLA DGIEOOTSOODITS_ FLUO MULTIFASICO EN TUBERIAS 6 ‘CAPITULO TT Sobre la base de la Ec. 3,2, el gradiente de presién en tuberias dP / dL. serd funcién del gradiente por: aceleracion 247 a. elevacion 2& Seno fe friccion L2H 2¢.4 Bajo condiciones de flujo multifésico, algunos términos de la Ec. 3.2 necesariamente deberdn ser ajustados, a fin de considerar una mezcla de gas y liquido. Basicamente, los términos a ajustar serian: densidad p por densidad de mezcla p,, velocidad 7 por velocidad de mezcla ¥, y factor de friccién f por factor de friccién de mezcla /,,. Figura 3.2. Seccién de Tuberia Inclinada. El efecto de cada una de las fases consideradas en la mezcla sera representado mediante la fraccién de la tuberia ocupada por el liquido, la cual es conocida como hold up liquide. La determinacién ue este hold up dependera a su vez si se considera 0 no el deslizamiento entre fases que ocurre a lo largo de la tuberia, En todo caso, sera necesario disponer de éste valor o de alguna correlacién para estimar el mismo. aN oTeswTEIE2 ‘ap. Lagat LROGI200TSOOROTSFLUIO MULTIFASICO EN TUBERIAS a ‘CAPITULO I 3.2 Curvas de Gradiente Es la representacién grafica de los cambios de presion que un fiuido tiene a lo largo de la tuberia que lo transporta, Esta representaci6n grafica puede clasificarse en dos tipos: Estatico y Dinamico. 3.2.1. Curvas de Gradiente Estatico Para una tuberia vertical (Seno 90" =1) y bajo condiciones estaticas, la Ec. 3.2 puede escribirse de la siguiente manera’ ap. dL 433 7, (33) Los gradientes por aceleracién y friccién son despreciados por cuanto el fluido se encuentra sin movimiento, La Ec. 3.3 asume también que el gradiente de presién es constante para fluidos compresibles y ligeramente incompresibles, como el agua y él petréleo, respectivamente. Sin embargo, para un fluido como el gas esta aseveracion puede no ser totalmente cierta por cuanto la densidad del mismo cambia a medida que aumenta la columna de fluido. Algunos graficos disponibles en Ia literatura, como el mostrado en Fig. 3.3, permiten estimar el gradiente de presién de: gas como una funcién de la presién y de la gravedad ‘especifica. Aun cuando exista una leve variacién del gradiente de presién para el gas, la tepresentacién grafica de Presién vs. Altura 0 columna de fuido cualquiera representara siempre una linea recta. Para una tuberia horizontal sin embargo, el gradiente de presion, de acuerdo a la Ec. 3.2, seria igual cero. Esto no significa que éste sea despreciable, sino que se asume que la presién de! fluido es constante en todos los puntos de la tuberia. 0 3.2.2. Curvas de Gradiente Dinai Como su nombre lo indica, el fluido se encuentra en movimiento y por 1o tanto se debe tomar en cuenta, adicionalmente de los efectos gravitacionales, los efectos debido a friccién y aceleracién. Para una mezcla multifésica, cambios de tasas de flujo, geometria de la tuberia y grado de inclinacién, y propiedades SEN STOTT ISTT Dap Legal WLP DOOTSOOROTSFLUO MULTIFASICO E ‘CAPITULO UI de los fluidos debido a presion y temperatura, contribuyen a que el gradiente de presién continuamente cambie a lo largo de la tuberia. En consecuencia, la representacion gréfica de Presién vs. Profundidad o largo de tuberia dejara de representar una linea recta Kemi Brown (1868) Figura 3.3. Curva Tipica de Gradiente de Gas. 3.3. Construccién de Curvas de Gradiente Para construir fa curva de gradiente de fluido que fluye a través de una tuberia de longitud L y presién de entrada P,, como la mostrada en Fig. 3.2, se debe dividir la tuberia en N intervalos de longitud AL. Se selecciona la correlacién y/o modelo mecanistico que corresponda, a fin de estimar el gradiente de presién en fa mencionada seccién, siguiendo paso a paso el procedimiento general que se detalla a continuacién: EN STENT IST Dep. agai ROS ROOTSOOROTSFLUJO MULTIFASICO EN TUBERIAS 9 ‘CAPITULOM a, Se selecciona el primer intervalo b+, Se estima un valor inicial (AP), (caida de presion en el intervalo considerado) c. Se determina la presion y temperatura promedio para el intervalo seleccionado d-. Mediante a correlacién y/o modelo mecanistico pre-seleccionado, se estima el nuevo gradiente de presion (AP/ AL) e. Se calcula el nuevo valor de (AP), = AL (AP/AL) fe Se compara el (AP), con el calculado en el paso anterior. Sino se satisface una tolerancia prefijada, se debera tomar entonces el (AP), como el nuevo valor asumido (AP), y se repetiré el procedimiento antes descrito desde el paso b hasta f. El procedimiento termina cuando la tolerancia es plenamente satisfecha g-. Una vez determinado el valor de AP, correspondiente a la seccién de tuberia AL, se procedera entonces a definir el valor de presion en el extremo de la seccién, denominada P, he. Se seleccionaré un nuevo intervalo, asumiendo P, como Fy se repetira el procedimiento desde el paso b, hasta cubrir la longitud total de la tuberia. Posteriormente si se desea, se puede construir la curva de gradiente de presién bajo cualquier condicién (vertical, horizontal o inclinada), simplemente graficando P vs. L. Las Figs. 3.4 y 3.5 muestran las curvas tipicas de gradiente de presi6n para una tuberia vertical y horizontal, propuestas por Brown (1984). 3.3.1 Curvas de Gradiente en Tuberia Vertical El primer estudio de flujo bifasico vertical fue realizado en 1914 por Davis y Weidner. Sin embargo, gracias al trabajo presentado por Verluys sobre la teoria basica del fiujo vertical, es a partir de 1930 cuando se comienza a desarrollar correlaciones para la caida de presion en tuberta vertical. aN Sa IRATE : ep Lage NAST OSOONOTSFLUO MULTIFASICO EN TUBERIAS ry Ket oan (184) Figura 3.5. Curva de Gradiente de Presion en Tuberia Horizontal 0 TSN 9980. Dap Legal LFDSIDOOTSOOIOTSFLUO MULTIFASICO EN TUBER oo CAPITULO En 1931, Moore y Wilde intentaron expresar las pérdidas de presion en flujo bifasico para tuberia vertical, como una combinacién de las pérdidas por elevacién y por friccién. En 1952, Poettman y Carpenter desarrollaron una correlacién basada en la ecuacion general de energia, donde la pérdida total resulta de la suma de las pérdidas por elevacién y friccién. La pérdida por friccién se calcula usando un factor de friccién, el cual esta relacionado con el numero adimensional de Reynolds que desprecia los efectos de la viscosidad. Los fluidos se consideraron como una mezcla homogénea de petréleo, gas y agua, para el calculo de la densidad del fluido y de la velocidad del flujo. Poettman y Carpenter propusieron una ecuacién para caloular la presién frente a la cara de la arena en pozos de gas, en los cuales la variacién del factor de compresibilidad se tomé en consideracion. Esta corretacién permitié calcular presiones de fondo con una buena aproximacién, cuando la taza de fiujo es alta y la relacién gas-liquido es baja. En 1954, Gilbert propuso un trabajo ante el Instituto Americano del Petroleo. API, en el cual se presenté por primera vez un conjunto de curvas de gradiente de presién dinamico. Dichas curvas son aplicables para diferentes didmetros de tuberias, tasas de produccién, relaciones gas-liquido, entre otros. En 1961, Tek incluy el nimero de Reynolds Bifasico con el fin de correlacionarlo con el factor de friccion f , obteniendo resultados satisfactorios. También en 1961, Ros realizé un estudio basado en el calculo del gradiente de presién, que requiere conocer el Hold up liquido 1, y el factor de fficcion. Un analisis adimensional indicé que ambos, tanto el Hold up liquido como el factor de friccién, esta relacionado a nueve grupos adimensionales. Mas tarde, se demostré que solamente cuatro de ellos son realmente importantes. Basado en estos cuatro grupos, se pudo seleccionar un programa experimental que trataria de cubrir la mayoria de las condiciones encontradas en pozos de petréleo. Este programa experimental fue instalado en un laboratorio, donde se TaN TOIT Dep Legal WF LFTEOOTSOOBO7SFLUO MULTIFASICO EN TunERIAS on CAPITULO TI determinaron tres patrones de flujo, divididos en tres regiones: baja, media y alta presencia de gas. Los gradientes de presién en dichas Tregiones fueron presentados en forma de correlaciones, las cuales se compararon la informacion disponible del campo, mostrando excelentes resultados. H, fue relacionado con la velocidad de deslizamiento del fluido, la cual es la diferencia promedio real entre las velocidades del gas y el liquido. En 1961, Baxendell y Thomas, basados en el método Propuesto por Poettman y Carpenter de 1952, desarrollaron una correlacion para calcular el gradiente de presion cuando existe flujo bifdsico. La correlacién para el factor de pérdidas por elevacién y friccién, propuesta por Poettman y Carpenter, era inoperable en condiciones de alto caudal de flujo. Por ello, Baxendell y Thomas se vieron en la necesidad de realizar una serie de experimentos con caudales Por encima de los 5000 BD, usando censores electrénicos de presién a lo largo de la tuberia experimental vertical. Como resultado de estos experimentos, se establecié una correlacién entre el factor de pérdidas por elevacién y pérdidas Por friccién, y la masa del caudal de flujo, la cual se creyé era aplicable para alto tango de tamafio de tuberia y tipos de crudo con alto caudal (por encima de 900 BD y para tuberia 27/8" de didmetro extern). Esto anticipé que los calculos de gradiente resultantes tendrian una precision en el orden de mas o menos 5%, Fancher y Brown (1983) utilizaron la correlacién de Poettman. y Carpenter, considerando la RGL, como parametro adicional en el calculo del factor de friccién. Los resultados obtenidos indican que existen ciertas desviaciones en los rangos de tasas de flujo y RGL . De igual forma, numerosas curvas demostraron desviaciones a bajas tasas de flujo y RGL por encima de 700 pen! Bn. En 1963, Duns y Ros desarrollaron una correlacién basados en datos de laboratorios obtenidos en tubos plasticos. Duns y Ros observaron la influencia de los patrones de flujo en el gradiente de presién. Presentaron ecuaciones para calcular la densidad de la mezcla, Hold up Liquido y el factor de friccion, de SENSES EIRTT ap. Legal WL FOSIROOTSOORTTSFLUJO MULTIFASICO EN TUBERIAS. 93 CAPITULO — acuerdo al patron de fiujo existente. También, derivaron una correlacién para predecir la velocidad de deslizamiento entre las fases. En 1964, Hagedorn y Brown realizaron pruebas en pozos con 1500 pie de tuberia y didmetros de 1", 14", 14". Los pozos de prueba fueron equipados con dos valvulas de Gas Lift, cuatro censores de presién electrénicos con medicion de tasa de produccién liquida, tasa de inyeccién de gas, temperatura y por supuesto las presiones. Las pruebas se realizaron variando tazas liquidas de produccién, relacion gas-liquido y viscosidades del liquido. A partir de esa informacion, se construyé una curva de presién-profundidad para cada prueba y diametro de tuberia, de los tres escogides para la prueba. Desarrollaron correlaciones, las cuales permitieron hacer una prediccién del gradiente de presién dindmico para los diferentes didmetros de tuberia, condiciones de flujo y propiedades del liquido. En 1967, Orkiszewski combiné el trabajo de Griffith para flujo burbuja, el de Griffith y Wallis para flujo tapén y el de Duns y Ros para flujo neblina. Para el caso de fiujo tapén, desarrollaron nuevas correlaciones para el célculo de la densidad de la mezcla y el factor de friccion, utilizando un parametro denominado coeficiente de distribucién de liquido, el cual fue correlacionado con el didmetro de la tuberia, la velocidad superficial y la viscosidad liquida, usando los datos de Hagedorn y Brown. En 1969, Holmes y Brow hicieron un andlisis de los trabajos presentados por Hagedorn y Brown, Duns y Ros, y Orkiszewski. El método fue evaluado con una base de datos de 44 pozos. Resultados mostraron que el método de Orkiszewski era el de mayor precision. En 1970, Acurero y Bohérquez aplicaron el método de Hagedorn y Brown para el estudio y prediccién de pozos bajo flujo natural y mediante levantamiento por inyeccién continua e intermitente de gas. En 1971, Cardozo propuso dos métodos para calcular pérdidas de presion en pozos perforados direccionalmente. El primer método consistio en el aN IHT ap. Lapa N*LFOGTZO0TSOOBOTSZ FLUIO MULTIPASICO EN TUBERIAS 9 ‘CAPITULO MT céloulo del gradiente por friccién. El segundo método consideré intraducir una funcién en la correlacién de Hagedorn y Brown, obteniendo resultados satisfactorios. En 1972, Aziz, Govier y Fogarasi utilizaron el mapa modificado por Govier, Radford y Duns (1957) para identificar los diferentes tipos de patrones de flujo. Los autores derivaron nuevas relaciones para estimar el factor de entrampamiento de liquido y la densidad de la mezcla, para patrones de flujo tapén y burbuja. Los resultados de Zuber y Findlay (1965), Neal (1963), Wallis (1969), Griffith y Wallis (1961) fueron utilizados para derivar las relaciones formuladas. Ademas, utiizaron el esquema de Duns y Ros (1963) para los patrones de flujo transicién y neblina, En 1974, Chierici y Colaboradores presentaron un mapa de identificacion de los patrones de flujo similar al de Orkiszewski (1967). La Unica diferencia radica en el valor de la constante para definir los limites entre el flujo tapon y burbuja. Utiizaron los patrones de flujo burbuja, transicién y neblina del esquema de Orkiszewski, Utilizaron nuevas relaciones para definir la velocidad de levantamiento de las burbujas y las pérdidas por friccion para el patron de flujo tapén. Lawson y Brill (1974) presentaron una evaluacién estadistica de las correlaciones de Poettman y Carpenter (1952), Baxendell y Thomas (1951), Fancher y Brown (1963), Hagedorn y Brown (1964), Duns y Ros (1964) y Orkiszewski (1967). Concluyeron que la correlacién de Hagedorn y Brown resulta ser la mas adecuada en la prediccién del gradiente depresién tomada a 276 pruebas de pozes. Vohra, Robinsén y Brill (1975) realizaron la evaluacién estadistica de las correlaciones de Beggs y Brill (1973), Aziz y Govier (1972), Chierici y colaboradores (1974). Determinaron que la correlacién de Beggs y Brill resulto la de menor desviacién estandar, y la de Azis y Govier la de menor error porcentual. ap aga FOGTEDOTSOHOT TROTFLUJO MULTIFASICO EN TUBERIAS 9s ‘CAPITULO TIE Parra y Gémez (1975), presentaron una evaluacién estadisticas de las correlaciones de Poettman y Carpenter (1952). Baxendell y Thomas (1961), Fancher y Brown (1963), Duns y Ros (1964), Hagedom y Brown (1964), Aziz y Govier (1972), Chierici y colaboradores (1974), Beggs y Brill (1973), para diferentes rangos de didmetros de tuberias, gravedad AP! del petréleo, relacion agua-petrdleo, relacién gas-liquido. Entre las variables que afectan el graciente de presién en tuberia vertical, se tienen: a. Efecto del didmetro de la tuberfa: A medida que aumenta el diametro de la tuberia, disminuyen las pérdidas de presién a lo largo de la misma, Sin embargo, si la tuberia es muy grande, el deslizamiento entre las fases incrementaria el gradiente. b- Efecto de tasa de flujo: A medida que aumenta la tasa de flujo, las pérdidas de presin en la tuberla son mayores. Sin embargo, cuando [a tase de fujo decrece, el deslizamiento entre las fases causa también un incremento en el gradiente de presion. ce. Efecto de la relacién gas-liquido RGL: A medida que aumenta la RGL., |e presién de fondo fluyente disminuye hasta llegar a un minimo, @ partir del cual un aumento de la RGL. provocaria un aumento en la presién de fondo fluyente de. Efecto de la densidad del Liquido: A medida que aumenta la densidad del liquido, aumenta las pérdidas de energia. Mientras mas pesada seré la columna del fuido, la presién de fondo fluyente aumentara, disminuyendo asf la tasa de produccién. e-. — Efecto de la relacién agua-petrdleo RAP: A medida que aumenta la proporcién de agua en la columna de fluidos, esta seré mas pesada produciéndose un incremento en las pérdidas de presion f. Efecto de la viscosidad liquida 4: A medida que aumenta 4,, mayor son tas pérdidas totales de energia por limitacion en el movimiento de los fluidos g- Efecto del desizamiento: A mayor desizamiento entre fases, mayores seran las pérdidas de energfa en la tuberia. Entre las correlaciones mas importantes para tuberla vertical, se tienen. BEN STROSS Dap Lag NT TFSI DDOTSOORTSCAPITULO TIE A-. Correlacién de Ros (1961) Se clasifica como una correlacién del “Tipo c’. Ros demostré que el gradiente de presién depende del hold up liquido y del factor de friccion. Un analisis adimensional indicé que ambos, tanto el hold up liquido como el factor de friccién, estén relacionados @ nueve grupos adimensionales. Mas tarde, se demostré que solo cuatro de ellos son realmente relevantes. Basado en estos cuatro grupos, se pudo seleccionar un programa experimental restringido que cubrié practicamente con todas las condiciones encontradas en los pozos de petréleo. Este programa experimental fue instalado en un laboratorio, donde se determinaron tres patrones de fiujo, los cuales se dividieron en tres regiones baja, media y alta presencia de gas. Los gradientes de presién en esas regiones fueron presentados en forma de correlaciones, las cuales se compararon con la informacién disponible de campo, mostrando excelentes resultados. De acuerdo a Ros, las pérdidas debido @ aceleracién son muy pequefias, por lo que pueden ser despreciadas. A continuacién, se presenta la correlacién propuesta por Ros para estimar el gradiente de presion: eH (3), : (34) donde el término (AP/AH), representa las pérdidas por friccién y su valor sera determinado como una funcién del patrén de flujo. El procedimiento a seguir para estimar el gradiente por fricci6n es el siguiente: a. Determinarién de los grupos adimensionales La prediccién del gradiente por friccién prevé el uso de cuatro grupos adimensionales propuestos por Ros, los cuales son: SBN DDR ISTT Dep La NL FOSTIOOTOONITSFLUO MULTICASICO EN TUBERIAS ‘CAPITULOTIT hs Numero de Velocidad Liquido 1, Ny, =1.938V,, 4|P2 (3.5) Numero de Velocidad Gas ,, Ny 19987, Nmero de Didmetro de Tuberia N, N, =120.872d & (3.7) e, Numero de Viscosidad Liquida NV’, (3.6) Ny, =0.1572 4, > (3.8) donde V, y V,, representa la velocidad superficial de las fases gas y liquido, respectivamente, en pie/seg. d representa el diémetro de la tuberia, en pie. p, yo; representa la densidad liquida y la tensién superficial, en Ibm/ pies’ y dina lem, respectivamente. SBN 9S THOITISTT Dep Lega LRTOTSOOROTSinet _____ Lo MuutiFAsico EN TUBERIAS CAPITULO b-. Determinacién de los parémetros L,, Los parametros Z,, L, y L, se determinan como una funcién de N,,, a partir de Figs. 3.6 y 3.7, respectivamente. El parametro L, se determina como una funcién de WV, , a partir de Fig. 3.8. c-.Determinacién del patrén de flujo Ros clasificé los tipos de patrones de flujo, basado en el siguiente criterio: Patron de flujo burbuja (Ly + Ly Nw) > Np Patron de flujo tapén (+L, Ny) < Nyy (50+36 N,,) d-. Determinacién del hold up liquide #7, como una funcién del pawén de flujo 7 Patron de flujo burbuja El hold up liquido #7, se determinara mediante la siguiente ecuacién: Vat {Vg Va #4, V, ay 2y, : BINDERS Dep Legal LPT OOPSFLUJO MULTIFASICO EN TUBERIAS. ” CAPITULO MI 4 29 0 ol w 1 ie Figura 3.7. Comrelacin para Estimar el Parémetros L, , Propuesto por Ros (1961). RENTED ITT Dep LAIN" LROTHOOTSOROTSFLUJO MULTIFASICO EN TUBERIAS 10 — “CAPITULO IT Figura 3.8. Correlacién para Estimar el Parémetro 1, Propuesto por Ros (1961) En Ec. 39, 7s, My ¥ My Fepresentan la velocidad de deslizamiento y Ia velocidad de las fases gas y liquido, respectivamente, expresadas en ple/ sex Yr, se determina mediante la siguiente ecuaci6n: Ss (3.10) 1,938 4 Vor pb Yo; se encuentran expresadas en Ibm/ pies’ y dina em, respectivamente. La velocidad de deslizamiento adimensional S puede ser estimada mediante la siguiente correlacion: (3.11) (3.12) [SGN FR OWATISTS ‘ep LIN LRORTROOTSOOROTS_ FLUO MULTIFASICO EN TUBERIAS 101 ‘CAPITULO Los parametros F,, F,, F, y F, se determinan como una funcién de N,, a partir de Fig. 3.9. Figura 3.9. Correlacion para Estimar los Pardmetros F,, F,, F, y F, , ropuesto por Ros (1961). Patron de flujo tapén El hold up liquido #7, y la velocidad de deslizamiento V, se determinaré mediante Ecs. 3.9 y 3.10. La velocidad de deslizamiento adimensional S debe ser estimada mediante la siguiente correlacién. Net +R san) | Aer] (3.13) R=RN th (3.14) SEN OTRSRD-IE ISTE Dep Legal LOOTDOOTSONOTSFLUJO MULTIFASICO EN TUBERIAS 02 CAPITULO Los pardmetros F,, F, y F, se determinan como una funcion de W,,, 2 partir de Fig, 3.10. Figura 3.10. Correlacion para Esti war los Parimetros F,, F, y F, , Propuesto por Ros (1961). De acuerdo a Ros, el parémetro F, es una constante, cuyo valor no fue definido inicialmente. Mas tarde, Duns y Ros (1963) sugirieron un valor de 0.0029 para el parametro F Patron de flujo neblina El hold up liquide 17, se determinaré mediante la siguiente ecuacién: Has (3.15) Ei ROTIFLUJO MULTIFASICO EN TUBERIAS 103 ‘CAPITULO UT €. — Daterminaci6n del gradiente de presion por friccion El gradiente de presion por friccién puede ser obtenido mediante: OP) gp Maly, Ma (3) -4% Y (-f) (3.16) La Ec. 3.16 asume que la friccién es causada por esfuerzo de corte en la fase liquida. Esta suposicién se muestra razonable para los patrones de flujo burbuja y tapén, donde el liquido representa la fase continua, Sin embargo, si el gas representa la fase continua, el gradiente de presi6n por friccién deberd ser estimada sobre la base de la fase gaseosa como: OP) gy Pelal), Va (35) =45. aa (#4) (3.17) La correlacién para estimar el factor de fricci6n /f, esta dado por: & 3.18 AG (3.18) donde jf, puede ser estimada mediante el diagrama de Moody, como una funcion del numero de Reynolds de la fase continua y la rugosidad relativa. La Fig. 3.11 presenta el diagrama de Moody, como una funcién del numero de Reynolds Ng, y la rugosidad relativa ¢/d. El factor f, es una correccién por efecto de la relacién V,,/V, y esta dado como una funcién del grupo (f Cog Vu) N32) y la Fig. 3.12. SEN Sr oOAT IST Dap. Lage TUS ROOTSDODOTSFLUO MULTIFASICO CAPTULO MT TUBERIAS. m7 Factor de Fricelin Figura 3.11. Diagrama de Moody. Figura 3.12. Correlacin para Estimarel Factor f,, Propuesto por Ros (1961). SEN OTE SRI IHTT - ‘Dep Legal WCF DGTZOOTSOOROTSFLUJO MULTIFASICO EN TUBERIAS tos, ‘CAPITULO tT Es de notar que la Fig. 3.12 posee dos curvas: una punteada y otra continua. Para flujo vertical, la curva ha utilizar es la continua. El factor f, es una correccién de segundo orden por efecto de Ia viscosidad liquida y la RGL Este factor puede ser estimado mediante la siguiente ecuacién: RGL (3.19) 1th Correlacién de Duns y Ros (1963) Se clasifica como una correlacion del “Tipo c”. El método de Duns y Ros es el resultado de un estudio de laboratorio, donde mas de 4000 pruebas de flujo bifasico fueron obtenidos de una instalacién vertical de 185 pies. Los diémetros de tuberias utilizados comprendieron un rango entre 1.26"~ 5.60", incluyendo configuraciones de flujo anular. La mayoria de las pruebas estuvieron bajo condiciones muy cereanas a la presién atmosférica y se utilizo como fiuidos experimentales, aire como la fase gaseosa y agua e hidrocarburo como la fase liquida, E! hold up liquido #7, fue medido mediante el uso de trazador radioactivo y su observacién fue posible debido al uso de una seccién transparente en la facilidad experimental. Duns y Ros propusieron una serie de correlaciones para estimar el factor de tricci6n y la velocidad de deslizamiento, como una funcién del hold up liquido y el patrén de flujo. Un anélisis adimensional elaborado por Duns y Ros indicd que 12 variables eran de particular importancia en la prediccién del gradiente de presién. Mediante un proceso de eliminacién, finalmente se demostré que solo cuatro de ello son realmente relevantes y los mismos fueron utilizados para seleccionar el rango de variables en el programa experimental. La ecuacién para estimar el gradiente de presién, propuesta por Duns y Ros, es la siguiente: TaN RRO Dep Lagan LFOGIDOGTIOOROTSLAO Mu (3.20) (4)- a) donde (AP/A#) representa el gradiente de presion en Lpc/ pie. p, representa la densidad de mezcla con deslizamiento, en /bni/ pies’. (AP/AH), representa el gradiente de presién por friccién, en Lpe/ pie. E, representa el término de energia cinética adimensional El procedimiento a seguir para estimar el gradiente por friccién es el siguiente a. _ Determinacion de los grupos adimensionales Referido a los cuatro grupos adimensionales propuestos por Ros (1961): N N,y N, (Eos. 3.5.4 3.8). Jape BY be. Determinacién de los parametros adimensionales Los pardmetros , y L, se estimaran mediante la Fig. 3.6, propuesta por Ros (1961). Los parémetros 1, y Ly, se estimaran mediante las siguientes ‘ecuaciones: (3.21) 275484. (3.22) c. Determinacién del patron de flujo Duns y Ros clasificaron los tipos de patrones de flujo, basado en el siguiente criterio: Patron de flujo burbuja NRIOL Dap Ca LA DOOTSORDTTEFLUO MULTIFASICO EN TUBERIAS rd “ERPITULO UT Patron de flujo {apon (hy Ny )ENq Sly Patron de flujo transitorio Ly £M py Slay Patron de flujo neblina La Fig. 3.13 muestra los distintos patrones de flujo propuesto por Duns y Ros, que ocurririan en una tuberia vertical d-. _Determinacion del hold up liquido #7, y el gradiente de presion (AP/AH), como una funcién del patron de flujo Patron de flujo burbuja El hold up liquido #, se estima mediante similar procedimiento propuesto por Ros (1961) para flujo burbuja (Ecs. 3.9, 3.10, 3.11 y 3.12). El gradiente de presin (AP/AH) se estima mediante la Ec. 3.20, despreciando el término de energia cinética adimensional £, . El gradiente de presion por fricci6n (AP/ AH), se determina utilizando similar procedimiento propuesto por Ros (1961). Patron de flujo tapén Bajo este patron de fiujo, el hold up liquido H, y la velocidad de deslizamiento V, se determina mediante Ecs. 3.9 y 3.10, respectivamente, propuestas por Ros (1961) para flujo tapén. La velocidad de deslizamiento adimensional $' debera ser estimada mediante Ecs. 3.13 y 3.14. TINO EIST Dep Legal LAOS ROOTSCOTOTFLO MULTIFASICO EN TUNER 10 ‘CAPITULO UL Duns y Ros sugirieron un valor de 0.0029 para el pardmetro F, el cual era desconocide su valor en el trabajo original presentado por Ros en 1961, De manera similar que en el patrén de flujo burbuja, el gradiente de presion (AP/AH) se estima mediante la Ec. 3.20, despreciando el término de energia cinética adimensional £, Duns y Res (1863) Figura 3.13. Mapa de Patrones de Flujo en Tuberia Vertical Patron de flujo neblina Bajo este patron de “ujo se asume que la fase continua es el gas. Adicionalmente, Duns y Ros asumieron que no existe deslizamiento entre fases. En consecuencia, el gradiente de presién puede estimarse como: aN TERT Dap LET LRUGTEROTSOOROTSLUO MULTIFASICO EN TUDERIAS at ws ‘CAPITULOMT a (3.23) donde, (2) aap Pete (3.24) y AM, 2d Debido a que no se considera deslizamiento entre fases, el factor de friccion f se obtiene del diagrama de Moody (Fig. 3.11), como una funcién del nlimero de Reynolds para la fase gaseosa, solamente. Esto es, Vig Hs Nee (3.25) Duns y Ros observaron durante sus experimentos que sobre la pared interna de la tuberia se formaba una delgada pelicula de liquido. Las ondas que se crean sobre esta pelicula por fa accién de la fase gaseosa, genera perdidas adicionales en el calculo del gradiente de presion debido al in-remento de los esfuerzos de corte entre el gas y la pelicula liquida. Variando la rugosidad de la tuberia ©, Duns y Ros determinaron que el proceso es afectado por la viscosidad liquida y también es gobemado por el numero de Weber, el cual es po Soe + cio ON 89 definido como: TEN S7STBTE ATT Dep LRN LOST ZOTSOODOTSFLUO Mu ‘CAPHTULO donde ¢ es la rugosidad de la tuberia. V, y o, representan la velocidad superficial del gas y la tensién superficial del liquido, respectivamente. Como se aprecia en Ec. 3.26, la viscosidad liquida no ejerce un efecto directo sobre el Ny. En consecuencia, el efecto de la viscosidad puede ser considerado haciendo el WV,,, una funcién de un némero adimensional que contiene el término de viscosidad liquida, y el cual es dado como: N, = (3.27) 206 Basados sobre data experimental, Duns y Ros establecieron una relacién funcional entre el nimero de Weber Ny, y W,. Esta relacién puede ser apreciada en Fig. 3.14. El valor de la rugosidad puede ser muy pequefa, pero la rugosidad relativa nunca podra ser menor al valor de la tuberia misma (107) Nel Duns y Ros (1863) Figura 3.14, Efecto de la Viscosidad Liquida sobre Ny, como una Funcién de NV, TON TTR IST ‘De Lagat LPOGTEDOT SOROSFLUJO MULTIFASICO EN TUBERIAS, — CAPITULO Sobre la base de la Fig. 3.14, la relacién s/d puede obtenerse de acuerdo a las siguientes condiciones: Para Ny, N,, $ 0.005 £_ 0.0749 6; fo 3.28 dp, Vad on Para Ny, N, > 0.005 (3.29) o, ¥ p, representa la tensién superficial y la densidad del gas, en dinalcm y Ibm! pies®, respectivamente. d representa el diémetro de la tuberia, pie. V,, representa la velocidad superficial de la fase gaseosa, Bajo la condicién de 107 <¢/d <0.05, el factor de friccién / puede ser estimado directamente del diagrama de Moody. Cuando la relacién (¢/d) > 0.05, los valores del factor de fricci6n _ bajo condiciones de flujo neblina pueden ser obtenidos mediante la siguiente extrapolacién del diagrama de Moody. fz4 Z = +0.067 way (3.30) [4 log (0.27 (e/a) Ya que bajo condiciones de flujo neblina, el término de aceleracién no puede ser despreciada, Beggs y Brill (1973) propusieron la siguiente expresién matemitica para estimar el término de energia cinética E, a aga EFOSTROOTSOCRTTSFLUO MULTIFASICO BN TUMRiAS m2 CAPITULO MT Vn Vow Pre aa (3.31) donde P representa la presién del segmento y puede ser estimada como el promedio aritmético entre la presién a la entrada y salida del segmento en estudio, respectivamente. Desafortunadamente, la Ec. 3.31 puede estimar valores incorrectos de £,, cuando su valor sea superior @ 1. El gradiente de presién total podra estimarse entonces mediante Ecs. 3.20 y 3.31 Patron de flujo transicion Bajo este patrén de flujo, el gradiente de presi6n total viene dado por: (Se)-«(e),_ (ae) ae donde el gradiente de presin (AP/ Af) es el estimado bajo condiciones de flujo Tapén y Neblina, discutidos en secciones anteriores del método de Duns y Ros. El coeficiente 4 puede determinarse mediante la siguiente ecuacién: Ly Np ae (3.33) donde N,., Ly y ly se determinan mediante Ecs. 3.6, 3.21 y 3.22, respectivamente. A fin de incrementar la exactitud en estimar el gradiente de presion bajo el patron de flujo de transicion, se ha recomendado corregir la densidad del gas, utilizando la siguiente ecuacion: ope PeNe., (3.34) RN ITT Dep La LFUTIOOTSOOEOTSFLUO MULTIFASICO EN TUBERIAS us — “CAPITULO donde p, representa la densidad del gas a condiciones operacionales de presién y temperatura. Esta modificacién tora en cuenta la presencia de una parte del liquido en la fase gaseosa. Por otra parte, el gradiente de presion por aceleracién es despreciado, bajo este patron de fujo. C-. — Correlacion de Hagedorn y Brown (1964) Ge clasifica como una cortelacién del “Tipo b”, Se baso en informacion obtenida de un pozo vertical de 1500 pies de profundidad. Aire y agua fue utlizada como fluids experimentales. Las pruebas fueron llevadas a cabo en tuberias de 1", 1%" y 1%, donde se variaron ampliamente los valores de 1a tasa de fujo, relacién gas liquide y viscosidad del fluido, También utlizaron la base de datos expuesta por Fancher y Brown (1963) para tuberias de 2" Esta correlacién ha sido modificada con el tiempo. La primera modificacién establece que sie! criterio de Grifith y Wallis (1961) predice la ‘ocurrencia de flujo burbuja, entonces el método de flujo burbuja propuesto por Gritfth (1962) deberd ser utilzado para estimar el gradiente de presion (esta aproximacion es parte del método de Orkiszewsk’). La segunda madificacion se reflere al célculo del hold up liquido. Los valores de #7, obtenido de las figuras propuestas por Hagedorn y Brown son algunas veces inferiores al comparatlas con los valores de la fracci6n vacia de gas 4,. Para flujo multfasico ascendente, cesta condicién resulta contradictoria debido a que el liquido no puede viajar mas rapido que el gas. La solucion propuesta a esta anormalidad prevé el uso de a por H,. La tercera modificacién recomienda despreciar los efectos Por aceleracién, ya que se ha demostrado que este término sobre predice los cdlculos de caida de presién. La ecuacién general para estimar las perdidas por presién esta dado por: fen te (3.35) “Dead! aN oTEDRTTE ITT Bap LAN TF UGLRODTSDTSTUBERIAS, 0 MULTIEASICO ~CAPTULOM donde (AP/AH) representa el gradiente de presion en Lpc/ pie. Py. YP. representa la densidad de mezcla sin y con deslizamiento, respectivamente, en Ibm! pies’. V,, y d representan la velocidad de mezcla y el didmetro de la tuberia, en pie/seg y pie, respectivamente Por otra parte, el valor numérico de! factor de friccién f puede ser obtenido del diagrama de Moody (Fig. 3.11), como una funcién del Nimero de Reynolds y la rugosidad relativa, asumiendo la condicion sin deslizamiento entre fases. El procedimiento a seguir para estimar el hold up liquido es el siguiente a-. _ Determinacién de los grupos adimensionales La prediccién del hold up liquide prevé el uso de cuatro grupos adimensionales propuestos por Ros (1961), los cuales son: Numero de Velocidad Liquido N,,; Numero de Velocidad Gas ,,; Numero de Didmetro de Tuberia N,; y Numero de Viscosidad Liquida W,. Eos. 35 a 3.8, respectivamente. b-. _ Estime la variable N,,. mediante la Fig. 3.15 y el grupo adimensional NV, c- —_Estime el valor de la relacién H, /‘¥’, mediante la Fig. 3.16. d-. — Estime el valor de la relacion ¥’, mediante la Fig. 3.17, e-. _ Estime el valor del hold liquido H/,, mediante la mediante la siguiente relaci6n matematica: (336) TON STEHT IST Dap: Lage TF USIDOOTSOOROTEFLUJO MULTIFASICO EN TUBERIAS Hosodam y Brown (1985) ‘CARTUCO — 0.05 0.0) 3 = ! 8 b07 0.07 0.70 1.0 N ogee Bm (1988) joura 3.15. Coreen de Hagedomy Brown para Estimar NY rr) Correlation based on ‘oak tubing sizes of 1 to 2in. viscosities of 0.86 to 110 cp fh o6 oa o2 ° 1077 10-8 10-5 10-¢ 1073 1072 wy, Wey" Me wllz), a Figura 3.16, Comelacién de Hagedorn y Brown para Estimar H7, / TNR RTT IST Dep-Legat *LFOSTROOTECORITS20 18 16 v 14 12 10 00.01 002 0.03 0.04 0.05 0.06 0.07 008 0.09 1.0 VN) Ni Hogedam y Brown (1685) Figura 3.17. Correlacin de Hagedorn y Brown para Estimar Y f. Verifique que se cumpla que H, 2h, G37 En caso contrario, se debe considerar #7, = 2, Finalmente, el gradiente de presién se obtiene mediante Ec. 3.35. D-. Correlacién de Orkiszewski (1967) Se clasifica como una correlacién del “Tipo c’. La correlacién de Orkiszewski es el resultado del andlisis de varios métodos publicados. El objetivo consistié en determinar si alguno de ellos proporcionaba una prediccién del gradiente de presién con mayor precisi6n, para un amplio rango de condiciones. El método de Orkiszewski fue dividido en tres categorias, cuya TAN SUAS a Dap Lag CFDSFLUIO MULTIFASICO EN TUBERIAS “CAPITULO MT discriminaci6n estuvo basada en la forma en que se consideré el hold up liquide i, para el célculo de la densidad, el factor de friccion y los patrones de flujo. Primera categoria: el hold up liquido 17, no es considerado en el calcul de la densidad. La densidad es simplemente la densidad de los fluidos producidos, corregidos por presion y temperatura, 17, y las pérdidas por friccion son expresadas mediante una correlacién empirica del factor de friccion que considera ambos efectos. No se realizan distinciones entre los patrones de flujo Segunda Categoria: el hold up liquido #, es considerado en el calculo de la densidad y es considerado separadamente 0 combinado en alguna forma con fas por fricién las cuales a su vez se basan sobre las propiedades las pé compuesta de las fases gas y liquide. No se realizan distinciones entre los patrones de flujo. Tercera Categoria: considera el hold up liquido H, en el célculo de la densidad y se determina a partir de la velocidad de deslizamiento. Las pérdidas por friccién se determinan a partir de las propiedades de la fase continua. Se distinguen cuatro patrones de flujo: burbuja, tapon. transicion y nebiina. La metodologia utilizada por Orkiszewski le permitié seleccionar aquellos parémetros que permitiesen mayor exactitud en el célculo del gradiente de presion, Propuso una nueva correlacién para la condicion de fio tapén, a partir de la data obtenida por Hagedorn y Brown (1964). Seleccioné el método de ara la condicién de fujo burbuja. Finalmente, el método de Duns y Ros (1963) fue utlizado para el patron de flujo nebiina. El Galculo del gradiente de presin dependeré del patron de flujo. El procedimiento para estimar este gradiente se describe detalladamente a continuacion: ae. Patron de fujo burbuja Este patrén de flujo existe si se cumple el siguiente erteio: (1-4) Le» 43) an G44) El gradiente de presién total resulta de la suma del gradients de presion por elevacién y friccién. El gradiente de presién por elevacién considera un procedimiento particular para estimar la densidad de la mezcla, Bajo el patron de flujo tapén, la densidad se estima mediante la siguiente ecuacion: Pi Vy + VAI 4 Pe Vor +a, 3.45) VV i (3.45) BON STAT Dep Legal LFDSIBOOTSOOIOTSFLUJO MULTIFASICO EN TUBERIAS a0 CAPITULO TT donde /, representa la velocidad de ascenso de una burbuja y Tel coeficiente de distribucion liquide. Orkiszewski desarrollo la Ec. 3.45 para tratar de estimar la densidad promedio, que considere simulténeamente la presencia de la burbuja de Taylor (1949) y e! tapén de liquide. Un oriterio muy similar tuvo Griffith y Wallis, solo que ellos despreciaron la presencia de una pelicula liquida alrededor de la burbuja de Taylor y la posibilidad de que gotas de liquide se encuentren dentro de esta, El Ultimo término de la Ec. 3.45, propuesto por Orkiszewski, toma en cuenia la distribucién del liquido en la burbuja y el tapén de liquide. Esta modificacién fue importante para extender el trabajo de Griffith y Wallis a condiciones de alta velocidad de flujo. De acuerdo a Griffith y Wallis, la velocidad de ascenso de una burbuja ¥, puede ser estimada mediante la siguiente ecuacion: Cr Ved. 3.46) donde d es el didmetro de la tuberia. Las constantes C, y C, son obtenidas de las Figs. 3.18 y 3.19, como una funcién de N,., ¥ Na,, . fespectivamente, Estos nuimeros de Reynolds son definidos como: Ny, =P, G47) 3.48) TSN 9789804 Dap Legal LF DGIOOTSOOSOTSFLUO MULTIFASICO EN TUBER! CAPITULO MT FLUJO TAPON (SLUG) oa 03 4 a 4 a 4 °o 10 20 30 cy 50 New, avd & tn y as 1808) Figura 3.18. Correlacién de Griffith y Wallis para EstimarC, m1 FLUJO TAPON (SLUG) ce y wats (1981) [SON SESE IT Figura 3.19, Coreacion de Griffith y Wallis para EstimarC ap La FSI 200TSOCATTSee cee Cuando Ia constante C, no pueda ser obtenida de Fig. 3.19, 1, puede ser calculado mediante el siguiente criterio: Cuando N,,, £3000 0.546 + 8.7410 Ny.) le d 6.49) Cuando N,,, > 8000 V, = (0.35+8.7410% Ny. Vga ~ (3.50) Cuando 3000 < N,,, £8000 B51) donde, (0.251+8.7410% Ny, Ved. (3.52) Debido a que V, ¥ Nye, Se encuentran interrelacionado, la determinacién de ¥,, requerira de un proceso iterative cuando se utilicen las Figs. 3.18 y 3.19 6 las Es. 3.49 a 3.52. El procedimiento a seguir, se lista a continuacién: 4-._ Estime un valor de ¥,. Una buena aproximacién seria: V,=05 Jed G53) TSUN SIEBOT TTT ‘a Legal LPOSTROOTROOEOTSFLUIO MULTIFASICO EN TUBERIAS ma za ‘CARTULO tt 2. Caloule N,,, ullizando el valor de 7, , obtenido en paso 1 3. Determine las constantes C, y C, y caleule nuevamente ¥,, utilizando Ec. 3.46 0 alguna de las Ecs. 3.49 a 3.52, dependiendo cual sea el caso. 4. Compare los valores de ¥,, obtenidos en paso 1 y 3. Si no satisfacen cl criterio de convergencia establecido, utilice el valor obtenido en el paso 3 como el nuevo valor supuesto y continue en paso 2. Repita el procedimiento hasta 10 se observa cambios entre los lograr convergencia, la cual sera aleanzada valores estimados y caloulados de ,,, Por otra parte, Orkiszewski utilizd la data de Hagedom y Brown para calcular y correlacionar el coeficiente de distribucién liquida I’. Sin embargo, la determinacion de este coeficiente demanda definir la fase continua liquida. Dependiendo de cual sea la fase continua y la velocidad de mezcia, el valor de Tr podrd ser estimado mediante el siguiente eriterio: Fase continua Agua: RAP >4 y V,, <10 pie/ see 0.013 log 1 _ 9,681 + 0.232 log V,, ~ 0.428 log d. (3.54) Fase continua Agua: RAP = 4 y V,, 210 pie /see 0.045 log 4 Fee a V,, - 0.888 log d 3.55) Fase continua Petréleo: RAP <4 y V,, <10 pie! sec 20127 Woes = o284 + 0.167 log V,, + 0.113 log d 56) TaN 759 ASTT Dep LER N RRSTROOTSOOROSFLUIO MULTIFASICO EN TUHERIAS ma Fase continua Petroleo: RAP <4 y V,, 210 pie/sec 0.0274 Nog (44, +) _ 7 MA Ng Hs +) 9161 + 0.569 ogd +X. G57) donde, +0 X=-logl,, [Sorte =0). 039700 log d (3.58) Las Eos. 3.54 a 3.58 consideran a viscosidad liquida j,, el diametro dy la velocidad de mezcla V, en cps, pie y pie/seg, respectivamente. Con el objeto de eliminar discontinuidades de presién entre los distintos patrones de fiujo, el valor de T’ se encuentra restringido a: SIV, <10 pie/sec, entonces T > — 0.065 ¥,, G59) i ¥, Ps Si V,, > 10 pie/sec, entonces 2 1-2 3.60) WatNs A El gradiente de presién por friccién, bajo flujo tapén, se determina mediante la siguiente ecuacion. ar) fave |(Yar%),r], ash AH 2d |\Wth donde / se obtiene del diagrama de Moody, como una funcién del numero de Reynolds definido por Ec. 3.48. Finalmente, la caida de presién total es la suma TaN OATS ap aga NP FGI OOTSOOSNTSLUO MULTIFASICO EN TURERIAS ns CAPITULO tt dol gradiente de presi6n por elevacién (asumiendo deslizamiento entre fases) ¥ el gradiente de presién por ficcién (Ec. 3.64). Bajo condiciones de fujo tapon el gradiente de presion por aceleracién es despreciado. Algunas discontinuidades de T pueden tener un efecto significativo, tal como puede apreciarse en Fig 3.20. BPs vase er é oe a = er a 3 4050 [SUPERFICIAL MIXTURE VELOCITY, tse a (1988) Figura 3.20. Discomtinuidades del Coeficiente de Distibucién Liquida T El método de Orkiszewski puede causar problemas de convergencia, Esto se debe a discontinuidades entre las Ecs. 3.54 y 3.55 y entre las Es. 3.56 y 3.57, para el agua y el petroleo como fase continua, respectivamente. Brill (1989) demostr6 que los limites establecidos entre las Ecs.3.59 y 3.60 no son suficientes para eliminar las discontinuidades de presion y suginio ullizar ta modificacién propuesta por Triggia (1984), la cual se encuentra dada por: Fase continua Agua: RAP >4 y V,, 10 pie/sec 0.013 log 4 Si — 0287 - 0.162 08 Vy — 0.428 log d (3.62) aN eT HTD ‘Bop aga LF OGTEOOTSONIOTS26 FLO MUL TIPASICY CAPITULO MT Fase continua Petréleo: RAP <4 y V,, >10 pie/ see 27 log Tie D oT + 0.113 logd + C(I log ¥,). donde 0.01 log (44, +1) Gra +0:397 + 0.63 log d (3.64) El coeficiente de distribucién liquido T puede que resulte negativo para altas tasas de flujo, debido a altos valores de V,. De ocurrir esta condicién, Orkiszewski propone reemplazar p, por p,., 10 que sugiere modificar la Ec. 3.61. La ecuacién resultante estaré dada por. AP) _f 2, Ve 3 (¥) yet (65) donde f se obtiene del diagrama de Moody, como una funcién del numero de Reynolds, el cual es definido por: (3.66) Bajo patrones de flujo Neblina y Transic'sn, Orkiszewski recomend6 estimar el gradiente de presién de la misma manera como se estiman en la correlacién de Duns y Ros (1963). TaN ORT ITT ‘Dep egal LOGI ROOTROTAOTSuorKias a ‘CAPITULO MT Gorrelacin de Aziz, Govier y Fogarasi (1972) Se clasifica como una correlacién del “Tipo c’. Aziz et al. propusieron un método dependiente del tipo de patron de flujo. Presentaron nuevas cortelaciones para estimar el hold up liquide #7, y la densidad de mezcla p,,. bajo condiciones de flujo burbuja y tapén. Bajo el patron de flujo neblina y transici6n, Aziz et al. recomendaron el método de Duns y Ros (1963). Para identifcar los patrones de flujo, Aziz et al. analizaron varios métodos disponibles en la literatura y encontraron que el método presentado por Govier, Radford y Dunn (1957) era el mas conveniente. La Fig. 3.21 muestra el mapa de patrones de flujo presentado por Govier et al 10 | My 10 ubbie sug ou 0.01 on 15 0 108 Ny aie tin (1972) Figura 3.21, Mapa de Patrones de Flujo Utlizado por la Coreeacidn de Aziz al. (1972) ‘A continuacién, se describe e! procedimiento a seguir para estimar el gradiente de presi6n: ISON 975950 RTD ‘ap. Legal N° LFOT2007SOODOTSFLUO MULTIFASICO IN TUBERIAS ns CAPITULO a a. _ Determinacién del patron de flujo El patron de flujo puede ser determinado graficamente, mediante el uso de Fig. 3.21, 0 numéricamente mediante la determinacién de las siguientes variables: a)" ((2\( ay fal ((E)@s] - Co &)] - (3.68) Ny, =0.51[100.N,P"", (3.69) 6+3.8.N,, (3.70) = 70[100n, J", (71) donde la velocidad superficial de las fases gas y liquido, V,, y Vy respectivamente, se expresan en pie/seg. p, Y PY o) S@ expresan en Ibm/ pies’ y dinalem, tespectivamente. El criterio para determinar el patron de flujo establece que: Patron de flujo burbuja ISIN FR-DHOTE ISTE Dep Legal LAOOTOTSOOOTSFLUJO MULTIFASICO EN TUBERIAS, 9 CAPITULO MT Patron de flujo tapén N,4 Patrén de flujo neblina Ny >N, para N, <4 Ny >265 para Ny >4 Patron de flujo transicién Ny 4, el patron de flujo no existe. b. _ Determinacién del gradiente de presién, dependiendo del patron de flujo Patron de flujo burbuja Bajo este patron de flujo, el gradiente de presion por friccién se determina como: LV oan 62) donde f se obtiene del diagrama de Moody, como una funcién del ntimero de Reynolds y la rugosidad relativa. N,, es definido por: SBN 97TH EIT Tap Legal NL FOSTIDEPSOOIOTSFLUO MULTIFASICO EN-TUBERAS 0. CAPITULO ME (3.73) La caida de presién total es la suma del gradiente de presién por elevacion (asumiendo deslizamiento entre fases) y el gradiente de presién por friccién (Ec. 3.73). Ya que el gradiente de presién por elevacion requiere conocer el valor del hold up liquido #,, bajo este patron de flujo, esta variable es calculado como: (3.74) donde V, es la velocidad de ascenso de las burbujas en un liquido en movimiento y puede estimarse mediante la siguiente ecuacién: 1.2 Py + Vig 6.75) El término ¥,, representa la velocidad terminal de las burbujas en un fluido estatico v puede estimarse mediante la siguiente ecuacién (3.76) donde, p, ¥ 9,.¥ 0, Se expresan en Ihm/ pies’ y Ihm/ seg”, respectivamente Patron de flujo tapén Requiere de un tratamiento muy similar al patrén de flujo burbuja. Bajo este patron de flujo, el gradiente de presién por friccién se determina como: TEN STRSTR Dap aga EROS TDOOTSOOIOTSZ a FLUIOMULTIEASICO EN TUBERIAS at CARTTULO G71) donde f se obtiene del diagrama de Moody, como una funcién del ntimero de Reynolds y la rugosidad relativa. V,, es definido mediante Eo. 3.73. En consecuencia, la caida de presi6n total es la suma del gradiente de presion por elevacién (asumiendo deslizamiento entre fases) y el gradiente de presién por friccién (Ec. 6.76). La caida de presién por aceleracién es despreciada bajo este patron de fiujo. El valor del hold up liquido H, es estimado mediante Ecs. 3.74 y 3.75, con la particularidad que la velocidad terminal de las burbujas en un fluido estatico V1, , se estima mediante la siguiente ecuacién (3.78) donde p, y p,, y d se expresan en /bm/ pies’ y pie, respectivamente. La constante C representa un factor de proporcionalidad, el cual seguin Wallis (1969) puede estimarse mediante la siguiente ecuaci6n: =0345 [1-2 G79) donde (3.80) ISSN TSOTSI Dap egal LOST OOTSORRTTSFLUO MULTIPASICO EN TUBE ua TE CRPTTULO TT G81 donde p, y 9. ¥ 9, se expresan en lbm/ pies’ y dinaem, respectivamente. El valor de m se determina sobre ta base del siguiente criterio: SiN, 250, entonces m=10, SI 250>N, >18, entonces m=69.N,*”, SIN, <18, entonces m=25 Patrén de flujo neblina ‘Aziz et al. recomendé utilizar el mismo método de Duns y Ros (1963), para estimar el gradiente de presion. Patron de flujo transicién ‘Aziz et al. recomendé utilizar el mismo procedimiento descrito por Duns y Ros (1963) para estimar el gradiente de presién bajo este patron de flujo, con la particularidad que el coeficiente 4 (Ec. 6.32) se determinaré de la siguiente manera 3.3.2. Curvas de Gradiente en Tuberia Horizontal El primer estudio de flujo bifasico horizontal fue realizado en 1949 por Lockhart y Martinelli. Este estudio asumié que la caida de presin generada por la fase gaseosa era igual a la generada por la fase liquida. Igualmente, se NOSIS ep. Legal LFOOTEOOTSOOROTSFLUJO MULTIFASICO EN TUBERIAS ey cari establecié que durante el flujo simultaneo de gas y liquid podian existir cuatro patrones de flujo, para las cuales propusieron 4 correlaciones: Ambas fases en flujo laminar; ambas fases en flujo turbulento, una fase en turbulento y otra en laminar, y una fase en laminar y otra en turbulento, Esta correlacién es considerada, sobre la base de sus resultados, muy buena para bajas tasas de flujo y didmetro pequefios de tuberia. En 1949, Bergelin y Gazeley presentaron un trabajo experimental en el que se describe la existencia de cinco patrones de flujo. Estos autores concluyen que la correlacién de Lockhart y Martinelli no es totalmente adecuada para estimar el gradiente de presion También en 1949, Kosterin realizo un trabajo teérico en el que introduce la idea de un factor de friccién bifésico, similar al desarrollado por Moody para flujo monofasico. Ademas, present6 una correlacién para estimar los patrones de flujo. En 1952, Jonson y Abou Sabe publicaron los resultados de un trabajo experimental en el cual construyen un grafico para predecir los patrones de flujo. En 1953, Schneider presenté los resultados de un trabajo experimental en el cual se desarrollé una correlacién para determinar el factor de friccién bifasico, basandose en 'a misma idea propuesta por kosterin. En 1954, Baker presenté un trabajo basado en data de campo, donde se Gescriben siete patrones de flujo y la ecuacién respectiva que permite determinar la caida de presién en cada uno de estos patrones. Esta correlacién dio muy buenos resultados bajo condicién de flujo tapén de gas. Esta correlacion es basicamente muy similar a la propuesta por Lockhart y Martinelli, con la unica diferencia en que se introduce el concepto de patrones de fiujo y se propone correlaciones para estimar el gradiente de presi6n en cada uno de estos. En 1955, Chenoweth y Martin llevaron a cabo un trabajo experimental para comprobar la correlacion propuesta por Lockhart y Martinelli, Se evaluaron 264 pruebas de laboratorio, en tuberia de gran diémetro y presiones promedio de 100 Zpe. Demostraron que la correlacién de Lockhart y Martinelli perdia SON OE DOA Da Legal LPOSROSONGOTSFiuuo Mun TUDERIAS_ a precision a medida que el didmetro aumenta. Ademas, obtuvieron una correlacion limitada para flujo turbulento, la cual puede ser aplicada en casos especiales. En 1955, White y Huntington, sobre la base de un trabajo experimental, propusieron seis correlaciones diferentes para seis diferentes tipos de fiujo, visualizados a través de tuberlas plasticas. Los gradientes de presién se correlacionaron por medio de un factor de pérdida y la tasa de flujo madsica. Finalmente, las correlaciones de White y Huntington solo son aplicables bajo presiones iguales a 1 am Bertuzzi, Tek y Poettman (1955) desarrollaron un balance de energia para flujo bifésico horizontal. Con datos experimentales calcularon las_pérdidas totales de energia debido a irreversibilidades. Introdujeron una funcién de Reynolds que puede correlacionarse con el factor de friccion o factor de pérdida de energia. Construyeron un grafico de factor de friccién contra Reynolds, obteniendo las curvas “K" como pardmetro; estas curvas presenten suposiciones y zonas transicién, muy similares al diagrama de Moody para flujo de una sola fase. Esta correlacién mostré pocos resultados satisfactorios. Chisholm y Laird (1958) llevaron a cabo un trabajo experimental, lo que les permitié proponer una modificacién al método de Lockhart y Martinelli (1949), para extenderla al uso de tuberias con diferentes grados de rugosidad. En un analisis de esta correlacién por Baker, hace notar que la mayoria de las lineas de conduccién de gas condensado caen en el rango excluido por Chisholm y Laird, lo cual limita el rango de aplicabilidad practica de esta correlaci6n. Chavez (1959) desarrollé una correlacién en la cual relaciona las diferentes variables que afectan al flujo bifasico en tuberias horizontales. Fue preparada con un amplio rango de datos de campo, por lo que su rango de aplicacién es mayor que el de las correlaciones anteriormente mencionadas. Chavez basé su estudio en las correlaciones de lockhart y Matinelli y Baker. Hoogendom (1959), mediante un trabajo experimental, estudio la caida de presion en flujo horizontal bifésico, utlizando tuberias de 24, 50, 91 y 140 TaN Om apap WEFT A00TSOOROTSmilimetros de didmetro. Midié el factor de entrampamiento de liquido "Hold Up", lo cual permitié calcular las verdaderas velocidades promedio de las fases. Determiné que la caida de presién por aceleracién puede tener un efecto de hasta el 15% sobre la caida de presi6n total, en tuberias de pequeito diametro. Dukler et al. (1964) publicaron dos trabajos que trataban sobre el fiujo multifasico en tuberias horizontales. En el primer trabajo, se hace una comparacién de las correlaciones presentadas por Baker, Bankoff, Chenoweth y Martin, Lockhart y Martinelli, y Yagi. Concluyeron que las correlaciones de Bankoff y Yagi son completamente inadecuadas, ademas se observé que en las correlaciones de Chenoweth y Martin y en la de Lockhart y Matinelli, hay una tendencia casi uniforme, destacandose ademas ciertas desviaciones en la medida en que el didmetro de la tuberia aumenta. Se demostré que las correlaciones que toman en cuenta el factor de entrampamiento son més exactas que las que dependen del conocimiento del tipo de flujo. En el segundo trabajo, se propone dos correlaciones: una aplicable cuando no existe deslizamiento entre fases, es decir, se supone flujo homogéneo, la otra aplicable cuando existe deslizamiento interfase. En ninguno de los casos estudiados, se considera los patrones de flujo. Eaton et al. (1966), basados en un trabajo experimental, propusieron una correlacién basada en un balance de energia para flujo multifasico. Propusieron correlaciones para determinar el hold up liquido y el factor de friccién Consideraron a las fases como una mezcla homogénea de propiedades promedias, evitando con ello considerar los diferentes patrones de flujo. Eaton afirmo que el cambio de un patron de flujo a otro era continuo, y no causaba una abrupta discontinuidad en las pérdidas totales de energia, con lo cual le quitaba importancia a la existencia de los factores o patrones de flujo. Se recomienda el uso de esta correlacién s6lo cuando el efecto de viscosidad sea despreciable, puesto que presenta fallas para el calculo de este factor o parémetro. Guzhov et al. (1967) desarroliaron una corretacién para calcular el factor de entrampamiento del liquido. El experimento consistié en hacer fiuir aire-agua ap. Legal WLFOGIDODT_ FLUO MULTIFASICO EN TUBERLAS 136 a través de tuberias de 15 pie de longitud, diametros que variaban entre 1" 2". y angulos de inclinacion de 0-10 grados. Las correlaciones presentadas para estimar el hold up liquido son vélidas sélo para flujo tapon, En 1967, Andrews lieve a cabo un estudio experimental y propuso una correlacién para lineas de 2" de didmetro, solamente. Degance y Atherton en 1970, realizaron una revisién de los métodos publicados hasta esa fecha y presentaron 7 articulos diferentes sobre todos los aspectos relevantes sobre flujo multifasico Trela en 1974 presenté una correlacién, en la cual se describe el uso de Un factor de friccién local y fa fraccién vacia de gas 2, para la determinacion de la caida de presién. Derivé dos relaciones analiticas simples para determinar sus valores en cualquier clase de flujo bifasico Armand (1976) correlacioné datos de caida de presién para flujo horizontal aire-agua, a presiones alrededor de} aim, cubriendo un amplio rango de calidad y de velocidad. Reconoce la validez de sus correlaciones solo a las condiciones de los datos en las cuales fueron desarrolladas En 1976, Friedel presente una correlacién para predecir la caida de presién en flujo bifésico. Esta correlacién esta basada principalmente en las propiedades del fluido y en los parémetros definidos por la geometria de la tuberia. También tomé en cuenta la tasa de flujo. Lombardi y Ceresa en 1978 presentaron una correlacién para el calculo de la caida de presién en flujo bifésico, basados en un balance de energii aplicado a un elemento infinitesimal de fluido, en una unidad de tiempo dado y atravesado por un flujo de masa bifésica. Ofiemans en 1979 publicé una correlacién en la cual describe el flujo de gasicondensado y de gas/petréleo en tuberia horizontal, Utilizé una tuberia de 30 pulgadas de diametro, con una longitud de 30 Km, operando en una presién igual o mayor de 1431 Lpca. Para determinar los patrones de flujo, se utiliz6 el mapa propuesto por Taitel y Dukler. La caida de presién se calculé utilizando un TaN Eom. aSHT Dep Lega NL FOGT ROTORSPLUJO MULTIFASICO EN TUBERIAS CAPITULO MT modelo simple, que esta sujeto a los limites correctos del flujo monofésico. Este modelo se encuentra en funcién del hold up liquide y hace uso de la expresion de Colebrook para estimar el factor de friccién bifasico. EI Método del Comité AGA-API permite predecir la caida de presion en un sistema bifasico, dentro de una aproximacién de mas o menos del 15%. Este método considera que la caida de presién total resulta de los gradientes de presion por fticci6n, elevacién y aceleracion. EI gradiente de presién por aceleracién se considera despreciable para tuberias de diametro superior a 4" Las variables que afectan las curvas de gradiente de presion en tuberia horizontal son esencialmente las mismas tomadas en cuenta en tuberia vertical. Exceptuando el hecho que las pérdidas de presién por efectos gravitacionales son despreciadas. Las caidas de presién en tuberias horizontales pueden llegar aserde 5 a 10 veces mayores que las ocurtidas en flujo monofasico, debido a que la fase gaseosa se desiiza, generaimente, a mayor velocidad sobre la fase liquida, incrementando las pérdidas. Entre las principales aplicaciones practices, el gradiente de presion en tuberia horizontal es utlizado en la industria petrolera bésicamente para determinar la contrapresién necesaria en el cabezal del pozo para llevar los fluidos producidos hasta el separador. La metodologia utilizada para estimar la caida de presién en una tuberia horizontal es basicamente el mismo al utlizado en tuberia vertical. Entre las correlaciones mas importantes para tuberia horizontal, se tienen’ A-. Correlacién de Dukler, Wicks y Cleveland (1964) El trabajo de Dukler et al. consistio de dos partes. La primera parte, clasificada como una correlacién del “Tipo a’, considera que no existe deslizamiento entre las fases liquide y gas, y no se toma en cuenta el patrén de flujo, lo que supone flujo homogéneo, De acuerdo a Dukler et al., el gradiente de presién total (AP/AL) puede ser estimado mediante la siguiente ecuacion: TERN FDR Dep Legal NLROSTROOTSOOROSFLUJO MULTIFASICO EN TUR ow uaa CATON (2) (AP/AH), se aH) I-a * ey donde (AP/AIN, representa el gradiente de presién por friccion y a el termine de aceleracién. (AP/AI/), puede ser estimada mediante la siguiente ecuacion: (at) -2Lae? vi SH), &Pmd * metro de la tuberia y la densidad de la mezcla donde d y p, representan el sin considerar deslizamiento entre fases, respectivamente. El factor de friccion {7 puede estimarse mediante la siguiente correlacién, propuesta por Dukler et al 7 =0.00140 2, (3.85) Re 6.86) donde 1, representa fa viscosidad de la mezcla, sin considerar deslizamiento entre fases. M; define la tas masica total y esta dada por: M;,=M,+M,, (3.87) 5.6152) Go Ba + 4u) 86400 G88 TON DIR TE TTT Dep LAIN LF OTSOODOTSFLUJO MULTIFASICO EN TUBERIAS, 1 ~ ‘CAPITULO MT at Be Pe RGLG~R Go) : $6400 3.89) El término de aceleracion @ puede ser estimado mediante: 16M; M, P wd RPP,’ foe donde P representa la presién promedia del segmento. 7, y P, representan la presién a la entrada y salida del segmento, respectivamente. La segunda parte de la correlacién de Dukler et al. es clasificada como una del “Tipo b’, ya que considera deslizamiento entre las fases liquide y gas, aun cuando no toma en cuenta el patrén de flujo. En este caso, se requiere determinar tanto el factor de entrampamiento del liquide, como el factor de friccion bifésico. De acuerdo a Dukler et al., el gradiente de presién total (AP/ AH) puede ser estimado mediante la siguiente ecuacién: AP AP (2)-(G) G91) donde (AP/AH), representa el gradiente de presién por friccion y £, las pérdidas por aceleracion. (AP/AH), puede ser estimada mediante la siguiente ‘ecuaci6n: Sf PaVe (i), od” (3.92) ENGST ap egal LPOG DOTSDOANTSFLUJO MULTIFASICO EN TUBERIAS 10 : ‘CAPTOLO MT donde 1, y d representan la velocidad de mezcla y el didmetro de la tuberia, respectivamente. La densidad de mezcla p, puede ser estimada, de acuerdo a Dukler et al., como: 3.93) Por otra parte, el factor de friccién f puede estimarse mediante la siguiente correlacién, propuesta por Dukler et al.: L=Sq[l-Ln(4/5)], (3.94) S=1.281-0.478y + 0.444 y? — 0.094 y° +0.00843y". (3.95) y=-Lnh, 3.96) El factor de friccién sin considerar deslizamiento entre las fases f,, puede ser estimado mediante la siguiente correlacién, propuesta por Dukler et al. 0.125 Jug = 0001404 os a7) ae Ode -(2] de® a 6.98) El hold up liquido H, se determina mediante el siguiente procedimiento de ensayo y error: a Se supone un valor inicial de #7, . Un valor inicial podria estimarse como: {SBN 97RD TTT Dep Lega LAG TOOTSOREOTSFLUO MULTIFASICO EN TUBERIAS a ‘CAPITULO MT H,=4,+0.001, 3.99) be Determine un nuevo valor de (N,.),, mediante Eo. 3.98 cc. Estime un nuevo valor de 11, mediante la Fig. 3.22 d-. Compare los valores obtenidos del hold up liquido 17, , en pasos “a y “co” Si la diferencia entre ambos es mayor al valor de la tolerancia considerada, asuma el valor del /1,, estimada en paso en "c’, como el asumido y repita los pasos “b” a “d’, hasta lograr convergencia. 1" » Dull etal (1964) Figura 3.22, Corrlacién para Estimar el Hold Up Liquido, Propuests por Dukler etal. (1964), Finalmente, Dukler et al. propuso Ia siguiente expresion para estimar las pérdidas por aceleracion F, ENDED AIST aeceaan Dep Lega LF 20OTSOCSOTS“Fuso MuLTIFASiCO ex Tune’ we 3.100) 1 1 a| ——— |-| —— 3 - (atm) (atm ( al: ae (eullan) an] 8 AH) (oH) \e Hib aa Be. Correlacién de Eaton, Andrews, Knowles, Silberberg y Brown (1966) Se clasifica como una correlacién del "Tipo b’. Eaton et al. realizaron una serie de estudios experimentales, utilizando lineas de pruebas de 1700 pies de 2"a" y 17" de didmetro. Como fluidos experimentales, se utiliz6 agua, destiado y petréleo de forma separada como la fase liquida y gas natural como la fase gaseosa. Se midié el hold up liquide por medio de valvulas de cierre rapido. La informacion de campo permitié desarrollar_correlaciones generalizadas de gradiente de presion, hold up liquido y patrones de flujo, mediante un balance mecénico de energia para ambas fases. Eaton et al propuso una ecuacién para estimar el gradiente de presion total para flujo horizontal, similar a la propuesta por Dukler et al. (1964) y la cual es dada por Ec. 3.91. El gradiente de presion por friccién (AP/AH), puede ser estimado: mediante la siguiente ecuacion: Pas Vn ad (3.103) aN va Deo Dep Legal NLR BOOTSOOETTSS100 EN TUBERIAS rs TULO FLO MUL ci donde V, y d representan la velocidad de mezcle y el diémetro de la tuberia, respectivamente. p, representa la densidad de la mezcla sin considerar deslizamiento entre las fases. Eaton et al. propuso la siguiente correlacién para estimar el factor de friccion f, Xx Me M; donde el valor de X es obtenido de la Fig. 3.23, como una funci6n del siguiente 04) término es (one) Md) ae M,\ d AH, } ce donde .M}. define la tasa médsica total y esta dada por la suma de las tasas on de la masicas del gas y liquide, M; y M;, respectivamente. La determina tasa masica esta definida por las Eos. 3.87 a 3.89. 4,, d y x, represertan el rea transversal y el didmetro de la tuberia, y la viscosidad de! gas, respectivamente. El gradiente de presién por energia cinética £, puede ser estimado mediante la siguiente ecuacién: (Mi, av? +M, A¥,) Eo . 2M" OH G.106) TENTHS ‘Bap -Legal N° LFOGITOOTSOOROTSAve =|07), - Ha), @.107) Ave =|), G.108) yt. t H, > (3.109) V, Fou * TH, G.110) Pee eel feet Eaton eta. (1966) Figura 3.23. Correlacidn para Estimar ef Factor de Friccién, Propuesta por Eaton et al (1966). TEIN FDRG ‘Dep Casal LFOGIOOTSONRGTSus is USO MULTIFASICO €* “CAPITULO Eaton et al. propuso el siguiente procedimiento para determinar a. Eaton et al. utiliza cuatro grupos adimensionales, similares a [os propuestos por Duns (1961), !os cuales son Numero de Velocidad Liquido ,, 938%, & _ valido para 0.0697 < N,, $13.246. Numero de Velocidad Gas N» 938 Vig & | valido para 1.5506 < N,, $140.537 o Namero de Dimetro de Tuberia Ny N, =120.8724 & svalido para 20.984 < N, £39.6277 o Numero de Viscosidad Liquida N, Ny, =0.1572 a 4 . Vac donde V y Yq representa la velocidad superficial de las fases gas y liquide, respectivamente, en pie! seg. d representa el diametro de la tuberia, en pie. 1 yo, representa la densidad liquida y la tensin superficial, en Ibn! pies? ¥ dina |em,, respectivamente. TENTED Dep LAAN LAIDOOTEOOOTSC0 EN TURERIAS i FLUO MULTHFASICO FITULO bs. El hold up liquido H, puede ser obtenido de la Fig. 3.24, como una funcién del siguiente termina: a i Mt : “\1463) 0.00226) Guy FACTOR OF ExTMAVEAMIENTO DEL LOuIBG Eaton etal (1966) Figura 3.24. Correlacién para Estimar el Hold Up Liquido, Propuesta por Eaton etal. (1966). 3.3.3. Curvas de Gradiente en Tuberia Inclinada En 1958, Flanigan propuso un procedimiento para estimar el gradiente de presién en tuberias inclinadas. Este gradiente resulta de la suma de los gradientes por friccién y elevacién. De acuerdo a Flanigan, el gradiente de presién por elevacion es independiente del dngulo de elevacion. Por otra parte y sobre la base de resultados experimentales, Flanigan propuso una correlacién Bp Legal N® LFOGLNTSOUDOTS TaN oRORLTEIGTTgeese FLUJO MULTIFASICO EN TUBERIAS 1 camaom para estimar el hold up liquide, como una funcién de la velocidad superficial del gas. El gradiente de presién por friccion se calcula mediante la ecuacién propuesta por Panhandle, a la cual se le incorporo el factor de eficiencia para tomar en cuenta la presencia de liquido en la tuberia. Baker (1960) sugirié (como Io hizo Flanigan) que a la caida de presion por fiiccion, estimada mediante el procedimiento propuesto por 6! para tuberias horizontales, se le adicionara la caida de presién por elevacion, de manera similar al propuesto por Flanigan. En 1971, Cardozo demostré que el método propuesto por Hagedom y Brown (1964) para estimar el hold up liquido era sobre estimado en el caso de pozos direccionales. Sobre la base de informacion de campo obtenida de 26 pozos, Cardozo sugirié un factor de ajuste para estimar el hold up y de esta manera poder extender el método de Hagedom y Brown en pozos con éngulos de desviacién por encima de 50° En 1973, Beggs y Brill publicaron una correlacién para estimar la caida de presi6n que oourre durante el flujo simultaneo de gas y liquide en una tuberia inclinada, Como fluidos experimentales, se consideré aire y agua fluyendo a través de una tuberia acrilica de 90 pie de longitud y I"-1.5" de diametro interno. Un total de 584 pruebas de flujo bifdsico se hicieron a diferentes angulos de inclinacién. Establecieron ecuaciones y generaron un mapa segiin los patrones de flujo. Beggs y Brill demostraron, valiéndose de parémetros adimensionales para el célculo el factor de friccién, que la inclinacién de la tuberia tiene efecto determinante en el factor de entrampamiento bifasico (hold up). Su trabajo se considera uno de los mas completes y su corretacion represent ser la primera eouacién en su tipo que permite predecir el gradiente de presién a cualquier grado de inclinacién de la tuberia En 1985, Mukherjee y Brill propusieron un nuevo método para predecir ef gradiente de presién en tuberias inclinadas, en un intento de superar algunas limitaciones presentes en el método inicialmente presentado por Beggs y Brill (1973), asi como también, tomar ventaja de una nueva instrumentacion TONED ATID ep. Legal NP LFDGIDOOTSOOROTSa FLuvo MULTIPASICO EN" —— ‘CAPITULO HI desarrollada para medir el hold up liquido. Utilizando aire y kerosén como fluidos experimentales, Beggs y Brill obtuvieron mediciones de gradiente de presion y hold up liquido, en un orden superior a los 1000 y 1500, respectivamente, y para un amplio rango de tasas de gas y liquido. Asheim (1986) formulé e implementé un nuevo modelo, mediante un programa computacional llamado MONA. El modelo tiene dos caracteristicas resaltantes: La primera, involucra el calculo del hold up y el factor de friccién descritos por tres parametros independientes, los cuales son relacionados a un fenémeno hidrodinamico y puede ser estimada a priori para una condicién de flujo dada; la segunda caracteristica, involucra un ajuste del modelo con data de campo, lo que permite ajustar los tres parémetros independientes, minimizando los errores de computacién. A: Correlacién de Beggs y Brill (1973) Se clasifica como una correlacién del “Tipo c”. Beggs y Brill investigaron el flujo multifésico en tuberias inclinadas, a fin de estimar el efecto que el angulo de inclinacién de la tuberia tiene sobre el hold up liquido y el gradiente de presin. Los fluidos utilizados en este estudio fueron aire y agua Beggs y Brill propusieron la siguiente ecuacién para predecir el gradiente de presion en tuberfas inclinadas: Lm Vn + pg seno 0 —E, G12) La predicoién del gr-diente de presién mediante Ec. 3.112, requiere inicialmente determinar el patron de flujo, para posteriormente predecir el hold up liquido. Prediccién del patron de flujo Beggs y Brill disefiaron un mapa empirico para predecir la ocurrencia de patrones de flujo en tuberia horizontal. Posteriormente, este mapa fue NTO ITE ‘Dep Legal NLFVGIDO0TSOODOTS__ FLUJO MULTIFASICO EN TUBERIAS w CAPITULO modificado ligeramente para incluir zonas de transicién entre los patrones de flujo segregado e intermitente. La Fig. 3.25 muestra el trabajo original y modificado, propuesto por Beggs y Brill a ‘rr Lagu CONTENT = 4, ‘Beggs y Brill (1973) Figura 3.25. Mapa de Patrones e Flujo en Tuberia Horizontal, Propuesto por Beggs y Bril (1973). Estos investigadores consideraron ademas correlacionar los limites de transicién de los patrones de flujo, como una funcién del hold up liquide sin deslizamiento entre fases 4, y el numero adimensional de Froude Nj» para la mezcla y definido mediante la siguiente ecuacion: 3.113) ap Legal NL FOGTROOTSONROTSFLUJO MULTIFASICO EN TUBERIAS CAPITULO UT donde ¥, y d representan la velocidad de mezcia y el diametro de Ia tuberia, respectivamente. Los limites de transicién para cada patron de flujo son estimados mediante las siguientes correlaciones L=3162, G.ll4) L, =0,0009254,°**, .115) L,=0.104)'", (3.116) 1,505 4,5 Gun Criterio utilzado para determinar el patron de flujo que existiria si la tuberia estuviese bajo una condicién horizontal. Patron de flujo Segregado Si 4, < 0.01 y Nex La (3.124) De acuerdo a Beggs y Brill, la prediccién del hold up liquide prevé iniciaimente, la determinacién de éste como si la tuberia estuviese bajo la condici6n horizontal, solamente, Posteriormente, éste valor es corregido por el Angulo de inclinacién de la tuberia. La Fig. 3.26 muestra a variacién del hold up liquido por efecto del angulo de inciinacién, para tres de sus pruebas. Bajo condiciones de flujo horizontal, el hold up liquide H, puede ser determinado mediante la siguiente ecuacion: G.125) La Ec, 3.126 tiene la restriccion que H,24. Los coeficientes a, b y ¢ dependerén del patrén de flujo existente en flujo horizontal, y los mismos se encuentran dados en Tabla 3.1. [SON 97ST ASTD ep Legal LPT OUT SOURIS7 FLUO MULTIFASICO EN TUBERIAS, CAPITULO It Tabla 3.1. Coeficientes para Estimar el Hold Up Liquido en Tuberia Horizontal. Beggs y Brill (1973). PATRON DE FLUJO poe CeNiES a b ‘ SEGREGADO 0.9800 0.4846 0.0868 INTERMITENTE 0.8450 0.5351 0.0173 DISTRIBUIDO 1.0650 0.5824 0.0609 BO 6 Ned64 40 + as.196 4 6 © 45.010 @ sof a % 20 2 2 4 00. =90 70-50-30 10 10 30 50 70 90 ANGLE OF PIPE FROM HORIZONTAL Beggs y Bril (1973) Figura 3.26, Efecto del Angulo de Inclinacién Sobre el Hold Up Liquido. El hold up liquido corregido por el efecto del Angulo de inclinacion My,., es determinado mediante la siguiente ecuacién: Hyg =H, (6.126) TaN TR DITA ‘ep. Legal WF LGTE0OTSOOROTS_ Fino MULTIFASICO EN TUBERIAS 1s EAPITULOU donde el factor para corregir M,.9, €8 dado por y=14C |Seno(1.8 0) -0.333(Seno(180))'] @.127) @ representa el Angulo de inclinacién de la tuberia y es medido desde la horizontal. La constante C es definida como C=C-A)Ln(e A Ni, New), 6.128) con la restriccion que C20. Los coeficientes ¢, f, ¢ y ft se encuentran disponibles en Tabla 3.2, como una funci6n del patron de flujo. ‘Tabla 3.2. Cooficientes para Estimar la Constante C . Beggs y Brill (1973). PATRON DE FLUO COEFICIENTES pseeneee ros easuaridaeeedenaa et ‘SEGREGADO (Fuso ASCENDENTE) oti 37680 35380-16140 INTERMITENTE (FLUO ASCENDENTE) 2.9000 03080 0.4873 00878 DisTRBUIDO NO HAY CORRECCION (FLUJO ASCENDENTE) coe ‘TODOS LOS PATRONES (Puno DESCENDENTE) 47000-05602 07244-05086 Nu, representa el ntimero de velocidad liquida y se encuentra definide como: (3.128) TEN TOTES ‘Da Legal WL FOGTD007SOTRN7SFLUIO MULTIFASICO EN TUBERIAS CAMTULO donde p,, ¥, ¥ &; representan la densidad, la velocidad superficial y la tension superficial de la fase liquida, respectivamente. Cuando el patron de flujo se ubica en la regién de transicion, el valor de H,», debe ser interpolado entre los valores correspondientes los patrones de flujo segregado @ intermitente, Hun uw ¥ Huordom Fespectivamente, utilzando la siguiente ecuacion (Hye, rane = A CH 0, see + A= A) Hoy aes (3.129) N, aoe 6.130) El factor de friccién /' es determinado mediante la siguiente ecuaci6n: (4), G.131) \ donde el factor de friccisn normalizado /, se determina mediante el diagrama de Moody (Fig. 3.11) para tuberias lisas, v como una funcién del_ ndmero de Reynolds, definido como: (3.132) Sobre la base de resultados experimentales, Beggs y Brill desarrollaron la siguiente correlacién para estimar la razon f/f, ONES ISTT ‘Bap Legal LRDSTRRTSDODNTSFLUJO MULTIFASICO EN TUBERIAS. 18 CAPITULO (3.133) ee EW REeee oo rains ay THs 3.182 LNG) -08725(Ln()) +0.01853LnOD” (3.135) (Arua La Ec. 3.134 presenta ciertas disconfinuidades para valores de y alrededor de 2.6310“ y 1.016, Cuando los valores de y se encuentren entre 1-12, Beggs y Brill propusieron la siguiente correlacién: s=Ln(2.2y-12). 6.136), Payne et al. (1979) encontré que la correlacién de Beggs y Brill sobre predice los factores de friccién y el hold up liquide #7... Payne et al recomendé utilizar los siguientes factores de correccion para mejorar los valores de Hi Para 0>0 Hyg) = 9.924 A sco) sesesy it~ (3.137) Para 6<0 Figg) = 0-685 H 0))beesy ait (3.138) TaN 7S OTT Bip CAINS LRUGTICOTSOODOTSFLUJO MULTIFASICO EN TUBERIAS. ~ ‘CAPITULO MT Se debe tener siempre presente que los valores de H,,,,, para 0 >0, no deberan exceder 4,. Finalmente, £, representa el término de energia cinética y es definide como: (3.139) donde P representa la presién del segmento y puede ser estimada como el promedio aritmético entre P, y P., las cuales representan la presion a la entrada y salida del segmento, respectivamente. TSUN STEDHOTEIOTS Dep Legal LPOGTSDUTSDOROTSCAPITULO IV ANALISIS NODAL PARA POZOS DE PETROLEO ‘élisis Nodal” representa una metodologia de diagnéstico, la cual ha sido aplicada por afos en diversas areas. Su principal objetivo se basa sobre el analisis de sistemas complejos, tales como: circuitos eléctricos, redes de lineas de fiujo, sistemas de bombeo centrifugo, entre otros. Ademés, esta técnica permite no solo definir los componentes o secciones principales que conforman estos sistemas, sino también la interrelacién que éstos poseen entre si. La aplicacion de este método a sistemas productores de hidrocarburos fue inicialmente propuesta por Gilbert en 1954, y discutida posteriormente por Nind (1964) y Brown (1977). El procedimiento consiste en elegir un nado solucion para dividir el sistema en este punto y poder analizar la interrelacion entre ambos cumponentes. En resumen, el andlisis nodal es un método muy flexible que puede ser utilizado para mejorar la eficiencia de produccién de un pozo de petréleo, En las proximas secciones, se explicara en detalle e! uso o aplicacion de esta metodologia no solo en pozos de petrdleo, sino también en pozos de gas. 4A Nodo Solucion Representa un punto ubicado en cualquier parte de un sistema de producci6n pozo-yacimiento, tal como se muestra en Fig. 4.1, y la ubicacion de éste dependerd del componente que se desee aislar para su evaluacion. TEN 9 DROTE I To Dep egal LFOGTROOTEDDROTSTanque Nodo Ubicacion Observacicn Separador Reductor Funcional Cabezal Val. Seguridad Funcional Restriocion Funcional Fondo del Pozo Perforacionos Yacimiento Figura 4.1, Ubicacién del Nodo en un Sistema de Produecién Pozo-Yacimiento, Todos los componentes aguas arriba del nodo solucién se encuentran representados por la curva de oferta, mientras que todos los componentes aguas abajo de éste se encuentran representados por la curva de demanda, A través de un nodo solucién, se debe cumplir ademas: 1) La tasa de flujo que fuye hacia el nodo solucion deberé ser la misma que sale de éste; 2) En este tipo de nodo existe evidentemente una y solo una presion. Es importante mencionar que hay dos puntos en el sistema de produccion pozo-yacimiento donde la presion no depende de la tasa de flujo, Una de estas presiones se encuentra representada por la presién promedio del yacimiento 7, La otra es la presién de salida del sistema, generalmente representada por la presién det separador P,,, aunque si el poz0 se encuentra controlado por un TaN FORTE ISTE ‘ep Taal N° LALGIEOOTSOOIOTSANALISIS NODAL PARA FOZ0S DE PETROLEO. 1 ‘CAPITULOIV estrangulador de flujo en superficie entonces la presién de salida podria ser la presién de cabezal del pozo P,, Una vez que el nodo solucién haya sido seleccionado, la presién en el mismo es determinada partiendo desde los puntos fijos de presién ubicados en ambos extremos, es decir, la presién de yacimiento y la presién de separacién o cabezal, P. y P.,, 0 P,,, respectivamente. En otras palabras, AP, ¥ Pa, representan todas las posibles pérdidas de presion que ocurren aguas arriba y aguas abajo, respectivamente, del nodo solucién. La representacién grafica de la P,,,. vs. tasa de flujo g definiré dos curvas, conovidas como las curvas de oferta y demanda, y cuya intercepcion establecera un balance entre la energia del fluido que aporta el yacimiento y la energia requerida por el fiuido que demanda la instalacion. La Fig. 4.2 presenta la representaci6n tipica de P..., VS. ¢ De manera generalizada, se recomienda seguir el siguiente procedimiento para aplicar la técnica de analisis nodal en pozos de petréleo 0 gas: 44 De acuerdo al sistema considerado, defina el componente que se desea aislar para su evaluacién. En ciertos casos, algunos cambios pueden resultar limitados, como por ejemplo: tamafio de hoyo perforado, tamafio de revestidor y tuberfa de produccién, entre otros 2. Seleccione el elemento que desea optimizar. 3. Seleccione Ia ubicacién del nodo que mejor representaria el efecto del cambio del elemento seleccionado. Este punto puede no ser critico ya que el mismo resultado se obtendra sin importar la ubicacién del nodo. SON TTRSROAT ASTD Dep, Lega LROGTROOTSOORITSDEPETROLO. ANALISIS NODAL PAI TAPITULOW Establezca las expresiones matematicas requeridas. Obtenga Ja informacion minima necesaria para calcular la caida de presién que ocurre a través de cada uno de los componentes, como una funcién de la tasa de flujo. Algunas veces, esta informacion puede resultar insuficiente y por lo tanto se debera simular algunas condiciones dentro de un rango estimado. Determine el efecto de cambiar algunos elementos del sistema seleccionado, mediante el calculo de la tasa de flujo que resulta de la intercepcién de las curvas de oferta y demanda generadas en cada condi Repita el mismo procedimiento para cada elemento que requiere ser optimizado. Establezca conclusiones y tome decisiones, siempre y cuando éstas se encuentren soportadas con un analisis econémico. Preston en el Nodo Poa Tasa de Filo 4 ‘Figura 4.2. Representaci6n Tipica dela Curva de Presién de un Nodo vs. Tasa de Flujo. Dag Lega FLAG SOODOTSANALISIS NODAL, PARA FOZ0S DE PETROLEO. st ‘CAPTULOTY 4.4.1 Procedimiento Solucién A-. Solucién en el Fondo del Pozo Probablemente, uno de los puntos del sistema de produccién donde preferencialmente se ubica un nodo solucién corresponde al fondo del pozo y es que el principal elemento a estudiar se encuentra representado por la tuberia 0 sarta de produccién, donde ocurre aproximadamente el 80% de las pérdidas totales de presién. Por otra parte y muy especialmente en pozos nuevos, resulta sumamente importante elegir el tamafio de tuberia adecuado para las condiciones de produccién de un pozo. Una tuberia muy pequefia restringiria la tasa de produccién debido a las excesivas pérdidas por fricci6n, mientras que una tuberia muy grande podria generar excesivas pérdidas por resbalamiento de liquido. Para considerar la solucién en el fondo del pozo, considérese un nodo ubicado en el punto 6 de la Fig. 4.1. En este caso, la presion en el nodo P,,4, S€ encuentra representada por P, y se asume que el pozo no se encuentra revestido y cafioneado (Pig =P.,). Ademés, el nodo divide el sistema pozo- yacimiento en dos componentes: el componente yacimiento y el componente sistema de tuberias, Se asume que el componente yacimiento no estara sometide a algin cambio y en consecuencia la curva de oferta deberé permanecer invariable. Sin embargo, el componente sistema de tuberias no solo se ha aislado, sino también algunos de los elementos que lo conforman podran ser cambiados y el efecto de este cambio podré ser analizado, recalculando la presién ?,,. Como se muestra en Fig. 4.3, la curva de oferta IPR podra ser determinada partiendo desde la presién de yacimiento P, y restandole todas las pérdidas ocasionadas por el flujo de fluidos a través del medio poroso. Estas pérdidas pueden ser estimadas mediante la ecuacién de flujo de fluides bajo una condicion de semi-estado estable o por cualquier otro método empirico, disponible en Ia literatura. La representacién de esta curva es posible mediante tal como se muestra en Fig. 4.4. un grafico de Pg, VS. g, donde Pag, = Pry ISON STE-98OAE IST Dep Lesa LFOGTROOTSONROTS__awALisis NODAL PARA N0708 DE RETROLEO a eno Pa Figura 43, Solucin en el Fondo det Pozo. ia IPR Figura 44. Curva Tipica de Comportamiento de AfIuenci TNR SANAE IT Dap Legal ROGETANALISIS NODAL PARA POZ0S DE PETROLEO CAPITULO TY 5 La curva de demanda se construye partiendo desde la presion de separacién ?,, tal como se muestra en Fig. 4.3, y adicionandole a la misma, todas y cada una de las pérdidas de presién que puedan ocurrir como consecuencia del flujo de fluidos a través del sistema de tuberias, hasta llegar al fondo de! pozo. Para la linea de flujo, las pérdidas de presién podran ser estimada mediante la seleccin de una adecuada correlacién 0 modelo mecanistico, disponible en la literatura, para estimar el gradiente de presin dinamico a través de tuberias. También, podra hacerse uso de alguna curva de gradiente dinamico, como la mostrada en Figs. 4.5 y 4.6. > & Pay Presion > " 3 4 R Ae, St % Longitud Figura 45, Curva Tipica de Gradiente Dindmico, Valida para Flujo en Tuberia Horizontal, En el caso de utilizar alguna curva de gradiente disponible en la literatura, el procedimiento a seguir es el siguiente: PNET ‘Dep Lagal WEFT BOOT SOORDTSANALISIS NODAL PAY tes ere arene Elja la curva de gradiente en tuberia horizontal y vertical que posea similar caracteristicas en cuanto a propiedades de los fluidos y variables geométricas, se refiere. Cada curva de este tipo se encuentra disponible para una serie de tasas de flujo. Seleccione una tasa de flujo y su correspondiente curva de gradiente, la cual y de acuerdo a Fig. 4.3, iniciaimente corresponderé para una tuberia horizontal Comience con la presién de separacién 7, (variable conocida) y descienda verticalmente sobre la curva hasta interceptar la curva de gradiente dindmico, correspondiente a la RGL tipica 0 caracteristica de! poze, tal como puede apreciarse en Fig. 4.5. Horizontalmente, intercepte el eje referido a longitud y descienda verticalmente una distancia equivalente 2 la existente entre el separador y el cabezal del pozo. Nuevamente, intercepte la curva de gradiente y ascienda verticalmente para determinar la presién de cabezal /, Manteniendo la misma tasa de fiujo, previamente elegida, seleccione la curva de gradiente dinémico pero para tuberia vertical, como la mostrada en Fig. 4.6. Con el valor de la presion de cabezal P., obtenida en paso 5, descienda verticalmente sobre la curva de gradiente e intercepte la curva de gradiente dinamico, correspondiente a la RGL tipica o caracteristica del ozo. Horizontalmente, intercepte el eje referido a profundidad y descienda verticalmente una distancia equivalente 2 la existente entre el cabezal y el fondo del pozo. Nuevamente, intercepte la curva de gradiente y ascienda verticalmente para determinar la presion de fondo fluyente P,, ep Legal LPT SOOT_ANALISIS NODAL PARA FOZOS DE PETROL ves CARTULOW 40. Grafique cada valor de 7, y su correspondiente 4- tal como se muestra en Fig. 47 11. Seleccione un nuevo valor de tasa de flujo y repta el paso 2 hasta 9 p, Presion Profundided Figura 4.6. Cana Tipies de Gradiente Dinmico, Valid pars Fluo en ‘Tuberia Vertical, La curva de demanda sera definida por ta union de cada uno de estos puntos, tal como puede apreciarse en Fig. 4.8 La intercepcién de ambas curvas define una tasa de equilibrio, la cual no representa la maxima, la minima, ni siquiera la éplima. Es simplemente la tasa de flujo a la cual éste pozo producira para el sistema de tuberia instalado, La ‘asa de flujo cambiaré solo si existe el cambio de algin elemento del sistema de produccién y en este caso posiblemente estarén referidos a cambios en él diémetro de la tuberia vertical u horizontal, estrangulador, presién del separador, entre otros, aN gTHTROTISTT aps Lagal N° LRDTEOOTSOOOOTSNODAL PARA POZOS DE PETROLEO i ‘CAPITULO TY curva de oferta u % Figura 4.8. Curvas de Oferta y Demand. Solucién en el Fondo del Pozo. TN SDE TTT pagal NLPOGT BOOTIESANALISIS NODAL PARA FOZOS DE PETROLEO. 16 CARTULOW Si se desea conocer el efecto de cambiar el didmetro de la tuberia de produccién por ejemplo, se repite el provedimiento anteriormente mencionado con la particularidad que las curvas de gradiente deberan ser ajustadas al nuevo valor det diémetro. La Fig. 4.9 presenta el efecto del cambio del diametro de la tuberia de producci6n. dy Blbmete Tusa de Prosuién Figura 49, Efeto del Cambio del Didmetro de la Tuberia de Produce, Solcién en Fondo del Pozo, La ubicacién de! nodo en el fondo del pozo puede también resultar estratégica en aquellos casos donde se requiera evaluar el efecto del cambio de presion del yacimiento 7 con tiempo. Asumiendo una relacion gas-petroleo ‘RGP constante, cualquier cambio de P, puede reflejarse perfectamente en a curva de oferta y por lo tanto su efecto sobre la tasa de produccién puede estimarse a partir de un grafico de P,, vs. q, como el mostrado en Fig. 4.10. Generaimente, cambios de P, conlleva necesariamente a cambios de I RGP y, RON TRTODITTT ep Legal N° LF OGTRDUTROAT