CINÉTICA DE LA GENERACIÓN DEL PETRÓLEO: PRINCIPIOS Y APLICACIÓN EN LAS CUENCAS COLOMBIANAS

GONÇÀLVES, F. T. T; GARCÍA, D. F; PENTEADO, H. L. B; GIRALDO, B. N; BEDREGAL, R. P; GÓMEZ, P

Services on Demand

Journal

Article

Indicators

Cited by SciELO
Access statistics

CT&F - Ciencia, Tecnología y Futuro

Print version ISSN 0122-5383On-line version ISSN 2382-4581

C.T.F Cienc. Tecnol. Futuro vol.2 no.2 Bucaramanga Jan./Dec. 2001

CINÉTICA DE LA GENERACIÓN DEL PETRÓLEO: PRINCIPIOS Y APLICACIÓN EN LAS CUENCAS COLOMBIANAS

F. T. T. GONÇÀLVES^*4, D. F. GARCÍA^*4, H. L. B. PENTEADO¹,
B. N. GIRALDO¹, R. P. BEDREGAL², P. GÓMEZ³

¹Petrobras - Cenpes, Cidade Universitária, Quadra 7, Rio de Janeiro, Brasil
²UFRJ - Coppe/Lamce, Cidade Universitária. CT Bloco I-2000, Río de Janeiro. Brasil
³Gems Ltda., Paraguitas I, Casa 13D. Floridablanta, Santander, Colombia
⁴Ecopetrol - Instituto Colombiano de Petróleo, A.A. 4185 Bucaramanga, Santander Colombia

e-mail: fgoncalv@ecopetrol.com.co - e-mail: dgarcia@ecopetrol.com.co

^*A quien debe ser enviada la correspondencia

RESUMEN

La mayoría de los modelos que describen el proceso de conversión de kerógeno en petróleo están basados en las formulaciones de la cinética química de primer orden. Aunque la aplicación de estos modelos requiera el conocimiento de los parámetros cinéticos (energía de activación y factor de frecuencia) del kerógeno de la roca madre, es práctica común en los estudios de modelamiento de cuenca usar los parámetros cinéticos y publicados de los kerógenos estándar (Tipos I, II y III) en ausencia de datos medidos. Por medio de análisis de Pirólisis Rack-Eval bajo diferentes tasas de calentamiento y usando técnicas de optimización numérica se determinaron los parámetros cinéticos de muestras de rocas madres de edades cretácica y terciaria de las cuencas del Valle Superior del Magdalena y los Llanos. Los resultados obtenidos revelaron una significativa variabilidad de los parámetros cinéticos, independiente de su clasificación, según los índices de hidrógeno y oxígeno, en kerógenos Tipos I, II y III. Ejercicios de moldeamiento con los parámetros medidos y una tasa de calentamiento constante (274.25K/M.y. 1,25°C/M.c.) revelaron que la conversión del potencial generador de faces orgánicos con parámetros cinéticos distintos puede estar desfasada hasta en 20-30 M.a., volviendo la historia de generación y expulsión del petróleo más larga y compleja de lo que se podría suponer, si consideramos el comportamiento cinético de estas rocas como homogéneo. Tales diferencias tienen particular importancia sobre la definición del sincronismo (Jímiug) entre la generación del petróleo y la formación/destrucción de trampas, más aún, si tenemos en cuenta la complejidad de la evolución tectónica de las cuencas colombianas.

Palabras clave: generación de petróleo, parámetros cinéticos, rocas madres, group villeta, formación mirador, Valle Superior del magdalena, cuenca de los llanos.

ABSTRACT

Most of the mathematical models that describe the conversion of kerogen into petroleum are based on the formulations of first-order kinetics. Although the application of such models requires the knowledge of the kinetic parameters (activation energies and frequency factor) of the kerogen, the usual practice in basin modeling studies is to use kinetic data of standard kerogen types (I. II or III) when measured data are not available. In this study, Rock-Eval pyrolysis under different heating rates and numerical optimization techniques were used to determine the kinetic parameters of Cretaceous and Tertiary source rocks of the Upper Magdalena and Llanos basins. The obtained kinetic parameters revealed a significant variability, which appears to be unrelated to the kerogen type classification based on hydrogen and oxygen indices Modeling exercises under a constant heating rate (1.25°C/M.y., 274.5K/M.y.) using the measured kinetic data indicates that kerogen conversion of organic facies with distinct kinetic parameters may be out of phase by 20-30M.y. Therefore. petroleum generation and expulsion history might be longer and more complex than if the kinetic behavior of these rocks was considered homogeneous. These differences are critical in defining the timing between petroleum generation an trap formation/destruction, particularly in the case of the Colombian sedimentary basins, characterized by a highly complex tectonic evolution.

