INGENIERIA DE YACIMIENTOS APLICADAS A ESTUDIOS INTEGRADOS
1. CARACTERISTICAS TERMODINAMICAS DE HIDROCARBUROS EN EL YACIMIENTO. PVT. INTERPRETACION Y VALIDACION.
2. FLUJO DE FLUIDOS EN MEDIOS POROSOS. PRUEBAS DE PRESIONES.
3. PROPIEDADES ROCAS FLUIDOS. Kr, Pc, CORRELACIONES. MANEJO EN SIMULADORES.
4. MECANISMOS DE PRODUCCION Y SU INTERPRETACION. MANEJO DE ACUIFERO Y CAPAS DE GAS.
5. CURVAS DE DECLINACION. RESERVACION. VALIDACION DEL POES.
6. PROGRAMAS COMPUTARIZADOS. MBAL. OTROS.
INGENIERIA DE YACIMIENTOSPROPIEDADES DEL MEDIO POROSO
POROSIDAD. DEFINICION. Ø = Volumen Poroso/Volumen Total TIPOS: ABSOLUTA(TOTAL) Y EFECTIVA GEOLOGICAMENTE: PRIMARIA Y SECUNDARIA SECUNDARIA:SOLUCIÓN, FRACTURAS Y
DOLOTIMIZACIÓN(CALIZAS EN DOLOMITAS) FACTORES QUE LA AFECTAN: EMPAQUE,
MATERIAL CEMENTANTE, DISTRIBUCIÓN DE LOS GRANOS, PRESENCIA DE FINOS(ARCILLA)
APLICACIONES
CALCULAR POROSIDAD DE UN EMPAQUE DE ESFERAS
Vt = Vg = 4/3 POROSIDAD =(Vt -Vg)/Vt 100%=47.6% INDEPENDIENTE TAMAÑO
ESFERAS
)2(3
rr 3
MEDICION DE POROSIDAD
MEDIDAS DEL VOLUMEN TOTAL, GRANOS, VOLUMEN POROSO
Vt SATURADA O CUBIERTA SUMERGIDA EN AGUA O EN MERCURIO
VOLUMEN DE LOS GRANOS: METODO DE MELCHER NUTTING. PESO DE LA MUESTRA SECA Y SATURADA. DETERMINA VOLUMEN
MEDICION DE POROSIDAD
MEDICION DE LOS GRANOS POR EL POROSIMETRO DE EXPANSION
MEDICION DEL VOLUMEN POROSO CON EL POROSIMETRO DE EXPANSION
METODO DE SATURACIONINYECCION DE MERCURIO
MEDICION DE POROSIDAD
METODOS ANTERIORES NO APLICAN PARA ROCAS DE CARBONATOS
MUESTRAS PEQUEÑAS NO INCLUYEN LAS FRACTURAS O CAVIDADES
REQUIEREN MUESTRAS MUY GRANDESLOS REGISTROS MIDEN POROSIDADES
QUE SE CORRELACION CON NUCLEOS
INGENIERIA DE YACIMIENTOSPROPIEDADES DEL MEDIO POROSO
POROSIDAD PROMEDIOSUMATORIA/NUMERO DE DATOSPONDERADO POR ESPESORPONDERADO POR AREAPONDERADO POR VOLUMENESTADISTICO: MEDIA, MODA-
DISTRIBUCIÓN ESTADISTICA
INGENIERIA DE YACIMIENTOSPROPIEDADES DEL MEDIO POROSO
LEY DE DARCY- PERMEABILIDADSISTEMA LINEAL Q = 1.127 K A(Pentrada -Psalida)/ LSISTEMA RADIALQ = 7.07 K h(Pe -Pwf)/ Ln(re/rw)SUPOSICIONES: MONOFASICO,
LAMINAR, 100% SATURADO
ANALOGIA DE LA LEY DE DARCY Y OTRAS LEYES FISICAS: OHM, FOURIER
LEY DE OHM: I= V/R, DONDE, R = L/A, = 1/ , LUEGO I = A V/L - ANALOGA A LA LEY DE DARCY
LEY DE FOURIER PARA LA TRANSMISION DE CALOR POR CONDUCCION
q = k´ A T/L ANALOGA A LA LEY DE DARCY
LAS ANTERIORES ANALOGIAS SON UTILES PUESTO QUE MUCHOS COMPLEJOS PROBLEMAS TANTO DEL FLUJO DE CALOR COMO ELECTRICIDAD HAN SIDO RESUELTOS ANALITICAMENTE Y SE PUEDEN EXTENDER AL FLUJO DE FLUIDOS EN MEDIOS POROSOS.
INGENIERIA DE YACIMIENTOSPROPIEDADES DEL MEDIO POROSO
TIPOS DE PERMEABILIDADABSOLUTA, EFECTIVA,
RELATIVA(EFECTIVA/ABSOLUTA)PROMEDIOS: PARALELO Y SERIE,
LINEAL Y RADIAL.P-L Y R: Kp = SUM(Kh)/SUM(h)S-L: Kp = SUM(L)/SUM(K/L)S-R: Kp = Ln(re/rw)/SUM((Ln Ri/Ri-1)/Ki)
APLICACION- COMBINACION CAPAS DIFERENTES PERMEABILIDADES
PERMEABILIDAD EQUIVALENTE DE CUATRO CAPAS PARALELAS CON IGUAL ANCHO Y LONGITUD QUE POSEEN LAS SIGUIENTES PROPIEDADES
CAPA ESPESOR, PIES PERMEABILIDAD, md
1 20 100 2 15 200 3 10 300 4 5 400 Kp = = 10000/50 = 200 md hkh /
APLICACION- COMBINACION CAPAS DIFERENTES PERMEABILIDADES
PERMEABILIDAD EQUIVALENTE DE CAPAS EN SERIE QUE TIENEN IGUALES ESPESORES PARA UN SISTEMA LINEAL Y RADIAL CON rw = 6 PULGS Y re = 2000 PIES CAPAS.. 1 2 3 4 LONGITUD, PIES 250 250 500 1000 PERMEABILIDAD,md 25 50 100 200
SISTEMA LINEAL Kp = = 2000/25 = 80 md SISTEMA RADIAL Kp = LN( re/rw)/
Kp = 30.4 md
kjLjLt
/
KjrjrjLN /)1/(
MEDICION DE PERMEABILIDAD
PERMEAMETROS - SE BASAN EN DARCY NO VALIDA PARA FLUJO TURBULENTO
EFECTO KLINKERBERG SI USA GAS, 1/P = 0SI USAN LIQUIDOS ASEGURARSE QUE NO
REACCIONAN CON LAS ROCASDIFICIL MEDIR ROCAS FRACTURADAS O
CON CAVIDADES
MEDIDA DE LA VARIACION DE PERMEABILIDAD
LA PERMEABILIDAD TIENE VARIACION GEOMETRICA
LA DEFINICION DE DYKSTRA-PARSONS SE RECOMIENDA
V = k A 50% Y 84.1%SUPONE DISTRIBUCION LOG NORMALDISTRIBUCION ACUMULADA DE kh vs
ACUM. h
kpkkp
CORRELACION ENTRE PERMEABILIDAD Y POROSIDADES
LA CORRELACION ENTRE POROSIDAD Y K, PERMEAB ES POBRE PERO CUANDO SE REALIZA USANDO LITOFACIES AUMENTA EL FACTOR DE CORRELACION
SE REQUIEREN MUCHOS NUCLEOS Y ANALISIS PARA LOGRAR UNA BUENA INFORMACION SOBRE ESTA RELACION
LOS MAPAS DE POROSIDAD ESPESOR, POROSIDAD VS PERMEABILIDAD, SECCIONES TRANSVERSALES, ENTRE MUCHOS OTROS, SE UTILIZAN PARA CARACTERIZAR LOS YACIMIENTO
EN UN YACIMIENTO SE USARON 11 H- AÑOS, 1.6 MM$
CORRELACIONES EN EL CAMPO ENTRE PERMEABILIDAD Y POROSIDAD
RELACION NUCLEO PERFIL ES NECESARIO PARA AJUSTAR LAS MEDIDAS EN LOS POZOS
MEDIDAS DE POROSIDAD Y PERMEABILIDAD SE CORRELACION DE MEDIDAS DE NUCLEOS
SE UTILIZAN DICHAS CORRELACIONES PARA ESTIMAR PERMEABILIDAD EN EL YACIMIENTO
CADA DIA SE UTILIZAN COMBINACIONES MAS COMPLEJAS QUE DEBEN USARSE CON RESERVA
INGENIERIA DE YACIMIENTOSDISTRIBUCION DE FLUIDOS
SATURACION Y DIST. DE FLUIDOSVOL FLUIDOS/VOL POROSO, o, w y gDETERMINADOS: REGISTROS Y LAB.DISTRIBUCIÓN: DENSIDADES- ROCASZONA DE TRANSICIÓN:CAPILARIDADCURVAS DE PRESION CAPILARTENSION INTERFACIAL: LABORATORIO
MEDIDAS DE SATURACIONES
METODO DE LA RETORTA: MIDE LOS VOLUMENES DE FLUIDOS. EXTRAIDOS. REQUIERE CALIBRACION
EXTRACTOR SOXHLETCENTRIFUGAREGISTROS - DIFERENTES TIPOS Y
ACTUALMENTE SE USAN MEDIDAS CONTINUAS EN LA GERENCIA DE YACIMIENTOS. VISUALIZACION.
CONDUCTIVIDAD ELECTRICA DE ROCAS SATURADAS CON FLUIDOS
FACTOR DE FORMACION: F = Ro/Rw, DONDE Ro ES LA RESISTIVIDAD DE LA ROCA Y Rw CUANDO ESTA SATURADA CON AGUA.
F ES FUNCION DE POROSIDAD Y GEOMETRIA F = C , C , CONSTANTE FUNCION DE
TORTUOSIDAD y m EN EL RANGO DE 1 A 2.SE MIDE EN LAB Y REGISTROS APLICANDO OHM Y EXISTEN VARIAS CORRELACIONES COMO ARCHIE C=1, m=1.3, HUMBLE C = 0.62, m = 2.15, OTROS
m
1
INGENIERIA DE YACIMIENTOSINTERACCION ROCA FLUIDOS
HUMECTABILIDAD: TENDENCIA- DRENAJE E IMBIBICION. APLICACIONES CASOS DE CAMPO.
EQUILIBRIO DE FUERZAS:W, O, ROCAANGULO DE CONTACTO. FIG. 2.9 MAYOR DE 90° HUMECTADO POR PET. PRESION CAPILAR FUERZAS DE RETENCIÓN DE o, w, g EN EL
YACIMIENTO
INGENIERIA DE YACIMIENTOSINTERACCION ROCA FLUIDOS
CURVAS DE PRESION CAPILAR RELACION Pc vs. Sw..Fig 2.17 FUNCION DEL TAMAÑO- DIST POROSMEDIDAS DE Pc- LAB: MEMB-MERCURIO CONVERSION DE LAB A CAMPO. Pcy = Pcl Ec. 2.38- Pc FUNCION DE K. Fig: 2-21
yl /
DATOS DE PRESION CAPILAR PROMEDIO
LA Pc DEPENDE DE LA PERMEABILIDAD Y SE MIDE EN NUCLEOS MUY PEQUEÑOS POR CONSIGUIENTE SE REQUIERE DETERMINAR CURVAS PROMEDIOS PARA LOS YACIMIENTOS
LA FUNCION J(Sw) = SE HA COMPROBADO QUE LA FUNCION J(Sw) MUESTRA
DISPERSION SIN EMBARGO SE PUEDE UTILIZAR PARA OBTENER CURVAS DE Pc PARA DIFERENTES NUCLEOS A TRAVES DE TODO EL YACIMIENTO
MANEJO ESTADISTICO PARA MANEJAR CORRELACIONES CON POROSIDAD Y PERMEABILIDAD COMO SIGUE:
Sw = a log K + C = a1 +a2 log k + C
)/(2/1
/ kPc
INGENIERIA DE YACIMIENTOSINTERACCION ROCA FLUIDOS
CURVAS DE PERMEABILIDAD RELATIVAS, KrKrf = Kef/Kabs…..VARIA DE 0 A 1. PETROLEO, AGUA Y GAS. PERM. RELATIVAS A 2 FASES. Fig 2.22DISTRIBUCION DE FLUIDOS: F(HUM.) Sor. FUNCION DE HUMECTABILIDADDETERMINACION:LAB, Pc, DATOS DE CAMPO,
ECUACIONES EMPIRICAS.
ECUACIONES EMPIRICAS - VALIDEZ. P61.WHAL Y ASOCIADOS: Kg/Ko ARENISCASCorey y Asoc: Kro y Krw, Arenas Cons. y no Cons.
Drenaje e Imbibición.Torcaso y Willie: Kg/Ko, drenaje en arenas. No
consolidadasPirson: Rocas con Porosidad Intergranular, dos
fases o, g y o, w para drenaje e imbibición.Willie y Gardner, Stone: Tres fases.
INGENIERIA DE YACIMIENTOSPERMEABILIDADES RELATIVAS
INGENIERIA DE YACIMIENTOS - PSEUDO CURVAS DE PERMEABILIDADES RELATIVAS
LAS PSEUDOS CURVAS DE PERMEABILIDADES RELATIVAS SON PROMEDIOS EN BASE A ESPESOR
Swp = Krw = Kro = GRAFICOS DE Krw y Kro vs Swp SE UTILIZAN
PARA TODO EL YACIMIENTO
Ni
i
Ni
i
ihiSwiihi11
/
N
i
N
i
hikihikikrwi1 1
/
N
i
N
i
hikikikroihi1 1
/
NUMERO CAPILAR - RELACION FUERZAS VISCOSAS A CAPILARES NUMERO CAPILAR SE DEFINE COMO LA RAZON DE LAS
FUERZAS VISCOSAS A CAPILARES EN UN PROCESO DE DESPLAZAMIENTO DE FLUIDOS EN UN MEDIO POROSO
Y ES IGUAL A O SEA LA RELACION ENTRE FUERZAS DESPLAZANTES Y RETENTIVAS EN EL MEDIO POROSO
EL PETROLEO RESIDUAL ES FUNCION DEL Nc VARIACION DEL Nc DE LA SATURACION
RESIDUAL DE PETROLEO VARIA DE 35 A 18% AL ALCANZAR Nc DEL ORDEN DE 10 A LAS (-2) SE OBTIENEN
REDUCCIONES IMPORTANTES DE LA Sor. UNA REDUCCION DEL ORDEN DE 1000 PARA UNA REDUCCION IMPORTANTE
/vNc LPkNc /
410710
DISTRIBUCION DE FLUIDOS EN MEDIOS POROSOS
SE HA COMPROBADO QUE LA FASE MOJANTE OCUPA LOS ESPACIOS MAS PEQUEÑOS
EL PETROLEO SE ENCUENTRA EN LOS CAPILARES MAS GRANDES
ESTUDIOS MICROSCOPICOS LO HAN DEMOSTRADO
EL PETROLEO RESIDUAL QUEDA ATRAPADO EN LOS ESPACIOS POROSOS MAS GRANDES
GRAFICO PRESION PROFUNDIDAD
GRADIENTE: g=.08, o=.33, w=.45PSI/P gr = gcn Vcn/Vcy=0.0763 g EE=1/(5.615 Bg), dPg/dD=.0763 gE/144 oy=( ocn5.615+Rs gcn)/5.615BodPo/Dd = oy/144. D(CGP)=(Po-Pb)/G,PSSUPONE GRADIENTE, G, CONSTANTE,
SIN EMBARGO VARIA
USO DE RFT PARA DETERMINAR CWP
LOS REGISTROS RFT SON LOS MEJORES MEDIOS PARA DETERMINAR CONTACTOS EN YACIMIENTOS
SE MIDEN LOS GRADIENTES EN LAS ZONAS DE PETROLEO Y AGUA, EN LOS POZOS
SE EXTRAPOLAN GRAFICOS DE P VS D PARA DETERMINAR LOS CONTACTOS
ES UNA TECNICA SENCILLA PERO DEBE APLICARSE CON SUMO CUIDADO.
MEDIDAS DE PRESIONES PARA APLICACIONES EN ING. DE YACIMIENTOS
PRESIONES EN LAS PERFORACIONES SE CALCULAN EN BASE A LA MEDIDA DE PRESION EN EL POZO A DETERMINADA PROFUNDIDAD, EL GRADIENTE Y LA DIF DE ALTURA
Pp = Pb + Gdh (Hp - Hb) PRESION AL DATUM DE REFERENCIA Pd = Pp + Gro(Hd - Hp) Pp = = PUEDE CALCULARSE CON OFM
N
i
N AiPiAi1
1/
N
i
N
i
ViPiVi1 1
/
GRAFICO PRESION PROFUNDIDAD
ADVENIMIENTO DE RFT Y MDTDETERMINACION DE LOS CAP Y CGPESTIMADO DE POES Y GOESPRESENCIA DE YACIMIENTOS
SEPARADOS POR LUTITAS COMPLICA EL ANALISIS
YACIMIENTOS CON CGP Y CAP
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
OOO
OOOO
OOO
to
dtto
tdtt
tdtt
B
1)()(
B
1)(
B
1
B
1)()(
B
1)(
B
1)(
B
1)()(
ASI ONCONSERVACI DE ECUACIONLA RESULTA d(OIP) OEXPANDIEND
B
1
B
1
B
1 ),(SS ¨),(VPVP
SIGUIENTES ESDEFINICION LAS USANDO
CALCULAN SE PRODUCCION DE MECANISMOS LOS POR RECOBRO EL
1VP
1VP d(OIP)
PETROLEO DEL ONCONSERVACI DE ECUACION
dSdPVddSdPVdSPV
SdPVdSdPVdSPVdSPVdOIPd
dSdVPd
BS
BS
oott
ot
t
o
t
ott
o
t
to
tdtt
o
t
o
to
dtt
o
dtto
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
ECUACIONES DE LOS DIF MECANISMOS
0DE Y DE,DGDG
RESULTA DG,A SUMA SE DE BURBUJEO, DE PUNTO DEL DEBAJO POR
1)(VPDG GAS, POR EMPUJE
1)(VPDW AGUA, POR EMPUJE
RESULTA DS, EN SUST LUEGO ),()()d(S ,
)(1
DE PETROLEO, DEL EXPANSION
1)(-VPDS ,SATURACION DE CAMBIO
1S d(VP)- DC ROCA,LA DE ONCOMPACTACI
dtt
dtt
o
dtt
to
t
og
t
ow
o
dtto
t
t
oo
dtt
o
BSd
BSd
SgdSwdADEMAS
VPdB
dSVP
BSd
B
MECANISMOS DE PRODUCCION POR SIMULACION DE YACIMIENTOS
EN EL CASO DE UNA CELDA SATURADA AL REDUCIR LA PRESION EL TERMINO (DE) ES NEGATIVO PUESTO QUE 1/Bo AUMENTA, LUEGO SE ACUMULA CON DG, POR CONSIGUIENTE
DG = DG - DE, Y DE = 0 LA PRODUCCION DE PETROLEO POR LOS
DIFERENTES MECANISMOS(DC,DW,DG,DE) SE SUMAN EN TODOS LOS BLOQUES E INTERVALOS DE TIEMPO.
EN ALGUNOS CASOS NO SE APLICA.
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
EL SIMULADOR ECLIPSE Y OTROS, LA COMPACTACION LA MODELAN COMO COMPRESIBILIDAD, DATOS TABULADOS EN FUNCION DE PRESION QUE PUEDE SER REVERSIBLES E IRREVERSIBLES, SE PUEDE INCLUIR LA POSIBILIDAD DE MODIFICAR LA TRANSMISIBILIDAD EN f(P), UN MODELO DE HISTERESIS.
CUANDO LA P DEL BLOQUE AUMENTA LA COMPACTACION PUEDE SER REVERSIBLE E IRREVERSIBLE
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
LA COMPACTACION MODIFICA LA TRANS. POR EL CAMBIO DEL VOLUMEN POROSO
LA PRESION DE SOBRECARGA MENOS LA PRESION DE FLUIDOS ES EL ESFUERZO EFECTIVO EN FUNCION DEL CUAL SE REPRESENTA LA COMPACTACION.
ALGUNAS ROCAS COMO LAS YESO SE PRODUCE UNA COMPACTACION INDUCIDA CUANDO SE CONTACTAN CON AGUA QUE PUEDEN SER SUMINISTRADAS POR TABLAS EN FUNCION DE P Y Sw
MECANISMOS DE PRODUCCIONEMPUJE POR GAS EN SOLUCION
EL PROCEDIMIENTO PRECEDENTE NO DIFERENCIA ENTRE LA PRODUCCION POR GAS LIBRE Y GAS EN SOLUCION
LA PRODUCCION PUEDE SUBDIVIDIRSE ENTRE LO CORRESPONDIENTE A GAS LIBRE Y GAS EN SOLUCION.
PUEDE UTILIZARSE UN INDICE Fs = 1 PARA GAS DISUELTO Y 0 PARA GAS LIBRE
LIBRE GAS POR PRODUCCION ,/)()1(
SOLUCION EN NTEORIGINALME Sg DE FRACCIONFs ,/)(**
BoSgdFsVPDF
BoSgdFsVPDS
MECANISMOS DE PRODUCCION - EMPUJE HIDRAULICO
UN PROCEDIMIENTO SIMILAR AL CASO PREVIO SE UTILIZA PARA DISTINGUIR ENTRE EL EMPUJE HIDRAULICO Y EL AGUA INYECTADA. ORIGINALMENTE Ft SE CONSIDERA 1.0, CERO AGUA INYECTADA.
ASI SE PUEDEN CALCULAR DWT Y DWR, LAS PRODUCCIONES POR We Y Wi
INICIALAGUA AL ECORRESPOND QUE FRACCIONLA ES
/)(*)1(*
/)(**
Ft
BoSwdFtVPDWR
BoSwdFtVPDWT
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
LA MALLA SE EXTIENDE AL ACUIFERO USANDO BLOQUES DIFERENTES CON MAYOR DEFINICION EN EL YACIMIENTO
REQUIERE + MEMORIA DEL COMPUTADOR PUEDE REPRESENTARSE EL ACUIFERO
MEDIANTE EL TERMINO PRODUCC/INYEC. SE UTILIZAN LAS DIVERSAS FORMAS DE
MANEJAR ACUIFERO PARA FLUJO CONT Y NO CONTINUO: HURST, VE&H, FEKOVITCH
UN EJEMPLO SE MUESTRA COMO SIGUE
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
TASA DE FLUJO DE AGUA EN BLOQUE m
lpca-Bw/día agua, deintrusión de constante C
acuífero el a referido bloque del área
contornoen prom.presión P /
)(/)()(
m
,
,11
mwtwmw
nmw
nwm
w
nenemw
PBCq
PPBCtB
tWtWq
m
m
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
EN ECLIPSE EL ACUIFERO SE MODELA POR UNA FILA DE CELDAS CONECTADAS AL YACIMIENTO
LA PALABRA CLAVE AQUNUM EN LA SECCION GRID SE USA PARA DEFINIR LAS PROPIEDADES COMO LONGITUD, AREA SECCIONAL, POROSIDAD, PERMEABILIDAD PRESION INICIAL PROFUNDIDAD, PVT Y Kr
EL ACUIFERO DEBE ESTAR CONECTADO A LA CARA DEL YACIMIENTO CON AQUCON
LAS CELDAS DEL ACUIFERO ESTAN CONECTADAS LAS CELDAS DEL ACUIFERO ESTAN AISLADAS DEL
YACIMIENTO
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
LA TRANSMISIBILIDAD ENTRE CELDAS DEL ACUIFERO VIENE DADO POR
ESTA EXPRESION SE USA TANTO PARA GEOMETRIAS CARTESIANAS Y CILINDRICAS
SE PRESENTAN LAS APROXIMACIONES DE FETKOVICH Y CARTER-TRACY
LONGi
XSECTi*PERMXi*2Ti DONDE
/1/1
TjTi
CDARCYTR
MANEJO DE ACUIFEROSFETKOVICH
EL MODELO DE FETKOVICH USA UNA APROXIMACION SIMPLIFICADA DE IP Y BALANCE DE MATERIALES ENTRE P Y We PARA EL ACUIFERO
LA INTRUSION DE AGUA DEL ACUIFERO SE MODELA POR LA ECUACION
acuífero elcon comunicado bloque del área :
acuífero-bloque del cara la de area el es m donde
)(
i
i
ii
iii
aiiai
A
Am
Am
donde
ddgppJQai
MANEJO DE ACUIFEROSFETKOVICH
LA PRESION DEL ACUIFERO SE OBTIENE POR BALANCE DE MATERIALES.
EL COMPORTAMIENTO DEL ACUIFERO DEPENDE DE DOS PARAMETROS, LA CONSTANTE DE TIEMPO DEL ACUIFERO Y EL INDICE DE PRODUCTIVIDAD
)( aaowota ppVCW
J
CtVwoTc
MANEJO DE ACUIFEROSFETKOVICH
BAJO LA SUPOSICION QUE LA PRESION DEL YACIMIENTO ES UNIFORME EN EL BLOQUE QUE CONECTA AL ACUIFERO, LA TASA DE INTRUSION DE AGUA PROMEDIO EN EL INTERVALO DE TIEMPO DELTA t ES
LA INTRUSION DE AGUA ACUMULADA SE CALCULA A CADA INTERVALO DE TIEMPO CUANDO SE ACTUALIZA LA PRESION, P
Tct
TctddgppJQ aiiaiai /
)/exp(1))((
MANEJO DE ACUIFEROSCARTER - TRACY
METODO SIMPLIFICADO QUE EVITA LA SUPERPOSICION, UTILIZANDO TABLAS Y SE APLICA A YACIMIENTOS DE FORMA ARBITRARIA
LOS PARAMETROS BASICOS SON Td, CONSTANTES TIEMPO E INTRUSION DE AGUA
cot
a
otwc T
trCh
Ck
rCT D
22
1
2
t, cy
MANEJO DE ACUIFEROSCARTER - TRACY
EL MODEL DE CARTER TRACY EXPRESA LA CAIDA DE PRESION EN EL CONTORNO EN TERMINOS DE PRESION ADIM, PID
LOS TERMINOS PDI Y SUS DERIVADAS SE OBTIENEN DE TABLAS SIMILARES A LAS CONOCIDAS DE VAN EVERDIGEN Y HURST
We SE CALCULA EXPLICITO AL FINAL Dt
conocidos términosdecalculan se by a donde )()(b(p-aQ
i bloque al acuífero del Qa flujo de tasala donde )(
iiai
I
tptt
tPIQ
pp
i
DDa
ao
MANEJO DE ACUIFEROSFLUJO CONSTANTE
EL FLUJO DEL ACUIFERO SE CONSIDERA CONSTANTE Y ES SUMINISTRADO AL SIMULADOR
LA TASA DE FLUJO EN EL BLOQUE DE UN ACUIFERO DE FLUJO CONSTANTE ES
LA CONSTANTE DEL ACUIFERO SE SUMINISTRA A CADA TIEMPO Y PUEDE VARIAR CON EL TIEMPO.
dato como acuífero del flujo Fa donde iiaai mAFQ
SIMULACION DE YACIMIENTOS DATOS DE PERMEABILIDADES
DEFINEN LAS CARACTERISTICAS FISICAS Y LA EXTENSION DEL YACIMIENTO PARA EVALUAR TRANSMISIBILIDADES
LAS FUENTES DE PERMEABILIDAD K SON LAS CURVAS DE PRESIONES TALES COMO RESTAURACION, DECLINACION, INTERFERENCIA, LABORATORIO, CORRELACIONES CON POROSIDAD
METODOS DE ANALISIS DE PRUEBAS DE PRESION: MUSKAT, HORNER, MDH, CURVAS TIPO, OTROS.
SIMULACION DE YACIMIENTOS DATOS DE PERMEABILIDADES
METODOS DE REGRESION: UTILIZA CORRELACIONES DE K EN FUNCION DE OTRAS VARIABLES COMO POROSIDAD, Sw,MEDIDAS DE REGISTROS
Swk
dSwcSwbak
3
22
250
SIGUIENTELA COMO NCORRELACIO OTRAS
SIMULACION DE YACIMIENTOS DATOS DE POROSIDADES
POROSIDAD PUEDE DETERMINARSE DE REGISTROS, PRUEBAS DE LABORATORIO Y CORRELACIONES
REGISTROS COMO DENSIDAD, NEUTRON Y SONICO. PREFERIDO EL PRIMERO.
LABORATORIO MEDIANTE MEDIDAS COMO SATURACIONES, POROSIMETRO DE EXPANSION E INYECCION DE MERCURIO
CORRELACIONES CON TIPO DE ROCA Y PROFUNDIDAD
SIMULACION DE YACIMIENTOS ESPESORES Y PROFUNDIDAES
ESPESOR SE OBTIENEN DE REGISTROS BIEN LA ARENA TOTAL Y NETA, TANTO PARA POTENCIAL COMO PARA POES.
EL ESPESOR TAMBIEN SE OBTIENE DE MAPAS ESTRUCTURALES.
LA PROFUNDIDAD SE OBTIENE DE REGISTROS Y RECORDS DE PERFORACION
LAS MEDIDAS SE CORRELACIONAN CON ESPACIO MEDIANTE GEOESTADISTICA.
OFM, ALMACENA LA BASE DE DATOS
DATOS DE SATURACIONES DE FLUIDOS Y Pc
LAS ZONAS DE INTERES SON EN LOS CONTACTOS DE FLUIDOS
POR ENCIMA DEL CONTACTO AGUA PETROLEO LA Sw ES CONSTANTE
PUEDEN OBTENERSE DE LOS REGISTROS, DATOS DE LABORATORIO Y CURVAS DE PRESION CAPILAR.
LAS PRESIONES CAPILARES PUEDEN DETERMINARSE DEL LABORATORIO
DATOS MUY IMPORTANTES
DATOS DE PERMEABILIDADES RELATIVAS PARA SIMULACION
LAS PERMEABILIDADES RELATIVAS SON LOS DATOS MAS DIFICILES DE EVALUAR
LAS CURVAS Kr QUE SE OBTIENEN SON Krow, Krog Y Krgw, SEGÚN EL CASO.
Kr SE DETERMINAN EN EL LABORATORIO CON MEDIDAS DIRECTAS, CURVAS DE Pc, INF DE CAMPO Y CORRELACIONES
LABORATORIO USA BUCKLEY LEVERETT, INTEGRANDO Pc, CAMPO Kg/Ko BALANCE DE MATERIALES Y CORRELACION DE STONE
SIMULACION DE YACIMIENTOSFENOMENO DE HISTERESIS
LA OPCION HISTERESIS EN SIMULACION REQUIERE ESPECIFICAR DIFERENTES FUNCIONES DE SATURACION PARA DRENAJE E IMBIBICION Y EN CADA CELDA SE SUMINISTRAN DOS TABLAS
LA Krd SE INICIA A LA MAXIMA SATURACION DE LA FASE MOJANTE, Swmaxd
EN FORMA SIMILAR, SI Sw AUMENTA SE USA LA CURVA DE IMBIBICION DESDE Swi
LA FIG SIGUIENTE MUESTRA DIFERENTES FASES.
SIMULACION DE YACIMIENTOSKr- HISTERESIS FASE NO MOJANTE
LA FASE NO MOJANTE ES PETROLEO EN O-W, GAS EN O-G, O-W, O-W-G
EN LA FIGURA SIGUIENTE LA CURVA 1-2 DRENAJE, Y 2-3 IMBIBICION
LA SATURACION CRITICA DE LA CURVA DE IMBIBICION ES MAYOR QUE PARA DRENAJE
LAS DOS CURVAS SE UNEN A Snwmáx. SI EL DRENAJE ES COMPLETO LA CURVA
ALCANZA 3 PERO SI NO SE REVIERTE EN 4 Y LAS Sncrt FUNCION Snw ALCANZADA
SIMULACION DE YACIMIENTOSKr- HISTERESIS FAS. NO MOJANTE
LA GENERACION DE LAS CURVAS PUEDE REALIZARSE POR LOS METODOS DE CARLSON-SPE 10157 Y KILLOUGH ACTAS AIME 1976
EL METODO DE CARLSON PRODUCE UNA CURVA PARALELA A LA CURVA DE IMBIBICION.
METODO KILLOUGH ES MAS ELABORADO Y GENERAL
SI EN LAS SIMULACIONES SE PRESENTAN PROBLEMAS DE CONVERGENCIA-REVISAR LAS CURVAS DE Kr
Krn curva de imbibición
Curva de drenaje
Saturación fas mojante
HISTERESIS EN LA FASE NO MOJANTE
S no mojante
HISTERESIS EN LA FASE NO MOJANTE
SATURACION CRITICA ATRAPADA, Sncrt
0.1 a de falta a Shy),-a(Snmax1A donde Sncrd)-C(ShyA
Sncrd-ShySncrdSncrt
USASE SncrdA TIENDESncrt SI ,))((
)(
)()()(
1
Sncrd-Sncri
1C
)(1
KILLOUGH DE METODO
SncrtShy
SncriSnmaxSncrtSnSncriSnorm
SnmaxKrnd
ShyKrndSnormKrniSnKrn
SncrdSnmaxdonde
SncrdShyC
SncrdShySncrdSncrt
SIMULACION DE YACIMIENTOSKr- HISTERESIS FASE MOJANTE
CURVAS TIPICAS DE Kr DE LA FASE MOJANTE VER EN LA FIGURA SIGUIENTE
CURVA 1-2 DRENAJE Y 2-3 IMBIBICION Y LAS DOS CURVAS SE UNEN A Swco
LA MAXIMA SATURACION DE IMBIBICION ES 1-Sncri. Swco = 1 - Snmax
SI EL PROCESO DE DRENAJE SE REVERSA EN 4 LA CURVA SE OBTIENE POR EL METODO DE KILLOUGH.
SI DRENAJE E IMBIBICION COINCIDEN SE SOLO SE UNEN EN 4 Y 5
SIMULACION DE YACIMIENTOSKr- HISTERESIS FASE MOJANTE
LA CURVA DE DRENAJE QUE SE REVERSA EN 4, Shy MAXIMA SATURACION NO MOJANTE Y LA Sw = 1 - Sncrt. SE USA EL METODO DE KILLOUGH Krnw
)1(
Sncri)- Krwi(1
Snorm)-1Shy))Krwi(-Krwd(1-Sncrt)-(Krw(1Shy)-Krwd(1Krw(Sw)
INTERMEDIACURVA EN Sw ,SATURACIONA DETERMINADA Kr LA
CURVATURA DE PARAMETRO
)))(1()1(()1()1(
SncriSnfSnorm
ASncrdSncri
SnrdScrtSncriKrwdSncriKrwiSncrtKrwdSncrtKrw A
curva de imbibición
Curva de drenaje
S fase mojante
HISTERESIS EN LA FASE MOJANTE
SIMULACION DE YACIMIENTOSMODELO SISTEMA MOJADO PETROLEO
SE APLICA EL MODELO DE CARLSON Y KILLOUGH A LA FASE NO MOJANTE AGUA
EL AGUA QUEDA ATRAPADO POR EL PETROLEO CURVA IMBIBICION SE TOMA AGUA
AUMENTANDO EN LUGAR DE LA CURVA DE PETROLEO AUMENTANDO. DRENAJE DEL PETROLEO
CURVA DE IMBIBICION SIEMPRE ES AGUA AUMENTANDO INDEPENDIENTE DEL MODELO
SIMULACION DE YACIMIENTOSHISTERESIS DE LA CURVA DE Pc
EN LA CURVA DE Pc, LA CURVA 1-2 DRENAJE Y LA 2-3 IMBIBICION
SI LA CURVA SE REVERSA EN 4 ALCANZA LA SATURACION CRITICA DE LA FASE NO MOJANTE EN 5 QUE ES UN PROMEDIO ENTRE CURVAS DRENAJE E IMBIBICION
LA ECUACION DE KILLOUGH. VER TRABAJO. F SE CALCULA COMO SIGUE
0.1)parametro f(S,F )( cdcicdc PPFPP
SIMULACION DE YACIMIENTOSHISTERESIS DE LA CURVA DE Pc
DONDE EN LA ECUACION ANTERIOR F ES
.histéresis igualmente ocurre gas,-agua ó fases tresde caso elEn
análogo es drenaje a imbibición de reversoun Para
Pci)-(PcdG Pci Pc
reverso, segundoun Para
Sncrt Sn para saturación Máxima Swma
.histéresis la de reverso elen Sw :Swhy
0.1 deorden del Parámetro :E
)11
/()1
ESwhy-Sw
1(
EESwhySwmaE
F
S fase
mojante
Curva de
Imbibición
Curva de
Drenaje
HISTERESIS DE LA CURVA DE Pc
SIMULACION DE YACIMIENTOSPERMEABILIDAD RELATIVA A 3 FASES
LAS PERMEABILIDAD RELATIVA A 3 FASES PUEDEN OBTENERSE DE FORMULAS O TABLAS PARA UTILIZARLA EN LOS SIMULADORES
ECLIPSE SUPONE EL MODELO DE LA FIGURA SIGUIENTE. EL AGUA Y EL GAS SE SUPONEN SEGREGADOS, MIENTRAS QUE EL PETROLEO SE SUPONE A LA SATURACION PROMEDIO DEL BLOQUE.
g-o w,-o SISTEMAS EN PET AL REL DADPERMEABILI
BLOQUE DEL PROMEDIOS ESSATURACION ,Sy ,S
GAS DELA ZONA ENAGUA DE SATURACION S , )(
ow
wco
rogrow
g
wcowg
rowwcowroggro
kyk
S
SSS
kSSkSk
PETROLEO
AGUA
GAS
SIMULACION DE YACIMIENTOSPERMEABILIDAD RELATIVA A 3 FASES
SIMULACION DE YACIMIENTO MODELO DE STONE-MODIFICADO
LA FORMULA BASICA VIENE DADA POR
MIN RESID PET ATS )1/(
)1/(
))1(/(
))1(/(
S DONDE )1/()(SS
CONNATAAGUA DE PRES EN PETA
RELATIVA PERM:k DONDE ,
om
oo
rocw
SSSSSS
SSSSS
DONDE
SSkkF
SSkkF
SSSSS
FFSSkk
omwcogg
omwcogw
grocwrogg
wrocwroww
omomwcoomo
gworocwro
SATURACION DE AGUA
SAT DE
G
A
S
DE AGUA Som SE INTERPOLA A PARTIR DE LA SATURACION Sw
SATURACION DE AGUA
SAT DE
G
A
S
DE AGUA Som SE INTERPOLA A PARTIR DE LA SATURACION Sg
SATURACION DE AGUA
SAT DE
G
A
S
PARA IMBIBICION Y DRENAJE SE CALCULAN DE LAS TABLAS CORRESPONDIENTES
SIMULACION DE YACIMIENTOSPERMEABIL. RELATIVAS A 3 FASES
SEGUNDO MODELO DE STONE MODIFICADO
Krog PERMEABILIDAD RELATIVA AL PETROLEO. EN UN SISTEMA PETROLEO GAS, Sw = Swc
Krow PERMEABILIDAD RELATIVA AL PETROLEO EN UN SISTEMA PETROLEO Y AGUA
rgrw
rgrocw
rogrw
rocw
rowrocwro
kk
kk
kk
k
kkk
-
)(( )
SIMULACION DE YACIMIENTOSESCALAMIENTO- TABLAS DE SATURACIONES
LAS TABLAS DE SATURACIONES DE PUNTOS EXTREMOS PERMITE DEFINIR SATURACIONES
CONNATAS, CRITICAS Y MAXIMAS EN LA DESCRIPCION DEL FLUJO DE FLUIDOS.
LA OPCION PERMITE SIMULAR YACIMIENTOS QUE POSEEN VARIACION INICIAL DE SATURACIONES CONNATAS O
CRITICAS EN UNA O MAS FASES PRESENTES. EL METODO TIENE APLICACIONES EN EL USO DE PSEUDO
FUNCIONES Y SATS VARIABLES.
SIMULACION DE YACIMIENTOSESCALAMIENTO- TABLAS DE SATURACIONES
MODELOS 3F SE DEFINEN 8 PUNTOS EXTREMOS SWL SATURACION DE AGUA CONNATA SWCR, SATURACION DE AGUA CRITICA SWU, SATURACION DE AGUA MAXIMA SGL, SATURACION DE GAS CONNATA SGCR, SATURACION DE GAS CRITICA SAGU, SATURACION DE GAS MAXIMA SOWCR,SATURACION CRITICA DE PETROLEO,O-W SOGCR,SATURACION CRITICA DE PETROLEO,O-G
SIMULACION DE YACIMIENTOSESCALAMIENTO- TABLAS DE SATURACIONES
LA OPCION ESCALAMIENTO PERMITE DEFINIR NUEVOS VALORES PARA CADA UNA CELDA MANTENIENDO DATOS CONSISTENTES EN LAS TABLAS DE SATURACION
EL CONJUNTO DE LOS 8 PUNTOS EXTREMOS SE APLICAN EN CORRIDAS DE 2 FASES
CUANDO LOS VALORES Kr Y Pc, SE REQUIERE CALCULAR DETERMINADAS SATURACIONES EQUIVALENTE PARA USAR LOS DATOS NO ESCALADOS. EJEMPLO SIGUIENTE
SIMULACION DE YACIMIENTOSESCALAMIENTO- TABLAS DE SATURACIONES
CUANDO SE REQUIERE Kr Y Pc A DETERMINDA SATURACION SE USA UNA TRANSFORMACION PARA DETERMINAR LA SATURACION EQUIVAL. PARA USAR LAS TABLAS NO ESCALADAS
UNA CELDA CON SATURACION DE AGUA S, CUYAS SATURACIONES CONNATAS Y MAXIMAS SON Sco Y Smax, DONDE CUYAS VALORES NO ESCALADOS SON S´co Y S´max, LAS Kr Y Pc SE EVALUAN A S´CALCULADAS COMO
)()´´)((´´ ScoSmax
coSmaxSSScoSS co
SIMULACION DE YACIMIENTOSESCALAMIENTO- TABLAS DE SATURACIONES
ADEMAS ES POSIBLE ESCALAR LOS VALORES DE Kr Y Pc USANDO PALABRAS CLAVES KRW Y PCW
LA OPCION HISTERESIS SE ACTIVA CUANDO SE USAN LAS PALABRAS CLAVES ANTERIORES PARA ESCALAR LAS CURVAS DE Kr PARA DRENAJE.
PARA LAS CURVAS DE IMBIBICION SE UTILIZAN LAS PALABRAS CLAVES ISWL, ISWCR, ISWU, ISGL, ISGCR, ISGU, ISOWCR, ISOGCR
CAMBIOS DE IMBIBICION DRENAJE: CARLSON
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Pc
LOS PUNTOS EXTREMOS DE LAS CURVAS DE Pc SON LAS SATURACIONES CONNATAS Y MAXIMAS, SWL Y SWU PARA W-O; SGL Y SGU PARA O-G.
ES POSIBLE MODIFICAR LOS PUNTOS EXTREMOS PARA LAS CURVAS DE Pc SIN MODIFICAR EL CORRESPONDIENTE ESCALAMIENTO PARA Kro
SWLPC Y SGLPC SE USAN PARA LAS SATURACS CONNATAS
DADAS SWL Y SWLPC SE ESCALAN Kr Y Pc.
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Pc
ESCALAMIENTO VERTICAL ES POSIBLE ESCALAR LA MAXIMA Pc EN UN
BLOQUE EN BASE A CADA BLOQUE SI SE ESPECIFICAN LOS MAXIMOS Pco-w,o-g COMO
PCW Y PCG. PARA EL CASO O-W Pc = Pct PCW/Pcm Pct: Pc DE LA TABLA Pcm: MAXIMO Pc DE LA TABLA A Sw= Swco PCW: MAXIMO Pc DE LOS DATOS PCW
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Kr
SE USAN DOS OPCIONES PARA ESCALAR Kr SI NO SE ESPECIFICA, EL ESCALAMIENTO
PRESERVA LA Kr EN DOS PUNTOS EXTREMOS SE SUPONEN LOS PUNTOS EXTREMOS DE Kr PARA CADA
FASE EN SISTEMAS O-W, O-G Krw : SWCR & SWU Krg : SGCR & SGU Krow : SOWCR & (1 - SWL -SGL) Krog : SOGCR & (1 - SWL - SGL)
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Kr
EN EL CASO DE 3 FASES SE USAN LOS PUNTOS EXTREMOS SIGUIENTES
Krw: SWCR, (1 -SWL-SGL) & SWU Krg: SGCR, (1-SOGCR -SWL) & SGU Krow: SOWCR, (1- SWCR-SGL) & (1 - SWL-SGL) Krog: SOGCR, (1 -SGCR-SWL) & (1. -SWL-SGL) EN EL CASO DE CORRIDAS EN SISTEMAS W-G LOS PUNTO
EXTREMAS SE TOMAN Krw : SWCR, (1 - SGCR) & SWU Krg : SGCR, (1 - SWCR) & SGU
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Kr
EL SEGUNDO METODO DEBE INTERPRETARSE COMO CORRIDAS EN DOS FASES MOVILES PRESERVANDO LAS Kr EN LOS EXTREMOS DE LA REGION DE 2 FASES
PUEDEN PRESENTARSE PROBLEMAS DE CONVERGENCIA CUANDO EL PUNTO MEDIO SE APROXIMA A LA SATURACION MAXIMA PUEDE ORIGINAR DISCONTINUIDADES EN Krw
ES NECESARIO TOMAR PRECAUCIONES PARA EVITAR LOS CAMBIOS BRUSCOS EN Kr
SIMULACION DE YACIMIENTOSESCALAMIENTO VERTICAL DE CURVAS DE Kr
ES POSIBLE ESCALAR LA Kr A LA MAXIMA Y CRITICA O RESIDUAL DE LA FASE ASOCIADA
SE UTILIZAN LOS PUNTOS EXTREMOS KRW, KRG Y KRO Y SUS DERIVADAS PARA LA SATURACION MAXIMA Y KRWR, KRGR, KRORW Y KRORG ASI COMO SUS DERIVADAS PARA LA SATURACION CRITICA O RESIDUAL DE LA FASE ASOCIADA.
EL ESCALAMIENTO DE Kr AL AGUA SE MUESTRA A CONTINUACION E IGUAL PARA Kro Y Krg
))máx(
)()((
tablaKrw
bloqueKRWtablaKrwKrw
SIMULACION DE YACIMIENTOSESCALAMIENTO DE CURVAS DE Kr
SI SE USA KRWR, EL ESCALAMIENTO DEBE HONRAR Kr A LA SATURACION CRITICA SR DE LA FASE ASOCIADA.
SR = 1 - SOWCR EN SISTEMAS O-W, O-W-G SR = 1 - SGCR EN SISTEMAS AGUA GAS LUEGO LOS DOS CASOS SON
SI Krwmax = Krw(SR) SE SUPONE UNA LINEA RECTA ENTRE KRW Y KRWR. EN CUALQUIER CASO LOS PUNTOS EXTREMOS SE HONRAN
))()((Krw(SR)-Krwmax
KRWR-KRWKRWRKrw
bla)Krw(SR)(ta
e)KRWR(bloquKrw(tabla) Krw
SRKrwtablaKrwSRSW
SRSW
VARIACIONES DEL PUNTO EXTREMO A TRAVES DE LA ZONA DE TRANSICION
EN YACIMIENTOS DONDE SE REQUIERE MODELAR LA VARIACION VERTICAL DE LA SATURACION CRITICA DE UNA FASE PARA QUE COMIENCE A FLUIR COMO SE PRESENTA EN LA FIGURA SIGUIENTE.
ES UNA FORMA CONVENIENTE ES ESPECIFICAR LAS SATURACIONES CON PROFUNDIDADES
ESTAS SATURACIONES TAMBIEN SE PUEDEN CALCULAR EN BASE A LAS MOVILIDADES DE LOS FLUIDOS
DURANTE LA SIMULACION LAS Kr SE RE-ESCALAN
SIMULACION DE YACIMIENTOSMOVILIDADES INICIALES DE FLUIDOS
EL MODELO DEL YACIMIENTO SE DEBE INICIAR CORRECTAMENTE EN EL VOLUMEN Y EN MOVILIDADES DE LOS FLUIDOS
SI LO ANTERIOR NO SE CUMPLE LOS FACTORES DE RECOBRO SE PRONOSTICAN CON ERROR.
SE CALCULA UNA SATURACIÓN PROMEDIO Y UNA MOVILIDAD PROMEDIO
LOS MODELOS TIENEN DIFERENTES MANERAS DE CALCULAR LAS VARIABLES ANTERIORES
ZONA DE
TRANSICION
ZONA DE
AGUA
CONTACTO
AGUA- PETROL.
PETROLEO
AGUA
VOLUMEN INICIAL DE PETROLEO = V(A+B)
VOLUMEN DE PETROLEO MOVIL = VA
SE USA UNA MALLA FINA PARA ESTIMAR Sprom y Mp
SE USA LA IDEA DE PSEUDO FUNCIONES
VARIACION DEL PUNTO EXTREMO DE EN LA ZONA TRANSICION
Pc
SIMULACION DE YACIMIENTOSFUNCION J DE LEVERETT
LA FUNCION J DE LEVERETT QUE SE CORRELACIONA CON PROPIEDADES DE ROCA. ECUACIONES BASICAS
)(
2/1
;
2/1
;
)(
unidadesconstwo
escalatablaswo
UK
CteF
FPcPcKPc
SwJ
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS.
PVT- MUESTRAS FLUIDOS DE YAC.ANALISIS COMPORTAMIENTO PVTPRUEBAS DE LABORATORIOCORRELACIONES EMPIRICASPRESION DE BURBUJEO(Pb): DEFINICIÓN.
MEDIDAS. CORRELACIONES DE STANDING, BORDEN Y RZASA, LASATER, VASQUEZ Y BEGGS, GLASE, OTRAS
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
PRESION DE ROCIO,Pd, DEFINICION. MEDIDA EN LAB., CORRELACIONES DE NEMETH Y KENEDY.
RELACION GAS EN SOL. PET. (Rs). DEFINICION. MEDICION. LABORATORIO Y CORRELACIONES. FIG. 6 PAG 16.
DEPENDE: P, T, °API, Ggas, TIP LIB.
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
CORRELACIONES ESTIMAR Rs.: BEAL, STANDING, LASATER, LASATER, VASQUEZ- BEGGS, OTROS
COMPRESIBILIDAD DEL GAS, Z.GASES REALES: PV = n Z R T. 0.8-1.2MEDIDA EN LAB Y CORRELACIONES
COMO STANDING Y KATZ. FUNCION DE GRAVEDAD, TEMPERATURA Y PRESION PSEUDOREDUCIDAS.
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
FACTOR VOLUMETRICO DEL GAS BgVOLUMEN EN EL YACIMIENTO DE LA
UNIDAD VOLUMEN EN SUPERFICIEBg = 0.00504 Z T/p, BY/PCN-FIG 14MEDIDO EN LAB- CORRRELACIONESFACTOR VOLUM. PET., Bo. F 16DEPENDE DE P Rs y T ORD 10EXP -3
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
Bo = (1 + Vp) (1 + Vt), FS.18-19, (1.2-1.7)EQUIVALENTE A VOLUMEN DE
PETROLEO EN SUPERFICIE MAS GAS EN SOLUCION. LAB Y CORR.
CORR.DE STANDING Y KATZ, STANDING, VASQUEZ Y BEGG, BORDEN Y RATZ, OTROS.
COMPRESIBILIDAD POR ENCIMA Pb.
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
FACTOR VOLUMETRICO TOTAL, BtRAZON DE VOL DE PETROLEO MAS
GAS DISUELTO Y LIBRE/ VOL. PET.FIG 26..POR ENCIMA Pb, Bo = Bt.MEDIDA LAB Y CORRELACIONES
COMO STANDING, GLASE, OTROS.CORRELACIONES Y REQUIERE
INFORMACION MAS DETALLADA.
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
FUNCION Y- AJUSTE DE DATOS PVTY = (Pb - P)/(P (Bt/Btb - 1)); Bt = vtY vs P, LINEA RECTA PAPEL NORMALDATOS PRESENTAN DISPERSIONCORRIGEN CON LA MEJOR RECTA YUNA FORMA DE VALIDAR DATOSDESVIACIONES CERCA DE Pb
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
FACTOR VOLUMETRICO AGUA, BwBw = VOL YAC/VOL CN-INCLUYE GASCORRELACION DODSON-STANDING,
NUMBERE-BRINHAM-STANDING, McCAIN. VALORES APROX A 1.O...
Rsw: GAS EN SOLUCION AGUA CORRELACIONES
CULBERSON .MACKETTA. VALORES 10-5O PCN/B
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
VISCOSIDAD. RESISTENCIA A FLUIRVISCOSIDAD DEL PETROLEO, DEPENDE DE P, T -GAS EN SOL.COMPORTAMIENTO FUNCION DE PbANALISIS FIG. 37. DETERMINA PVTCORRELACIONES: BEAL, BEGGS Y
ROBINSON, ASTM, GENERALIZADA (PETROLEOS MUERTOS)
O
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
ERRORES CORR. DE VISCOSIDADASTM MENOR DE 25%, BEAL 4.64%, BEGGS
Y ROBINSON O.64%, GENERALIZADA 20%.UNIDADES CP, POISE, DINAS SEG/CM
SEGUNDOS SAYBOLT, FUROL, ENGLER, REDWOOD
CONVERSION FIG 44 - VCS PET. ENTRE (O.3-MILLONES)CPS
2
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
VISCOSIDAD DEL GAS MENORES QUE PETROLEO Y AGUAORDEN ENTRE 0.01 Y O.05 CPSDETERMINADO EN LABORATORIO Y
CORRELACIONES: GPSA, CALHOUN, CARR, KOBAYASHI Y BURROWS, LEE Y ASOCIADOS.
VISCOSIDAD AGUA..0.2 - 2 CP. CORR. VAN WINGEN-OTRO, MCCAIN.
g
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
COMPRESIBILIDAD DEL PETROLEO CAMBIO DE VOLUMEN CON P-T CTE
Co = - 1/V (dV/dP) = 1/Bo(dBo/dP), T=C Bo=Bob EXP(- Co(P - Pb)) Ec 141. P PbCo = -((Bo - Bob)/(P - Pb))/Bob, P PbCo = - 1/Bo((dBo/dP)-Bg(dRs/dP)) a T
CTE A PRESIONES DEBAJO Pb
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
CORRELACIONES PARA CoCALHOUN, TRUBE, VASQUEZ Y BEGGS,
RAMEY, MACAIN Y ASOC.COMPRESIBILIDAD DEL GASCg = 1/P - 1/Z dZ/dP A T CTE.Cg = 1/P A T CTE, PARA GAS IDEAL.CORRELACIONES TRUBE, MATTAR Y
ASOC., ORDEN 5 * 10 EXP(-4)
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
COMPRESIBILIDAD DEL AGUA, CwCORRELACIONES DE DODSON Y STANDING,
MEEHAN, OSIF, OTROSVARIA ENTRE 2 - 4 10 EXP(-6)Cw=-(dBw/dP -Bg dRsw/Dp)Bw, P PbCOMPR. PROMEDIO Y EFECTIVACt = So Co + Sw Cw + Sg Cg +CfCoe=Ct/So, Cwe=Ct/Sw,Cge=Ct/SgK/ = Ko/ o + Kw/ w + Kg/ g
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
COMPRESIBILIDAD DE ROCAS-FORM.CORRELACIONES HALL, VAN DER
KNAAP, RATT Y NEWMAN.CALIZAS Y ARENISCAS FUNCIONPOROSIDAD VIENE EN FRACCION ORDEN 1-1OO EXP(-6) EN LPC(-1)MUY IMPORTANTE EN VARIOS YAC.
DENSIDAD DEL PETROLEOFUNCION DE LAS DENSIDADES DEL
PETROLEO Y DEL GAS, Rs o = ( ost + 0.0136 g Rs)/ Bo, lbs/PCCORRELACIONES STANDING, METODO
DE STANDING, ALANI Y KENNEDY Y GOTTFRIED
g = P M/ZRT, VER UNID. TODAS
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
DENSIDAD DEL GAS g = P M / Z R T. UNIDADESDENSIDAD AGUA, w = wst /BwFACTOR DE EXPANSION DEL PET. o = op/ or (T - Tr)..Ec 218FAROUQ ALI: SIN DAT 5 10(-4)°F(-1)EXP GAS. g = 1/T+1/Z(dZ/dT) P CTE
INGENIERIA DE YACIMIENTOPROPIEDADES DE LOS FLUIDOS
EXPANSION TERMICA DEL AGUA w = Bw/ (Bw T) PVT vs PRUEBAS DE LAB.SOFTWARE/PVT. MABAL, ECLIPSEUSO DE EXCELFACTORES DE CONVERSION..p 159UNIDADES…p…162
BALANCE DE MATERIALESSUPOSICIONESVOLUMEN POROSO DEL YACIMIENTO
SE CONSIDERA CONSTANTELOS CALCULOS SE EFECTUAN A UNA
PROFUNDIDAD DE REFERENCIA (DATUM)
PVT REPRESENTAN LOS FLUIDOS EN EL YACIMIENTO
LA EXPANSION DE ROCA Y AGUA CONNATA SON DESPRECIABLES
BALANCE DE MATERIALESSUPOSICIONESLOS FLUIDOS EN EQUILIBRIO . NO HAY
DIRECION AL FLUJOLA TEMPERATURA SE CONSIDERA
CONSTANTEBw, Rsw, PROPIEDADES DE LA ROCA SON
CONSTANTES Y SE CONOCEN.SE CONOCEN LOS DATOS PRODUCCION
Y PRESION Y SE PUEDE ESTIMAR EL FUTURO
BALANCE DE MATERIALESECUACION GENERAL
BALANCE DE MATERIALES EN EL YACIMIENTO ES IGUAL A:
FLUIDOS PRODUCIDOS, BY = (EXP. DEL(PETROLEO + GAS EN SOLUCION) BY + (EXPANSION DE LA CAPA DE GAS) BY + (REDUCCION DEL VPHC) BY, TODO A CONDICIONES DEL YACIMIENTO
BALANCE DE MATERIALES:DEFINICIONES
N: PETROLEO INICIAL, BNG: GAS LIBRE INICIAL EN LA CAPA
DE GAS, BYm = G/ NBoi, BY/BYGp: GAS PRODUCIDO ACUM, PCNNp: PETROLEO PROD ACUM, BNRp = Gp/Np, PCN/BN
BALANCE DE MATERIALES:CALCULOS BASICOSEXPANS. DEL PETROLEO= N(Bo - Boi), BYEXPANS. GAS LIBERADO=N(Rsi-Rs)Bg, BYEXPANS CAPA GAS=mNBoi(Bg/Bgi-1), BYREDUCCION VPHC= CAMBIOS DE
VOLUMENES DE AGUA, ROCA, CAPA GAS Y ENTRADA DE AGUA(ACUIFERO)
d(VPHC) = -dVw + dVporos - dVacuífero - dVcapa de gas
BALANCE DE MATERIALESECUACIONES BASICAS -d(VPHC) = (1+m)NBoi (Cw Swc+ Cf) P /(1-
Swc) + We = Efw + WePROD FLUIDOS=Np(Bo+(Rp-Rs)Bg) + WpBw
- WiBw - Gi Bg = F, BY IGUALANDO SE OBTIENE EGBMVARIANTES DE LA EC DE LA BAL MAT.SE CONOCE LA INFORMACION DE Np,
Gp, PVT, WP, m; EXCEPTO N Y We
BALANCE DE MATERIALESEGBM COMO LINEA RECTA
LA EGMB COMO UNA LINEA RECTAF = N( Eo + mEg + Efw) + WeEo = EXP DEL PET + GAS INIC EN SOLEo = Bo - Boi + (Rsi -Rs) Bg, BY/BNEo = Bt - BtiEg = Boi(Bg/Bgi - 1), BY/BNWe = C f(P,t), SOLUCIONES DE HURST,
V.E&H, FEKOVITCH, CARTER Y TRACY
BALANCE DE MATERIALESECS SIMPLIFICADAS
CASO We =0, Wi, Gi=0, EGBM REDUCE ANp(Bt+(Rp-Rsi)Bg) +WpBw = N((Bt-Bti) +
Bo(Cw Swi+ Cf)/(1-Swi) P)CASO P ENCIMA DE Pb, PETROLEO NO
SATURADO, Rp = Rs = Rsi, Bt = BoNp/N = Boi/Bo Ce PCe = (Co So + Cw Sw + Cf)/(1 - Swi)
BALANCE DE MATERIALESECS SIMPLIFICADAS
DEBAJO Pb, CASO We = 0, Cf = 0Np/N = (Bt - Bti)/(Bt + (Rp - Rsi) Bg)PREDICCIONES REQUIEREN LAS ECS.So = (1 - Np/N)(Bo/Boi)(1 - Swi)R = Rs+(Bo/Bg)(krg/kro) ( o/ g)Rp = Gp/Np = R Np/ Npkr: PERMEABILIDADES RELATIVAS
BALANCE DE MATERIALESYACIMIENTOS DE GAS
EGBM COMO UNA LINEA RECTA ESF = G( Eg + Efw) + WeF: FLUIDOS PRODUCIDOS, BYF= Gwgp Bg + Wp BwGwp : GAS HUMEDO ACUM. PROD,
PCN = Gp + Npc FcGp: GAS SECO ACUMULADO PROD.
BALANCE DE MATERIALESYACIMIENTOS DE GAS
Npc: PROD ACUM DE CONDENSADOS, BNFc : FACTOR CONV CONDENS. PCN/BNFc = 132.79 c /Mc c= GRAVEDAD ESP CONDENSADO( w=1)Mc = PESO MOLECUL DEL CONDENSADO
Mc = 6084/(°API - 5.9)G : GAS HUMEDO EN EL YAC., PCN
BALANCE DE MATERIALESYACIMIENTOS DE GAS
Eg = Bg - BgiEfw = EXPANSION DEL AGUA Y
REDUCCION DEL VP, BY/PCNEfw = Bgi Ce (Pi - P)Efw = Bgi(Cw Swi + Cf)/(1-Swi) PYAC. GAS EMP. POR AGOTAMIENTO, We
y Efw SON CERO, EBM GAS SECOP/Z = Pi/Zi(1 - Gp/G), P/Z vs Gp RECTA
BAL. DE MATERIALESEBM COMO LINEA RECTA
YACIMIENTOS DE PETROLEOF = N(Eo + mEg + Efw) + WeF = N Et + We, We =0, GAS EN SOLUCIONF = N Et, F vs Et. N = PENDIENTEYACIMIENTO CON CAPA DE GAS, m 0SE SUPONE m y SE CALCULA F vs EtPENDIENTE AUMENTA, m PEQUEÑO,PENDIENTE DISMINUYE, m GRANDE
BAL. DE MATERIALESEBM COMO LINEA RECTA
YACIM. CON CAPA DE GAS, We = 0F/Eo = N + m N (Eg/Eo), F/Eo vs Eg/EoN = INTERCEPTO, mN = PENDIENTEMETODO DE HAVLENA Y ODEH YAC.
CON EMPUJE POR AGUA.F/Et = N + We/Et = N + C f(P,t)/EtREPRESENTAR F/Et vs f(P,t)/EtN = INTERCEPTO, C = PENDIENTE
BALANCE DE MATERIALESEBM - LINEA RECTA
We REQUIERE ANALISIS ESPECIALGRAFICO F/Et vs We/Et, LINEA RECTA, We
OKPENDIENTE DISMINUYE, We GRANDEPENDIENTE AUMENTA, We PEQUEÑOPENDIENTE CAMBIA DE DIRECCION, LA
GEOMETRIA ES INCORRECTA: LINEAL, RADIAL, ANGULAR, OTRA
BALANCE DE MATERIALESEBM LINEA RECTA
METODO DE CAMPBELL-INT. AGUAF/Et vs F, N = INTERCEPTOHORIZONTAL, We = 0, INCLIN. We 0METODO DE AJUSTE DE PRESIONESP vs Np, N, m y We, SE OBTIENEN POR
MINIMOS CUADRADOS U OTROS.EL PROGRAMA MBAL INCLUYE ESTOS
ANALISIS.
BALANCE DE MATERIALESEBM - LINEA RECTA
YACIMIENTOS DE GASECUACION GENERAL, F = G Et + WeEt = Eg + EfwYACIMIENTOS SIN We, AGOTAMIENTOF = G Et, F vs Et, G = PENDIENTEYACIMIENTO DE GAS CON We, Efw = 0F/Eg = G + We /Eg, G = INTERCEPTOPENDIENTE = C, INTRUSION - AGUA
BALANCE DE MATERIALESEBM -LINEA RECTA
HAVLENA Y ODEH, YACIMIENTO DE GAS CON We
F/Et = G + C f(P,t)/Et, F/Et vs f(P,t)/EtG = INTERCEPTO, C = PENDIENTEMETODO DE COLEF/Et vs F, G= INTERCEPTO,
HORIZONTAL We=0, INCLIN. We 0
BALANCE DE MATERIALESEBM LINEA RECTA
YACIMIENTOS DE GASMETODO DE AJUSTE DE P vs GpAJUSTAR PARA OBTENER G Y We CON
LA MININA DESVIACIONP/Z PARA CUALQUIER YACIMIENTOP/Z = (1- Gp/G)Pi/Zi, EBMLR, PEND=1/GRECTA AGOTAMIENTO, PENDIENTE
AUMENTA We PENDIENTE VARIA, P
ANALISIS DE LA EGBM
LA EGBM SE PUEDE ESCRIBIR COMO SIGUE (N(Eo + Boi Ce P) + mN(Eg + CeBoi P)+We)/F = 1, LUEGO LA EGBM SE PUEDE SEPARAR EN VARIOS
TERMINOS Y DEFINIENDO LOS INDICES DE EMPUJES POR GAS EN SOLUCION, CAPA DE GAS E HIDRAULICO COMO IGS, ICG Y IEH
IGS = N(Eo + Ce Boi P)/F, IEH = We/F ICG = mN(Eg + Boi Ce P) RESULTA QUE IGS + ICG+ IEH = 1 PUEDEN REPRESENTARSE GRAFICAMENTE
INTRUSION DE AGUA- METODOS DE CALCULO
VAN EVERDIGEN Y HURST 1949SOL ECUACION DIFUSIVIDAD, Plim CTEHISTORIA DE PRESION SE SUSTITUYE POR
ESCALONADA DISCRETAUSA SUPERPOSICION- SUMACARTER-TRACY Y FETKOVICH
SIMPLIFICAN Y FACILITAN CALCULO
CR-TR, USA TASA TERMINOS CTE
INTRUSION DE AGUA- METODOS DE CARTER-TRACY
Wej = Wej-1+((C PDj - Wej-1 P´j)/ Pj - tDj-1 P´j)) (tDj - tDj-1)
We: INT. AGUA ACUMULADA, BY C = 1.119 f h c , BY/psi, ro, RADIO DEL YACIMIENTO, pies, c COMPRES= Cw + Cf, f= ANG INT
tD : TIEMPO ADIMENSIONAL, tD = 0.00634 kt/( c ), t DIAS
Pj =ao+ a1 Td +a2 LTd+a3LtDE(2)
2ro
ro2
ro2
INTRUSION DE AGUA- METODOS DE CARTER-TRACY
LAS CONSTANTES HAN SIDO CALCULADAS POR FANCHI PARA DIF. rD = re/ro, re: R ACUIFERO
SIMILAR A VE&H- MAS FACILMODELO QUE PUEDE USARSE CON
HAVLENA Y ODEH.f(P,t) DE EGBM SE CALCULA CON
ECUACIONES DE CARTER Y TRACY.
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
OOO
OOOO
OOO
to
dtto
tdtt
tdtt
B
1)()(
B
1)(
B
1
B
1)()(
B
1)(
B
1)(
B
1)()(
ASI ONCONSERVACI DE ECUACIONLA RESULTA d(OIP) OEXPANDIEND
B
1
B
1
B
1 ),(SS ¨),(VPVP
SIGUIENTES ESDEFINICION LAS USANDO
CALCULAN SE PRODUCCION DE MECANISMOS LOS POR RECOBRO EL
1VP
1VP d(OIP)
PETROLEO DEL ONCONSERVACI DE ECUACION
dSdPVddSdPVdSPV
SdPVdSdPVdSPVdSPVdOIPd
dSdVPd
BS
BS
oott
ot
t
o
t
ott
o
t
to
tdtt
o
t
o
to
dtt
o
dtto
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
ECUACIONES DE LOS DIF MECANISMOS
0DE Y DE,DGDG
RESULTA DG,A SUMA SE DE BURBUJEO, DE PUNTO DEL DEBAJO POR
1)(VPDG GAS, POR EMPUJE
1)(VPDW AGUA, POR EMPUJE
RESULTA DS, EN SUST LUEGO ),()()d(S ,
)(1
DE PETROLEO, DEL EXPANSION
1)(-VPDS ,SATURACION DE CAMBIO
1S d(VP)- DC ROCA,LA DE ONCOMPACTACI
dtt
dtt
o
dtt
to
t
og
t
ow
o
dtto
t
t
oo
dtt
o
BSd
BSd
SgdSwdADEMAS
VPdB
dSVP
BSd
B
MECANISMOS DE PRODUCCION POR SIMULACION DE YACIMIENTOS
EN EL CASO DE UNA CELDA SATURADA AL REDUCIR LA PRESION EL TERMINO (DE) ES NEGATIVO PUESTO QUE 1/Bo AUMENTA, LUEGO SE ACUMULA CON DG, POR CONSIGUIENTE
DG = DG - DE, Y DE = 0 LA PRODUCCION DE PETROLEO POR LOS
DIFERENTES MECANISMOS(DC,DW,DG,DE) SE SUMAN EN TODOS LOS BLOQUES E INTERVALOS DE TIEMPO.
EN ALGUNOS CASOS NO SE APLICA.
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
EL SIMULADOR ECLIPSE Y OTROS, LA COMPACTACION LA MODELAN COMO COMPRESIBILIDAD, DATOS TABULADOS EN FUNCION DE PRESION QUE PUEDE SER REVERSIBLES E IRREVERSIBLES, SE PUEDE INCLUIR LA POSIBILIDAD DE MODIFICAR LA TRANSMISIBILIDAD EN f(P), UN MODELO DE HISTERESIS.
CUANDO LA P DEL BLOQUE AUMENTA LA COMPACTACION PUEDE SER REVERSIBLE E IRREVERSIBLE
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
LA COMPACTACION MODIFICA LA TRANS. POR EL CAMBIO DEL VOLUMEN POROSO
LA PRESION DE SOBRECARGA MENOS LA PRESION DE FLUIDOS ES EL ESFUERZO EFECTIVO EN FUNCION DEL CUAL SE REPRESENTA LA COMPACTACION.
ALGUNAS ROCAS COMO LAS YESO SE PRODUCE UNA COMPACTACION INDUCIDA CUANDO SE CONTACTAN CON AGUA QUE PUEDEN SER SUMINISTRADAS POR TABLAS EN FUNCION DE P Y Sw
MECANISMOS DE PRODUCCIONEMPUJE POR GAS EN SOLUCION
EL PROCEDIMIENTO PRECEDENTE NO DIFERENCIA ENTRE LA PRODUCCION POR GAS LIBRE Y GAS EN SOLUCION
LA PRODUCCION PUEDE SUBDIVIDIRSE ENTRE LO CORRESPONDIENTE A GAS LIBRE Y GAS EN SOLUCION.
PUEDE UTILIZARSE UN INDICE Fs = 1 PARA GAS DISUELTO Y 0 PARA GAS LIBRE
LIBRE GAS POR PRODUCCION ,/)()1(
SOLUCION EN NTEORIGINALME Sg DE FRACCIONFs ,/)(**
BoSgdFsVPDF
BoSgdFsVPDS
MECANISMOS DE PRODUCCION - EMPUJE HIDRAULICO
UN PROCEDIMIENTO SIMILAR AL CASO PREVIO SE UTILIZA PARA DISTINGUIR ENTRE EL EMPUJE HIDRAULICO Y EL AGUA INYECTADA. ORIGINALMENTE Ft SE CONSIDERA 1.0, CERO AGUA INYECTADA.
ASI SE PUEDEN CALCULAR DWT Y DWR, LAS PRODUCCIONES POR We Y Wi
INICIALAGUA AL ECORRESPOND QUE FRACCIONLA ES
/)(*)1(*
/)(**
Ft
BoSwdFtVPDWR
BoSwdFtVPDWT
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
OOO
OOOO
OOO
to
dtto
tdtt
tdtt
B
1)()(
B
1)(
B
1
B
1)()(
B
1)(
B
1)(
B
1)()(
ASI ONCONSERVACI DE ECUACIONLA RESULTA d(OIP) OEXPANDIEND
B
1
B
1
B
1 ),(SS ¨),(VPVP
SIGUIENTES ESDEFINICION LAS USANDO
CALCULAN SE PRODUCCION DE MECANISMOS LOS POR RECOBRO EL
1VP
1VP d(OIP)
PETROLEO DEL ONCONSERVACI DE ECUACION
dSdPVddSdPVdSPV
SdPVdSdPVdSPVdSPVdOIPd
dSdVPd
BS
BS
oott
ot
t
o
t
ott
o
t
to
tdtt
o
t
o
to
dtt
o
dtto
DETALLES DEL CALCULO DE MECANISMOS RECOBRO DE PETROLEO
ECUACIONES DE LOS DIF MECANISMOS
0DE Y DE,DGDG
RESULTA DG,A SUMA SE DE BURBUJEO, DE PUNTO DEL DEBAJO POR
1)(VPDG GAS, POR EMPUJE
1)(VPDW AGUA, POR EMPUJE
RESULTA DS, EN SUST LUEGO ),()()d(S ,
)(1
DE PETROLEO, DEL EXPANSION
1)(-VPDS ,SATURACION DE CAMBIO
1S d(VP)- DC ROCA,LA DE ONCOMPACTACI
dtt
dtt
o
dtt
to
t
og
t
ow
o
dtto
t
t
oo
dtt
o
BSd
BSd
SgdSwdADEMAS
VPdB
dSVP
BSd
B
MECANISMOS DE PRODUCCION POR SIMULACION DE YACIMIENTOS
EN EL CASO DE UNA CELDA SATURADA AL REDUCIR LA PRESION EL TERMINO (DE) ES NEGATIVO PUESTO QUE 1/Bo AUMENTA, LUEGO SE ACUMULA CON DG, POR CONSIGUIENTE
DG = DG - DE, Y DE = 0 LA PRODUCCION DE PETROLEO POR LOS
DIFERENTES MECANISMOS(DC,DW,DG,DE) SE SUMAN EN TODOS LOS BLOQUES E INTERVALOS DE TIEMPO.
EN ALGUNOS CASOS NO SE APLICA.
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
EL SIMULADOR ECLIPSE Y OTROS, LA COMPACTACION LA MODELAN COMO COMPRESIBILIDAD, DATOS TABULADOS EN FUNCION DE PRESION QUE PUEDE SER REVERSIBLES E IRREVERSIBLES, SE PUEDE INCLUIR LA POSIBILIDAD DE MODIFICAR LA TRANSMISIBILIDAD EN f(P), UN MODELO DE HISTERESIS.
CUANDO LA P DEL BLOQUE AUMENTA LA COMPACTACION PUEDE SER REVERSIBLE E IRREVERSIBLE
MECANISMOS DE PRODUCCION COMPACTACION DE LA ROCA
LA COMPACTACION MODIFICA LA TRANS. POR EL CAMBIO DEL VOLUMEN POROSO
LA PRESION DE SOBRECARGA MENOS LA PRESION DE FLUIDOS ES EL ESFUERZO EFECTIVO EN FUNCION DEL CUAL SE REPRESENTA LA COMPACTACION.
ALGUNAS ROCAS COMO LAS YESO SE PRODUCE UNA COMPACTACION INDUCIDA CUANDO SE CONTACTAN CON AGUA QUE PUEDEN SER SUMINISTRADAS POR TABLAS EN FUNCION DE P Y Sw
MECANISMOS DE PRODUCCIONEMPUJE POR GAS EN SOLUCION
EL PROCEDIMIENTO PRECEDENTE NO DIFERENCIA ENTRE LA PRODUCCION POR GAS LIBRE Y GAS EN SOLUCION
LA PRODUCCION PUEDE SUBDIVIDIRSE ENTRE LO CORRESPONDIENTE A GAS LIBRE Y GAS EN SOLUCION.
PUEDE UTILIZARSE UN INDICE Fs = 1 PARA GAS DISUELTO Y 0 PARA GAS LIBRE
LIBRE GAS POR PRODUCCION ,/)()1(
SOLUCION EN NTEORIGINALME Sg DE FRACCIONFs ,/)(**
BoSgdFsVPDF
BoSgdFsVPDS
MECANISMOS DE PRODUCCION - EMPUJE HIDRAULICO
UN PROCEDIMIENTO SIMILAR AL CASO PREVIO SE UTILIZA PARA DISTINGUIR ENTRE EL EMPUJE HIDRAULICO Y EL AGUA INYECTADA. ORIGINALMENTE Ft SE CONSIDERA 1.0, CERO AGUA INYECTADA.
ASI SE PUEDEN CALCULAR DWT Y DWR, LAS PRODUCCIONES POR We Y Wi
INICIALAGUA AL ECORRESPOND QUE FRACCIONLA ES
/)(*)1(*
/)(**
Ft
BoSwdFtVPDWR
BoSwdFtVPDWT
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
LA MALLA SE EXTIENDE AL ACUIFERO USANDO BLOQUES DIFERENTES CON MAYOR DEFINICION EN EL YACIMIENTO
REQUIERE + MEMORIA DEL COMPUTADOR PUEDE REPRESENTARSE EL ACUIFERO
MEDIANTE EL TERMINO PRODUCC/INYEC. SE UTILIZAN LAS DIVERSAS FORMAS DE
MANEJAR ACUIFERO PARA FLUJO CONT Y NO CONTINUO: HURST, VE&H, FEKOVITCH
UN EJEMPLO SE MUESTRA COMO SIGUE
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
TASA DE FLUJO DE AGUA EN BLOQUE m
lpca-Bw/día agua, deintrusión de constante C
acuífero el a referido bloque del área
contornoen prom.presión P /
)(/)()(
m
,
,11
mwtwmw
nmw
nwm
w
nenemw
PBCq
PPBCtB
tWtWq
m
m
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
EN ECLIPSE EL ACUIFERO SE MODELA POR UNA FILA DE CELDAS CONECTADAS AL YACIMIENTO
LA PALABRA CLAVE AQUNUM EN LA SECCION GRID SE USA PARA DEFINIR LAS PROPIEDADES COMO LONGITUD, AREA SECCIONAL, POROSIDAD, PERMEABILIDAD PRESION INICIAL PROFUNDIDAD, PVT Y Kr
EL ACUIFERO DEBE ESTAR CONECTADO A LA CARA DEL YACIMIENTO CON AQUCON
LAS CELDAS DEL ACUIFERO ESTAN CONECTADAS LAS CELDAS DEL ACUIFERO ESTAN AISLADAS DEL
YACIMIENTO
SIMULACION DE YACIMIENTOS-MANEJO DE ACUIFEROS
LA TRANSMISIBILIDAD ENTRE CELDAS DEL ACUIFERO VIENE DADO POR
ESTA EXPRESION SE USA TANTO PARA GEOMETRIAS CARTESIANAS Y CILINDRICAS
SE PRESENTAN LAS APROXIMACIONES DE FETKOVICH Y CARTER-TRACY
LONGi
XSECTi*PERMXi*2Ti DONDE
/1/1
TjTi
CDARCYTR
MANEJO DE ACUIFEROSFETKOVICH
EL MODELO DE FETKOVICH USA UNA APROXIMACION SIMPLIFICADA DE IP Y BALANCE DE MATERIALES ENTRE P Y We PARA EL ACUIFERO
LA INTRUSION DE AGUA DEL ACUIFERO SE MODELA POR LA ECUACION
acuífero elcon comunicado bloque del área :
acuífero-bloque del cara la de area el es m donde
)(
i
i
ii
iii
aiiai
A
Am
Am
donde
ddgppJQai
MANEJO DE ACUIFEROSFETKOVICH
LA PRESION DEL ACUIFERO SE OBTIENE POR BALANCE DE MATERIALES.
EL COMPORTAMIENTO DEL ACUIFERO DEPENDE DE DOS PARAMETROS, LA CONSTANTE DE TIEMPO DEL ACUIFERO Y EL INDICE DE PRODUCTIVIDAD
)( aaowota ppVCW
J
CtVwoTc
MANEJO DE ACUIFEROSFETKOVICH
BAJO LA SUPOSICION QUE LA PRESION DEL YACIMIENTO ES UNIFORME EN EL BLOQUE QUE CONECTA AL ACUIFERO, LA TASA DE INTRUSION DE AGUA PROMEDIO EN EL INTERVALO DE TIEMPO DELTA t ES
LA INTRUSION DE AGUA ACUMULADA SE CALCULA A CADA INTERVALO DE TIEMPO CUANDO SE ACTUALIZA LA PRESION, P
Tct
TctddgppJQ aiiaiai /
)/exp(1))((
MANEJO DE ACUIFEROSCARTER - TRACY
METODO SIMPLIFICADO QUE EVITA LA SUPERPOSICION, UTILIZANDO TABLAS Y SE APLICA A YACIMIENTOS DE FORMA ARBITRARIA
LOS PARAMETROS BASICOS SON Td, CONSTANTES TIEMPO E INTRUSION DE AGUA
cot
a
otwc T
trCh
Ck
rCT D
22
1
2
t, cy
MANEJO DE ACUIFEROSCARTER - TRACY
EL MODEL DE CARTER TRACY EXPRESA LA CAIDA DE PRESION EN EL CONTORNO EN TERMINOS DE PRESION ADIM, PID
LOS TERMINOS PDI Y SUS DERIVADAS SE OBTIENEN DE TABLAS SIMILARES A LAS CONOCIDAS DE VAN EVERDIGEN Y HURST
We SE CALCULA EXPLICITO AL FINAL Dt
conocidos términosdecalculan se by a donde )()(b(p-aQ
i bloque al acuífero del Qa flujo de tasala donde )(
iiai
I
tptt
tPIQ
pp
i
DDa
ao
MANEJO DE ACUIFEROSFLUJO CONSTANTE
EL FLUJO DEL ACUIFERO SE CONSIDERA CONSTANTE Y ES SUMINISTRADO AL SIMULADOR
LA TASA DE FLUJO EN EL BLOQUE DE UN ACUIFERO DE FLUJO CONSTANTE ES
LA CONSTANTE DEL ACUIFERO SE SUMINISTRA A CADA TIEMPO Y PUEDE VARIAR CON EL TIEMPO.
dato como acuífero del flujo Fa donde iiaai mAFQ