Schlumberger
Sistema de Producción de Petróleo y Gas
Diseño y Manejo con
PIPESIM
Resumen de Historia de PIPESIM
1984 - PIPESIM desarrollado en Plataforma Unix
1985 - PIPESIM cambiado a DOS y formado en Baker Jardine
1990 - Desarrollo de PIPESIM GOAL (Gas Lift Optimization & Allocation)
1993 – Interface Gráfica de Usuario (IGU) en Windows
1994 – Lanzamiento de PIPESIM Net
1996 – Lanzamiento de PIPESIM FPT
1997 – PIPESIM FPT conectado a ECLIPSE
2000 – Desarrollo de PIPESIM 2000 (Nueva IGU de 32-bit)
2001 – Schlumberger compra Baker Jardine
2003 – Lanzamiento de PIPESIM 2003
2006 – Lanzamiento de PIPESIM 2006 y Avocet IAM
Clientes Principales - Operadores
ExxonMobil British Gas
BP Petronas
Shell SUCO
Pemex PDVSA
ChevronTexaco UNOCAL
ConocoPhillips ADCO/ADMA
Saudi Aramco BHP
Amerada Hess JNOC
Clientes Principales - Contratistas
AMEC Halliburton
Bechtel Intec Engineering
Brown & Root Kvaerner
CSO Aker McDermott
Fluor Daniel MW Kellogg
Foster Wheeler Paragon
Genesis Snamprogetti
Granherne Worley
PIPESIM: Integración total del Sistema
Network Analysis
• Modelaje de Redes de Producción
Pipeline & Facilities Analysis
• Estudios detallados de líneas de producción
Field Planning
• Integración Yacimiento-
Facilidades de Superficie
Well Performance Analysis
• Modelo de Influjo, Análisis Nodal
Production Optimization
• Optimización del sistema y
manejo de limitaciones
Multilateral Well Simulator
• Modelos de pozos Horizontales y Multilaterales
PIPESIM: Integración total del Sistema
Network Analysis
• Modelaje de Redes de Producción
Pipeline & Facilities Analysis
• Estudios detallados de líneas de producción
Field Planning
• Integración Yacimiento-
Facilidades de Superficie
Well Performance Analysis
• Modelo de Influjo, Análisis Nodal
Production Optimization
• Optimización del sistema y
manejo de limitaciones
Multilateral Well Simulator
• Modelos de pozos Horizontales y Multilaterales
Sistema Total de Compresor
Producción Separador
Riser Bomba
Estrangulador gas
Línea de Flujo
petróle
o
Válvula de Líneas de
Seguridad Exportación
Tubería
Yacimiento
Completación
P8 = Psep – PCD
Pérdidas de Presión en el
Sistema
P5 = Ptf – PDSC P7 = PRB – Psep
P9 = Psep – PPD
P6 = PDSC – PRB
P10 = PCD –
P3 = PUSV – PDSV
PGD
P11 = PPD – PLD
P4 = Pwf – Ptf
P1 = PR – Pwfs
P2 = Pwfs – Pwf
Relaciones de Comportamiento de Afluencia
Yacimientos de Petróleo: Gas and Gas Condensado
IP del Pozo Reservoirs:
Vogel IP del Pozo
Fetkovich Back Pressure
Jones Jones
Estado Pseudo-Estable Estado Pseudo-Estable
Índice de Productividad (IP)
Pwf > Pb :
Q = PI x (Pws - Pwf)
Pwf < Pb :
Para pozos de Petróleo:
Q = PI x (Pws - Pb) + (PI / 2Pb)(Pb2 - Pwf2)
Para yacimientos de gas compresible:
Q = PI x (Pws2 - Pwf2)
donde, Pws = presión estática del yacimiento
Pwf = presión de fondo fluyente
Pb = presión de burbujeo, Q = tasa de flujo
Ecuación de Vogel
Relación empírica para fluidos por debajo de la presión de
burbujeo:
q/qmax = 1 - (1 - C)(Pwf/Pws) - C(Pwf/Pws)2
donde, C = Coeficiente IP, valor normal es 0.8
qmax = Potencial absoluto a hoyo abierto
Pws = Presión estática del yacimiento
Pwf = Presión de fondo fluyente
Ecuación de Fetkovich
Ecuación alternativa a Vogel
Correlación Empírica
q / qmax = [ 1 - ( Pwf / Pr )2 ] n
Mientras menor sea el valor de n, mayor será el grado de turbulencia
También conocida como “Normalised Backpressure Equation”
Ecuación de Jones
Yacimientos de Gas y Petróleo saturados
Ecuaciones:
Gas: (P2) = AQ + BQ2
Oil: (P) = AQ + BQ2
donde
A : Coeficiente de flujo laminar (Darcy)
B : Coeficiente de flujo turbulento (Non Darcy)
También conocida como “Ecuación de Forcheimer”
Back Pressure
Para Pozos de Gas
Q = C (Pws2 - Pwf2)n
Ecuación empírica de Schellhardt & Rawlins
Normalmente, 0.5 < n < 1.0
Ecuación de Estado Pseudo-Estable
Yacimientos de Petróleo y Gas
Ecuación de Darcy Modelo radial de yacimiento
Parámetros usados en las ecuaciones:
Permeabilidad
Espesor
Radio (drenaje externo del yacimiento)
Daño (factor de daño adimensional)
Diámetro del hoyo
Pozo de Gas: flujo laminar y turbulento
Pozo de Petróleo: flujo laminar
Opciones de Completación
SOLO válido cuando se utiliza la ecuación de Estado Pseudo-
Estable como modelo de comportamiento de afluencia
Para calcular el factor de daño y el coeficiente de turbulencia
(para pozos de gas)
Opciones de completación:
Ninguna (es decir, no daño que cause resistencia a influjo)
Hoyo Abierto (pozo no cementado o revestido)
Cañoneado (Modelo McLeod)
Empacado con Grava (Modelo Jones)
Modelos de Completación Horizontal
IP distribuido (conductividad finita):
IP Distributivo: IP por unidad de longitud
Estado Estable (Joshi)
Estado Pseudo Estable (Babu & Odeh)
IP Puntual / Single Point PI (conductividad infinita):
Estado Estable (Joshi)
Estado Pseudo-Estable (Babu & Odeh)
Flujo Multifásico
Descripción de Orkiszewski (Flujo Vertical):
Patrones de flujo
Bubble Flow
Slug Flow
Transition – Churn Flow
Annular – Mist Flow
Mapa de régimen de Flujo Horizontal - Taitel
Dukler
Los patrones de flujo son determinados con mapas de velocidad
superficial del líquido vs. velocidad superficial del gas
75.0
Bubbly
10.0
Intermittent
UsL 1.00 Annular
(ft/s)
0.10
Stratified Stratified
Smooth Wavy
0.01
0.1 1.0 10.0 100.0 900.0
UsG
(ft/s)
Factor de colgamiento (Hold-up):
AL AG
HL
APIPE AL
Líquido:
Hold-up = 0 Todo el gas fluye
Hold-up = 1 Todo el líquido fluye
Hold-Up Gas: 1-Hold Up Líquido
No-Slip Liquid Holdup (No deslizante):
Colgamiento de líquido, si el gas y líquido estuviesen viajando a la misma
velocidad.
qL
L
qL q g
Flujos deslizante y no deslizante (slip & no slip)
VL=VG
VL<VG
Liquid Holdup, HL
Función del grado de deslizamiento de la fase de gas
Q g = V gA g Vg Ag
AG
AL
QL = VLAL VL AL
Efectos de deslizamiento de fluidos
Ascendente Descendente
vG vG
vL qG = vG AG vL
qL = vL AL
v L < vG v L > vG
Por qué es importante el colgamiento del líquido?
Se requiere para:
Estimación precisa de las pérdidas de presión entre dos fases
Pronóstico de formación de hidratos y deposición de parafinas
Predicción del volumen de tapón presente
Método de Cálculo de Cambio de Presión
Determinar las fases presentes
Determinar el ángulo de inclinación
Determinar el patrón de flujo
Calcular el colgamiento de líquido
Calcular las pérdidas o aumentos de presión por elevación,
fricción y aceleración
Fases Presentes
.00001 < fracción de vol. de Liq. < .99
Gas Multifásico Líquido
Monofásico Monofásico
Ángulo de Inclinación
Si el ángulo de inclinación > 45° ó < -45° entonces se
aplican los patrones de flujo vertical y correlaciones de
cambios de presión
En caso contrario aplican los patrones de flujo
horizontal y correlaciones de cambio de presión
Cambios de Presión – Flujo Monofásico
(dp/dl) = (dp/dl)elev.+ (dp/dl)fric.+ (dp/dl) acc.
= -gsin (por elevación)
+ -2/2d (por fricción)
+ - x d/dl (por aceleración)
Donde,
= factor de fricción
= densidad del fluido
= velocidad del fluido
g = constante gravitacional a la elevación actual
= ángulo de flujo
d = diámetro de tubería
Correlaciones de Flujo Monofásico
Disponibles
Moody (predeterminada)
AGA – Ecuación de Gas Seco
Panhandle A
Panhandle B
Hazen-Williams
Weymouth
Correlaciones de Flujo Multifásico
Correlaciones estándar publicadas en la industria:
Duns & Ros
Orkiszewski
Hagedorn & Brown
Beggs & Brill (original & revised)
Mukherjee & Brill
Govier, Aziz & Fogarasi
AGA & Flanigan
Oliemans
Gray
Noslip
Correlaciones de Flujo Multifásico (cont.)
Correlaciones desarrolladas por BJA
BJA para líneas de condensado
Correlaciones mecanísticas
Ansari
Xiao
OLGA-S
Correlaciones Recomendadas:
(Hallazgos de Baker Jardine JIP, 1995)
Sistemas Monofásicos - Moody
Pozos de Petróleo - Hagedorn & Brown
Pozos de Gas/Condensado - Hagedorn & Brown
Líneas de Flujo (pipelines) de Petróleo - Oliemans
Líneas de Flujo de Gas/Condensado – Correlación BJA
Proyectos de Flujo de Fluidos de la
Universidad de Tulsa
Líder en investigación de flujo
multifásico por más de 30 años
Schlumberger volvió a asociarse en
Agosto 2002
Recientemente ha desarrollado el Modelo
Mecanístico Unificado, válido para todos
los ángulos de inclinación
Modelado de Fluidos - Alternativas
Método Black oil:
Utiliza las correlaciones publicadas de Black Oil para
determinar las propiedades físicas de los fluidos
Método Composicional:
Ecuación de estado, modelando componentes puros y
fracciones pesadas de hidrocarburos.
Modelo Black Oil
Modelo Black oil
Válido para fluidos con relaciones gas petróleo (RGP)
menores a 1500pcn/bn
Se utilizan correlaciones empíricas para determinar las
propiedades derivadas (gas en solución -Rs- y factor
volumétrico de formación -Bo-)
Se utilizan las propiedades derivadas para calcular las
propiedades físicas de los fluidos, como densidad y
viscosidad del petróleo
Datos de entrada básicos para Black Oil
1. Corte de agua (o RGA para sistemas de gas)
2. Relación Gas Petróleo en el tanque (o RPG para
sistemas de gas)
3. Gravedad específica en el tanque
4. Gravedad específica del agua en el tanque
5. Gravedad API (o densidad del fluido muerto)
Datos opcionales de calibración
Gas en solución medido a presión y temperatura de burbujeo
Viscosidad del petróleo muerto medida a una o dos
temperaturas
Factor volumétrico del petróleo a presión y temperatura por
encima de la presión de burbujeo
Valores de factor volumétrico del petróleo, viscosidad del
petróleo vivo, viscosidad del gas, y/o compresibilidad del gas
medidos a condiciones de presión y temperatura igual o
menor que el punto de burbujeo
Viscosidad de mezclas petróleo/agua
Woelfin
Mixture
emulsion
viscosity
correlation
µ (oil)
Volume ratio µ (water)
Inversion
method
method
(default)
0% watercut 100%
Modelo Composicional
Modelo Composicional
Usa “Multiflash”
Modelado con Ecuación de Estado
Se utiliza para sistemas de crudos volátiles o
condensado
Incluye hidratos
Envolventes de fase
Procedimientos de Cálculo
Datos de entrada:
selección/especificación de componentes
fracciones de los componentes (en %Mol)
Fracciones de Petróleo
medidas / propiedades críticas
Tablas PVT generadas M L , hL , L , L ,
Parámetros calculados: M G , hG , G , G
Fracciones de Petróleo
Propiedades medidas:
Punto de ebullición (BP)
Peso Molecular (PM)
Gravedad específica (Gesp) o densidad
[Nota: Se requiere al menos dos propiedades físicas]
Propiedades críticas (información adicional):
Temperatura crítica (Tc)
Presión crítica (Pc)
Factor acéntrico (Omega)
Ecuación de Estado (EOS)
Tres opciones:
Soave-Redlich-Kwong (predeterminada)
p = (NRT/(V-b)) + (a/(V2+2bV-b2))
Peng-Robinson
p = (NRT/(V-b)) + (a/(V(V+b)))
BWRS Benedict-Webb-Rubin-Starling
Nota: En las ecuaciones anteriores, “a” y “b” son derivadas como
funciones de temperatura crítica, presiones y factores acéntricos de los
componentes puros
Operaciones con una Tubería
Operaciones con una tubería
Análisis del Sistema
Perfíl Presión/Temperatura
Ajuste de la correlación de flujo multifásico
Análisis Nodal
Optima longitud de un pozo horizontal
Tablas hidráulicas
Gasto (flujo) de BN-LAG vs Presión de Tub. de Rev.
Comportamiento del Sistema Artificial
Ajuste de correlación de Flujo Multifásico
Para determinar la correlación de flujo que mejor se
ajusta
Seleccionar las correlaciones de flujo requeridas
Introducir presiones y temperaturas conocidas a
diferentes profundidades medidas (MD)
Introducir condiciones de borde
Los resultados muestran cada correlación y los datos
introducidos
Ajuste de correlaciones de flujo
Ansari
Hagedorn Brown
Orkezewski
Beggs Brill Measured
Perfíl Presión/Temperatura
Calcula el perfil de presión y temperatura para un
sistema y también varía algunos otros parámetros
dentro del sistema
Se introducen las variables de sensibilidad
Se introducen las condiciones de borde
Como resultado se obtiene la gráfica de presión y
temperatura contra profundidad (del pozo) o elevación
(de tubería)
Perfíl Presión/Temperatura
Corrida de
Simulación Sencilla
Análisis Nodal
Análisis nodal clásico en cualquier punto
Dividir el sistema en dos y calcular el flujo que entra y
sale alrededor de ese punto
El gráfico resultante muestra las curvas clásicas de
flujo que entra y sale ( Demanda / Oferta )
Análisis Nodal
Psep
Nodo Tomado en el fondo del Hoyo
Modelado hasta el separador
Capacidad de fluir(Oferta) = PR –
Pwf
Capacidad de transportar (Demanda) = Pwf – Psep
Pwf PR
Análisis Nodal
Psep
PR Inflow
Outflow
Pwf
Pwf PR Psep
Flow rate
Análisis Nodal
Tubing ID =
3”
Tubing ID = 3 1/2”
Tubing ID =
4”
Skin = -2
Skin = 0
Skin = 2
Tablas Hidráulicas
Generar una tabla de presión de fondo que pueda ser
utilizada por el simulador de yacimientos
Interfase a simuladores de yacimiento comunes
como:
ECLIPSE
VIP
PORES
COMP4
MoRes
Tablas Hidráulicas
Comportamiento de Levantamiento Artificial
Permite que las curvas de comportamiento de
levantamiento artificial (gas / BES ) sean generadas y
también varía algunos otros parámetros dentro del sistema
Generar curvas de comportamiento de entrada para GOAL
El gráfico resultante es cantidad de gas lift (o poder de la
BES) versus tasa de producción de petróleo
Comportamiento de Levantamiento Artificial
Psep = 75 psig
125 psig
150 psig
175 psig
Comportamiento de Pozo Horizontal
Un pozo horizontal con múltiples fuentes a lo largo del
pozo, el cual calcula el drenaje del yacimiento y las caídas
de presión a lo largo del pozo
Predice en forma precisa el comportamiento hidráulico del
pozo en la completación
Evaluación de la productividad mediante la opción
“Longitud optima de completación horizontal”
Optimum Horizontal Well Length
Specify:
- reservoir pressure
- outlet pressure
- completion length
Diseño de Levantamiento
Artificial
Sistema de Levantamiento Artificial
Gas Lift
*Un punto de inyección
*Puntos de inyección múltiples (GLV’s)
BES (Bombas Electro Sumergibles)
Rod Pump Bombeo Mecánico
* Diagnóstico
* Diseño
Gas Lift Design
• Espaciamiento de Mandriles
• Selección de Válvulas
• Simulación de diseños existentes
Válvulas controladas por presión de Revestimiento y tubería
(válvulas IPO / PPO)
Espaciamiento de válvulas, cálculos de presiones de
apertura /Cierre y tamaño de válvulas
• Opciones de múltiples gradientes estáticos
Diseño de Gas Lift
Dos etapas distintas:
(1) Cálculo del punto de inyección más profundo
(2) Determinación del espaciamiento de Mandriles
Simulación de Gas Lift
Simular un diseño existente
Calcular estatus de las válvulas (abiertas o cerradas)
Determinar la tasa de inyección de gas en las válvulas
Realizar revisiones al diseño o reacondicionamiento
propuesto
Determinar la posibilidad de un nuevo diseño
Profundidad de las válvulas de Gas Lift
Bomba Electro Sumergible (BES)
Base de datos con una lista de fabricantes de BES y modelos
(es decir, Reda, Centrilift, etc) disponible
Datos básicos: diámetro de revestimiento, tasas de flujo
mínima/máxima y velocidad base
Datos de diseño: velocidad de la bomba, número de etapas,
factor de carga
Selección de BES
Tabla de Comportamiento BES
Curva de Comportamiento BES
RP Bombeo Mecánico
BM Diagnóstico (para pozos existentes) :
• Analiza & Diagnostica las condiciones de la unidad
de bombeo, cabillas, condiciones de descarga y
bombeo.
• Completo entendimiento del estado de la producción
del pozo, i.e automáticamente calcula el
corportamiento del bombeo de fondo (vía carta
dinagráfica) con red inteligente neural.
RP Bombeo Mecánico
BM Diseño (para pozos nuevos / re-diseño de pozos
viejos) :
• Determina el mejor diseño para obtener alta
producción y ahorro de enegía
• Diseño de parámetros swab y re-diseña la sarta de
cabillas actual (para pozos viejos)
Módulo RP
BM Diagnóstico : BM Diseño :
•Datos PVT* y Datos de Diseño* •Datos PVT* y Datos de Diseño*
•Datos de Producción •Datos de Productividad del Pozo(IPR)*
•Datos de la sarta de cabillas •Selección del Equipo y Bomba
•Carta Dinagráfica •Diseño Optimizado Específico
Nota: * Datos transferidos del Modelo de Pozo PIPESIM
RP Diagnóstico – Resultados
RP Diseño –Resultados
Gracias por su atención