Keywords: petroleum generation, kinetic parameters, source rocks, group villeta, mirador formation, Upper Magdalena and Llanos basins, llanos Basin.

INTRODUCCIÓN

El petróleo es producto de la degradación térmica del kerógeno debido al progresivo incremento de temperatura asociado al enterramiento de las rocas madres (Tissot y Welte. 1984). Diversos estudios han demostrado que este proceso de degradación puede ser descrito en gran parte de los casos por modelos basados en las formulaciones de la cinética química. En estos modelos cinéticos, la conversión del kerógeno es representada a través de una serie de reacciones paralelas. irreversibles e independientes de primer orden (Tissol et al, 1987). La evolución de estas reacciones es función de la temperatura y depende de los parámetros cinéticos del kerógeno, los cuales reflejan su grado de "reactividad".

En la actualidad, el modelamiento cinético es aplicado de forma rutinaria en los estudios de sistemas petrolíferos para la definición del sincronismo entre la generación del petróleo y la formación de las trampas. La aplicación de estos modelos requiere el conocimiento de los parámetros cinéticos (energía de activación y factor de frecuencia) del kerógeno de la roca madre. Estos parámetros pueden ser determinados a través de análisis de pirólisis Rock-Eval bajo diferentes tasas de calentamiento, combinados con técnicas de optimización numérica (Ungerer, 1993). Aunque sea práctica común en los estudios de modelamiento usar parámetros cinéticos publicados de los kerógenos estándar (Tipos I. II y III; Tissot et al, 1987) en ausencia de datos medidos, algunos autores (Tegelaar y Noble, 1994) han demostrado que dependiendo de la estructura molecular y de la composición química del kerógeno los parámetros cinéticos pueden variar significativamente. afectando el sincronismo de generación y expulsión de petróleo del sistema petrolífero.

Son pocos los estudios de caracterización de las propiedades cinéticas de las rocas madre colombianas (Giraldo et al, 1996. Jarvie y Wavrek, 1996). El presente trabajo pretende proveer los principios de la cinética química aplicada al estudio de la generación de hidrocarburos y discutir los métodos analíticos usados en la determinación de los parámetros cinéticos y la relación entre el tipo del kerógeno y su comportamiento cinético. Se presentan también los primeros parámetros cinéticos adquiridos después de implementada la metodología en el Laboratorio de Geoquímica del Instituto Colombiano del Petróleo (ICP) y la aplicación de estos parámetros en la simulación de las historias de generación para rocas generadoras del Cretáceo y Terciario en cuencas de Colombia.

MODELOS CINÉTICOS DE GENERACIÓN DEL PETRÓLEO

Con el objeto de facilitar la discusión de los modelos cinéticos de la transformación del kerógeno en petróleo, se presenta de forma resumida los conceptos básicos de la cinética química de las reacciones de primer orden. aquellas en las cuales un único reactivo (A) genera los productos (B. C. etc.):

En estas reacciones, la disminución relativa de la cantidad de reactivo es directamente proporcional al tiempo en que ocurre la reacción, conforme es expresado por la relación:

Donde x es la concentración del reactivo, k es la constante de tasa y t es el tiempo. La tasa de conversión (v) del reactivo en los productos también es proporcional a su concentración:

Por su parte, la constante de tasa (k) es función de la temperatura, conforme lo expresa la Ley de Arrhenius:

Donde. A es el factor de frecuencia o pre-exponencial, E es la energía de activación, R es la constante universal de los gases y T es la temperatura. El factor de frecuencia representa la frecuencia de choques entre las moléculas mientras la energía de activación puede ser comprendida como el umbral de energía por encima del cual las reacciones químicas comienzan. Estas dos variables (factor de frecuencia y energía de activación) son comúnmente conocidas como parámetros cinéticos y son características de cada reacción.

De acuerdo con el esquema cinético propuesto por Tissot y Espitalié (1975). la conversión del kerógeno en petróleo y del aceite en gas. procesos conocidos respectivamente como "craqueo" primario y secundario. pueden ser simuladas como una serie de reacciones químicas paralelas, independientes e irreversibles que tienen, cada una. su propia velocidad de reacción. A cada reacción corresponde un porcentaje del potencial generador de la roca, expresado bajo la forma del índice de Hidrógeno (Figura 1). A causa del poco conocimiento del reactivo (kerógeno). los productos, los coeficientes estequiométricos y los mecanismos de reacción. lo que llamamos "reacción", constituye en realidad la suma de pseudo-reacciones con energías de activación aparentes.

Partiendo de modelos cinéticos globales (Tissot y Espitalié. 1975), desarrollos posteriores han permitido describir los procesos de "craqueo" primario y secundario dividiendo los productos generados en distintas clases de compuestos (ej. C₁, C₂ -C₅, C₆-C₁₅, etc.). cada cual con potencial generador y distribución de energías de activación características (Espitalié et al, 1988: Ungerer, 1990; Behar et al., 1997). Estos desarrollos más recientes, conocidos como modelos cinéticos composicionales, permiten predecir la composición y propiedades físicas de los hidrocarburos formados.

