caracterización petrofísica del yacimiento bach-01,...

47
102 ( ) 384 , 0 * 009 , 0 384 , 0 * 0314 , 0 * 009 , 0 + = + = RQI PHI Ka S or (39) Donde: RQI : Índice de calidad de roca en función de permeabilidad (md) y porosidad (fracción). Figura 73. Gráfico y correlación de saturación de petróleo residual “Sor” Vs Índice de calidad de roca “RQI” con la data de núcleo del pozo LL-3689 y LS-5169 en las arenas de edad mioceno del yacimiento Bach-01. A continuación se muestra un resumen de los valores resaltantes para la saturación de los fluidos (Swir y Sor) para el tipo de roca en el yacimiento Bach-01 (Ver sección 4.3.4, tabla 7). Tabla 6. Rangos de valores de la saturación de fluidos en función de la data de núcleo (LL-3689 y LS-5169) de los puntos de “Swir” y “Sor”, Yacimiento Bach-01. Saturación de petroelo residual "Sor" Vs Indice de Calidad de roca "RQI" Yacimiento Bach-01 y = -0,009x + 0,384 R 2 = 0,134 0,00 0,05 0,10 0,15 0,20 0,25 0,30 0,35 0,40 0,45 0,50 0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50 5,00 RQI (micrones) Sor (fracción) Krwo_LL3689 Krw_LS5169 Lineal (Krwo_all) Tipo de Roca Rangos Porosidad Rangos Permeabilidad Swir Sor cod % md % % 1 34 -_40 2400-7300 10 -_15 33 -_34 2 30 -_34 1400-2400 15 -_24 34 -_36 3 25 -_30 700-1400 24 -_30 36 -_37 4 20 -_25 50-700 30-_35 37 -_39 YACIMIENTO FORMACIÓN MIEMBRO Bach-01 Lagunillas Bachaquero

Upload: dodien

Post on 24-Sep-2018

234 views

Category:

Documents


5 download

TRANSCRIPT

Page 1: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

102

( ) 384,0*009,0384,0*0314,0*009,0 +−=+

−= RQIPHIKaSor (39)

Donde:

RQI : Índice de calidad de roca en función de permeabilidad (md) y porosidad

(fracción).

Figura 73. Gráfico y correlación de saturación de petróleo residual “Sor” Vs Índice de calidad de roca “RQI” con la data de núcleo del pozo LL-3689 y LS-5169 en las arenas de edad mioceno del yacimiento Bach-01.

A continuación se muestra un resumen de los valores resaltantes para la

saturación de los fluidos (Swir y Sor) para el tipo de roca en el yacimiento Bach-01

(Ver sección 4.3.4, tabla 7).

Tabla 6. Rangos de valores de la saturación de fluidos en función de la data de núcleo (LL-3689 y LS-5169) de los puntos de “Swir” y “Sor”, Yacimiento Bach-01.

Saturación de petroelo residual "Sor" Vs Indice de Calidad de roca "RQI" Yacimiento Bach-01

y = -0,009x + 0,384R2 = 0,134

0,00

0,05

0,10

0,15

0,20

0,25

0,30

0,35

0,40

0,45

0,50

0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50 5,00

RQI (micrones)

Sor

(fra

cció

n)

Krwo_LL3689Krw_LS5169Lineal (Krwo_all)

Tipo de Roca Rangos

PorosidadRangos

PermeabilidadSwir Sor

cod % md % %

1 34 -_40 2400-7300 10 -_15 33 -_34

2 30 -_34 1400-2400 15 -_24 34 -_36

3 25 -_30 700-1400 24 -_30 36 -_37

4 20 -_25 50-700 30-_35 37 -_39

YACIMIENTOFORMACIÓN

MIEMBRO

Bach-01Lagunillas

Bachaquero

Page 2: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

103

La distribución de fluidos original en las arenas de edad mioceno del yacimiento

Bach-01 han variado considerablemente en algunas áreas o parcelas, debido al

histórico de explotación de producción desde el año 1934, lo cual se observa

durante el monitoreo de los pozos el control de la dirección y avance pasivo tanto

del agua original del acuífero como la inyectada mediante el IAV (Inyección

alternada de Vapor) a partir de 1971. Finalmente, se podría utilizar la saturación de

agua calculada por perfiles para determinar el petróleo original en sitio “POES”,

considerando estimar la saturación de agua inicial sobre el contacto agua petróleo

original “CAPO”, y a partir del cálculo de la saturación de agua irreducible “Swir” (EC

38), para así validar que la saturación de agua estimada por registro no presente un

avance o barrido considerable, sin embargo; el acuífero hacia la parte baja de la

estructura (Sur) del yacimiento no es activo, lo cual presenta un barrido pasivo.

4.3.3.3 Determinación del modelo matemático de saturación de agua a partir de

modelos conocidos y validado con data de núcleo.

La saturación de agua de las areniscas de edad mioceno del yacimiento Bach-01,

se calculó para los 904 del estudio, se selecciono del modelo matemático de

saturación de agua que ajuste con la saturaciones del reservorio en las arenas

arcillosa del yacimiento Bach-01, se probaron varios métodos conocidos,

determinándose y validándose el modelo matemático de “Simandoux” mediante los

pozos claves, ya que el mismo representa las característica de saturación del agua

del yacimiento Bach-01 para arenas arcillosas. Se obtuvo un buen cotejo con la data

de núcleos (LL-3689 y LS-5169) para el modelo matemático de saturación de agua

“Simandoux”, mediante los análisis convencionales (Dean Stark) y los análisis

especiales de presión capilar y permeabilidad relativa.

A continuación se describe en detalle cada modelo matemático de saturación de

agua conocido, con su respectiva justificación y ecuación para cada tipo de litología

con el tipo de arcilla predominante y distribución (dispersa, laminar y estructural).

Modelos matemáticos de saturación de agua según literatura; Archie “SW_archie”,

Indonesia (Poupon) “SW_indonesia”, Waxman&Smits “SW_waxman&smits”, y

Simandoux “SW_SIMe” (calculo a partir PHIE) o “SW_SIMt” (calculo a partir PHIT),

Page 3: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

104

se efectuaron comparaciones, seleccionando así el modelo matemático de

saturación de agua “Simandoux”, ya que ajusta con los datos de saturación de agua

de los análisis convencionales (SW_DnStrk), análisis especiales de presión capilar

(SW_Pc) y permeabilidad relativa agua petróleo (SW_Krel) de los núcleos cortados

en los pozos LL-3689 y LS-5169. Ver figura 73.

Modelo matemático de saturación de agua “Archie”: En el año 1942 archie

deriva este modelo matemático para calcular la saturación de agua en

areniscas limpias, donde no es necesario efectuar corrección por arcillosidad,

por lo cual este modelo no honra las características de las arenas arcillosas

del yacimiento Bach-01. La ecuación básica de este modelo matemático esta

representada de la siguiente manera:

)*