DETERMINACIÓN DE LOS PARÁMETROS CINÉTICOS DE LAS ROCAS GENERADORAS

La determinación de los parámetros cinéticos ha sido efectuada mediante el uso de diversos métodos de pirólisis, en sistema cerrado o abierto, con o sin agua (Espitalié et al, 1988: Castelli et al., 1990: Schenk y Horsfield, 1993). La pirólisis anhidra en medio abierto es la técnica experimental más usada por la industria para la determinación de los parámetros de las reacciones de "craqueo" primario (Ungerer. 1984, 1990 y 1993: Ungerer et al. 1986: Braun y Burnham, 1987; Burnham et al. 1987: Ungerer y Pelel. 1987: Tissot et al, 1987; Espitalié et al, 1993). Sin embargo, este tipo de pirólisis no reproduce precisamente la forma en la cual las reacciones tienen lugar en la naturaleza, ya que el producto generado no presenta exactamente la misma composición que el petróleo (Lewan. 1985). A pesar de esta limitación, la aplicación de los parámetros cinéticos asi obtenidos en estudios de modelamiento de cuencas reproduce bien las observaciones hechas sobre series de muestras naturales (cantidad de extractos orgánicos, potenciales residuales, etc.; Ungerer et al. 1990; Forbes et al. 1991; Burros et al, 1996).

A partir de los datos experimentales de análisis de pirólisis realizadas a diferentes tasas de calentamiento, los parámetros cinéticos son rutinariamente obtenidos por medio de la optimización de funciones matemáticas que siguen las formulaciones de la cinética de primer orden y la Ley de Arrhenius. Dos aproximaciones que describen matemáticamente las reacciones de "craqueo" primario son comunes en la literatura. De una parte, algunos de los autores mencionados anteriormente usan una distribución discreta de energías de activación y un único factor de frecuencia. Por otro lado. (Qüigley et al, 1987). Braun & Burnham (1987) y (Buniham el al., 1987) usan una distribución gausiana de energías de activación y un factor de frecuencia. Además de no existir argumento teórico que indique que la distribución de energías de activación es estrictamente gausiana, Burnham et al. (1987 y 1988) demostraron que una distribución discreta explica mejor los resultados de laboratorio, Es importante mencionar que trabajos más recientes (Bumham et al. 1995 y 1996) han demostrado que el modelo cinético de reacciones paralelas de primer orden no describe de forma adecuada la conversión de algunos casos específicos (carbones, kerógenos algales).

Los procedimientos de adquisición y tratamiento de datos deben ser rigurosos, ya que el proceso de optimización numérica y las condiciones analíticas (tipo de gas de arrastre, precisión de las medidas de temperatura. tasas de calentamiento, etc.) pueden afectar la determinación de los parámetros cinéticos (Jarvie, 1991; Sundararaman et al., 1992: Espitalié et al., 1993). Otro factor que influye en los resultados es el grado de madurez termal de las rocas, ya que en muestras muy inmaduras (con valores de reflectancia de vitrinita menores que 0,4% Ro) la producción de no-hidrocarburos (H₂0, CO₂, NSO, etc.) afecta la medición del potencial generador y en muestras maduras (>0.6%Ro) una parte del potencial original ya fue convertido (Jarvie, 1991; Espitalié et al., 1993).

La conveniencia de usar muestras de roca total o de kerógeno aislado también es motivo de controversias. Si de una parte el uso de la muestra total preserva las condiciones naturales de la roca, el empleo del kerógeno aislado evita efectos de adsorción y/o catalíticos de la matriz mineral sobre los productos generados y el kerógeno respectivamente (Dembicki, 1992). Sin embargo, algunos estudios han demostrado que para rocas generadoras con altos contenidos de carbono orgánico (>2%). las diferencias entre los parámetros cinéticos obtenidos en roca total o kerógeno son mínimas (Pelet. 1994; Reynolds y Burnham, 1995).

COMPORTAMIENTO CINÉTICO DE LOS KERÓGENOS CLÁSICOS Y LA VARIABILIDAD DE LAS ROCAS MADRES

Actualmente se reconoce que los parámetros cinéticos son función de la composición química y de los tipos de enlaces presentes en las moléculas de un kerógeno. Por lo tanto la distribución de energías de activación y los factores de frecuencia que describen el "craqueo térmico" del kerógeno son también útiles en su caracterización (Ungerer y Pelet, 1987; Tissot et al, 1987; Espitalié et al, 1988; Behar et al., 1997). Esta dependencia entre la composición química y el comportamiento cinético es claramente demostrada en los tipos "clásicos" de kerógenos (Tissot y Espitalié. 1975; Tissot et al. 1987). En el Kerógeno Tipo I (Figura 2). el carácter prácticamente unimodal de la distribución de energías de activación (54 kcal/mol para un valor de A de 1.026 -10¹⁴s^-1) debe reflejar la relativa homogeneidad de los tipos de enlaces químicos, esencialmente C-C. El Kerógeno Tipo II. por contener mayor cantidad de heteroátomos (O.S.etc.) tiene una mayor diversidad de tipos de enlaces y por tanto presenta una distribución más amplia de energías de activación (variando entre 44 y 60 kcal/mol para un valor de A de 1.11-10¹⁵s^-1). La energía de activación principal, aquella a la cual se relaciona la mayor parte del potencial generador, es de 52kcal/mol. valor inferior a la del kerógeno tipo I. Finalmente, en el Kerógeno Tipo III (Figura 2), la presencia de diversos tipos de enlaces con predominio de aquellos con mayor estabilidad térmica, resulta en una distribución amplia pero con valores de energías de activación más elevados (entre 50 y 74kcal/mol. con un valor de A de 5,460 · 10¹⁶s^-1).

Las diferencias del comportamiento cinético entre los distintos tipos de kerógenos. tienen importantes efectos sobre el sincronismo de generación, ya que frente a la misma historia térmica ellos reaccionan de manera distinta (Figura 3). Debido a la distribución de energías de activación, el kerógeno Tipo II empieza a transformarse más temprano y su historia de generación es más larga comparada con la del kerógeno Tipo I. cuya distribución de energías unimodal resulta en una historia de generación prácticamente "instantánea". Por su parte, los altos valores de energías de activación del kerógeno Tipo III conducen a una historia de generación más tardía.

A pesar de utilizarse rutinariamente en el modelamíento de cuencas, los parámetros cinéticos de los kerógenos clásicos pueden no describir satisfactoriamente la variedad composicional de las rocas madres. En algunos casos, para la misma unidad de roca madre han sido observadas diferencias significativas en los parámetros cinéticos (Jarvie, 1991; Sundararaman et al. 1992; Bunvood el al. 1992: Gonçalves et al. 1994; Soldán et al. 1995). Con base en el análisis de un conjunto de 71 rocas madres. Tegelaar y Noble (1994) reconocieron siete patrones de distribuciones de energías de activación con comportamientos independientes de su clasificación según sus indices de hidrógeno en kerógenos tipos I, II y III. Estos autores sugieren que la distribución y magnitud de energías de activación refleja la proporción relativa de los diferentes tipos de bíomacromoléculas preservadas durante la diagénesis y el contenido de azufre orgánico. Otros trabajos (Tissot et al., 1987: Baskin & Peters 1992; Araújo & Tríguis. 1994. Tomic et al. 1995) también han demostrado que altos contenidos de azufre orgánico resultan en valores de energías de activación más bajos, posiblemente debido al hecho de que los enlaces C-S son más débiles que los enlaces C-C.

DETERMINACIÓN DE PARÁMETROS CINÉTICOS DE ROCAS MADRES COLOMBIANAS

Como una primera aproximación a la caracterización del comportamiento cinético de las rocas madres colombianas, fueron seleccionadas 11 muestras del Grupo Villeta en la Cuenca del Valle Superior del Magdalena (VSM), siendo ocho muestras de una sección aflorante en la Quebrada Bambuca (Kairuz et al., 2000) y tres muestras del pozo Toldado-1 (Figura 4 y Tabla 1). Las rocas pelíticas marinas de este grupo y sus equivalentes en edad en otras cuencas sedimentarias son responsables de la generación de más del 80% de las reservas de petróleo colombianas (Mora, 2000). Otras dos muestras de la Formación Mirador en el pozo Floreña-N2F fueran seleccionadas como posibles representantes de las rocas fuentes terciarias que contribuyen para la carga de los yacimientos del Piedemonte Llanero (Range) et al., 1999).

Todas las muestras analizadas poseen un elevado contenido orgánico, con tenores de carbono orgánico total (COT) variando entre 2.66% y 8.80% (Tabla 1). Las muestras del Grupo Villeta (VSM) están inmaduras (Tmax entre 415°C y 431 °C, 688K y 704K) y las de la Formación Mirador (Llanos) se encuentran justo en el inicio de la ventana de generación aceite (Tmax alrededor de 435°C,708K). Los índices de Hidrógeno (IH) variando entre 41() y 560mgHO/gCOT indican un predominio de kerógeno del tipo II en las muestras del Grupo Villeta. mientras que los valores de III alrededor de 300mgHC/gCOT presentados por las muestras de la Formación Mirador sugieren una mayor proporción de kerógeno del tipo III (Tabla 1).