*(/1

_

RPHITRw

Tm

na

S archieW = (40)

Donde:

a: es la constante de Archie o factor de tortuosidad.

m: es el exponente de cementación.

n: es el exponente de saturación depende de la humectabilidad.

PHIT: representa la porosidad total (fracción).

RT (RD_0): es la resistividad verdadera de la formación (Ohm-m).

Rw: es la resistividad del agua de formación a la temperatura del

Yacimiento (Ohm-m).

Modelo matemático de saturación de agua “Indonesia”: En el año 1954 el

modelo matemático de Indonesia creado por Poupon, fue diseñado para

calcular y corregir los valores de saturación de agua en areniscas arcillosas,

pero con distribución de arcilla laminar, donde la presencia de componentes

arcillosos dispersos en el medio poroso incrementan la conductividad o reducen

la resistividad verdadera de la formación. Este modelo no representa la

saturación de agua en las arenas arcillosas con un tipo de arcilla dispersa

Page 4: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

105

predominante, lo cual si lo hace el modelo de Simandoux. La ecuación que

representa al modelo de Indonesia tiene la siguiente forma:

( )

−= − VCL

RshVCL

RTRwPHIEaindonesiaSW m 1/1*)**(_ (41)

Donde:

a: es la constante de Archie o factor de tortuosidad.

m: es el exponente de cementación.

n: es el exponente de saturación depende de la humectabilidad.

PHIE: representa la porosidad efectiva (fracción).

RT (RD_0): es la resistividad verdadera de la formación (Ohm-m).

Rsh: es la resistividad de la lutita (Ohm-m).

Rw: es la resistividad del agua de formación a la temperatura del

Yacimiento (Ohm-m).

VCL: Volumen de arcilla del área evaluada (fracción).

Modelo matemático de saturación de agua “Waxman&Smits”: En el año 1968,

Waxman y Smits publicaron su renombrado trabajo el cual más tarde fue

conocido con el nombre de ecuación de Waxman&Smits. Desde ese entonces la

ecuación fue modificada por Waxman y Thomas (1974) y más tarde por Juhasz

(1981). Este fue y sigue siendo catalogado como un refinado método para

calcular la saturación de agua en arenas sucias a partir de información

aportada por los registros de resistividad, sin embargo; este modelo no

representa las saturaciones reales de agua del yacimiento Bach-01, y con muy

poca data de análisis de intercambio cationico para corroborar el modelo. La

ecuación se expresa de la manera siguiente:

+=

nm

SwQvBRw

aRwPHIERTwSW

1

**1*

* smits&axman_

(42)

Page 5: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

106

Donde:

Sw: Saturación de agua en la zona virgen, (fracción)

RT (RD_0): Resistividad verdadera en la zona virgen, (ohm-m)

Rw: Resistividad del agua de formación, (ohm-m)

PHIE: porosidad efectiva, (fracción)

a: es la constante de Archie o factor de tortuosidad.

m: es el exponente de cementación.

n: es el exponente de saturación depende de la humectabilidad.

Qv: Capacidad de intercambio catiónico en la zona virgen, (meq/cc)

B: Llamado conductividad específica de los cationes y es un índice de la movilidad

de los cationes absorbidos sobre la superficie arcillosa, dada por la siguiente

ecuación:

( )27.0*045.0128.1*000406.0*225.0

23.1

2

−+−−

=TformRw

TformTformB (43)

Donde:

formT : Temperatura calculada de la formación de interés. (ºF)

Rw: Resistividad del agua de formación, ohm-m.

Modelo matemático de saturación de agua “Simandoux”: En 1963 Simanduox

reportó experimentos en mezclas homogéneas de arena y arcilla, para calcular

y corregir los valores de saturación de agua en areniscas arcillosas, donde la

presencia de arcillas dispersa el medio poroso incrementan la conductividad y

reduce la resistividad verdadera de la formación, lo cual es el caso de las

arenas arcillosas, con arcillas dispersas del yacimiento Bach-01, este modelo

además de honrar las características litológicas y de fluido del yacimiento,

ajusta muy bien con las saturaciones de los análisis de los pozos núcleos LL-

3689 y LS-5169, como se muestra en la figura 74. A continuación la ecuación

del modelo matemático de saturación de agua “Simandoux”, “SW_SIMe” o

“SW_PHIE” del yacimiento Bach-01, expresada de la manera siguiente:

Page 6: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

107

+=

RshPHIEVCLRwa

RshPHIEVCLRwa

RTPHIERwaSIMeSW m

n

mm **2**

**2**

**_

12

(44)

Donde:

RT (RD_0): Resistividad verdadera en la zona virgen, (ohm-m)

Rw: Resistividad del agua de formación, (ohm-m).

Rsh: es la resistividad de la lutita (Ohm-m).

PHIE: porosidad efectiva, fracción (fracción)

a: es la constante de Archie o factor de tortuosidad.

m: es el exponente de cementación.

n: es el exponente de saturación depende de la humectabilidad.

VCL: Volumen de arcilla del área evaluada (fracción).

Figura 74. Cotejo del modelo matemático de simandoux en los pozos núcleos LL-3689 y LS-5169 para las arenas de edad mioceno del yacimiento Bach-01. Se observa claramente un solapamiento de las curvas del modelo “SW_SIMe y SW_SIMt” en las zonas de arenas arcillosas (circulo azul) y limpias (circulo rojo) con las curvas de la data de núcleos “SW_Pc, SW_Krel y SW_DnStrK”.

0 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

0 SW_Krel 10 SW_Pc 1

LL-36890 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

LS-5169

0 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

0 SW_Krel 10 SW_Pc 1

LL-36890 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

0 SW_Krel 10 SW_Pc 1

LL-36890 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

LS-51690 VCL_LRJ 1

0 GRN1 1500 RT 1000

0 RS 1000 0 SW_SIMe 1

0 SW_SIMt 10 SW_DnStrk 1

LS-5169

Page 7: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

108

4.3.4 Tipos de rocas mediante propiedades petrofísicas, facies, petrofacies y

unidades de flujo, litofacies, descripción de núcleos y ambiente sedimentario.

Los rangos de las propiedades petrofísicas de data de núcleo, petrofacies y

unidades de flujo, facies sedimentarias o depofacies, litofacies, descripción del

núcleo y ambiente sedimentario son vitales para el desarrollo y calibración del

modelo de tipos de roca y así generar una curva “TIPO_ROCA”, que puede derivarse

de la evaluación petrofísica de los perfiles de pozos. Se caracterizó siete litofacies

para la sección estratigráfica del yacimiento Bach-01, mediante la descripción de los

pozos núcleos LL-2318 y LS-5169, definidas mediante porcentaje de arena, eventos

depositacionles, estructuras sedimentarias y/o cementación. Estas siete litofacies

representan una progresión desde un yacimiento excelente (SS a SSa), a un

yacimiento de de menor calidad (SSva a SScg), hasta una roca no-yacimiento sin

prospectividad (Het, M y SH), ver tabla 7. El bajo espesor de las litofacies

representadas en escala de pulgadas y en algunos casos sobrepasa levemente la

escala de pies, lo cual se dificulta un poco la calibración de los tipos de roca,

derivándose mediante perfiles de pozos. Los informes reportados de los análisis

geológicos de los pozos núcleos LL-2318 y LS-5169 presenta una descripción

litológica detallada de cada una de las litofacies y ambientes sedimentarios

establecidos, Ver tabla 7. Los mapas de facies proporcionan una síntesis de todos

los aspectos del modelo sedimentológico los cuales fueron divididos en cuatro

categorías de depofacies (I Canal, II Barra, III Abanico y IV Llanura), ver tabla 7.

Los rangos de las propiedades petrofísicas, análisis de los resultados de las pruebas

convencionales y especiales de núcleos (porosidad, permeabilidad absoluta/relativa

y presión capilar) de los pozos LL-3689 y LS-5169 permitió determinar los tipo de

petrofacies y unidades de flujo, comprobando así el grado de heterogeneidad de las

arenas de edad mioceno del yacimiento Bach-01, ver tabla 7.

Para desarrollar el modelo de tipos de roca, se codifico en formato digital para

cargarse en Openworks (Petroworks) las litofacies a cada una de las profundidades,

ya que el formato de las hojas sedimentologicas solo se disponía en Imagen Tiff

entregadas por las compañías de laboratorio. Luego se volvió a chequear la

correlación núcleo perfil hecha para los modelo matemáticos petrofísicos donde se

llevaron los datos de núcleo a profundidad de registro mediante el core gamma.

Page 8: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

109

También se corroboró el control de base de datos de las curvas crudas y

procesadas. A partir de lo anterior, se genero un análisis interpretativo confiable

para el tipo de roca con los pozos claves (LL-3689, LL-2318, LS-3803 y LS-5169).

Inicialmente a partir de las hojas sedimentologicas de los pozo claves (núcleos)

se clasificaron rangos de propiedades petrofísicas tanto de núcleo como registros

(VCL, PHI, Ka y NGT) para cada una de las siete Litofacies, identificando así tipo de

rocas que se caracterizan por sus propiedades petrofísicas mediante varios gráficos

cruzados “crossplot”, integrándose con una nueva descripción litológica que

involucra tipo de grano, afinamiento y engrosamiento, material arcilloso o no, y

calidad de yacimiento, además de integrar al tipo de roca las petrofacies (mega,

macro, meso, micro y nano), depofacies (I Canal, II Barra, III Abanico y IV

Llanura) de los mapas facies y el ambiente sedimentario (desde fluvial en su base

hasta deltaico en el tope), ver tabla 7. Según la metodología expuesta, se comenzó

analizando la relación de la porosidad efectiva “PHIE_D” en función de volumen de

arcilla “VCL_LRJ” de los pozos claves (núcleos), mediante el crossplot de la figura

75, se observa una buena correlación con muy poca dispersión en los valores de

“PHIE_D” con valores iguales de “VCL_LRJ”, por lo cual esto evidencia pocos

procesos diagenéticos que afectaron el yacimiento Bach-01, indicando que presenta

bajo grado de complejidad (heterogeneidad) y variabilidad entre los tipos de roca.

Figura 75. Gráfico de crossplot de la relación de porosidad efectiva “PHIE_D” Vs volumen de arcilla “VCL_LRJ”, con los pozos claves del yacimiento Bach-01.

VCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ

VCL_LRJVCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ

Page 9: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

110

Continuando con la metodología, se analizó la relación entre la porosidad efectiva

“PHIE_D” y volumen de arcilla “VCL_LRJ” en función de la resistividad profundad

“RD_0”, lográndose visualizar en el crossplot de la figura 76, una buena relación con

la calidad de la roca al incrementar los valores de “RD_0” (>20 ohm-m), con esta

interpretación se podría inferior zonas prospectivas y un parámetro de corte “Cut-

off” con la curva de RD_0. Esta interpretación aplica en areniscas saturadas por

petróleo y agua irreducible como es el caso del yacimiento Bach-01, por encima del

contacto agua petróleo (CAP).

Figura 76. Gráfico de crossplot de la relación de “PHIE_D” Vs “VCL_LRJ” en función de la resistividad profundad “RD_0”, pozos claves del yacimiento Bach-01.

Posteriormente se realizó el análisis mediante el gráfico crossplot con la relación

de porosidad efectiva “PHIE_D” y volumen de arcilla “VCL_LRJ” en función de la

permeabilidad absoluta del modelo “K_PHIe” y los rangos de propiedades

petrofísicas de las litofacies en los pozos claves (LL-2318 y LS-5169) que

representan una mejor calidad de roca (SS a SSa), a una menor calidad (SSva a

SScg) para la roca reservorio (Bach-01). Se puede observar una buena

correspondencia con las litofacies entre el incremento de “PHIE_D” y “K_PHIe” con

la disminución del “VCL_LRJ”. Esta integración permitió obtener zonas o valores de

propiedades petrofísicas para cada tipo de roca (ver tabla 7), representando la

buena calidad de la roca a regular para el yacimiento Bach-01, como puede

observarse en la figura 77.

Crossplot PHIE_D Vs VCL_LRJ en Función RD_0

VCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ en Función RD_0

VCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ en Función RD_0

VCL_LRJ

Page 10: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

111

Figura 77. Gráfico de crossplot de la relación de “PHIE_D” Vs “VCL_LRJ” en función de la permeabilidad absoluta “K_PHIe” y los rangos de propiedades petrofísicas de las litofacies con los pozos claves del yacimiento Bach-01.

Continuando con el análisis para determinar el tipo de roca con varias fuentes de

data y metodologías variadas, permitiendo así un tipo de roca integrado con mayor

certidumbre, por lo cual se determino las petrofacies (Radio de garganta poral “RI”

por Winland y Pitman) y unidades de flujo para cada tipo de roca “TIPO_ROCA” del

yacimiento Bach-01, utilizando los análisis de las pruebas especiales de núcleos

(Presión capilar por inyección de mercurio y unidades de flujo “FZI” e índice de

calidad de yacimiento “RQI”) de los pozos LL-3689 y LS-5169. En los gráficos de la

figura 78, se observa del lado izquierdo el perfil de garganta de poro, el mismo se

realiza con la finalidad de obtener el radio de garganta poral “RI” a diferente niveles

de saturación de mercurio, pero con la utilización de la metodología de Winland (El

sistema poroso que domina el flujo a través de la roca yacimiento a una saturación

de mercurio “Shg” de RI35%) se interpretaron 2 tipo de petrofacies de buena

calidad de roca (mega y macro). Mientras que se observa del lado derecho el gráfico

de Ápices (Apex Plots) en la figura 78, el mismo permite identificar rangos de

saturación de mercurio donde la roca alcanza el mayor desplazamiento de radio de

garganta poral dominante “RI”. Posteriormente en los gráficos de la figura 79, se

observa del lado izquierdo la saturación incremental de mercurio “Shg” permitiendo

observar en cual tipo de roca ocurre el mayor desplazamiento de la fase no mojante

(mercurio), obteniéndose en el yacimiento Bach-01 de mayor a menor calidad de

Crossplot PHIE_D Vs VCL_LRJ en Función Litofacies de Núcleo y K_PHIe

VCL_LRJ

SS

SSa

SSva

LITOFACIES NÚCLEO

SScg

Crossplot PHIE_D Vs VCL_LRJ en Función Litofacies de Núcleo y K_PHIe

VCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ en Función Litofacies de Núcleo y K_PHIe

VCL_LRJ

Crossplot PHIE_D Vs VCL_LRJ en Función Litofacies de Núcleo y K_PHIe

VCL_LRJ

SS

SSa

SSva

LITOFACIES NÚCLEO

SScg

SS

SSa

SSva

LITOFACIES NÚCLEO

SScg

Page 11: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

112

roca “3 petrofacies predominante”, Roca 1 “Mega” (RI: 17,75), Roca 2 “Macro” (RI:

8,87) y Roca 3 “Meso” (RI: 1,77). Continuando con la figura 79, se observa del lado

derecho el gráfico uno-uno “one-one” (con Línea 45º de correlación), el mismo

permite comparar los valores de “RI” obtenidos a partir del Perfil de Garganta de

Poro con los calculados usando las correlaciones empíricas de Winland y Pittman,

para cada nivel de saturación de mercurio perteneciente al rango previamente

establecido a través del gráfico de ápices. Se concluye que para las arenas de edad

mioceno del yacimiento Bach-01, en el sistema poroso la saturación que domina el

flujo a través de la roca y radio de garganta poral es “RI35”.

Figura 78. Gráfico del perfil de garganta de poro (del lado izquierdo) y gráfico de Ápice (del lado derecho), con los pozo claves del yacimiento Bach-01.

Figura 79. Gráfico de Shg incremental (del lado izquierdo) y gráfico uno-uno (del lado derecho), con los pozo claves del yacimiento Bach-01.

PERFIL DE GARGANTA POROInyección Mercurio Yac Bach 01

1

10

100

1000

10000

020406080100 % Volumen de Poro Ocupado

Pc (p

si)

Muestra_302_A Muestra_324_A Muestra_606_A Muestra_722_A Muestra_806_A Muestra_812_AMuestra_821_A Muestra_1013_A Muestra_1016_A Muestra_1215_A Muestra_1045_A

100

25

.025

.05

.01

.1

.5

.25

1

2.5

5

10

TIPO PORO

2

MA

CR

ON

AN

OM

ICR

OM

ESO

MEG

A

Shg :35 % Winland

Yac Bach 01 SHg/Pc vs SHg Grafico Apices

0.00

5.00

10.00

0 10 20 30 40 50 60 70 80 90 100

SHg

SHg/

Pc

302A_Macro

324A_Mega

606A_Macro

722A__Mega

806A_Mega

812A_Macro

821A_Macro

1013A_Macro

1016A_Mega

1215A_Mega

1405A_Mega

Shg :35 % Winland

PERFIL DE GARGANTA POROInyección Mercurio Yac Bach 01

1

10

100

1000

10000

020406080100 % Volumen de Poro Ocupado

Pc (p

si)

Muestra_302_A Muestra_324_A Muestra_606_A Muestra_722_A Muestra_806_A Muestra_812_AMuestra_821_A Muestra_1013_A Muestra_1016_A Muestra_1215_A Muestra_1045_A

100

25

.025

.05

.01

.1

.5

.25

1

2.5

5

10

TIPO PORO

2

MA

CR

ON

AN

OM

ICR

OM

ESO

MEG

A

Shg :35 % Winland

Yac Bach 01 SHg/Pc vs SHg Grafico Apices

0.00

5.00

10.00

0 10 20 30 40 50 60 70 80 90 100

SHg

SHg/

Pc

302A_Macro

324A_Mega

606A_Macro

722A__Mega

806A_Mega

812A_Macro

821A_Macro

1013A_Macro

1016A_Mega

1215A_Mega

1405A_Mega

Shg :35 % Winland

SHg Incremental Yac BACH 01 Vol hg vs Radio (micrones)

0

5

10

15

20

25

30

35

0.01 0.10 1.00 10.00 100.00

Radio de garganta (micras)

Volu

men

Hg

Incr

emen

tal

302A_Meso

324A_Mega

606A_Mega

722A_Mega

806A_Mega

812A_Mega

821A_Macro

1013A_Meso

1016A_Mega

1215A_Mega

1405A_Mega

NANO MICRO MESO MACRO MEGA2.000.5

Roca 3 Meso R: 1.77

Roca 2 Macro R: 8.87

Roca 1 Mega R: 17.75

R35_Pc vs. R35_Pitmann Yac Bach 01y = 1.0126xR2 = 0.6757

05

1015202530354045505560657075

0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75

R35_Pc

R35

_Pitm

ann

SHg Incremental Yac BACH 01 Vol hg vs Radio (micrones)

0

5

10

15

20

25

30

35

0.01 0.10 1.00 10.00 100.00

Radio de garganta (micras)

Volu

men

Hg

Incr

emen

tal

302A_Meso

324A_Mega

606A_Mega

722A_Mega

806A_Mega

812A_Mega

821A_Macro

1013A_Meso

1016A_Mega

1215A_Mega

1405A_Mega

NANO MICRO MESO MACRO MEGA2.000.5

Roca 3 Meso R: 1.77

Roca 2 Macro R: 8.87

Roca 1 Mega R: 17.75

R35_Pc vs. R35_Pitmann Yac Bach 01y = 1.0126xR2 = 0.6757

05

1015202530354045505560657075

0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75

R35_Pc

R35

_Pitm

ann

Page 12: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

113

Luego de obtener 3 tipo de rocas por los gráficos anteriores mediante la

metodología de petrofacies, se valido nuevamente con los gráficos de la figura 80,

determinándose los tipos de roca “RI Petrofacies” con sus unidades de flujo “FZI” e

índice de calidad de yacimiento o roca “RQI”. Se observa del lado izquierdo el

gráfico de la relación de permeabilidad Vs porosidad de núcleo en función de radio

de garganta poral “RI” (clasificándose según los colores, Mega RI>10 micrones,

Macro 2<RI<10 y Meso 0,5<RI<2) obteniéndose para el yacimiento “Bach-01” 3

tipo de roca de mayor a menor calidad (Mega, Macro y Meso). Posteriormente del

lado derecho de la figura 80, se realizo el gráfico Log-Log de “RQI” Vs porosidad

normalizada “PHIE_Z” (grupo de familia de porosidad), obteniéndose que para cada

tipo de roca caracterizada por petrofacies, se aprecia 5 tendencias de indicadores de

zonas de flujo “FZI” variando el índice de calidad de yacimiento “RQI” para Bach-01,

de una característica de roca de mayor calidad a menor calidad, clasificándose que

para la petrofacies “Mega” existe 3 unidades de flujo “UF” (ver gráfico, rojo FZI_1,

verde FZI_2, y azul FZI_3), “Macro” existe 1 “UF” (ver gráfico, amarillo FZI_4), y

por ultimo “Meso” existe 1 “UF” (ver gráfico, morado FZI_5).

Figura 80. Gráfico de permeabilidad Vs porosidad de núcleo (del lado izquierdo) y gráfico Índice de calidad de roca “RQI” y unidades de flujo “UF” (del lado derecho), con los pozo claves del yacimiento Bach-01.

En seguida se realizo el gráfico de Lorenz modificado (Gunter) y el diagrama de

almacenamiento y flujo estratigráfico (Miller), mediante estas nuevas metodologías

se identificaron y validaron 5 UF para el yacimiento Bach-01, detallándose mediante

Petrofacies K vs Phie Yac Bach 01 "FZI (Unidades de Flujo)"

0.001

0.01

0.1

1.

10.

100.

1000.

10000.

100000.

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45Porosidad

Perm

eabi

lidad

(md)

Mega_Fzi_1 Mega_Fzi_2 Mega_Fzi_3 Mega_Fzi_4

Macro_Fzi_4 Macro_Fzi_5 Meso_Fzi_5

Rad

io d

e G

arga

ntas

de

Poro

s (m

icro

nes)

MA

CRO

ME

GA

ME

SOM

ICR

ONA

NO

0.5

2

10

0.1

100

0.01

Indice de Calidad Yacimiemto (RQI) vs Grupos de Porosidad (Phie_z)Unidades de Flujo BACH 01

y = 16.891x - 5.828R2 = 0.87

y = 8.3131x - 1.6544R2 = 0.6792

y = 4.3156x - 0.8681R2 = 0.6605

y = 4.1969x - 1.2143R2 = 0.8561

y = 2.4922x + 0.5461R2 = 0.3325

0.0

0.1

1.0

10.0

100.0

0.01 0.10 1.00PHIE_Z

RQ

I,mic

rone

s

Mega_FZI _1 Mega_FZI _2 Mega_FZI _3 Mega_Macro_FZI _4 Macro_Meso_FZI _5

Petrofacies K vs Phie Yac Bach 01 "FZI (Unidades de Flujo)"

0.001

0.01

0.1

1.

10.

100.

1000.

10000.

100000.

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45Porosidad

Perm

eabi

lidad

(md)

Mega_Fzi_1 Mega_Fzi_2 Mega_Fzi_3 Mega_Fzi_4

Macro_Fzi_4 Macro_Fzi_5 Meso_Fzi_5

Rad

io d

e G

arga

ntas

de

Poro

s (m

icro

nes)

MA

CRO

ME

GA

ME

SOM

ICR

ONA

NO

0.5

2

10

0.1

100

0.01

Indice de Calidad Yacimiemto (RQI) vs Grupos de Porosidad (Phie_z)Unidades de Flujo BACH 01

y = 16.891x - 5.828R2 = 0.87

y = 8.3131x - 1.6544R2 = 0.6792

y = 4.3156x - 0.8681R2 = 0.6605

y = 4.1969x - 1.2143R2 = 0.8561

y = 2.4922x + 0.5461R2 = 0.3325

0.0

0.1

1.0

10.0

100.0

0.01 0.10 1.00PHIE_Z

RQ

I,mic

rone

s

Mega_FZI _1 Mega_FZI _2 Mega_FZI _3 Mega_Macro_FZI _4 Macro_Meso_FZI _5

Page 13: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

114

su distribución vertical que la UF 1 es la que presenta el mayor aporte de capacidad

de flujo y mejor característica de la mejor calidad de roca (Ver Figura 81).

Figura 81. Gráfico de Lorenz modificado (arriba), y gráfico de diagrama de almacenamiento y flujo estratigráfico (abajo), identificando 5 UF con el pozo núcleo LL-3689 para las arenas del mioceno del yacimiento Bach-01.

Luego de tener caracterizado los tipos de roca por petrofacies y unidades de flujo,

se integro esta data con el tipo de roca determinado por litofacies, propiedades

petrofísicas, facies, descripción de núcleos y ambiente sedimentario. (Ver tabla 7)

Page 14: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

115

Tabla 7. Caracterización del rango de las propiedades petrofísicas para determinar el tipo de roca integrado “TIPO DE ROCA” para el yacimiento Bach-01, mediante las propiedades petrofísicas, facies, petrofacies y unidades de flujo, litofacies, descripción de núcleos y ambiente sedimentario de la data de los pozos claves o núcleos.

Tipos de

Roca

Representación

gráfica Litologia

Rango de valores VCL [%] PHIE [%] Ka [md] NTG [%] RI35 (u)

RQI(adim) FZI (adim)

Tipo de roca

“Petrofacies”

(RI35), (UF)

Facies Sedimenta

ria (Depofacie

s en mapas)

“Litofacies”equivalentes

del estudi

o”

Descripción del núcleo

interpretada

Ambiente sediment

ario

1

VCL <= 10 PHIE >= 34 Ka > 2400 NTG>50 RI35>23 RQI>2,64 FZI>5,12

Mega FZI_1

Facies I Canal II Barra

SS

Arena de grano medio a fino, relativamente

libre de material arcilloso.

Excelente Roca Yacimiento.

Deltaico/ Fluvial

2

10< VCL <= 14 30 <= PHIE < 34

1400 < Ka <= 2400 NTG>50

18<RI35<=23 2,15<RQI<=2,64 5,01<FZI<=5,12

Mega FZI_2

Facies I Canal II Barra

SSa

Arena de grano fino a muy fino,

con poco material arcilloso

interlaminado y/o agregado. Buena Roca Yacimiento.

Deltaico/ Fluvial

3

14 < VCL <= 18 25 <= PHIE < 30 700 < Ka <= 1400

NTG>50 10<RI35<=18

1,66<RQI<=2,15 4,98<FZI<=5,01

Mega FZI_3

Facies I Canal II Barra

SSva

Arena de grano muy fino, con

abundante material arcilloso interlaminado y/o

agregado. Regular Roca Yacimiento.

Deltaico/ Fluvial

4

18 < VCL <=22 20 <= PHIE < 25 50 < Ka <= 700

5>NTG>50 2<RI35<=10

0,50<RQI<=1,66 1,99<FZI<=4,98

Macro FZI_4

Facies II Barra

III Abanico SScg

Arena conglomerática con clastos de

rocas. Regular a pobre Roca Yacimiento.

Deltaico/ Fluvial

5

0.5<RI35<=2 VCL <=40 Ka <=50

Meso FZI_5

Facies III Abanico

Het

Heterolítica con mezclas de arena,

limo y/o lodo. Roca no

Yacimiento.

Deltaico/ Fluvial

6

40 < VCL <= 70 Facies IV Llanura

M Arcillas

Deltaico/ Fluvial

7

VCL > 70 Facies IV Llanura

SH Lutitas, sellantes Deltaico/ Fluvial

Page 15: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

116

Una vez caracterizado el rango de las propiedades petrofísicas para cada uno de

los 7 tipo de roca integrado mediante la correlación núcleo-perfil, y a pesar de la

diferencia de escala entre la data de núcleos y perfiles de pozos, se pudo desarrollar

un modelo confiable de tipos de roca integrado para el yacimiento Bach-01 (Ver

Figura 82), a partir de las relaciones de los gráficos cruzados o crossplot entre

volumen de arcilla, porosidad efectiva y resistividad, facies, petrofacies, unidades de

flujo, litofacies, descripción de núcleos y ambiente sedimentario.

Figura 82. Modelo de tipo de roca integrado establecido para las arenas de edad mioceno formación lagunillas del yacimiento Bach-01.

El modelo de tipos de roca integrado, es un modelo que puede ser usado con

gran confiabilidad en los 904 pozos del estudio, utilizando las curvas procesadas de

la evaluación petrofísica. Este modelo presenta mayor confianza en los pozos con

datos de porosidad efectiva derivada de los registros de densidad y volumen de

arcilla, sin embargo; su incertidumbre se incrementa en los pozos donde la

porosidad total y efectiva se determine por algoritmos o correlación de propiedades.

Litofacies Núcleo

Petro- Facies

SS Mega

SSa Mega

SSva Mega

SScg Macro

Het Meso

M

SH

Crossplot PHIE_D Vs VCL_LRJ en Función Tipo de roca integrado

VCL_LRJ

Litofacies Núcleo

Petro- Facies

SS Mega

SSa Mega

SSva Mega

SScg Macro

Het Meso

M

SH

Crossplot PHIE_D Vs VCL_LRJ en Función Tipo de roca integrado

VCL_LRJ

Page 16: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

117

4.3.5 Modelo de Permeabilidad

En las arenas edad mioceno del miembro Bachaquero, las tendencias de

permeabilidad con la porosidad de núcleos y modelos matemáticos muestran una

clara relación con el volumen de arcilla y tipo de roca integrado, obteniéndose un

modelo matemático de permeabilidad cotejado con data de núcleos. En la figura 83

y 84 se puede observar una buena correspondencia con la data de núcleo y modelo

matemáticos respectivamente. Donde se observa el incremento de permeabilidad y

porosidad con la reducción del volumen de arcilla, lo cual esta representado por el

tipo de roca de mayor calidad “Tipo roca 1” en la arenas del yacimiento Bach-01.

Figura 83. Gráfico de permeabilidad de núcleo “Core K” Vs porosidad de núcleo “Core PHI” en función del volumen de arcilla (del lado izquierdo) y Tipo de Roca (del lado derecho), con los pozo claves del yacimiento Bach-01.

Figura 84. Gráfico de permeabilidad del modelo “K_PHIe” Vs porosidad del modelo “PHIE_D” en función del Tipo de Roca (del lado izquierdo) y “K_PHIe” (del lado derecho), con los pozo claves del yacimiento Bach-01.

Page 17: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

118

4.3.5.1 Determinación de permeabilidad absoluta a partir de modelos conocidos y

validado con data de núcleo

La permeabilidad absoluta “K_PHIe” de las areniscas de edad mioceno del

yacimiento Bach-01 se representa con un máximo valor según data de núcleo de

7300 md, se determinó para los 904 del estudio, seleccionándose el modelo

matemático de permeabilidad absoluta “K_PHIe” que ajuste con los datos de

permeabilidad de los análisis de pozos núcleos (claves) para yacimiento Bach-01. Se

utilizaron varios métodos conocidos según literatura, determinándose y validándose

el modelo matemático compuesto de “Timur-Coates” mediante los pozos claves, ya

que el mismo representa la característica de saturación del agua irreducible “Swir”,

tipo de roca y la permeabilidad absoluta para el yacimiento Bach-01. Se obtuvo un

buen cotejo con la data de núcleos (pozos claves) para el modelo matemático

compuesto “Timur-Coates”, mediante los análisis convencionales (Dean Stark), ver

figura 86 y87.

Para el modelo de permeabilidad se procedió a determinar los modelos

matemáticos conocidos según literatura, timur, timur modificado, Coates, Tixier, y

Willie–Rose, con la finalidad de efectuar comparaciones y seleccionar el que ajuste

con la permeabilidad medida en los análisis de núcleos de los pozos claves.

Posterior a esta comparación se tuvo que generan un modelo matemático

compuesto “Timur-Coates” de permeabilidad absoluta “K_PHIe” a partir de 2 modelo

conocidos, el cual represento las característica del tipo de roca, Swir y la

permeabilidad absoluta de las arenas de edad mioceno del yacimiento Bach-01.

Modelo matemático de Timur: 225.2 ))/()*100((_ SwirPHIEtimurK = (45)

Modelo matemático de Timur modificado:

22.2 ))/()*93((mod_ SwirPHIEtimurK = (46)

Page 18: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

119

Modelo matemático de coates:

22 ))/))1(**70(((_ SwirSwirPHIEcoatesK −= (47)

Modelo matemático de Tixier: 23 ))/()*250((_ SwirPHIEtixierK = (48)

Modelo matemático de Willie y Rose:

2))/()*79((__ SwirPHIERosewillieK = (49)

Modelo matemático compuesto Timur-Coates “K_PHIe”: 24 )1))06.0_/(06.1((*))12.0/((_ −+= LRJVCLPHIEcompuestoK (50)

A continuación se muestra una buena relación en la figura 85, donde se compara

los valores de la permeabilidad calculada mediante el modelo compuesto “K_PHIe”

con las medidas de permeabilidad de núcleos “Core K” de los pozos claves en

función del tipo de roca. Adicional los valores de la curva de permeabilidad “K_PHIe”

derivados a partir de la correlación establecida se ajustan a los datos medidos en

los núcleos de los pozos LL-3689 y LS-3803, observándose en la evaluación

Petrofísica de las figura 86 y 87.

Figura 85. Comparación entre la permeabilidad calculada del modelo compuesto “K_PHIe” y la obtenida de núcleos “Core K” de los pozos claves, en función de los tipos de roca para el yacimiento Bach-01.

Cor

e K

Crossplot Core K Vs K_PHIe en Función del Tipo Roca

Cor

e K

Cor

e K

Crossplot Core K Vs K_PHIe en Función del Tipo Roca

Page 19: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

120

Figura 86. Correlación núcleo-perfil de la evaluación petrofísica pozo clave (núcleo) LL-3689 de los Lentes superiores del yacimiento Bach-01. Cotejo de los modelos matemáticos de las propiedades petrofísicas y tipo de roca, incluyendo la permeabilidad absoluta “K_PHIe” de la curva del modelo compuesto de permeabilidad “Timur-Coates” solapando con la obtenida de núcleos “Core K o NU_PRM” del track N-6 en el pozo clave.

Modelo PHI Modelo K Modelo Sw

Modelo VCL

Curvas crudas (litología, Resistividad y porosidad) TIPO ROCA

ANT ANR ANP

Page 20: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

121

Figura 87. Correlación núcleo-perfil de la evaluación petrofísica pozo clave (núcleo) LS-3803 de los Lentes superiores del yacimiento Bach-01. Cotejo de los modelos matemáticos de las propiedades petrofísicas y tipo de roca, incluyendo la permeabilidad absoluta “K_PHIe” de la curva del modelo compuesto de permeabilidad “Timur-Coates” solapando con la obtenida de núcleos “Core K o NU_PRM” del track N-6 en el pozo clave.

Modelo PHI Modelo K Modelo Sw

Modelo VCL

Curvas crudas (litología, Resistividad y porosidad) TIPO ROCA

ANT ANR ANP

Page 21: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

122

En forma general para el modelo matemático de permeabilidad se derivo la

permeabilidad directamente de la porosidad utilizando facies y estratigrafía en los

pozos a los cuales se les había cortado núcleo (pozos claves). Las facies se

definieron a partir de las electrofacies (registro de resistividad y/o de rayos gama) y

según el ambiente predominante (Deltaico/ Fluvial) validado con la sedimentologia.

Las secuencias de bloque o de adelgazamiento hacia arriba se interpretaron como

canales y secuencias de engrosamiento hacia arriba podían ser en general barras o

como segunda opción abanicos de rotura, y por último secuencia completamente

arcillosas se interpreto como una facies de Llanura. En la figura 88, se observa en la

sección superior de los pozos claves que contiene más arenas limpia y tienen una

distribución similar de K-PHI en las arenas inferiores, ya que a lo largo de toda la

sección de la columna, tienen un nivel de compactación, tamaño del grano, y una

mineralogía consistente, sin embargo; la influencia principal sobre la relación baja

de K-PHI es el contenido de volumen de arcilla.

Figura 88. Correlación núcleo-perfil e interpretación de facies a partir de electrofacies (registros de gamma ray espectral y resistividad), y su comparación de la relación de permeabilidad y porosidad “K-PHI” en toda la columna de los pozos claves validado con los datos de núcleos del yacimiento Bach-01.

Coarse Up Dom. (Lo K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Coarse Up Dom. (Lo K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

LL-2318 LL-3689

Resultados de las Facies con relación K-PHI

Coarse Up Dom. (Lo K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

Coarse Up Dom. (Lo K-PHI)

Channel Dom. (Hi K-PHI)

Channel Dom. (Hi K-PHI)

LL-2318 LL-3689

Resultados de las Facies con relación K-PHI

Page 22: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

123

4.3.5.2 Determinación de permeabilidad efectiva mediante análisis especiales de los

pozos núcleos LS-5169 y LL-3689.

La permeabilidad determina el comportamiento del yacimiento y del pozo, pero el

término puede referirse a muchos tipos de mediciones. Por ejemplo, la

permeabilidad puede ser absoluta o efectiva, horizontal o vertical. Cuando un solo

fluido fluye a través de la formación, se puede medir una permeabilidad absoluta

“Ka”, mientras cuando dos o mas fluidos están presentes, cada uno reduce la

capacidad de fluir del otro a una saturación irreducible, lo cual es la permeabilidad

efectiva “Ke”, y la permeabilidad relativa es la relación entre la Ke y Ka.

Se determino una permeabilidad efectiva aproximada “Ke” en función de las

curvas de permeabilidad relativa de los análisis especiales de los pozos claves,

integrándose con los puntos del capitulo (4.3.3.2 y 4.3.4), y la caracterización del

tipo de roca del tabla 7, obteniéndose rangos de unas permeabilidad efectiva

aproximada al petróleo “Ko” a una saturación de agua irreducible “Swir” y unas

permeabilidad efectiva al agua “Kw” a una saturación de petróleo residual “Sor”

para los principales tipo de roca que integran lo prospectivo como roca yacimiento

(Bach-01), permitiendo determinar una relación de permeabilidad relativa para la

fase no mojante (Ko/Ka) y mojante (Kw/Ka) para el tipo de roca del yacimiento

Bach-01, lo cual se podrá determinar una permeabilidad efectiva aproximada en

una fase agua-petróleo para cualquiera de los 4 tipo de roca descritos en base a la

permeabilidad absoluta “Ka”. Ver tabla 8.

Tabla 8. Estimación de la permeabilidad efectiva al petróleo “Ko a Swir” y permeabilidad efectiva al agua “Kw a Sor” en base a la relación de permeabilidad relativa determinada de la fase agua-petróleo de los análisis especiales de los pozos claves (LS-5169 y LL-3689), integrada con la caracterización de los 4 tipos de roca que representan la prospectivo del yacimiento Bach-01.

Tipo de Roca Integrada

Rangos PHIE

Rangos Ka

Rangos Swir

Rangos Sor

Rangos Ko

Rangos Kw

Relación Ko/Ka

Relación Kw/Ka

cod % md % % md md % %

1 34 -_40 2400-7300 15 -_10 34 -_33 1359-4133 47-140 56,616 1,918

2 30 -_34 1400-2400 24 -_15 36 -_34 963-1359 40-47 68,786 2,857

3 25 -_30 700-1400 30 -_24 37 -_36 373-963 29-40 53,286 4,143

4 20 -_25 50-700 35 -_30 39-_37 26-373 3_-29 52,000 6,000

YACIMIENTOFORMACIÓN

MIEMBRO

Bach-01Lagunillas

Bachaquero

Page 23: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

124

4.4 Determinación de valores de corte (Cut-off).

Se integraron las evaluaciones petrofísicas de las curvas “VCL_LRJ, PHIE_D y

SW_SIMe” con los gráficos “crossplot” de los datos de núcleos de los pozos claves,

impregnaciones de petróleo en los núcleos e intervalos abiertos de producción, para

así determinar los valores de corte “Cut-off” de las propiedades petrofísicas.

Mediante estos valores “Cut-off” se permitió determinar la arena neta total “ANT”,

arena neta reservorio “ANR” y arena neta petrolífera “ANP”. Mediante el análisis

integrado de tipo de roca en los pozos claves (núcleos) y las impregnaciones de las

manchas de petróleo en los núcleos y datas de producción, se permitió identificar el

“TIPO DE ROCA” 1, 2, 3 y 4 como roca yacimiento o reservorio del Bach-01 (ANR) y

productora (ANP). Para estimar el volumen de arcilla de corte (VCL.cut-off<=22 %)

y determinar así las arenas neta total “ANT” del yacimiento Bach-01, se visualizaron

las frecuencias de la mayores manchas de petróleo <=22 % del volumen de arcilla

indicando zonas de rocas yacimiento prospectivas en la figura 89, como también en

la figura 37 se obtienen un valor de 22 % del modelo “VCL_LRJ” validado con data

de difracción de rayos X “VCL_XRD”.

Figura 89. Gráfico de frecuencias de manchas de petróleo identificando el volumen de arcilla de corte “VCL.cut-off” <=22 % para las arenas prospectivas del yacimiento Bach-01, mediante la data de los pozos claves (núcleos).

Para determinar el valor de corte o “Cut-off” de la porosidad (“PHIE cut-off”

=>20%), y obtener así las arenas neta reservorio “ANR” aplicando el valor de corte

de porosidad en las “ANT” del yacimiento Bach-01, se visualizaron en la figura 90

(de lado izquierdo) las frecuencias de la mayores manchas de petróleo =>20% de

Pequeñas manchas de petroleo >22 % VCL

Mayor frecuencia de Mancha de petroleo <22% VCL.cutoff

Poco recobro de núcleo

VCL_LRJ

Pequeñas manchas de petroleo >22 % VCL

Mayor frecuencia de Mancha de petroleo <22% VCL.cutoff

Poco recobro de núcleo

VCL_LRJ

Page 24: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

125

porosidad, indicando zonas de rocas yacimiento prospectivas, mientras que del lado

derecho de la figura 90, se realizo un gráfico cruzado (crossplot) de los modelos

matemáticos, volumen de arcilla “VCL_LRJ” con la porosidad “PHIE_D” versus la

saturación de agua, obteniéndose que lo mayoría de los puntos bajos de saturación

de agua (Sw<=50%, color rojo y rosado) caen por de bajo del 22% del VCL_LRJ y

por encima de 20% de porosidad, indicando zonas prospectivas de arena neta

petrolífera “ANP”.

Figura 90. Gráfico de frecuencias (del lado derecho) de manchas de petróleo identificando la porosidad de corte “PHIE_D.cut-off” =>20 %, y gráfico cruzado “crossplot” (del lado derecho) de los 3 parámetros de corte (VCL, PHIE_D y SW) para las arenas prospectivas del yacimiento Bach-01, mediante la data de los pozos claves (núcleos).

Adicionalmente según lo comentado, se determino y valido el valor de corte o

“Cut-off” de Sw (“SW cut-off” <= 50%) y “Cut-off” resistividad profundad

(RD_0>20 ohm-m), mediante los análisis especiales de permeabilidad relativa

agua-petróleo de los lentes superiores e inferiores del yacimiento Bach-01, se pudo

observar en el cruce de las curvas Kro-Krw en la figura 91 que la “SW cut-off” esta

en el orden de 50-60%, mientras que el gráfico cruzado (crossplot) en la figura 92

del modelo matemáticos de volumen de arcilla “VCL_LRJ” con la resistividad

profundad “RD_0” versus la saturación de agua, se obtuvo nuevamente (igual a la

figura 90) que lo mayoría de los puntos bajos de saturación de agua (Sw<=50%,

color rojo y rosado) caen por de bajo del 22% del VCL_LRJ, pero con resistividades

profundad => 20 ohm-m como “Cut-off”, indicando zonas prospectivas de arena

neta petrolífera “ANP”.

20%

PHIE_D

Mayor frecuencia de Mancha de petroleo >20% PHIE_D.cutoff

20%

PHIE_D

Mayor frecuencia de Mancha de petroleo >20% PHIE_D.cutoff

VCL cutoff <= 0.22

PHI=

>0.

20

SW

PHIE_D

VCL_

LRJ

VCL cutoff <= 0.22

PHI=

>0.

20

SW

PHIE_D

VCL_

LRJ

Page 25: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

126

Figura 91. Curvas de permeabilidades relativa agua-petróleo del núcleo LL-3689, de los lentes superiores e inferiores del yacimiento Bach-01.

Figura 92. Gráfico cruzado “crossplot” de los 3 parámetros de corte (VCL, RD y SW) para las arenas prospectivas “ANP” del yacimiento Bach-01, mediante la data de los pozos claves (núcleos).

Los parámetros de corte “Cut-off” se validaron además de impregnaciones del

núcleo, con producción del yacimiento Bach-01. A continuación se muestra un

resumen de los valores de corte “Cut-off” para determinar la evaluación Petrofísica

resultante (ANT, ANR y ANP) para el yacimiento Bach-01 (Ver tabla 9).

Tabla 9. Valores de corte “Cut-off” para la evaluación Petrofísica (ANT, ANR y ANP) del yacimiento Bach-01, validado con núcleos y producción.

Curva de permeabilidad relativa agua-petroleo Bach-01 " Lentes Inferior EE"

0,001

0,01

0,1

1

0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0Sw

Kr

Imbibición Kro (muestra 1902) Drenaje Krw (muestra 1902)

Curva de permeabilidad relativa agua-petroleo Bach-01 "Lentes Superior AP"

0,001

0,01

0,1

1

0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0Sw

Kr

Imbibición Kro (muestra 720) Drenaje Krw (muestra 720)

RD

_0=>

20 o

hm-m

Sw

VCL cutoff <= 0.22

RD_0

VCL_

LRJ

RD

_0=>

20 o

hm-m

Sw

VCL cutoff <= 0.22

RD_0

VCL_

LRJ

ANTVCL PHIE SW RD

Bach-01Lagunillas

Bachaquero<=22 % _=>20 % <=50 % _=>20 ohm-m

YACIMIENTOFORMACIÓN

MIEMBRO

PARÁMETROS DE CORTE “CUT-OFF”

ANRANP

Page 26: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

CAPITULO V

PRESENTACIÓN E INTEGRACIÓN DE RESULTADOS

5.1 Mapas de propiedades petrofísicas.

Las areniscas de edad mioceno del miembro Bachaquero formación lagunillas,

presentan un gran potencial petrolífero, por lo cual se generaron mapas de

propiedades petrofísicas con una buena correspondencia con lo mapas de facies y

tipo de roca para los 904 pozos evaluados, con el objetivo de analizar y validar la

distribución espacial de las propiedades petrofísicas de la roca, y así visualizar zonas

prospectivas del yacimiento Bach-01. Luego de obtener la curvas procesadas de la

evaluación petrofísica mediante los software Interactive Petrophysics “IP” y

Petroworks, se realizaron mapas 2D-3D de iso-propiedades en el software “PETREL”

para las 19 zonas que involucran toda la columna del yacimiento Bach-01, arena

neta total “ANT”, arena neta reservorio “ANR”, arena neta petrolífera “ANP”,

volumen de arcilla “VCL_LRJ”, permeabilidad absoluta “K_PHIe”, porosidad efectiva

“PHIE_D”, columna equivalente de petróleo “CEP”, saturación de agua “SW_SIMe”,

y Tipo de roca. Para la elaboración de los mapas iso-propiedades en 2D y 3D se

cargaron en el software “PETREL” las curvas procesadas con un step de 0.5 pies de

resolución vertical, lo cual presenta mayor certidumbre que cargar los sumarios con

un valor constante para cada zona. Mediante los valores de corte se permitió

establecer las rocas reservorios “ANR” y las areniscas prospectivas “ANP”,

validándose con pruebas de producción, y con la saturación de agua irreducible

“Swir” se lleva a condiciones iniciales del yacimiento Bach-01. Para la elaboración de

los mapas se utilizaron los sistemas de fallas, y limites del yacimiento,

determinados por las disciplinas de geología, geofísica e ingeniería de yacimientos,

orientándose en los mapas iso-propiedades los cuerpos arenosos con sus

propiedades petrofísicas según la dirección de la sedimentación establecida en el

mapa de facies. El control de calidad de los mapas se dirigió en la identificación de

valores inconsistentes y ojos de buey, por lo cual se corrigieron y suavizaron los

mismos. En la figura 93 se presenta los mapas iso-propiedades del yacimiento

Bach-01, “ANP”, “VCL_LRJ”, “PHIE_D”, “K_PHIe”, “CEP” y “SW_SIMe” con la

Page 27: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

128

superficie del Contacto agua petróleo actual CAPA a 3240 pies. Estos mapas

concuerdan con el modelo geológico y sedimentológico del yacimiento para los 19

Lentes, en la figura 93 se observa el ejemplo de la zona AP50 con una buena

distribución de sus propiedades con los mayores espesores de ANP y CEP para el

área del central del yacimiento Bach-01.

Figura 93. Mapas de propiedades petrofísicas en 2D del yacimiento Bach-01 AP50, con la superficie de CAPA a 3240 pies.

CEP ( acre-pies) SW_SIMe (fracción)

ANP (pies)

PHIE_D (fracción) K_PHIe (mD)

VCL_LRJ (fracción)

7300

CEP ( acre-pies) SW_SIMe (fracción)

ANP (pies)

PHIE_D (fracción) K_PHIe (mD)

VCL_LRJ (fracción)

73007300

Page 28: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

129

5.2 Metodología del Modelo petrofísico 3D en Petrel del yacimiento Bach-01,

integrando los mapas isopropiedades con los mapas del tipo de rocas y facies

sedimentarias.

El yacimiento Bach-01 corresponde a una acumulación petrolífera bastante

compleja, por lo tanto es de suma importancia la identificación de las zonas con

mejores propiedades petrofísica para la producción, por lo cual la integración del

modelo geológico-sedimentologico y petrofísico a través de la herramienta de

modelado PETREL, es de suma importancia para la caracterización Petrofísica 3D del

yacimiento con los 904 pozos evaluados. Por lo cual se observa en la figura 94, que

a partir de los resultados obtenidos del modelo estructural y estratigráfico se

construyó una malla 3D geológica de 396, 334 y 145 celdas en las direcciones X, Y

y Z respectivamente para un total de 19.178.280 celdas, con una resolución areal

de 210 mts aproximadamente y la resolución vertical del modelo estratigráfico la

cual es de aproximadamente 40 pies. Esta resolución se considera consistente para

generar los mapas isopropiedades en 3D, facies y Tipo de Roca para el yacimiento

Bach-01.

Figura 94. Representación del mallado geológico en 3D del yacimiento Bach-01, integrando con las 19 zona del modelo estratigráfico y con el sistema de fallas mas importantes del modelos estructural.

Page 29: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

130

Luego de tener el grillado 3D del modelo geológico, se cargaron los archivos

(.MFD) de los mapas de facies generados en el software ZMAP, en la herramienta de

modelado PETREL, desarrollándose un modelado de facies 3D determinístico (Ver

figura 7, y 100 en la parte de abajo), el cual captura desde el punto de vista

conceptual, la configuración de las facies predominantes en cada una de las

unidades estratigráficas del yacimiento, estas facies fueron agrupadas en 4 grupos,

que representan, en un ambiente Fluvial en su base de la columna y hacia el tope

Deltáico, los siguientes tipos de depósitos de: Canales fluviales o distributarios,

barras de desembocadura, abanicos de rotura y depósitos de llanura de inundación.

Estas facies fueron definidas en función de la interpretación de las electrofacies en

pozos (engrosamiento y afinamiento), tendencia en la relación del espesor bruto vs

la arena neta total “ANT” y arcillosidad.

Mediante la herramienta de modelado PETREL como se observa en la figura 100

(la parte de arriba) el modelo Tipo de Roca en 3D presentando un buen cotejo con

el modelo 3D de facies y según la metodología descrita a continuación se genero de

la siguiente manera: se uso el grillado mencionado con anterioridad, cargándose las

curvas de tipo de roca generadas en el software Petroworks, según la metodología

expuesta en el capitulo 4 (punto 4.3.4), luego se cambio la curvas de general a

variables discretas (Roca 1, 2, 3, 4, 5, 6 y 7, Ver figura 95), y posteriormente se

realizo el escalamiento (Ver figura 96), promediando la data a una escala micro (0,5

pies) a una mayor (+/-5 pies) a través de promedios aritméticos, condicionadas a

las interpretación de facies del modelo sedimentologico. La figura 97 muestra el

histograma de la variable discreta del tipo de roca, donde se puede observar como

la data escalada (Upscaled cells) a la malla geológica mantiene una distribución

estadística similar a la data de registro (Well logs). Luego de tener un control de

calidad de un buen cotejo de la data escalada a la de registro, se realizo análisis

geoestadístico convencional mediante variogramas y curva de proporción vertical

“CPV” para analizar la variabilidad o anisotropía natural de los datos y su

distribución en el espacio, como también ver el porcentaje que aporta cada tipo de

roca en la columna por la 19 zonas interpretadas del modelo estratigráfico del

yacimiento Bach-01. Ver figura 98.

Page 30: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

131

Figura 95. Representación del cambio de una curva general a discreta mediante cálculo en el software PETREL, caso ejemplo “Tipo de Roca” del yacimiento Bach-01.

Figura 96. Representación del escalamiento del tipo de roca en función de las facies sedimentarias, mediante cálculo en el software PETREL, yacimiento Bach-01.

Figura 97. Histograma de comparación de la data escalada del tipo de roca con un buen cotejo con la data de registro, mediante en el software PETREL, zona AP60.

Page 31: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

132

En lo que respecta al modelo de variograma, se seleccionó un modelo tipo

esférico, esto fue determinado con el buen ajuste de los variogramas verticales, los

cuales por lo general tienden a ser más estables, tal como muestra la siguiente

figura 98 (del lado izquierdo), donde se puede observar el efecto pepita o efecto

nugget (valor del variograma a distancia 0), se definió igualmente con base en el

variograma vertical, estimándose con un valor de 0.054 (<20% buena correlación y

mejora acercándose al valor cero), mientras que efectos nugget mayores a 40% por

lo general se consideran atribuibles a errores en la data y no representan la física de

la correlación espacial (efectos nugget de 100% indican cero correlación). En la

dirección horizontal, para generar un análisis variográfico representativo, se utilizó

lo que se denomina data suave en lo que respecta a la dirección de máxima

continuidad horizontal, correlacionando el variograma como se observa del lado

izquierdo en la figura 98, para un valor máximo (Major range) de 500 y mínimo

(Minor range) de 300 con un rango vertical de 7, para este caso, siguiendo la

tendencia de los mapas de facies y el conocimiento geológico previo se estimo una

tendencia de mayor continuidad de N 45° E.

Se realizó un análisis geoestadístico de tipo de roca a través de las variaciones de

proporción vertical (Curvas de Proporción Vertical “CPV”) y probabilidad. La

distribución de aporte de cada tipo de roca por unidad geológica y distribución areal

a partir de los mapas de facies del modelo determinístico. Estas CPV fueron

definidas para cada intervalo del modelo estratigráfico, lo cual permite al algoritmo

de “Sequential indicador simualation” honrar la distribución vertical de cada roca

con respecto a cada facies. En cuanto a la a distribución areal de la proporción de

cada tipo de roca, se determinó en función de la fracción que presente cada facies

del modelo determinístico. Por ejemplo, la proporción del tipo de roca 1 en la Zona

AP50 como se observa en la figura 98 (del lado derecho), presenta las mejores

propiedades petrofísicas representando por facies predominantes de depósitos de

Canal y Barra.

Una vez modelado los variogramas y CPV para todos los tipos de rocas en las 19

unidades estratigráficas, se procedió a aplicar el algoritmo de “Sequential indicador

simualation” para generar la realización del volumen y la distribución del tipo de

roca del yacimiento en 3D, condicionada con la data de pozos y la distribución de

Page 32: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

133

facies del modelo 3D. Ver figura 99 ejemplos del volumen del tipo de roca en las

zonas AP60 y AP50.

Figura 98. Análisis geoestadísticos mediante el software PETREL, ejemplo Tipo de roca 1 zona AP50, del lado izquierdo modelado de variograma vertical, y del lado derecho la curva de proporción vertical “CPV” y la probabilidad, integrándose el modelo 3D de Facies sedimentarías. Figura 99. Generación del volumen en 3D del tipo de Roca zona AP60 y AP50 en Facies Modeling (PETREL), integrando la dirección del modelo 3D de Facies del yacimiento Bach-01 y los valores obtenidos del variograma y CPV.

Page 33: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

134

Figura 100. Modelo 3D del Tipo de Roca (arriba), con un buen cotejo con el modelo 3D de Facies sedimentarías (abajo), para las 19 zonas del yacimiento Bach-01.

Luego de generar el modelo de tipo de roca 3D integrado con el modelo de

facies, se elaboro los mapas en 3D de las isopropiedades (VCL_LRJ, PHIE_D,

SW_SIMe, K_PHIe, ANT, ANR, ANP, CEP y Net_ANR: ANR/ANT “% de Arena

Reservorio”), siguiendo la metodología del tipo de roca a partir del escalamiento,

pero como las propiedades petrofísicas, ejemplo porosidad son variables continuas

se utiliza el modulo de la herramienta PETREL “Petrophysical modeling” para

generación del volumen tridimensional de las propiedades petrofísicas siguiendo la

dirección de las tendencia del tipo de roca y las depofacies sedimentarias (canal,

barra, abanico y llanura). Los mapas de tipo de roca se realizaron con el objetivo de

caracterizar las propiedades petrofísicas y visualizar las tendencias de las iso-

Page 34: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

135

propiedades correspondientes a las facies sedimentarias. Los tipos de rocas

corresponde a los rangos de propiedades interpretadas en la sección 4.3.4 del

capitulo 4. El modelo petrofísico 3D (Net_ANR, VCL_LRJ, PHIE_D y K_PHIe)

integrado con los modelos 3D del Tipo de roca y facies, se puede observar en la

figura 101 ejemplo en la zona AP50 (hacia el tope de la columna del yacimiento

Bach-01) validando el modelo sedimentologico con un ambiente deltáico, y en la

figura 102 ejemplo zona EInf (hacia la base de la columna del yacimiento Bach-01)

como un ambiente fluvial, donde los mapas de las propiedades petrofísicas en 3D de

las figuras 101 y 102, se observa que los mejores Tipo de roca 1, 2, 3 y 4

representa la roca yacimiento, honrando con las mejores propiedades petrofísicas el

modelo 3D de facies siendo estos canales, barras principalmente y abanicos de

manera secundaria, mientras que las rocas no yacimiento (desde la roca 5 hasta la

7), depofacies de abanico y llanura, con propiedades petrofísicas bajas de mala

calidad, por lo cual se agrupan por razones prácticas, denominándose como roca 5,

para un modelo de “Tipo de Roca simplificado en 3D" comparándose de manera mas

fácil con los mapas en 3D de las propiedades petrofísicas, y correspondiéndose así

con un buen cotejo con las facies sedimentarías.

Como se hizo referencia con anterioridad, los mapas de tipos de roca nos pueden

ayudar a identificar ambientes geológicos y la correspondencia de las mejores y

malas propiedades petrofísicas con las facies sedimentarias para las 19 zonas del

yacimiento Bach-01. Los mapas del % de arena reservorio e isopropiedades son

útiles para confirmar que la interpretación geológica y sedimentológica está

representada con el modelo 3D del Tipo de roca. Las figuras 101 y 102 presentan

los mapas 3D de isopropiedades, % de arena reservorio, Tipos de roca simplificado

y facies sedimentarias interpretados para las 19 zonas, pero en las figuras

mencionadas se refleja el ejemplo con las zonas AP50 (Hacia el tope del yacimiento)

y EInf (base del yacimiento). Se muestran en la figuras 101, zona AP50 que los

mapas de propiedades petrofísicas y los tipos de rocas (1, 2 y 3, colores claros)

muestran una alta calidad de roca, con porosidades (marrón) y permeabilidades

(amarillo y rojo) altas, volumen arcilla bajos (amarillo con gris) con alto % de arena

reservorio (amarillo claro), representando depofacies de canales (amarillo claro) y

barras (amarillo), mostrando una correspondencia entre el modelo Petrofísico 3D

con el modelo de Tipo de roca y Facies sedimentarias. Se distinguen en la figuras

Page 35: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

136

102 que la parte inferior del yacimiento Bach-01 (ejemplo zona EInf,) es de menor

calidad que la parte superior del yacimiento (ejemplo zona AP50, Ver figura 101),

ya que en cuanto a propiedades petrofísicas presenta valores altos de volumen de

arcilla (verde oscuro), baja permeabilidad (rosado) y porosidad (amarillo claro) con

bajo % de arena reservorio (gris oscuro), representándose depofacies de abanico

(naranja) y Llanura (verde) con tipos de roca 4 (celeste) y 5 (gris) predominantes.

Figura 101. Comparación de mapas de isopropiedades del modelo Petrofísico 3D (arriba y el centro), con el modelo Tipo de Roca en 3D simplificado (abajo derecha) y el modelo de facies sedimentarias en 3D (abajo izquierda). Las facies interpretadas por el modelo sedimentologico, son canales (amarillo claro), barra (amarillo), abanico (naranja) y llanura (verde). Zona AP50, yacimiento Bach-01.

Net_ANR “% Arena Reservorio” VCL_LRJ

PHIE_D

FACIES TIPO DE ROCA SIMPLIFICADO

K_PHIe 7300

Net_ANR “% Arena Reservorio” VCL_LRJ

PHIE_D

FACIES TIPO DE ROCA SIMPLIFICADO

K_PHIe 7300

Page 36: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

137

Figura 102. Comparación de mapas de isopropiedades del modelo Petrofísico 3D (arriba y el centro), con el modelo Tipo de Roca en 3D simplificado (abajo derecha) y el modelo de facies sedimentarias en 3D (abajo izquierda). Las facies interpretadas por el modelo sedimentologico, son canales (amarillo claro), barra (amarillo), abanico (naranja) y llanura (verde). Zona EInf, yacimiento Bach-01.

K_PHIe7300

Net_ANR “% Arena Reservorio” VCL_LRJ

PHIE_D

FACIES TIPO DE ROCA SIMPLIFICADO

K_PHIe7300

K_PHIe7300

Net_ANR “% Arena Reservorio” VCL_LRJ

PHIE_D

FACIES TIPO DE ROCA SIMPLIFICADO

Page 37: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

138

5.3 Determinación de contactos originales de los fluidos “CAPO” y su distribución

actual “CAPA”.

Para obtener una mayor certidumbre a la hora de determinar el petróleo original

en sitio “POES”, es necesario predecir con la mayor exactitud posible el contacto

agua petróleo original “CAPO” y su avance.

El yacimiento Bach-01 se limita hacia el suroeste y noroeste por un acuífero de

irregular a pasivo, en donde existe un contacto agua-petróleo original “CAPO”

común para todos los intervalos estratigráficos, con una inclinación del mismo hacia

la zona norte y noroeste de la estructura de mayor buzamiento que el yacimiento, y

una inclinación del CAPO hacia la zona sur y suroeste menor al buzamiento del

yacimiento. El avance de este acuífero a lo largo de la vida productiva del

yacimiento ha sido irregular, evidenciándose en distintas regiones de la zona sur

una producción de agua en los pozos. Los pozos que han mostrado irrupción de

agua, según el análisis petrofísico provienen de las arenas inferiores a producción

como consecuencia principal del mencionado avance irregular del acuífero por

fenómeno de conificación y/o adedamiento, y el aporte de la inyección alternada de

vapor “IAV” (1971), resultantes de la relación adversa de movilidad entre el crudo

pesado y el agua del acuífero, los efectos de las fuerzas capilares en el yacimiento y

las características locales de la sedimentación. Mientras que en la parte alta de la

estructura el acuífero presenta un avance pasivo y el mismo tiene procedencia de

agua fresca, donde CAPO se aprecia hacia los lentes superiores. Además, la

acumulación de petróleo en las capas más profundas parece seguir las tendencias

de las facies y está atrapada localmente por cambios en las facies. Esto da por

resultado “agua sobre petróleo” tanto estructuralmente como estratigráficamente.

Con el objetivo de representar en forma detallada el comportamiento

hidrodinámico del acuífero pasivo e inclinado presente en el yacimiento Bach-01, y

así identificar el CAPO y el avance del contacto agua petróleo actual “CAPA”, se

realizo la interpretación de los registros integrando con la evaluación petrofísica y la

producción. Primero se identificaron 2 épocas donde el cambio de los fluidos es

evidente, la cesta de pozo perforados al inicio de la producción del yacimiento

donde no fue sometido a ningún proceso de inyección térmica (IAV) hasta el año

Page 38: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

139

1971, por lo que se pudiese utilizar este periodo inicial (37 años), para cotejar el

tiempo de irrupción del agua en algunos pozos, y los pozos perforados hasta la

actualidad como segundo periodo. Luego se cargo las interpretaciones del CAPO del

al inicio de la producción del yacimiento, y del CAPA de los pozos perforados

después del IAV, con pozos nuevos perforados hasta la actualidad, en el software de

modelado PETREL, para así crea una superfiecie de los contactos agua-petróleo

mediante profundidades interpretadas de CAPO y CAPA según la metodología

mencionada con anterioridad (Ver figura 103). Posteriormente en la figura 103 se

observa, la selección 2 pozos con interpretaciones cercanos al CAP, uno hacia la

zona sur (LL-2268) y otro en la zona norte (LL-1791), con el objetivo atravesar la

mayor parte de la columna, visualizando el tope “AP60” y la base “AB” de

yacimiento con los mapa de CEP para detallar zonas mojadas, obteniendo como

resultado que el CAPO inclinado pudiese tener una variación desde profundidades de

3830´en la parte suroeste del campo hasta 3425´en la parte alta del yacimiento

Bach-01 (zona noroeste), mientras que el avance del CAPA varia desde

profundidades de 3210´ al suroeste del yacimiento hasta 2810´en la parte noroeste

y alta de la estructura, evidenciando el avance del acuífero en forma pasiva, debido

a que el yacimiento tiene 75 años de producción.

A pesar que esta metodología expuesta con anterioridad puede presentar algún

detalle, debido a la interpretación puntuales de CAP en los pozos evaluados, se

valido la misma con los análisis SCAL (Permeabilidad relativa y Presión capilar)

aplicados en la inicialización del modelo de simulación del Yacimiento Bach-01, Fase

III PDVSA-CMG, obteniéndose que los valores de los 2 pozos referencias

mencionados con anterioridad entran dentro del rango con valores muy similares

del CAPO para la zona suroeste y la zona alta de estructura (noroeste del

yacimiento). La figura 104 (fuente Informe modelo de simulación del Yacimiento

Bach-01, Fase III PDVSA-CMG) muestra una sección transversal del modelo en 3D

realizado en el software STARS, donde se puede observar claramente la inclinación

del contacto agua-petróleo, con los pozos LL231 y UNA35 del año 1934 que dan

inició con la producción del yacimiento Bach-01, con la saturación de agua inicial

validando el CAPO.

Page 39: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

140

Figura 103. Comparación de las superficies del CAPO y CAPA, ubicando el tope y base del yacimiento mediante las superficies de la CEP en diferentes formas trasversales en 3D en PETREL. Pozos con interpretación de CAP hacia la zona noroeste y suroeste del yacimiento, evidencia el avance del acuífero en forma pasiva en el yacimiento Bach-01. Figura 104. Representación del sección transversal del modelo en 3D del contacto agua-petróleo en el software STARS, con la saturación de agua inicial y los pozos que dieron inicio a la producción del yacimiento Bach-01.

Zona Suroeste CAPO @3830´ “AB”

Zona Suroeste CAPA @3210´ “AP30”

Zona Noroeste CAPA @2810´ “AP40”

Zona Noroeste CAPO @3425´ “AB”

LL-1791

LL-2268

Superficie CAPO

Superficie CAPA

Superficie CEP AB

Superficie CEP AP60

LL-2268 LL-1791

Zona Suroeste CAPO @3830´ “AB”

Zona Suroeste CAPA @3210´ “AP30”

Zona Noroeste CAPA @2810´ “AP40”

Zona Noroeste CAPO @3425´ “AB”

LL-1791

LL-2268

Zona Suroeste CAPO @3830´ “AB”

Zona Suroeste CAPA @3210´ “AP30”

Zona Noroeste CAPA @2810´ “AP40”

Zona Noroeste CAPO @3425´ “AB”

LL-1791

LL-2268

Superficie CAPO

Superficie CAPA

Superficie CEP AB

Superficie CEP AP60

LL-2268 LL-1791

Superficie CAPO

Superficie CAPA

Superficie CEP AB

Superficie CEP AP60

LL-2268 LL-1791

Page 40: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

141

5.4 Presentación de plantilla (Template) de evaluación petrofísica integrada con

Tipos de rocas y facies litológicas (Litofacies).

Las figura 105 y 106 presenta los ejemplos de la representación gráfica de los

resultados obtenidos de la evaluación petrofísica de los pozos LL-2318 y LS-5169

integrando el tipo de roca, validándose así el modelo de TIPO DE ROCA con las

Litofacies, antes de extrapolarlo el mismo a los 904 pozos del estudio. La tabla 6,

muestra la descripción de la información contenida en cada grillado o “track”

presente en los gráficos que muestra los resultados de la evaluación petrofísica

realizada en cada pozo.

Tabla 10. Información de los Tracks del Template de la evaluación petrofísica.

Tracks Descripción

1

Curva de diámetro o calibrador del hoyo “CAL”. Rayos gamma

normalizado en el pozo “GRN1”, y por último presenta la curva de

potencial espontáneo “SP_0”.

2 Curvas de resistividad somera y profunda “RD_0 y RS_0”.

3 Curvas de porosidad corridas en el pozo, y procesadas o sintéticas

generadas por el modelos matemáticos “PHIT_SP y PHIT_GR”.

4 Cotejo de las curvas procesadas de porosidad total y efectiva con data

de núcleo “NU_PHIpcs”, y la distribución de los volumen de fluidos.

5 Modelo de permeabilidad “K_PHIe” generado para el yacimiento Bach-

01, validado con data de núcleo “NU_PRMpcs”.

6 Modelo de saturación de agua “SW_SIMe” generado para el yacimiento

Bach-01, cotejado con data de núcleo “NU_SW”. Escala de 0 a 1.

7 Modelo de volumen de arcilla, cotejado con data de núcleo “UN_Vcl”

8 Representación Litológica del los Tipo de Roca Integrado.

9 Litofacies de la descripción sedimentológica de núcleos. Este Template

es utilizado para los pozos con interpretación de litofacies en núcleo.

10 Topes y Lentes estratigráficos de la formación Lagunillas. Yacimiento

Bach-01

Page 41: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

142

Figura 105. Correlación núcleo-perfil de la evaluación petrofísica pozo clave (núcleo) LL-2318 de los Lentes superiores del yacimiento Bach-01. Ejemplo de evaluación petrofísica con registros de porosidad y análisis de núcleos, validando los modelos matemáticos de las propiedades petrofísicas y tipo de roca integrado con Litofacies sedimentarias.

Modelo PHI Modelo K Modelo Sw Modelo VCL Curvas (litología, Resistividad y porosidad) TIPO ROCA

ANT ANR ANP

LITOFACIES Topes “FM”

Page 42: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

143

Figura 106. Correlación núcleo-perfil de la evaluación petrofísica pozo clave (núcleo) LS-5169 del Lente superior AP del yacimiento Bach-01. Ejemplo de evaluación petrofísica con registros de porosidad y análisis de núcleos, validando los modelos matemáticos de las propiedades petrofísicas y tipo de roca integrado con Litofacies sedimentarias.

Modelo PHI Modelo K Modelo Sw Modelo VCL Curvas (litología, Resistividad y porosidad) TIPO ROCA

ANT ANR ANP

LITOFACIES Topes “FM”

Page 43: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

CAPITULO VI

CONCLUSIONES

Se obtuvo un 100 % del control de calidad de las curvas de los pozos claves,

pozos control, y pozos no control del presente estudio “Según procesos de

diagrama, espalme de certificación de Registros”.

Se determinaron los parámetros petrofísicos para el yacimiento Bach-01,

mediante los análisis convencionales y especiales de núcleos de los pozos LL-

2318, LL-3689, LS-3803 y LS-5169, a una presión de sobrecarga de 1570 Lpc,

teniéndose como resultado 1,66 de exponente de cementación “m”, 1,70 para el

exponente de saturación “n”, 1,0 de coeficiente de tortuosidad “a”, 2,65 gr/cc

densidad de grano “ρm” representando el tipo de litología de arenas clásticas, y

6 ohm-m de resistividad de las arcilla “Rsh”.

Se caracterizo el agua de formación del yacimiento Bach-01, a partir de las

sensibilidades de los análisis físicos químicos de agua antes y después de la

inyección de vapor, validándose mediante las sensibilidades del método de

Pickett Plot. Obteniéndose un patrón de agua fresca con una resistividad de

agua de formación “Rw” de 2,4 ohm-m a 130 ºF con equivalente de 1275 ppm

NaCl, y un patrón de agua salobre con un “Rw” 1,53 ohm-m a 130 ºF con un

equivalente NaCl de 2012 ppm para el yacimiento Bach-01.

Se hizo una caracterización del modelo de arcillosidad desde todo punto de

vista, estos resultados evidenciaron que para el yacimiento Bach-01, el tipo de

arcilla predominante es la “Caolinita” de origen detrítico, y su distribución con

mayor proporción es como “arcilla dispersa” y menor proporción “arcilla

laminar”. Adicional el modelo matemático de volumen de arcilla que representa

la arcillosidad y la litología de edad mioceno del yacimiento, es el modelo de

Larionov rocas terciarias o jóvenes “VCL_LRJ”, debido a que fue el que mas

ajustó con las pruebas de núcleos de difracción de rayos X, SF y SEM.

Page 44: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

145

El modelo de porosidad se cotejo con data de pozos núcleo o claves para el

yacimiento Bach-01, determinándose en los pozos de control petrofísico la

porosidad total “PHIT” con la curva del registro de densidad “RHOB”, mientras

que en los pozos no petrofísicos (sin registros de porosidad) se determino

mediante correlaciones “PHI&GR” y algoritmos matemáticos de redes neuronal,

validándose los resultados en los pozos de control petrofísicos o claves. Se

cálculo la porosidad efectiva condicionando para arenas limpias y arcillosas,

siendo el modelo que mas ajusto.

El modelo de saturación de los fluidos “Swir y Sor” se cálculo por la tendencia

con el índice de calidad de reservorio “RQI” de los análisis de especiales de

núcleos de los pozos claves. La saturación de agua de las areniscas de edad

mioceno del yacimiento Bach-01, se determinó mediante el modelo matemático

de “Simandoux” siendo el mismo el mas representativo de las característica de

saturación del agua para arenas arcillosas, con resultados similares a las

saturación de agua de los análisis convencionales y saturaciones de agua

irreducible “Swir” de los análisis especiales de núcleos.

Se genero un modelo confiable de TIPO DE ROCA INTEGRADA para el

yacimiento Bach-01, a partir de la correlación núcleo-perfil de gráficos cruzados

“crossplot” entre volumen de arcilla, porosidad efectiva y permeabilidad,

depofacies o facies, petrofacies, unidades de flujo, litofacies, descripción de

núcleos y ambiente sedimentario. Los tipos de roca representan la litología y las

propiedades petrofísicas del yacimiento, presentando la ROCA 1, 2, 3, 4, 5, 6 y

7 una buena correlación con las litofacies sedimentarias descritas en los núcleos.

El modelo matemático de permeabilidad absoluta “K_PHIe” es un modelo

compuesto de permeabilidad “Timur-Coates” integrado al tipo de roca, los

resultados obtenidos cotejan con los valores de las medidas en los pozos

núcleos. Se estimo una permeabilidad efectiva al petróleo aproximada “Ko a

Swir” y al agua “Kw a Sor” en base a la relación de permeabilidad relativa de los

análisis especiales de los pozos claves, integrada con la caracterización de los 4

tipos de roca que representan la roca reservorio prospectiva del Bach-01.

Page 45: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

146

Mediante en el software de modelado PETREL, se desarrollo un Modelo

Petrofísico y Tipo de Roca en 3D del yacimiento Bach-01, donde los mapas

isopropiedades presentan buena correspondencia con los mapas del modelo de

Tipo de Roca y Facies sedimentarias en 3D, representando el modelo geológico.

Los mapas isopropiedades muestran que las mejores propiedades petrofísicas

integrado al Tipo de Roca 1, 2, 3 y 4 se encuentran hacia la zona central del

yacimiento, representando la zona mas atractivas de roca-fluido, honrando las

mejores facies de canal y barra del modelo 3D, donde la parte inferior de la

columna estratigráfica del yacimiento es de menor calidad que la parte superior.

El yacimiento Bach-01 presenta un contacto de agua petróleo “CAP” inclinado

desde el suroeste del campo hasta la parte norte de la estructura. Mediante el

software de modelado PETREL, se represento el comportamiento hidrodinámico

del acuífero pasivo e inclinado, realizando dos superficies de CAP de épocas

diferentes de cestas de pozos perforadas separadas en 1971 (IAV), obteniendo

como resultado que el CAPO y CAPA inclinado tiene una variación de 3830´ y

3210´al suroeste del campo hasta 3425´ y 2810´ en la parte alta del

yacimiento Bach-01 respectivamente, evidenciando el avance del acuífero en

forma pasiva, debido a que el yacimiento tiene 75 años de producción.

Page 46: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

147

CAPITULO VI

RECOMENDACIONES

Se recomienda efectuar la evaluación petrofísica al resto de los pozos

adyacentes o no incluidos dentro de este estudio mediante la metodología,

modelos y parámetros petrofísicos derivados en el presente trabajo. En vista a

los resultados se recomienda utilizar el modelo petrofísico en 3D para integrarlo

con las demás disciplinas y generar un modelo dinámico para el yacimiento.

En vista del avance leve del agua debido al acuífero pasivo en la zona suroeste

del yacimiento, se recomienda la perforación y/o rehabilitación de pozos en la

zona central. Se recomienda incorporar nuevas pruebas de producción para

validar la prospectividad de la zona norte, debido al acuífero de agua fresca.

Se recomienda incorporar el último núcleo cortado en el yacimiento Bach-01 con

el modelo petrofísico 3D, integrándose con nuevas innovaciones de registros

que permitan medir a hoyo entubado en pozos viejos de la franja del Kilómetro.

Page 47: Caracterización Petrofísica del Yacimiento Bach-01, …tesis.luz.edu.ve/tde_arquivos/150/TDE-2014-10-07T14:36:45Z-5269/... · Los rangos de las propiedades petrofísicas, análisis

148

BIBLIOGRAFÍA

Asquith G. B, (1991) Log Evaluation of Shaly Sandstones: A Practical Guide for SPWLA. Continuing Education Series#31. Tulsa, Oklahoma USA. Archie, G. E. The Electrical Resistivity Log as an Aid in Determining Some Reservoir Characteristics, Transactions. AIME 146, 1942, pp. 54-67. Bassiouni, Z. (1994) Theory, Measurement, and Interpretation of Well Logs. SPE Textbook Series Vol. 4 Louisiana USA. Corelab, (2005) Estudio de las propiedades rutinarias de la roca y Estudio de las propiedades avanzadas de la roca. Pozo LL-3689. Maracaibo, Venezuela. CMG-PDVSA (2006) ASISTENCIA TECNICA Bachaquero-01 Fase III Reporte de Progreso. Dresser Atlas, (1982) Well Logging and Interpretation Techniques. The Course For Home Study. Third Edition. EXXON (1982), P. W. Holbrook; Facies and environmental interpretation of Bachaquero-1 core from LL 2318, Lagunillas field, Lake Maracaibo. EXXON-LAGOVEN (1981), P. W. Holbrook and F. J. Moretti; Reservoir geology of the Bachaquero-1 reservoir, Bolivar Coastal Field, Lake Maracaibo, Venezuela. IRT-PDVSA (2005); Estudio Integrado del Yacimiento Bach-01, Campo Lagunillas Fase II, Descripción del yacimiento. Larinov, V. V Borehole Radiometry. Nedra, M, 1969, p.327. Malandrino, G. (2004) Petrofísica de Yacimientos. Materia de Maestría en Geología Petrolera de Grado. División de Postgrado. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela. Schlumberger, (1991) Principios/Aplicaciones de la Interpretación de Registros. Schlumberger Education. Salazar, V. (2006) Caracterización Petrofísica del Yacimiento Lagunillas Inferior 05. Trabajo de Grado. División de Postgrado. Facultad de Ingeniería. Universidad del Zulia. Maracaibo, Venezuela. Schlumberger Petrel (2007), Manual Property Modeling Course. 2nd Edition.