Las muestras seleccionadas fueron pulverizadas y extraídas por 24h con diclorometano en un aparato tipo Soxhlet para la remoción del bitumen. Fracciones de 40-50mg de cada muestra fueron sometidas a la pirólisis en un equipo Rock-Eval 6 (Lafargue et al., 1998) bajo 6 tasas de calentamiento diferentes: 1, 2. 5. 10, 15 y 25°C/min (298K/min). La optimización numérica fue realizada con la ayuda del software Optkin (Beicip-Franlab), lo cual adopta el modelo de "craqueo* del kerógeno a través de una serie de reacciones paralelas de primera orden como base para calcular la distribución de energías de activación y el factor de frecuencia que mejor se ajustan con los resultados de la pirólisis. El proceso de optimización fué exitoso para todas las muestras analizadas, permitiendo reproducir muy bien las curvas de pirólisis con una distribución de energías espaciada a cada 2kcal/mol (Figura 5).

Los parámetros cinéticos obtenidos en este estudio revelaron la existencia de significativas diferencias entre las muestras analizadas, con factores de frecuencia variando entre 9.9-10¹² y 2.3-10¹⁵s^-1 y valores de energía activación principal (a la cual está asociada la mayor parte de la conversión del kerógeno) entre 50 y 58kcal/mol (Tabla 2). En el pozo Floreña-N2F. la muestra a 13.869 pies posee una distribución de energías de activación similar al kerógeno Tipo III (comparar Figuras 6 y 2) mientras la muestra a 13.440 pies presenta valores de E más bajos, aunque ambas muestras presenten valores de IH compatibles con un predominio de kerógeno Tipo III (Tabla 1).

Entre las muestras analizadas del Grupo Villeta se encontró también una significativa diversidad de patrones de energías de activación a pesar de la relativa homogeneidad de los valores de IH indicativa de un predominio de kerógeno Tipo II (Figura 7 y Tabla 1). Las distribuciones de energías de activación varían desde algo muy similar a la del kerógeno Tipo II clásico (muestra del pozo Toldado-1, 4.740 pies) hasta algo próximo a la distribución clásica de un kerógeno Tipo III (muestra BV-137 de la Quebrada Bambuca; comparar Figuras 7 y 2).

Los resultados obtenidos para las muestras analizadas indican que el empleo de la clasificación de la materia orgánica en tipos I, II o III a través de los índices de hidrógeno y oxigeno como base para definir su comportamiento cinético puede inducir a errores. La variabilidad cinética encontrada probablemente refleja diferencias en la estructura molecular, tipos de enlaces químicos y/o en la composición de los kerógenos de las muestras analizadas, aparentemente sin una relación directa clara con el valor del IH. En esta fase del estudio aún no es posible definir los factores que controlan la variabilidad cinética observada ni tampoco establecer los comportamientos cinéticos típicos de cada facies orgánica de las rocas madres estudiadas.

IMPLICACIONES DE LA VARIABILIDAD CINÉTICA DE LAS ROCAS MADRES SOBRE EL SINCRONISMO DE GENERACIÓN

Con el objetivo de evaluar el impacto de la variabilidad cinética encontrada sobre la historia de generación del petróleo, se realizaron simulaciones de la conversión del kerógeno con el software GENEX (Beicip-Franlab) considerando cada conjunto de parámetros cinéticos obtenido. Se adoptaron para la simulación temperaturas aumentando desde 20°C (293K) hasta 260°C (533K.) bajo una tasa de calentamiento de 1,25°C/M.a. (274.25°K/M.y.). permitiendo una completa conversión de los kerógenos analizados. Aunque una tasa de calentamiento única no permite describir la diversidad del régimen térmico de las cuencas abordadas en este estudio (Llanos y VSM), el valor adoptado en esta simulación preliminar se encuentra dentro del rango de variación esperado para estas cuencas.

Las gráficas de las Figuras 8 y 9 presentan los resultados de la simulación para las muestras de la formación Mirador y el Grupo Villeta respectivamente. Estas gráficas muestran la evolución de la tasa de transformación (TT) del kerógeno en petróleo en función del aumento de la temperatura y del tiempo de calentamiento para cada muestra analizada y para los kerógenos clásicos I, II y III.

Para la Formación Mirador, se observa que la conversión de las dos muestras se encuentra desfasada en aproximadamente 20°C-30°C (293K-303K) y 20-25M.a. para la tasa de calentamiento adoptada (Figura 8). La muestra de 13.869 pies tiene un comportamiento bastante similar al del kerógeno Tipo III. mientras la muestra de 13.440 pies se aproxima a la del kerógeno Tipo I. aunque ambos presenten valores de IH y 10 bastante similares (Tabla 1). En términos de sincronismo de generación y expulsión las diferencias observadas son significativas, por ejemplo alrededor de 110 M.a. después del inicio del calentamiento, el kerógeno de la muestra 13.440 pies estaría prácticamente agotado (TT>80%) mientras la muestra 13.869 pies apenas estaría alcanzando condiciones de expulsión (TT alrededor de 20%).

Con relación al Grupo Villeta (Figura 9), aunque la mayoría de la muestras sean compatibles con el kerógeno Tipo II. hay tres muestras con comportamiento bien distinto. La dispersión observada bajo esta tasa de calentamiento représenta un desfase de aproximadamente 20°C-25°C (293K-298K) y 20-25M.a. Como en el caso anterior, mientras los kerógenos similares al Tipo II clásico se encuentran en un nivel avanzado de conversión (TT>90% a 100M.a.), el kerógeno menos reactivo (Q. Bambuca. BV-137) apenas está empezando la conversión (TT> 10%). Estas diferencias no serían normalmente esperadas para un conjunto de muestras con valores de IH y IO relativamente similares y compatibles con un predominio de materia orgánica del tipo II (Tabla 1).

Las diferencias observadas entre las muestras deben ser analizadas bajo dos aspectos. De una parte, se observa que los parámetros cinéticos de los kerógenos clásicos no describen satisfactoriamente la variabilidad inética de las rocas madres estudiadas, ya que muestras que por sus valores de IH y IO serian clasificadas como un determinado tipo de kerógeno, no necesariamente presentan el comportamiento cinético típico de este kerógeno. De otra parte, tales diferencias sugieren que la historia de generación y expulsión a partir de las rocas madres estudiadas puede haber sido más larga y compleja de lo que se podría suponer, considerando el comportamiento cinético de estas rocas como homogéneo. En cuencas con una historia de evolución tectónica compleja como las que abordamos en este estudio, diferencias de 20 M.a. en la historia de conversión del kerógeno pueden ser criticas para la definición del sincronismo entre las edades de generación y de formación/destrucción de las trampas. Vale mencionar que bajo tasas de calentamiento más rápidas las diferencias en tiempo y temperatura disminuyen.

Los resultados de este estudio demuestran que existe una significativa variabilidad cinética. Sin embargo antes de definir la expresión volumétrica (en términos de potencial generador) de las facies orgánicas con comportamiento cinético diferentes del convencional, es prematuro evaluar la importancia de las variaciones encontradas sobre la historia de carga de hidrocarburos.

CONCLUSIONES

Con la implantación de la metodología de pirólisis aplicada al análisis cinético en el ICP, se pudieron determinar los parámetros cinéticos de algunas muestras de rocas generadoras colombianas. Se encontró una significativa variabilidad de comportamiento cinético para muestras del Terciario (Fm. Mirador) y del Cretáceo (Gr. Villeta). Los resultados obtenidos demuestran que el empleo de la clasificación de la materia orgánica en tipos I, II o III como base para definir su comportamiento cinético puede inducir a errores. La variabilidad cinética encontrada en nuestro estudio probablemente refleja diferencias en la estructura molecular, tipos de enlaces químicos y/o en la composición de los kerógenos analizados, sin una relación directa clara con el valor del IH.

La simulación de la conversión del kerógeno con los parámetros cinéticos obtenidos bajo una tasa de calentamiento de l.25°C/M.a. (274.25K/M.y.) reveló que las historias de generación pueden estar significativamente desfasadas en tiempo (20M.a.) y temperatura 20³C-,30°C (293K-303K). De esta forma, dentro de una única roca madre, mientras una facies orgánica se encuentra próxima de agotar su potencial generador, otra puede estar apenas empezando la conversión del kerógeno, volviendo la historia de generación y expulsión del petróleo más larga y compleja de lo que se podría suponer considerando el comportamiento cinético de estas rocas como homogéneo.

Los resultados aquí presentados demuestran la necesidad de estudios adicionales para definir de modo sistemático el comportamiento cinético de lass diversas facies orgánicas de las rocas madres colombianas, averiguar los factores que controlan la variabilidad encontrada y determinar la importancia década facies orgánica para la carga de hidrocarburos. Los estudios de modelamiento de cuenca deben tener en cuenta los diferentes comportamientos cinéticos para evaluar su impacto sobre el sincronismo de generación y entrampamiento.

REFERENCIAS

Araújo, L.M., and Trigüis. J.A., 1994. "Evaluation of hydrocarbon generation processes of the irati formation source rock, paraná basin, brazil". IV Latin American Congress on Organic Geochem., Bucaramanga, Colombia, Bolétin de resumos, 92. [ Links ]

Baskin, D.K., and Peters, K.F., 1992. "Early generation characteristics of a sulfur rich monterey kerogen". AAPG Bulletin. 76 (1): 1-13. [ Links ]

Behar, F., Vandenbroucke. M., Tang, Y., Marquis. F., and Espitaliè. I., 1997. "Thermal cracking of kerogen in open and closed systems: determination of kinetic parameters and stoichiometrice coefficients for oil and gas generation". Organic Geochem., 26 (5-6): 321 - 339). [ Links ]

Braun, R.L., and Burnham, A.K., 1987. "Analysis of chemical reaction kinetics using a distribution of activation energies and simpler models". Energy & Fuels, 1: 153 - 161. [ Links ]

Burnham, A.K., Braun. R.L., and Gregg, H.R., 1987. "Comparison of methods for measuring kerogen pyrolysis rales and fitting kinetic parameters". Energy & Fuels, 1:452 - 158. [ Links ]

Burham, A.K., Braun. R L., and Samoun, A.M., 1988. "Further comparison of methods for measuring kerogen pyrolysis rates and fitting kinetic parameters". Org, Geochem., 13 (4-6): 839-845. [ Links ]

Bumham, A.K., Schmidt. BJ. and Braun, R.L., 1995. "A test of parallel reaction model using kinetic measurements on hydrous pyrolysis residues", Org. Geochem., 23 (10): 931-939. [ Links ]

Bumham, A.K., Braun, R.L..Cobum, T.T., Sandvik. E.I., Curry E.J., Schmidt, BJ. and Noble, R.A., 1996. "An appropriate kinetic model for well-preserved algal kerogens". Energy & Fuels, 10:49-59. [ Links ]

Burwood. R., Leplat, P., Mycke, B. and Paulet. J., 1992. "Rifted margin source deposition: A carbon isotope and biomarket study of a West African Lower Cretaceous lacustrine section". Org. Geochem., 19:41-52. [ Links ]

Burros, .J., Osadetz, K., Wolf. S., Doligez, B., Visser, K. and Dearborn, D., 1996. "A two-dimensional regional basin model of Williston Basin hydrocarbon systems". AAPC Bulletin, 80(2): 265-291. [ Links ]

Castelli. A., Chiaramonte. M.A., Beltirame, P.L., Carniti, P., Del Bianco, A. and Stroppa. F., 1990. "Thermal degradation of kerogen by hydrous pyrolysis: a kinetic study". Org. Geochem., 16 (1-3):75-82. [ Links ]

Dembicki, H., 1992. "The effects of the mineral matrix on the determination of kinetic parameters using modified Rock-Eval pyrolysis". Org, Geochem,, 18:531-539. [ Links ]

Espilalié, J., Ungerer, P., Irwin. I. and Marquis, F., 1988. "Primary cracking of kerogens: experimenting and modeling C₁, C₂, C₅, C₆-C₁₅, and C₁₅, classes of hydrocarbons formed". Org. Geochem., 13 (4-6): 893 - 899. [ Links ] [ Links ]

Forbes, P.L., Ungerer. P.M., Kuhfuss, A.B., Riis. F. and Eggen. S., 1991. "Compositional modeling of petroleum generation and expulsion: trial application to a local mass balance in the Smørhukk Sør Field. Haltenbanken Area, Norway". AAPG Bulletin 75 (5): 873-893. [ Links ]

Giraldo, B.N., Rangel, A., Talukdar, S.C. and Burnham. A.K., 1996. "Oil generation pyrolysis kinetic of upper cretaceous source rocks, upper Magdalena Valley Basin. Colombia". Memorias del IV Congreso Latinoamericano de Geoquin. Org., Bucaramanga. 199 - 202. [ Links ]

Gonçalves, F.T.T., Trigüis. J.A., Gonzaga. F.G. and Coutinho, L.F.C. 1994. "Kerogen composition vs. kinetic variability in marine source rocks from Amazonas Basin. Brazil". IV Latin American on Org. Geochem., Bucaramanga, Colombia. [ Links ]

Jarvie. D.M., 1991. "Factors affecting Rock-Eval derived kinetics parameters". Chemical Geol., 93: 79 - 99. [ Links ]

Jarvie. D.M. and Wavreck. D. 1996. "Comparison of the timing of hydrocarbon generation for major petroleum source rocks in North and South America". II AAPG/SVG International congress and exhibition. Caracas, Venezuela. AAPG Bulletin. 80: 1303-1304. [ Links ]

Kairuz, E.C., Ferretra, P., Ortiz, A. and Blanco. Y., 2000. "Distribución e interpretación facial de las unidades calcáreas Tetuan, Bambuca y la Luna Grupo Villeta. Cuenca del Valles Superior del Magdalena". Informe Técnico Ecopetrol/AEX-ICP, 89. [ Links ]

Lafargue, E., Marquis. F. and Pillot, D., 1998. "Rock-Eval 6 applications in hydrocarbon exploration, production, and soil contamination studies". Revue de L'Institut Français du Petróle. 53(4). [ Links ]

Lewan, M.D., 1985. "Evaluation of petroleum generation by hydrous pyrolysis experimentation". Philos. Trans. R. Soc. London. Ser. A. 315:123-134. [ Links ]

Mora, C. A., 2000. "Evaluación del potencial de los sistemas petrolíferos en las cuencas Cretácicas con producción comercial en Colombia". Tesis de Maestría, Universidad Federal do Rio de Janeiro. Brazil, 250pp. [ Links ]

Pelet, R., 1994. "Comments on the paper The effects of the mineral matrix on the determination of kinetic parameters using modified Rock-Eval pyrolysis by H. Dembicki Jr.". Org, Geochem., 21 (8-9): 979-981. [ Links ]

Quigley, T.M., Mackenzie. A.S. and Gray. J.R., 1987. "Kinetic theory of petroleum generation". In: Migration of Hydrocarbons in Sedimentary Basins, Doligez, B. ed., Technip, 649-666. [ Links ]

Rangel, A., Brooks. P.W. and Giraldo, B., 1999. "Petroleum geochemistry of oils in the foothills of the Llanos Basin, Colombia". Rev. lat. Amer. de Geochem. Org., 5:5-23. [ Links ]

Reynolds, J.G. and Bumham, A.K., 1995. "Comparison of kinetic analysis of source rocks and kerogen Concentrates". Org. Geochem., 23(1): 11-19. [ Links ]

Schenk. H.J. and Horsfíeld, B., 1993. "Kinetics of petroleum generation by programmed-temperature closed- versus open-system pyrolysis". Geochim. et Cosmochim. Acta. 57:623-630. [ Links ]

Soldán, A.L., Gonçalves, F.T.T., Araújo, C.V, Martins, M.L., Menezes. S.M.C., Aquino. N., F.R. and Simoneit. B.R.T. 1995. "Kinetics of some Tasmanacean-rich kerogens". 17th Meeting on Org. Ceochem., San Sebastián. Extended Abstract. [ Links ]

Sundararaman, P., Merz, P,H. and Mann, R.G.. 1992. "Determination of kerogen aclivalion energy distributions". Energy & Fuels, 6: 793 - 803. [ Links ]

Tegelaar, E.W. and Noble. R.A., 1994. "Kinetics of hydrocarbon generation as a function of the molecular structure of kerogen as revealed by pyrolysis-gas chromatography". Org. Geochem., 22 (3-5); 543 - 574. [ Links ]

Tissot, B.P., and Espitalié, J., 1975. "L'évolution thermique de la matiére organique des sediments: applications d'une simulation mathématique". Revue de l'Institut Français du Pétrole, 30:743-777. [ Links ]

Tissot, B.P., Pelet, R. and Ungerer P., 1987. "Thermal history of sedimentary basins, maturation indices, and kinetics of oil and gas generation". AAPG Bulletin, 71(12): 1445-1466. [ Links ]

Tissot, B.P., and Welte, D.H., 1984. "Petroleum formation and occurrence". Springer-Verlag, Berlin. 2ª ed. [ Links ]

Tomic, J., Behar, F., Vandenbroucke, M. and Tang, Y., 1995. "Artificial maturation of Monterey kerogen (Type II-S) in a closed system and comparison with Type II kerogen: implications on the fare of sulfur". Org. Geochem., 23 (7): 647-660. [ Links ]

Ungerer, P., 1984. "Models of petroleum formation: how to take into account geology and chemical kinetics", In: Thermal Phenomena in Sedimentary Basins. Durand B. ed., Paris, Editions Technip, 235 - 246. [ Links ]

Ungerer, P., 1990. "State of the art of research in kinetic modeling of oil formation and expulsion". Org. Geochem,. 16 (1-3): 1-25. [ Links ]

Ungerer, P., 1993. "Modeling of petroleum generation and migration". In: Appl. Petroleum Geochem., Bordenave M.L. ed., Paris, Editions Technip, 395 - 442. [ Links ]

Ungerer, P., Espitalié, J., Marquis, E. and Durand, B., 1986. "Use of kinetic models of organic matter evolution for the reconstruction of palcotemperatures". Application to the Gironvillie Well (France). In: Thermal Modeling in Sedimentary Basins, Bnrrus J. ed., Paris, Editions Technip. 531-546. [ Links ]

Ungerer, P. and Pelet, R., 1987." Extrapolation of the kinetics of oil and gas formation from laboratory experiments to sedimentary basins", Nature, 327:52 - 54. [ Links ]

Ungerer, P., Burrus, J., Doligez, B., Chenet, P.Y. and Bessis. F., 1990. "Basin evaluation by integrated two-dimensional modeling of heat transfer, fluid How, hydrocarbon generation, and migration". AAPG Bulletin, 74 (3): 309-335. [ Links ]