CT&F Ciencia, Tecnología y Futuro
ISSN: 0122-5383
ECOPETROL S.A.
Colombia
González, Santiago; Idrobo, Eduardo A.
Caracterización dinámica de yacimientos estratigráficamente complejos usando algoritmos genéticos
CT&F Ciencia, Tecnología y Futuro, vol. 2, núm. 5, diciembre, 2004, pp. 23-51
ECOPETROL S.A.
Bucaramanga, Colombia
Disponible en: http://www.redalyc.org/articulo.oa?id=46520502
Cómo citar el artículo
Número completo
Más información del artículo
Página de la revista en redalyc.org
Sistema de Información Científica
Red de Revistas Científicas de América Latina, el Caribe, España y Portugal
Proyecto académico sin fines de lucro, desarrollado bajo la iniciativa de acceso abierto
23CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS
ESTRATIGRÁFICAMENTE COMPLEJOS USANDO
ALGORITMOS GENÉTICOSSantiago González*1 y Eduardo A. Idrobo*2
1Ecopetrol – Gerencia Centro Oriente, El Centro, B/bermeja, Santander, Colombia2Ecopetrol – Instituto Colombiano del Petróleo, A.A. 4185 Bucaramanga, Santander, Colombia
e-mail: [email protected] e-mail: [email protected]
(Recibido Junio 11 de 2003; Aceptado Septiembre 20 de 2004)
En este trabajo se presenta la aplicación a un caso real de una metodología novedosa para la caracte-rización de yacimientos altamente heterogéneos mediante la integración de la información dinámica del yacimiento al modelo estático actualizado. La metodología fue desarrollada por el Dr. Carlos
Romero, quien la aplicó a un yacimiento sintético.
El proceso de caracterización de yacimientos tiene como objetivo construir un modelo del yacimiento, lo mas realista posible, mediante la incorporación de toda la información disponible. El enfoque clásico consiste en elaborar un modelo que se fundamenta en la información estática del yacimiento, teniendo como etapa final del proceso la validación del modelo con la información dinámica disponible. Es importante aclarar que el término validación implica un proceso por naturaleza puntual, generalmente asociado a asegurar solamente la coherencia requerida entre zonas productoras y propiedades petrofísicas.
El objetivo de la metodología propuesta es mejorar la capacidad de predicción del modelo del yacimiento mediante la integración a priori de los parámetros intrínsecos a la dinamicidad de los fluidos del yacimiento mediante un proceso de inversión dinámica de datos a través de un procedimiento de optimización basado en computación evolutiva.
La metodología propuesta parte de la construcción del modelo estático del yacimiento de alta resolución, el cual es escalado mediante técnicas híbridas buscando preservar la heterogeneidad del yacimiento. Pos-teriormente, usando como marco de referencia un simulador analítico, el modelo escalado es metódica-mente modificado mediante un proceso de optimización que usa algoritmos genéticos y como información condicional los datos de producción.
Palabras claves: caracterización de yacimientos, modelamiento geoestadístico, simulación streamline, algoritmos genéticos, optimización global, ajuste histórico automático.
* A quien debe ser enviada la correspondencia
24 CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
25CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
El Producto final de este proceso es un modelo que respeta las condiciones estáticas y dinámicas del yaci-miento con capacidad de minimizar el impacto económico que genera el ajuste histórico de producción en las tareas de simulación. Este modelo final, entonces, presenta unas propiedades petrofísicas (porosidad, permeabilidad y saturación de agua) modificadas para lograr un mejor ajuste de la historia de producción simulada con la real (History Matching). Adicional a lo anterior, el proceso involucra una liguera modifica-ción de las permeabilidades relativas, las cuales se han cambiado buscando calibrar estas propiedades que igualmente presentan un grado alto de incertidumbre.
A novel methodology is presented in this paper for the characterization of highly heterogeneous oilfields by integration of the oilfields dynamic information to the static updated model.
The objective of the oilfield’s characterization process is to build an oilfield model, as realistic as possible, through the incorporation of all the available information. The classical approach consists in producing a model based in the oilfield’s static information, having as the process final stage the validation model with the dynamic information available. It is important to clarify that the term validation implies a punctual process by nature, generally intended to secure the required coherence between productive zones and petrophysical properties.
The objective of the proposed methodology is to enhance the prediction capacity of the oilfield’s model by previously integrating, parameters inherent to the oilfield’s fluid dynamics by a process of dynamic data inversion through an optimization procedure based on evolutionary computation.
The proposed methodology relies on the construction of the oilfield’s high resolution static model, escalated by means of hybrid techniques while aiming to preserve the oilfield’s heterogeneity. Afterwards, using an analytic simulator as reference, the scaled model is methodically modified by means of an optimization process that uses genetic algorithms and production data as conditional information.
The process’s final product is a model that observes the static and dynamic conditions of the oilfield with the capacity to minimize the economic impact that generates production historical adjustments to the simulation tasks. This final model features some petrophysical properties (porosity, permeability and water saturation), as modified to achieve a better adjustment of the simulated production’s history versus the real one (History Matching). Additionally, the process involves a slight modification of relative permeabilities, which have been changed to allow calibration of these properties that also feature a high degree of uncertainty.
Keywords: reservoir characterization, geostatistical modeling, streamline simulation, genetic algorithms, global optimization, automatic history matching
24 CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
25CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
INTRODUCCIÓN
El proceso de caracterización de yacimientos tiene como objetivo construir un modelo del yacimiento, lo más realista posible, mediante la incorporación de toda la información disponible. El enfoque clásico consiste en elaborar un modelo que se fundamenta en la información estática del yacimiento, teniendo como etapa final del proceso la validación del modelo con la información dinámica disponible.
La metodología aplicada fue, inicialmente desarrolla-da por el Dr. Carlos Romero (Romero, 2000) como tesis de doctorado; sin embargo, aquí fue ajustada mediante la incorporación de las siguientes modificaciones:
• Se trabajó en la definición de la distribución volu-métrica de las propiedades del yacimiento median-te técnicas geoestadísticas detalladas, utilizando definiciones litológicas (estratigrafía de detalle).
• Se mejoró la velocidad del método utilizando el concepto de simulación streamline.
• Se redefinieron:
- los cromosomas (variables) a optimizar (hubo una disminución, ya no se contemplaron las transmisi-bilidades de las fallas ni factor Skin en los pozos).
- la Función Objetivo.
- Calculo de la Aptitud.
- Estrategias de Selección.
- La forma como se aplican los operadores genéticos.
- El procedimiento de actualización de la población.
• Todo el proceso se programó en plataforma PC que utiliza como sistema operativo Microsoft Windows©. Se debe resaltar que el programa que permite la ejecución de la metodología se realizó exclusivamente para esta investigación (INVER-DIMANICA©), el cual se encarga de controlar (manejar) todas las aplicaciones adicionales, necesarias para automatizar el proceso de ajuste histórico. Estas aplicaciones son: Microsoft Ex-cel, S3D© (Simulador de Streamline) y SGSIM© (Simulador de Geoestadística). Para comprobar
los resultados obtenidos se utilizó el streamline (Simulador Numérico).
El objetivo de la metodología propuesta es mejorar la capacidad de predicción del modelo del yacimiento mediante la integración a priori de los parámetros in-trínsecos a la dinamicidad de los fluidos del yacimiento usando un proceso de inversión dinámica de datos a través de un procedimiento de optimización basado en computación evolutiva. Esta tesis se enfoca en la aplicación de los Algoritmos Genéticos (AGs) a la ca-racterización de yacimientos de hidrocarburos.
La metodología utiliza un AG especialmente di-señado para buscar descripciones del yacimiento que hacen honor tanto a los datos de los pozos como a sus historias de producción, siempre y cuando se tengan los modelos geológico, estratigráfico y petrofísico previamente definidos. Este AG se combina con el modelamiento geoestadístico para generar múltiples descripciones del yacimiento. El resultado es un conjunto de posibles soluciones (correspondientes a diferentes realizaciones del modelo del yacimiento), a partir de las cuales se puede seleccionar la más apta para análisis posteriores.
Los yacimientos de hidrocarburos están compuestos, por lo general, de rocas sedimentarias heterogéneas que contienen dependiendo del caso: aceite, gas y agua. Para conocer cómo se comportan estos fluidos en el yacimiento, se utilizan los simuladores. Para su aplicación, el yacimiento debe estar dividido en celdas discretas que conforman una gran malla, por lo gene-ral, las propiedades en estas celdas son desconocidas. Posteriormente, se debe llevar a cabo un proceso de calibración en el cual los resultados son comparados con la información histórica del yacimiento. Esta cali-bración la hace el ingeniero mediante la variación de algunas de las propiedades de las celdas hasta obtener un ajuste adecuado.
Esta metodología automatiza este proceso, mini-mizando el impacto económico que genera el ajuste histórico de la producción en las tareas de simulación y el modelo resultante incorpora toda la información disponible (estática y dinámica), por lo tanto será el más realista posible. Esta automatización se logra al consi-derar el ajuste histórico como un proceso de optimiza-ción multivariable (el cual puede tener varios óptimos). Teniendo en cuenta esta última reflexión, se impone
26
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
27CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
como requerimiento un método de optimización global, dentro de los cuales, los AGs han demostrado una alta eficiencia en la solución de este tipo de problemas.
Los AGs usan ideas de la teoría de evolución de las especies de Darwing, donde sólo los organismos (individuos) más fuertes sobreviven adaptándose a las diferentes condiciones que el medio les imponga (más aptos). Estos organismos tienen una información genética que los hace únicos y, por lo tanto, cambios pequeños en estos códigos los hacen más fuertes o por el contrario más vulnerables.
Los resultados obtenidos superan las expectativas que se tenían en relación con esta metodología, la cual, sólo había sido probada en un yacimiento sintético, proveyendo un procedimiento que se puede aplicar aún cuando se tenga una alta incertidumbre en cuanto a parámetros tales como: propiedades petrofísicas de-terminadas con registros “viejos”, datos de producción e inyección de agua. Otra ventaja que se tiene es que se puede utilizar en yacimientos tanto en producción primaria como en recobro secundario; sin embargo, los tiempos de cómputo serán mayores para los casos de campos en producción primaria.
Todo el programa corre en plataforma PC ambiente Microsoft Windows©, lo cual le permite sacar provecho de: uso en cualquier equipo, utilización de la última tecnología en PC para acortar los tiempos de cómputo y menores costos en relación con los equipos en pla-taforma UNIX.
DESCRIPCIÓN PROCESO CARACTERIZACIÓN
A continuación se hace una presentación y descrip-ción general de la secuencia que realiza el programa INVERDINAMICA ©. El software fue desarrollado, en particular, para automatizar el proceso de ajuste históri-co de un modelo de yacimiento en plataforma PC.
Cada corrida involucra: cargar los datos básicos del yacimiento, hacer simulaciones de flujo de fluidos (streamline), definir los puntos piloto, aplicar los AGs (optimización), hacer simulaciones Geoestadísticas (SGSIM) para cada propiedad y actualizar la pobla-ción de individuos. La Figura 1 presenta un diagrama que resume toda la secuencia que hace el programa
Figura 1. Esquema general metodología propuesta
FIN
algorítmo genético
permebilidadesrelativas
propiedadespetrofísicas
parámetrosgeoestadística
parámetrossimulación
procesosimulación
simulcióngeoestadística NO
evaluarfunciónobjetivo
puntospiloto
SI
26
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
27CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
durante su ejecución (Romero, 2000; Romero etal., 2000).
Simulación StreamlinePara definir la “estratificación social” de cada in-
dividuo y, por ende, los operadores genéticos a aplicar, es necesario determinar la Función Objetivo (Ecuación 1) (Baker, 1985) y para esto, es necesario tener una producción de aceite simulada, esta curva se obtiene utilizando un simulador streamline, el cual, es muy conveniente para esta metodología debido a su gran velocidad de convergencia (Batycky, 1997; Batycky et al., 1997; Caers, et al., 2001; Wang et al., 2000).
Los archivos “Grid_Pozos.prn”, “Field.prn” y “Production.prn”, necesarios para la ejecución del simulador, no cambian durante el proceso de ajuste histórico; por el contrario, los archivos “Files.prn” y
“Fluid.prn” deben ser actualizados para cada individuo (realización), ya que las propiedades roca - fluido (per-meabilidades relativas), son cromosomas dinámicos y en cada generación van siendo modificados por el AG (Deutsch et al., 1998).
Los datos de las propiedades se encuentran organi-zados en archivos texto, existe un archivo de este tipo por propiedad y por individuo, en total se tienen 60 archivos para las tres propiedades petrofísicas conside-radas y los 20 individuos. Adicionalmente, los puntos finales de las permeabilidades relativas se encuentran en otro archivo plano.
El programa diseñado hace uso de la dependencia que existe entre las permeabilidades y sólo se alma-cenan los datos de permeabilidad horizontal en X, los datos de permeabilidad en Y y las permeabilidades en Z, son calculados por el simulador, utilizando la infor-mación de anisotropías (Ky/Kx y Kz/Kx) suministrada (Dake, 1978 ; Renard et al., 1997).
Definición puntos pilotoLa metodología utilizada para disminuir la cantidad
de datos a manejar en el AG (upscaling), se denomina Método de los puntos piloto. Los puntos piloto deben considerarse como pseudopozos y se utilizan como los parámetros a variar durante el desarrollo del proceso de optimización con AGs, tal como se describe en la Figura 1. Estos pseudopozos y los pozos son utilizados como datos de entrada al Simulador de Geoestadística,
este a su vez se encarga de asignarle valores a todas las celdas que conforman el grid de simulación numérico. Estos puntos serán ubicados por coordenadas (X, Y), de forma análoga a las localizaciones de los pozos ver-daderos dentro de un grid de simulación. De acuerdo con el conocimiento heurístico de la geoestadística, se determinó que la cantidad de pseudopozos a utilizar debía ser equivalente al 5% de las celdas que conforman el enmallado (Bissell et al., 1987)
Una vez determinadas las coordenadas (X, Y) para todos los puntos piloto, se deben formar arreglos bi-dimensionales –uno por cada propiedad petrofísica en consideración– (Figura 2), donde se almacena tanto los pseudopozos como los pozos y su respectiva informa-ción petrofísica disponible para las celdas en dirección Z que componen la coordenada (X, Y) seleccionada (Cuypers et al., 1998)
1
2
3
4
5
6
7
8
9
10
11 12
Capas
puntos piloto
1 2 3 4 5 6 7 8 9 10
Figura 2. Arreglo bidimensional puntos piloto
Resumiendo, durante la definición de los puntos piloto el programa se encarga de: calcular la cantidad de celdas que equivalen al 5% del enmallado, seleccionar aleatoriamente las celdas que serán utilizadas, crear arreglos bidimensionales –uno por cada propiedad petrofísica– (Figura 2) e ir a los archivos planos para cargar los datos de cada celda.
Formulación del algoritmo genético modificadoLa selección de esta metodología de optimización
se hizo teniendo en cuenta las características que hacen de los AGs una técnica conveniente para realizar el ajuste histórico automático en yacimientos de hidro-carburos, una vez este se formula como un problema de optimización (Chen et al., 1974; Davis, 1991; Yang et al., 1987).
28
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
29CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
La primera característica es que los AGs pueden ser usados eficazmente para la optimización de fun-ciones “ruidosas”, porque no requieren información adicional (tales como derivadas) de la función que va a ser maximizada o minimizada (Mitchell, 1996; Stoffa et al., 1991; Thomas et al., 1971).
La segunda característica esta muy ligada a la pri-mera. Porque los AGs usan búsquedas aleatorias, por lo tanto no siguen reglas determinísticas y pueden escapar a mínimos locales. Esta característica les da el potencial para alcanzar óptimos globales (Deschamps et al., 1998; Gomez et al., 1999; Goldberg, 1989).
La capacidad para manejar muchos parámetros es otra característica importante de los AGs. Dependiendo de la complejidad y nivel de detalle necesario, una si-mulación completa de un yacimiento puede requerir el cálculo de muchos parámetros, generalmente del orden de miles o algunas veces de cientos de miles (Bush et al., 1996; Zalzala et al., 1997).
La cuarta característica de los AGs son sus atributos de paralelismo. Por una parte, los AGs tienen un para-lelismo implícito porque ellos trabajan con una pobla-ción de soluciones y no con una sola solución. De esta manera, varias soluciones propuestas son evaluadas y combinadas durante el proceso de optimización. Adi-cionalmente, las soluciones propuestas son basadas en la evaluación de soluciones previas, en otros términos, los AGs hacen uso de todas las experiencias exitosas anteriores, gracias a la llamada “memoria implícita”. Esto es particularmente importante cuando una función es muy costosa para ser evaluada, como es el caso del ajuste histórico en los yacimientos. De otro lado los AGs tienen el potencial de retornar un conjunto de posibles soluciones, es decir, la población final puede ser usada para evaluar parámetros y predecir la incer-tidumbre (Ounes, 1992; Sen et al., 1995).
Estructura del genoma y el genotipoEl primer paso en la definición de la estructura de
un genoma es la selección de las variables que van a ser incluidas en la optimización. Estrictamente, el genoma debe ser dividido en partes: el genoma dinámico, con las estructuras que cambiarán durante el proceso de optimización y el genoma estático, con la información que contribuirá a la conformación de los modelos nu-méricos pero permanecerá inalterable a lo largo de la
optimización y puede evolucionar después. (Goldberg et al., 1991; Holland, 1975).
Los modelos numéricos de yacimientos requieren especificar un gran número de variables. Para este caso en particular, se han seleccionado las siguientes variables para ser incluidas en la optimización: propie-dades petrofísicas (permeabilidad horizontal, porosidad y saturación de agua) tanto en los puntos piloto como en los pozos; parámetros que controlan el proceso geoes-tadístico (rango en la dirección de máxima correlación, rango en la dirección de mínima correlación y rango de correlación vertical) y permeabilidad relativa en los puntos finales.
La información restante no incluida en el genoma dinámico, tal como: propiedades de los fluidos (PVT), estructura geológica, contacto de fluidos, etc., que es necesaria para una descripción completa del modelo de simulación del yacimiento, esta contenido en el genoma estático y permanece ahí a lo largo de todo el proceso. De aquí en adelante, el término “genoma” hará referencia al genoma dinámico.
El genotipo solo provee la información a partir de la cual el fenotipo es construido. Este tipo de mapeo es conocido como “Embriogenia” porque el genotipo solo establece las reglas para el desarrollo de las carac-terísticas fenotípicas, las cuales son producidas por la interacción de varios genes. En este caso, el genotipo contiene las propiedades petrofísicas de los bloques de los puntos piloto y de pozos, además de parámetros para el control del proceso geoestadístico y permeabi-lidades relativas en los puntos finales. Cada modelo se construye bajo la interacción de estos parámetros. Un solo genotipo puede producir una variedad de fenotipos, dependiendo de los elementos aleatorios de los métodos geoestadísticos (Bäck, 1996).
En este trabajo, un genoma con estructura no es-tándar fue seleccionado por las razones ya descritas anteriormente. Las variables han sido divididas en ocho grupos, y a cada grupo de variables tridimensionales, se le asignó un cromosoma separado (porosidad, per-meabilidad y saturación de agua). Cromosomas indi-viduales, para cada una de las variables que componen los grupos de parámetros restantes, fueron construidos: parámetros geoestadísticos y permeabilidades relativas para los puntos finales (Figura 3).
28
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
29CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Cromosomas de propiedades petrofísicasLos tres cromosomas de propiedades petrofísicas
tienen la estructura que se presenta en la Figura 4. Estos arreglos tridimensionales presentan inserciones para los pozos. Con esta estructura, los elementos vecinos en los arreglos tienen la misma relación espacial a los puntos piloto equivalentes en el fenotipo (Palazzari et al., 1999).
Los valores de las propiedades petrofísicas son re-presentados por números reales (Surry et al., 1997), en lugar de formas binarias. Esto se hace para reducir la complejidad del manejo de los datos (Wright, 1991).
Cromosomas misceláneosLas variables misceláneas están contenidas en
cinco cromosomas unidimensionales. Cada variable es codificada dentro de una cadena de bits binarios con una longitud determinada (Calvete et al., 2001; Martínez, 1990).
Para este trabajo, el número de bits osciló entre cuatro y siete por parámetro, de esta manera las va-riables fueron discretizadas entre 16 y 128 intervalos, dependiendo de la precisión que se les quisiera dar. Las cadenas binarias son estándar durante todo el desarrollo de la metodología (Whitley, 1999).
Propiedades Petrofísicas
Geoestadística
Permeabilidades Relativas
K SwF
P P P
P P P
P P P
W
W
WP P P
P P P
P P P
W
W
WP P P
P P P
P P P
W
W
W
Figura 3. Estructura del genoma no-estándar (Romero, 2000)
1
2
3
4
5
6
7
8
9
10
11 12
Capas
puntos piloto
1 2 3 4 5 6 7 8 9 10
Figura 4. Arreglo tridimensional para los puntos piloto (Romero, 2000)
30
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
31CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
forma, este modelo corregido podrá ser utilizado en la predicción de la producción.
El ajuste del modelo se hace optimizando una fun-ción objetivo, la cual fue definida como:
El cromosoma para los datos geoestadísticos contiene tres parámetros: rango en la dirección de máxima correlación, rango en la dirección de mínima correlación y rango de correlación vertical. Para este trabajo, un conjunto de sencillo de parámetros contro-la las simulaciones geoestadísticas en todas las capas (Jensen et al., 1997).
El cromosoma que incluye los parámetros de las permeabilidades relativas en los puntos finales, fue codificado usando números binarios. Este cromosoma se incluyó en el estudio para dar integridad al modelo, aunque su variabilidad se mantuvo al mínimo.
Habiendo definido la estructura del genoma y el procedimiento para codificar las variables dentro de los cromosomas, otro conjunto de estructuras deben ser definidas: las estrategias de selección y los opera-dores genéticos. La Figura 3 muestra una ilustración simplificada del genoma no-estándar, implementado en este trabajo.
Estrategias de la poblaciónLos AGs trabajan con una población de soluciones
potenciales individuales del problema que se está con-siderando. A medida que el proceso de optimización avanza, la población actual va siendo reemplazada progresivamente por nuevas soluciones. La manera exacta como se logra esto depende de la formulación de las estrategias de la población (Baker, 1985).
Las estrategias de población requieren la definición de tres aspectos principales: la evaluación de la aptitud de cada individuo es decir, obtener una medida de que tan conveniente es un modelo particular (realización) como solución del problema que se está considerando; el procedimiento usado para asignarle descendencia a cada individuo, según la última aptitud calculada; y cómo los individuos entran o salen de la población (actualización) (De Jong, 1975).
Función objetivo
La metodología busca ajustar las propiedades petrofísicas del yacimiento para que este, al ser si-mulado, reproduzca la curva de producción real del campo, obteniendo de esta forma un modelo que ha sido ajustado automáticamente (sin intervención del ingeniero de yacimientos) y que reproduce lo más pa-recido posible el comportamiento del campo. De esta
��
��
�� ��������å -=�� (1)
Donde:
FO : es el valor de la Función Objetivo.
Qosi : es la producción simulada de aceite prome-dio para un mes i.
Qori : es la producción real de aceite promedio para un mes i.
NT : número total de meses considerados.
El AG busca disminuir hasta su mínima expre-sión el valor de FO, haciendo ajustes a las variables representadas por los cromosomas, para este caso en particular son: porosidad, permeabilidad, satura-ción de agua, permeabilidades relativas y paráme-tros geoestadísticos.
Los cromosomas de porosidad, permeabilidad y saturación de agua, se les aplica una metodología denominada puntos piloto mediante la cual se hace un muestreo de valores aleatorio y estos son los datos que se van a ajustar. Posteriormente, las propiedades petrofísicas son nuevamente distribuidas en todo el modelo mediante una simulación geoestadística, el modelo resultante nuevamente es simulado y se repite el proceso.
Cálculo de la aptitud
En los AGs, la aptitud de un individuo es un indicativo de qué tan bueno es, en particular, para resolver el problema, en comparación con los otros individuos de la población. Para aplicar los AGs en un proceso de ajuste histórico en un yacimiento de hidrocarburos con datos de producción, valores de aptitud deben ser asignados a cada realización (mo-delo). Estos valores deben reflejar el grado, en que una realización particular puede reproducir la historia de producción. Esto se hace comparando los datos de producción reales con los datos obtenidos por medio de una simulación numérica.
30
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
31CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Los datos generados por la simulación geoesta-dística (H3D) son convertidos (arreglados), por una subrutina del programa, al formato de streamline, gene-rando los diferentes archivos de entrada que contienen los datos de porosidad, permeabilidad y saturación de agua de cada celda del grid utilizado por el programa de simulación. Este paso es de conversión de archivos y se requiere porque los sistemas de coordenadas de los programas de simulación geoestadística y de flujo son diferentes.
El paso final del proceso es el cálculo de la aptitud, el cual se realiza después de que los resultados de la producción simulada están disponibles. Este proceso se realiza internamente en el programa y utiliza como base la Función Objetivo definida por la Ecuación 1. Una vez, se han obtenido los valores de la función objetivo para los 20 individuos, se comparan entre sí y se selec-ciona el valor más bajo, luego mediante la Ecuación 2 se calcula la aptitud para cada individuo.
baja para que puedan encontrarse nuevos individuos con buenas aptitudes, mientras que la fase de explota-ción requiere que la presión de selección sea alta para que los individuos con buenas aptitudes se acentúen en espera que ellos produzcan descendencia con apti-tudes mejores. Debido a la complejidad del espacio de búsqueda a la exploración se le da la prioridad aunque no debe descuidarse la explotación.
La estrategia de selección implementada en el desarrollo de este proyecto, después de evaluar todas las estrategias propuestas, fue el elitismo combinado con el clasismo; esta estrategia mixta que consiste en clasificar los individuos de acuerdo a su aptitud, creando, por decirlo así, una realeza, una clase media y una clase baja. De acuerdo a su clasificación se aplican los diferentes operadores genéticos (Mitchell et al., 1992).
Esta clasificación se hace mediante el cálculo de una variable que indica la diferencia (amplitud) máxi-ma que existe entre los valores de aptitud calculados para una misma generación, para tal fin se utiliza la Ecuación 3:
Donde:
Aptitud (ni, ng): es el valor de la aptitud de un indi-viduo ni en la generación ng
FO (ni, ng): es el valor de la Función Objetivo de un individuo ni en la generación ng
FOmin (ng): es el valor de la Función Objetivo míni-mo en la generación ng
Estrategias de selección
La estrategia de la selección es uno de los factores principales que controlan el desempeño de la optimiza-ción. La selección apunta a dar las preferencias a aque-llas realizaciones con aptitudes cercanas a la unidad (1), para pasar sus atributos a otras realizaciones.
Una estrategia que proporcione una óptima presión de selección, es crítica para el buen desempeño del AG propuesto. La selección debe favorecer la realización con mejor aptitud, sin causar una convergencia prema-tura. Una presión de selección óptima debe mantener el equilibrio entre la exploración y la explotación. La fase de exploración requiere que la presión de selección sea
Donde:
Dif Aptitud (ng): es valor de máxima diferencia (amplitud) en la generación ng
Aptitudmax (ng): es el valor máximo de la Aptitud en la generación ng
Al valor máximo obtenido se le aplican los porcentajes de elitismo, los cuales son definidos por el usuario depen-diendo del comportamiento del proceso; en el desarrollo de este trabajo se utilizaron los siguientes valores:
Reproducción : 10%
Cruzamiento : 70%
Mutación : 20%
Los valores de evaluación del clasismo de los indi-viduos se calculan de utilizando la Ecuación 4:
���
���������� ������������������������� =���������� (4)
������
������������
���
=����������� (2)
��������������� ���=������������ (3)
32
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
33CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Donde:
Clasismo (ng): Determina el límite para una deter-minada clase social en una generación ng
Porcentaje Elitismo: Valor máximo permisible para una clase social
Al valor obtenido para cada una de las clases sociales se le debe sumar uno (1) y de esta forma se determinan los límites superiores. Esto significa que los individuos con aptitudes que se encuentren entre 1 y el 10% de desviación se seleccionarán para reproducción; los individuos con aptitudes que se encuentren entre el 10% y el 80% de desviación, se seleccionarán para cruzamiento; y los individuos con aptitudes superiores al 80% de desviación serán mutados.
De esta forma, todos los individuos que com-ponen una generación se estratifican y, simultánea-mente, quedan definidos los operadores genéticos que se les aplicarán dependiendo de la clase social a la cual pertenezcan.
Operadores genéticos
Se han creado operadores de cruzamiento no están-dar para trabajar con los cromosomas tridimensionales de las propiedades petrofísicas. Existen pocas pautas acerca de cómo estos operadores deberían recombinar sus genes, sin embargo se optó por un sistema sencillo que más adelante se explicará.
Los cromosomas normales fueron combinados apli-cando las normas de los operadores genéticos binarios.
Reproducción
En este trabajo se definió la reproducción como la acción de “clonar” un individuo de la generación actual, en la siguiente, logrando de esta forma preservar una memoria implícita respecto a la mejor solución lograda y, de paso, disminuir el tiempo de cómputo requerido, ya que a este individuo no se le debe hacer ni simulación geoestadística ni numérica.
Solo el individuo con aptitud igual a uno (1) será reproducido en cada generación, es decir, solo un indi-viduo tendrá esta condición y será denominado como
“rey” dentro de la estratificación clasista diseñada para aplicar el AG.
Sin embargo, en una misma generación uno o más individuos pueden ser clasificados para ser reproducidos, estos individuos cuya aptitud es buena pero no mejor que la del individuo “rey” serán denominados realeza y estos se cruzarán solo con el individuo rey buscando mejorar la aptitud de los individuos descendientes. Este cruza-miento se hace siguiendo las mismas reglas descritas en el párrafo Cruzamiento (Mitchell et al., 1992).
Cruzamiento
Dentro del grupo de individuos seleccionados para cruzamiento se verifica inicialmente que sea un grupo par, es decir, que todos los individuos tengan un com-pañero de cruce. En caso, de no ser así, el individuo con mejor aptitud dentro de los clasificados para mutación se convierte a individuo de combinación o cruzamiento, completando de esta forma el número par requerido.
El programa se encarga, aleatoriamente, de formar parejas con los individuos disponibles para cruza-miento, los individuos tendrán una y sólo una pareja, y su descendencia serán otros dos individuos, donde se combinan las características genéticas de los dos individuos padres.
En los cromosomas tridimensionales (propiedades petrofísicos), una vez escogidas las parejas y habiendo seleccionado los puntos piloto previamente, se forman dos arreglos bidimensionales; posteriormente, se selec-ciona aleatoriamente un punto de corte (capa), punto por donde se combinarán los dos padres, tal como se muestra en la Figura 5, de esta forma la descendencia estará compuesta nuevamente por individuos que tienen información genética de ambos padres.
En los cromosomas binarios, el cruzamiento se hará de la forma tradicional en los AGs, haciendo referencia explícitamente al método del punto sencillo; método finalmente seleccionado para aplicar en este trabajo.
Mutación
En este trabajo solo una clase de mutación fue se-leccionada y es aplicada sobre el individuo padre y su resultado es un individuo en el cual se modifica uno de sus genes; el tipo de mutación utilizado se denomina de salto.
En los cromosomas binarios el procedimiento es muy sencillo y consiste en cambiar, aleatoriamente un
32
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
33CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
bit de la cadena que compone el cromosoma (0 por 1 o 1 por 0).
En los cromosomas tridimensionales se aplica un procedimiento similar, pero se tiene el cuidado de no cambiar (alterar) las celdas que representan las arcillas, la Figura 6 muestra lo que ocurre en un individuo selec-cionado para mutación, en una generación dada.
Las celdas correspondientes a facies de arcillas pueden ser fácilmente reconocidas, ya que los valo-res de porosidad y permeabilidad fueron ajustados a 0.0001, y en el caso de la saturación de agua, su valor siempre se encuentra por encima del 90%. Cuando se esta ejecutando la subrutina de mutación y el programa encuentra una de estas celdas, inmediatamente salta ese punto piloto y continua con el siguiente.
Figura 6. Mutación cromosomas tridimensionales
1
2
3
4
5
6
7
8
10
11 12
Capas
puntos piloto
1 2 3 4 5 6 7 8 9 10
Celdas Seleccionadas Mutación
INDIVIDUO "A"
Punto PilotoNo Modificado
1
2
3
4
5
6
7
8
9
10
11 12
Capas
puntos piloto
1 2 3 4 5 6 7 8 9 10
1
2
3
4
5
6
7
8
9
10
11 12
Capas
puntos piloto
1 2 3 4 5 6 7 8 9 10
Capa Corte Seleccionada
INDIVIDUO "A" INDIVIDUO "B"
Figura 5. Combinación cromosomas tridimensionales
34
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
35CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
milla” a utilizar en la simulación, rango en la máxima dirección de correlación, rango en la mínima dirección de correlación, rango vertical de correlación, ángulos de anisotropía, tipo de estructura del variograma, tipo de Kriging, nugget y sill.
Tal como se puede evidenciar, los nombres de los archivos de entrada y salida, y la información prove-niente del cromosoma de parámetros geoestadísticos, son los únicos datos que se actualizan constantemente, los demás se mantienen constantes.
De acuerdo con el análisis de anisotropía hecho, el ángulo de la elipse de búsqueda es de 92° (358° en azimut); el tipo de estructura del variograma seleccio-nado fue el modelo exponencial; el tipo de Kriging utilizado fue Ordinary Kriging (OK); el nugget se definió como cero (0) y el sill es uno (1) por defecto, ya que el programa normaliza las variables de entrada (Jensen et al., 1997).
Una vez se han generado los archivos texto de geoestadística, estos deben ser convertidos a formato streamline para que puedan ser utilizados por el si-mulador en la siguiente generación. Adicionalmente, se debe tener presente que algunas celdas del modelo total se encuentran inactivas y por lo tanto estas de-ben ser nuevamente deshabilitadas. Los archivos acá generados reemplazan a los existentes, lo que facilita que el proceso pueda ser completamente cíclico (WU et al., 1998).
CASO DE CAMPO: AREA 07, ZONA C, CAMPO LA CIRA (COLOMBIA)
El campo La Cira-Infantas esta localizado en la cuenca del Valle Medio del Magdalena (VMM) en Colombia (Figura 7). Se encuentra ubicado en la parte central de la antigua Concesión De Mares, al Este del río Magdalena y al Sur del río Sogamoso, abarcando un área de aproximadamente 160 km²; a una distancia aproximada de 22 km al SE de la Ciudad de Barran-cabermeja y a aproximadamente 250 km al NNE de Bogotá (Figura 7). El VMM esta limitado al Este por la Cordillera Oriental, al Oeste por la Cordillera Central, y se extiende desde las localidades del Banco (Magdalena) en el Norte hasta Honda (Tolima) en el Sur (Dengo et al., 1993).
En las celdas que se pueden mutar, se selecciona un número aleatorio, que se encuentra entre el mínimo y máximo definido para la propiedad petrofísica en evalua-ción (porosidad, permeabilidad y saturación de agua).
Actualización de la poblaciónDentro de un AG, una población de soluciones se
actualiza constantemente, por norma, las nuevas solu-ciones reemplazan aquellas cuya aptitud tenga que ser evaluada. Una estrategia de actualización ideal debe poder mantener la diversidad en la población, mien-tras se va acercando a un óptimo global. El número de generaciones debe ser suficientemente grande para permitir que la población evolucione hasta acercarse al nivel de aptitud requerido. En este trabajo, sin embargo, el número de evaluaciones de la función objetivo se ha restringido a menos de 600 debido a la limitación de recursos computacionales.
Un enfoque generacional fue usado, es decir, la población se reemplazaba casi en su totalidad en cada nueva generación, con excepción del individuo con mejor aptitud, el cual se “clonaba” en la siguiente ge-neración. Los parámetros de AG típicos en este trabajo son: el tamaño de la población es de 20 individuos y el número máximo de generaciones fue de 20; aunque se hicieron pruebas con 25 y 30 generaciones sin encon-trar mejorías significativas, mientras que el consumo de tiempo máquina si aumento considerablemente. Al preservar el mejor individuo de una generación a la otra, se hizo una reducción del 5% en el tiempo de cómputo, ya que no fue necesario hacer 400 evalua-ciones (en el caso de 20 individuos y 20 generaciones) sino solamente 380.
Simulación gaussiana secuencialUna vez generados los archivos de entrada al
Simulador Geoestadístico (SGSIM) (Deutsch et al., 1998), es necesario leer el archivo texto que contiene los parámetros geoestadísticos considerados y generar un número aleatorio impar de cinco cifras.
Estos datos se utilizan para generar el archivo llamado “SGSIM.prn”, el cual contiene la siguiente información: nombres y rutas para los archivos de en-trada de datos y salida de resultados, valores mínimo y máximo permitidos, número de celdas en todas las direcciones (X, Y y Z) que componen el modelo, “se-
34
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
35CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Campo La Cira: Esta estructura consiste de un an-ticlinal en forma de domo alargado en dirección N-S al Noroeste del Anticlinal de Infantas, con cabeceo (plun-ge) hacia el Norte y hacia el Sur y está cortada, al Este del eje del Anticlinal, por la falla La Cira, que es una falla inversa de alto ángulo (70°-80°) que buza hacia el Occidente y cuyo desplazamiento oscila entre 200 y 300 pies en el tope de la estructura. Hacia el Este de la falla La Cira la estructura se hunde un poco en forma de “U” a manera de sinclinal entre esta falla y la falla de Infantas. De acuerdo con los datos de producción del campo se puede considerar que la falla La Cira actúa como una barrera efectiva entre el flanco Occidental y Oriental de la estructura. Existe en este campo un gran número de fallas transversales al eje de la estructura las cuales son normales con buzamientos entre 80° y 90° y desplazamientos de hasta 250 pies y que dividen al campo en varios bloques. Esta estructura tiene una longitud aproximada de nueve km y un ancho de seis km . Los buzamientos de la estructura hacia el Norte y el Oeste son de aproximadamente 10° a 15°.
El Campo La Cira produce de tres zonas “A”, “B” y “C” (Figura 8). La zona “C” es un yacimiento de arenas fluviales de permeabilidad baja a moderada. Esta zona contiene aproximadamente el 80% del OOIP en el Campo La Cira. Las zonas “C” y “B” producen de la Formación Mugrosa de edad Eoceno – Oligoceno y la zona “A” produce de la formación Colorado de edad Oligoceno.
En la actualidad se han perforado en el área de la Cira-Infantas 1716 pozos, de los cuales hay aproximadamente 580 pozos productores activos, 81 pozos inyectores de agua activos, 592 pozos productores inactivos, 106 po-zos inyectores inactivos y 357 pozos abandonados. El OOIP del campo la Cira-Infantas es de 3594 Mbls de 24° API. La producción promedio de aceite actual es de 5700 BOPD (Diciembre de 2002), con una producción de aceite acumulada de 798,15 Mbls y una inyección acumulada de 957,4 Mbls (Diciembre de 2001).
En el área 07, la zona C inició su explotación en 1933 con el pozo LC-745 con 1132 bopd. El principal
Pavas
BonanzaCristalinaLas Garzas
Cantagallo - Yarigui
Provincia
SalinasAguas ClarasSogamoso
LlanitoGala Lisama
Tesoro
Nutria
Peroles
GalánCasabe
La CiraInfantas
San Luis
Tenerife
Aguas Blancas
Colorado
Peñas Blancas
Convenciones
CamposPoblaciones
San Pablo
Puerto Wilches
Sabana de Torres
Barrancabermeja
Yondo
K 0 10 20 K
1000000
10500001000000
12500001300000
12500001300000
1050000
Fl or enci a
Villavicencio
Nei va
Tumaco
BuenaventuraIbagué
M aniz ales
Cúcu t a
Valledupar
Pasto
Cali
Mede llín
Santa MartaBarranquillaCartagena
Bogota
Caribbean Sea
Paci f ic
Ocean
BRAZIL
VENEZUELA
PERU
ECUADOR
Per eir a
M ont er í a
R í ohacha
Figura 7. Mapa ubicación campo La Cira (ECP, 1999)
36
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
37CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
desarrollo de esta área se llevó a cabo entre los años 1933 y 1941 cuando se perforaron 124 pozos, de los cuales 99 penetraron la zona C. La inyección se inició en 1971 con modelos regulares de seis y siete puntos, modificados posteriormente a modelos de cuatro a siete puntos, por problemas mecánicos de los pozos.
Para realizar el modelo estratigráfico y el mode-lo petrofísico se seleccionaron 65 pozos, los cuales se pueden observar en la Figura 9. Algunos de los criterios que se tuvieron a la hora de realizar esta selección fueron:
• Pozos Corazonados (Carrillo et al., 1990).
Figura 8. Columna estratigráfica general para el VMM y el campo La Cira (ECP, 1999)
Esmeraldas
ZonaD
400
500
600
700
800
900
1000
1100
1200
1300
1400
1500
1600
1700
1800
1900
2000
2100
2200
2300
2400
2500
2600
2700
2800
2900
3000
3100
3200
3300
3400
3500
3600
3700
300
3800
Arenas 58
Arenas 116
Arenas 244
Arenas 13
Arenas 13 L
Arenas La Cira
Arenas 47
Arenas C1
Arenas C2
Arenas C3
Arenas C4
Arenas D
AguaArenas 157
UA3
UA4a
UA4b
A1
AA
A1a
A2
A2a
A3A3a
A4
A/B
B1
B1aB1bB1cB1dB2
B2aB2b
B3B3a
B3bB3c
B4B4a
B4b
C1C1BtC1AbC2C2DtGtbC3
C3cbC4C4AbC4CbC4EtC5D1
D2D3D4D5
FormaciónColorado
FormaciónMugrosa
ZonaA
ZonaB
ZonaC
Período
Cuatern.
Terciario
Cretáceo
Jurásico
Premesozoico
Época
Pleistoceno
Plioceno
Mioceno
Oligoceno
Eoceno
Paleoceno
Maestrichtiano
Campaniano
Santoniano
ConiacianoTuroniano
Cenomaniano
Albiano
Barremiano
Valanginiano
Formación
Mesa
Real
Colorado
Mugrosa
Esmeraldas
La Paz
DiscordanciaEoceno-Paleoceno
Lisama
Umir
LaLuna
Galembo
Pujamana
Salada
Simití
TablazoPaja
Rosablanca
Tambor
Girón
Basamento
Litología
Aptiano
• Pozos perforados más recientes (calidad registros eléctricos).
• Que tuvieran una buena distribución areal (cober-tura total del área).
• Que estuviera representada la geomorfología de las unidades definidas.
Modelo sedimentarioPara la definición del modelo sedimentario y dia-
genético se contó con la información litológica (cora-zones) del pozo La Cira 1884, perforado y corazonado en 1988 como reemplazo del pozo LC–1133. En la
36
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
37CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
vertical entre las diferentes unidades de la zona “C” (Bernal et al., 1992).
Rellenos de canales abandonados
Estas facies limo-arcillosas se caracterizan por pre-sentar la base del canal, normalmente de poco espesor, suprayacida por una secuencia de material limoso o arcilloso que representa el gradual o súbito abandono del canal (Bernal et al., 1992).
Lodolitas de llanura de inundación
Conformada por arcillolitas varicoloreadas (gris verdosa, rojiza, amarillenta, etc.) y por limolitas gris rojizas con restos de raíces, indicios de exposición subaérea y sectores bioturbados rellenos de arenas calcáreas (Bernal et al., 1992).
Arenas finas y limos de borde de canal (Natural Levee)
Constituyen los umbrales de los canales y consisten de material clástico fino con espesores que varían entre 2 y 10 pies aproximadamente (Bernal et al., 1992).
Arenas finas y limos de desborde de canal (Cre-vasse Splay)
Consiste de arenas finas, normalmente con tenden-cia granocreciente al tope, espesores de 4 a 10 pies aproximadamente (Bernal et al., 1992).
Modelo EstratigráficoLos cambios en las propiedades petrofísicas con
respecto a la posición estratigráfica son proporcionales al grado de preservación de facies idénticas deposita-das bajo condiciones variables de A/S (acomodación y suministro de sedimentos). En areniscas con estratifica-ción cruzada, los cambios en porosidad y permeabilidad son exactamente coincidentes con variaciones sutiles en los atributos que son función de la posición estrati-gráfica (Ramón et al., 1997).
Teniendo en cuenta las anteriores consideraciones, se procedió a clasificar cada una de las arenas presen-tes en la zona “C” (objeto de este estudio) de acuerdo con su posición estratigráfica, clasificando estas de acuerdo a su unidad genética (depósitos de Canales o depósitos de Llanuras de Inundación). Las “barras de arenas” y los apilamientos de “barras de arenas” se
Figura 9. Área descripción estratigráfica
descripción de los 350 pies de corazón del pozo LC-1884, que incluye la mayoría de las arenas y arcillas que conforman la unidad “C”, se pudieron identificar facies correspondientes a: depósitos de canales flu-viales y depósitos de llanura de inundación (Reyes et al., 1989).
Para la identificación de estas facies, se utilizó como guía los modelos y conceptos resultantes del estudio realizado por la División de Yacimientos, Geología de Desarrollo y Asociados y el Distrito de Producción El Centro de Ecopetrol en 1992.
Arenas de barras de meandro
Estas facies arenosas son las más características dentro de la secuencia estudiada en muestras de cora-zón, con espesores que oscilan entre 10 y 35 pies ca-racterizados por presentar hacia la base de las arenas de grano grueso a conglomerático, intraclastos arcillosos y fragmentos leñosos, descansando sobre superficies de erosión, con estratificación cruzada de alto ángulo en artesa que gradualmente cambia a estratificación plano paralela a subparalela discontinua y a ondulitas, acompañado por disminución gradual del tamaño de grano hasta llegar a grano fino pasando finalmente a facies arcillosas de planicie de inundación, (Bernal et al., 1992). Normalmente esta secuencia se ve interrum-pida debido al desarrollo de apilamientos verticales (“stacking point bars”), estableciendo comunicación
38
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
39CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
clasificaron bajo la denominación de Point Bars (PB); las arenas Crevasse Splay y Natural Levee se clasifi-caron bajo sus respectivos nombres y se encontrarán en adelante con la siguiente nomenclatura: CS y NL respectivamente. Adicionalmente, las arcillas que ac-túan como sellos tanto en los depósitos de canal (clay plugs) como en los depósitos de llanuras de inundación fueron clasificadas, quedando denominadas como CP y FP respectivamente.
Esta clasificación permitió que se generaran va-riables discretas, así: 1 para Point Bars (PB), 2 para Crevasse Splay (CS), 3 para Natural Levee (NL) y 4 para Clay Plugs (CP) y arcillas de Flood Plain (FP); las cuales fueron utilizadas para generar el modelo
Figura 10. Columna estratigráfica zona “C”
C1D
C1C
C1B2
C1C2
C2A
C1B1
C1A
C1A1
C2B
C2C
C2D
C2E
C2F
C2G
C2H
C2I
Clay1
Clay2
Clay3
Clay4
Clay5
Clay6
Clay7
Clay9
Clay15
Clay13
Clay10
Clay12
Clay14
Clay11
Clay8
CS
CS
PB
PB
PB
NL
NL
NL
PB
NL
PB
PB
PB
PB
PB
PB
FP
FP
FP
FP
FP
CP
FP
CP
CP
FP
CP
CP
CP
CP
FP
Litología Unidades Facies petrofísico que se utilizó para determinar la porosidad y permeabilidad en todos los pozos involucrados, es importante resaltar que las arcillas se unieron bajo una misma clasificación por: simplificación de los datos manejados y por que las arcillas, en este caso particular, no tienen importancia económica.
En el desarrollo de este trabajo se conservó la nomenclatura que se utilizó en el estudio “Modelos de Caracterización del Yacimiento: Geológico, de Distribución de Fluidos y Petrofísicos” de Bernal et al. (1992), para zonificar la zona “C” superior (C1 y C2); la subzona C3 no fue tenida en cuenta por su carácter acuífero.
El espesor promedio para el intervalo seleccionado es de 347 ft aproximadamente, en donde se identifi-caron 16 arenas de barras de arena, “natural levee” y
“crevasse splay”, las cuales fueron asociadas a 16 uni-dades genéticas, cada una de ellas representando ciclos de sedimentación coetáneos (Figura 10). Se destaca la relativa mayor continuidad de las fajas de meandros en dirección Norte – Sur, dirección en la cual se creen corrían las corrientes principales.
En términos generales se pueden dividir las dos subzonas presentes de la siguiente manera:
En la subzona C1 se identificaron 6 unidades genéti-cas llamadas de base a tope: C1A1, C1A, C1B1, C1B2, C1C y C1D. En la secuencia que separa las subzonas C1 y C2 se presenta una unidad no continua denominada C1C2 que registra condiciones de alta energía (tipo canal) dentro de una “llanura de inundación” (FP), la cual se comunica verticalmente con la unidad C2A de la subzona C2 (Bernal et al., 1992).
En C2 se han identificado 9 unidades que presen-tan un mayor acoplamiento vertical que el observado durante la depositación de las unidades del C1. Estas unidades se denominan de base a tope: C2I, C2H, C2G, C2F, C2E, C2D, C2C, C2B y C2A (Bernal et al., 1992).
Modelo estructuralLa profundidad de tope y base de cada una de las
unidades genéticas fue definida utilizando los regis-tros eléctricos y en particular el registro litológico disponible para cada pozo (SP o GR). Adicional a las
38
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
39CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
unidades descritas en el párrafo anterior, se relaciona-ron dos marcadores adicionales denominados TOPE, BASE y C3CB, las profundidades de estos, fueron determinadas utilizando la sísmica 3D y se tomaron del estudio “Evaluación Integrada Yacimientos del Campo La Cira-Infantas” (Lobo et al., 1999).
Los marcadores TOPE, C1D y C3CB fueron carto-grafiados y se utilizaron para definir la geomorfología del área en estudio. El marcador BASE fue utilizado solo para definir la base de la unidad genética C2I y no fue posible cartografiarlo debido a la escasez de datos, ya que fueron pocos los pozos que llegaron a esta profundidad.
Para definir la geomorfología en la parte superior de la zona “C” se seleccionó la unidad genética C1D (Figura 11(a)), este patrón estructural fue utilizado en la simulación geoestadística. Para la parte inferior exis-
Figura 11. (a) Mapa estructural marcador C1D
885-2283.50
886
-2268.50
908
-2265.00
1040
-2284.50
1041
-2301.50
1042
-2325.50
1043
-2329.00
1047
-2485.00
1049
-2524.00
1056
-2492.00
1057
-2543.50
1112
-2518.50
1113
-2605.001117
11181119
-2504.50
1121
-2547.50
1125
-2354.00
1126
-2354.00
1129
-2235.50
1130
-2242.00
1131
1134
-2242.00
1135
1136-2343.00
1137
1138
1140
1181
1182
1183
-2488.50
1184
-2575.50
1206
-2584.00
1207
-2538.00
1208
-2441.501209
1210
-2390.00
1211
1212
1213
1214
1218
1227
1228
1229
1233
1244
1245
-2155.00
1246
-2182.50
1247
1254
-2441.00
1256-2212.00
1257
1258
1483
-2300.00
1819
-2357.50
1820
-2286.00
1821
1822
1823
18241825
-2070.50
1884
-2415.00
1886
-2428.50
250 0 250 500 750 m 750 m250 0 250 500
885
-2737.00
886-2627.00
908
-2626.00
1040
-2646.00
1041
-2657.00
1042
-2652.00
1043
-2666.00
1047
-2798.00
1049
-2793.00
1056
-2705.00
1057
-2892.00
1112
-2871.00
1113
-2942.001117
-2882.00
1118
-2839.00
1119
-2840.00
1121
-2824.00
1125
-2675.00
1126
-2601.00
1129
-2606.00
1130
-2611.00
1131
-2641.00
1134
-2599.00
1135
-2600.00
1136-2644.00
1137
-2696.00
1138
-2729.00
1140
-2812.00
1181
-2687.00
1182
-2816.00
1183
-2831.00
1184
-2913.00-2918.00
1206
-2909.00
1207
-2858.00
1208
-2795.001209
-2780.001210
-2743.00
1211
-2780.00
1212
-2721.00
1213
-2662.00
1214
-2615.00
1218
-2536.00
1227
-2903.00 1228
-2893.00
1229
-2636.00
1233
-2552.00
1244
-2479.00
1245
-2512.00
1246
-2547.00
1247
-2566.00
1254
-2852.00
1256
-2529.00
1257
-2456.00
1258
-2326.00
1483
-2694.00
1819
-2711.00
1820
-2651.00
1821
-2582.00
1822-2530.00
1823
-2444.00
1824
-2389.00
1825
-2426.00
1884-2705.00
1886
-2748.50
tían dos posibles marcadores a utilizar: BASE y C3CB, siendo C3CB el más profundo; sin embargo, el mapa hecho con los puntos de control de BASE no era muy representativo de la geomorfología existente, por tal motivo se seleccionó el marcador C3CB (Figura 11(b)) para definir la geomorfología de la parte inferior de la zona “C” del área seleccionada.
SELECCIÓN DEL GRID SIMULACION
Con los modelos sedimentológico, estratigráfico y estructural definidos, se procedió a definir el modelo litológico más probable del área. Para tal fin, se gene-raron 51 realizaciones litológicas, las cuales fueron pos-teriormente evaluadas bajo el criterio de conectividad de los geo-objetos presentes en el yacimiento (Deutsch et al., 1998; Idrobo, 1999).
40
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
41CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
885-2283.50
886
-2268.50
908
-2265.00
1040
-2284.50
1041
-2301.50
1042
-2325.50
1043
-2329.00
1047
-2485.00
1049
-2524.00
1056
-2492.00
1057
-2543.50
1112
-2518.50
1113
-2605.001117
11181119
-2504.50
1121
-2547.50
1125
-2354.00
1126
-2354.00
1129
-2235.50
1130
-2242.00
1131
1134
-2242.00
1135
1136-2343.00
1137
1138
1140
1181
1182
1183
-2488.50
1184
-2575.50
1206
-2584.00
1207
-2538.00
1208
-2441.501209
1210
-2390.00
1211
1212
1213
1214
1218
1227
1228
1229
1233
1244
1245
-2155.00
1246
-2182.50
1247
1254
-2441.00
1256-2212.00
1257
1258
1483
-2300.00
1819
-2357.50
1820
-2286.00
1821
1822
1823
18241825
-2070.50
1884
-2415.00
1886
-2428.50
250 0 250 500 750 m 750 m250 0 250 500
885
-2737.00
886-2627.00
908
-2626.00
1040
-2646.00
1041
-2657.00
1042
-2652.00
1043
-2666.00
1047
-2798.00
1049
-2793.00
1056
-2705.00
1057
-2892.00
1112
-2871.00
1113
-2942.001117
-2882.00
1118
-2839.00
1119
-2840.00
1121
-2824.00
1125
-2675.00
1126
-2601.00
1129
-2606.00
1130
-2611.00
1131
-2641.00
1134
-2599.00
1135
-2600.00
1136-2644.00
1137
-2696.00
1138
-2729.00
1140
-2812.00
1181
-2687.00
1182
-2816.00
1183
-2831.00
1184
-2913.00-2918.00
1206
-2909.00
1207
-2858.00
1208
-2795.001209
-2780.001210
-2743.00
1211
-2780.00
1212
-2721.00
1213
-2662.00
1214
-2615.00
1218
-2536.00
1227
-2903.00 1228
-2893.00
1229
-2636.00
1233
-2552.00
1244
-2479.00
1245
-2512.00
1246
-2547.00
1247
-2566.00
1254
-2852.00
1256
-2529.00
1257
-2456.00
1258
-2326.00
1483
-2694.00
1819
-2711.00
1820
-2651.00
1821
-2582.00
1822-2530.00
1823
-2444.00
1824
-2389.00
1825
-2426.00
1884-2705.00
1886
-2748.50
Figura 11. (b) Mapa estructural marcador C3CB
Una de las aplicaciones de evaluar la conectividad de los geo-objetos es el ranqueo de múltiples reali-zaciones geoestadísticas. Los volúmenes de celdas conectadas pueden ser usados para clasificar múl-tiples realizaciones geoestadísticas, por ejemplo de alta, media y baja conectividad. El ranqueo solo puede ser desarrollado a partir de un número escalar, pero generalmente el carácter del modelo tridimensional no esta referenciado por escalares. Sin embargo, las medidas promedio de conectividad consideradas aquí, suministran una forma de ranquear realizaciones. Para este propósito se dispone del programa RANK_OBJ, desarrollado por Deutsch (1998).
Una vez seleccionado el modelo litológico más representativo para el área utilizada en la definición del modelo estratigráfico, se hicieron 20 descripcio-nes petrofísicas de dicho modelo. Posteriormente y por razones prácticas, se hizo una reducción del modelo a
caracterizar, para seleccionar esta área se escogieron tres patrones de inyección/producción, los cuales deli-mitaron un nuevo grid.
Modelo de simulaciónEl área seleccionada para llevar a cabo el ajuste
histórico automático, corresponde a tres modelos (patrones) de inyección/producción, su ubicación dentro del área de modelamiento estratigráfico se puede observar en la Figura 12 (Bradley, 1987). Estos tres patrones fueron seleccionados debido: a su buena respuesta al proceso de recobro secundario por inyección de agua, y a la no presencia de fallas, lo que simplificaba el modelo numérico (Buckley et al., 1942; Dykstra et al., 1950).
El área seleccionada tiene una superficie de aproxi-madamente 376 acres y esta compuesta por 12 pozos: LC-1131, LC-1137, LC-1181, LC-1208, LC-1209, LC-
40
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
41CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Figura 12. Área descripción petrofísica
12091208
1210
1212
11811137
1213
1214
1247
1229
1319
1131
1262850
1262100
1261350
1260600
1259850
1262850
1262100
1261350
1260600
1259850
1033280 1034080 1034880 1035680 1036480
Figura 13. Representación gráfica del yacimiento
12345678910
282726252423222120191817161514131211
1 2 3 4 5 6 7 8 9 10 191817161514131211
J
i=
Celdas Activas
xxxxxxxx
x
xxx
xx
x
xxxxxxxx
x
xxx
xx
x
xxxxxxxx
x
xxx
xx
x
xxxxxxxx
x
xxx
xx
x
xxxxxxxx
x
xxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xxxxxxxx
xx
x
xxxxxxxxxxxxxxxx
x
xx
xx
x
xxxxxxxxxxxxxxxx
x
xx
xx
x
xxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxx
xxxxxxxxxxxxxxxx
42
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
43CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
1210, LC-1212, LC-1213, LC-1214, LC-1229 LC-1247 y LC-1819, todos (excepto el pozo LC-1819) tuvieron producción primaria.
Las coordenadas de esta área son:
Xo = 1’035,040 m : Yo = 1’260,994 m
Xf = 1’036,085 m : Yf = 1’262,450 m
Una vez definido este nuevo grid, se optó por deshabilitar algunas de las celdas que componen este enmallado; estas, no se tuvieron en cuenta porque ha-cen parte del área de influencia de otros pozos que no fueron considerados dentro de este estudio. En la Figura 13 se pueden observar marcadas con X y sombreadas, las celdas que finalmente quedaron activas dentro del desarrollo del presente trabajo (Satter et al., 1994).
Para resumir los ajustes que se le han hecho al enmallado, se presentan los siguientes datos: grid de referencia, 60 x 60 x 1000; grid de geoestadística, 60 x 60 x 347; y grid de simulación escalado teniendo en cuenta tanto celdas activas como inactivas, 19 x 28 x 25 (total 13300 celdas), el tamaño de cada bloque es de 55 m x 52 m y con tamaño variable en z.
RESULTADOS
Una vez finaliza la ejecución del programa, este automáticamente genera un libro de Microsoft Excel compuesto por cinco hojas de cálculo donde se relacio-nan los resultados obtenidos después de haber hecho el ajuste histórico automático.
En la hoja llamada “Presentacion_1” se encuentran graficados los datos de producción de aceite (bopd) reales, simulados y promedio utilizado en streamline (Figura 14).
En la hoja llamada “Presentacion_2”, se presenta una tabla con los valores de permeabilidades relativas finalmente utilizados durante el proceso de simulación streamline, en la parte inferior de esta se encuentran los puntos finales de las permeabilidades relativas, así como también los valores de saturación de agua connata, la saturación de aceite residual y los factores de Corey (Corey, 1954). La Tabla 1 se presenta a continuación y allí se resumen estos parámetros. En la Figura 15 presenta de manera gráfica los valores relacionados en la Tabla 1.
Figura 14. Comparación curvas de producción
50,00
0,00
100,00
150,00
200,00
250,00
300,00
350,00
400,00
450,00
500,00
0 50 100 150 200 250 300
BOPD
Meses
Prod_Real
Prod_Sim
Prod_Sim_Ajustada
42
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
43CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Figura 15. Curvas de permeabilidades relativas
0,000
0,100
0,200
0,300
0,400
0,500
0,600
0,700
0,800
0,900
1,000
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1
Krw
Kro
Tabla 1. Permeabilidades relativas (modificadas), área 7, zona C – Campo La Cira
Sw Krw Kro So
0,3133 0,000000 0,925000 0,6868
0,3200 0,000069 0,811657 0,6800
0,3500 0,001386 0,442780 0,6500
0,4000 0,006354 0,145251 0,6000
0,4500 0,014240 0,040428 0,5500
0,5000 0,024743 0,009018 0,5000
0,5500 0,037682 0,001469 0,4500
0,6000 0,052926 0,000148 0,4000
0,6300 0,063138 0,000026 0,3700
0,6500 0,070378 0,000006 0,3500
0,6650 0,076032 0,000002 0,3350
0,6750 0,079906 0,000001 0,3250
0,6850 0,083864 0,000000 0,3150
0,6868 0,084586 0,000000 0,3132
0,7500 0,111598 0,000000 0,2500
0,7590 0,115708 0,000000 0,2410
0,7653 0,118623 0,000000 0,2347
endkro endkrw sor swc real_no real_nw
0,9250 0,1186 0,23475 0,31325 8,6875954 1,7730235
44
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
45CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Tabl
a 2.
Rel
ació
n ev
oluc
ión
indi
vidu
os
Indi
vidu
oSu
ma
Apt
itud
Ope
raci
on4
77
40
,03
SEG
79
5,6
7M
IN.
13
,26
HO
RAS
12
34
56
7
11.
675.
477,
751,
007
Repr
oduc
ir1.
543.
891,
381,
941
Com
bina
r1.
234.
673,
751,
552
Com
bina
r1.
494.
070,
881,
878
Com
bina
r1.
711.
925,
882,
370
Mut
ació
n1.
184.
561,
881,
640
Com
bina
r1.
645.
723,
882,
278
Mut
ació
n
22.
055.
567,
131,
236
Com
bina
r1.
118.
020,
751,
405
Com
bina
r1.
471.
882,
751,
850
Com
bina
r1.
056.
153,
631,
327
Com
bina
r1.
298.
298,
381,
797
Com
bina
r97
8.87
1,94
1,35
5C
ombi
nar
992.
734,
131,
374
Com
bina
r
31.
778.
239,
381,
069
Com
bina
r1.
034.
780,
811,
301
Com
bina
r1.
087.
800,
751,
367
Com
bina
r1.
252.
017,
131,
574
Com
bina
r1.
257.
096,
501,
740
Com
bina
r1.
188.
440,
251,
645
Com
bina
r1.
280.
410,
881,
772
Com
bina
r
41.
847.
969,
501,
111
Com
bina
r1.
404.
770,
501,
766
Com
bina
r93
6.72
4,94
1,17
7C
ombi
nar
1.06
8.34
3,50
1,34
3C
ombi
nar
1.02
4.90
1,88
1,41
9C
ombi
nar
1.27
5.10
1,88
1,76
5C
ombi
nar
1.26
1.66
9,00
1,74
6C
ombi
nar
51.
724.
279,
001,
037
Repr
oduc
ir95
8.07
5,75
1,20
4C
ombi
nar
1.69
9.53
7,25
2,13
6M
utac
ión
1.56
4.70
0,25
1,96
7C
ombi
nar
1.18
1.12
9,38
1,63
5C
ombi
nar
1.59
3.63
5,63
2,20
6M
utac
ión
1.39
6.97
0,25
1,93
4C
ombi
nar
61.
803.
142,
251,
084
Com
bina
r1.
094.
220,
631,
375
Com
bina
r1.
133.
733,
251,
425
Com
bina
r1.
071.
939,
131,
347
Com
bina
r1.
232.
981,
131,
707
Com
bina
r1.
160.
621,
881,
607
Com
bina
r1.
155.
333,
251,
599
Com
bina
r
72.
791.
061,
501,
678
Com
bina
r1.
543.
046,
381,
939
Com
bina
r1.
199.
021,
131,
507
Com
bina
r1.
518.
493,
131,
909
Com
bina
r1.
324.
940,
631,
834
Com
bina
r1.
620.
727,
132,
243
Mut
ació
n1.
294.
648,
381,
792
Com
bina
r
81.
673.
598,
751,
006
Repr
oduc
ir98
2.38
9,44
1,23
5C
ombi
nar
1.04
6.09
4,31
1,31
5C
ombi
nar
814.
629,
501,
024
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
91.
861.
659,
251,
119
Com
bina
r99
9.66
9,56
1,25
6C
ombi
nar
1.14
7.98
9,00
1,44
3C
ombi
nar
1.54
2.56
5,63
1,93
9C
ombi
nar
1.19
5.61
7,50
1,65
5C
ombi
nar
1.21
5.60
5,38
1,68
3C
ombi
nar
1.25
6.51
7,00
1,73
9C
ombi
nar
102.
086.
833,
131,
255
Com
bina
r1.
208.
426,
251,
519
Com
bina
r1.
315.
966,
001,
654
Com
bina
r1.
160.
329,
751,
458
Com
bina
r1.
320.
352,
631,
828
Com
bina
r1.
259.
173,
501,
743
Com
bina
r1.
170.
489,
751,
620
Com
bina
r
112.
300.
965,
751,
383
Com
bina
r1.
125.
493,
501,
415
Com
bina
r1.
340.
329,
251,
685
Com
bina
r1.
288.
546,
881,
620
Com
bina
r1.
280.
494,
131,
773
Com
bina
r1.
007.
133,
881,
394
Com
bina
r97
4.12
2,25
1,34
8C
ombi
nar
121.
802.
924,
251,
084
Com
bina
r1.
222.
624,
631,
537
Com
bina
r1.
369.
141,
131,
721
Com
bina
r95
1.55
3,50
1,19
6C
ombi
nar
867.
745,
311,
201
Com
bina
r1.
481.
005,
252,
050
Mut
ació
n99
8.20
0,81
1,38
2C
ombi
nar
131.
906.
713,
251,
146
Com
bina
r79
5.61
2,31
1,00
0Re
prod
ucir
795.
612,
311,
000
Repr
oduc
ir79
5.61
2,31
1,00
0Re
prod
ucir
795.
612,
311,
101
Repr
oduc
ir1.
034.
162,
631,
432
Com
bina
r1.
155.
815,
631,
600
Com
bina
r
142.
129.
620,
251,
280
Com
bina
r88
2.13
7,19
1,10
9Re
prod
ucir
1.59
9.24
4,63
2,01
0C
ombi
nar
1.46
9.14
0,75
1,84
7C
ombi
nar
1.73
6.79
6,25
2,40
4M
utac
ión
1.62
2.42
8,13
2,24
6M
utac
ión
1.32
1.48
0,75
1,82
9C
ombi
nar
152.
123.
554,
001,
277
Com
bina
r89
4.63
9,50
1,12
4Re
prod
ucir
1.00
6.98
3,25
1,26
6C
ombi
nar
1.83
5.54
3,75
2,30
7C
ombi
nar
1.00
9.59
9,06
1,39
8C
ombi
nar
1.04
3.54
0,94
1,44
5C
ombi
nar
1.26
1.79
2,13
1,74
7C
ombi
nar
161.
663.
364,
381,
000
Repr
oduc
ir1.
663.
364,
382,
091
Com
bina
r1.
295.
207,
131,
628
Com
bina
r1.
128.
353,
381,
418
Com
bina
r1.
129.
455,
251,
563
Com
bina
r1.
484.
827,
752,
055
Mut
ació
n1.
123.
822,
251,
556
Com
bina
r
171.
845.
265,
751,
109
Com
bina
r1.
080.
298,
751,
358
Com
bina
r1.
277.
129,
381,
605
Com
bina
r1.
322.
550,
001,
662
Com
bina
r1.
575.
207,
502,
180
Com
bina
r1.
387.
021,
131,
920
Com
bina
r1.
198.
519,
381,
659
Com
bina
r
181.
946.
947,
001,
170
Com
bina
r1.
100.
980,
881,
384
Com
bina
r97
8.68
2,63
1,23
0C
ombi
nar
1.02
0.70
2,81
1,28
3C
ombi
nar
1.18
8.98
9,38
1,64
6C
ombi
nar
1.05
7.76
9,75
1,46
4C
ombi
nar
1.47
1.89
3,00
2,03
7C
ombi
nar
192.
140.
972,
751,
287
Com
bina
r94
1.17
7,69
1,18
3C
ombi
nar
1.32
7.74
4,00
1,66
9C
ombi
nar
1.19
8.80
8,88
1,50
7C
ombi
nar
1.12
5.73
7,00
1,55
8C
ombi
nar
1.13
5.48
3,75
1,57
2C
ombi
nar
1.63
0.72
5,63
2,25
7M
utac
ión
201.
857.
158,
381,
117
Com
bina
r3.
077.
224,
002,
500
Mut
ació
n1.
307.
456,
631,
643
Com
bina
r1.
536.
861,
881,
932
Com
bina
r1.
203.
909,
501,
666
Com
bina
r1.
297.
800,
751,
796
Com
bina
r1.
489.
462,
752,
062
Mut
ació
n
1.9
50
.76
5,6
71
.23
3.5
42
,21
1.2
28
.54
7,6
71
.25
4.5
45
,83
1.2
09
.16
0,3
91
.23
7.5
16
,58
1.2
40
.13
7,4
7
89
10
11
12
13
14
1.34
4.44
7,13
1,86
1C
ombi
nar
1.26
7.37
5,50
1,75
4C
ombi
nar
1.28
7.16
0,75
1,78
2C
ombi
nar
1.17
9.85
5,50
1,63
3C
ombi
nar
1.33
1.84
6,50
1,84
4C
ombi
nar
1.11
6.64
2,63
1,54
6C
ombi
nar
1.07
8.85
6,88
1,49
3C
ombi
nar
903.
725,
631,
251
Com
bina
r1.
181.
630,
381,
636
Com
bina
r1.
275.
385,
381,
765
Com
bina
r1.
001.
502,
631,
386
Com
bina
r1.
282.
585,
751,
775
Com
bina
r1.
341.
818,
001,
857
Com
bina
r1.
381.
294,
131,
912
Com
bina
r
1.29
7.72
1,88
1,79
6C
ombi
nar
1.08
9.39
8,88
1,50
8C
ombi
nar
1.07
8.32
7,50
1,49
3C
ombi
nar
1.99
1.11
6,50
2,50
0M
utac
ión
1.51
4.99
9,50
2,09
7M
utac
ión
1.56
6.75
9,38
2,16
9M
utac
ión
1.31
1.98
3,75
1,81
6C
ombi
nar
1.07
5.11
8,38
1,48
8C
ombi
nar
1.23
8.69
1,13
1,71
5C
ombi
nar
1.10
8.78
6,25
1,53
5C
ombi
nar
1.13
9.51
0,75
1,57
7C
ombi
nar
864.
401,
311,
197
Com
bina
r87
6.31
2,63
1,21
3C
ombi
nar
1.10
4.81
2,63
1,52
9C
ombi
nar
1.00
9.37
1,38
1,39
7C
ombi
nar
1.17
1.26
3,50
1,62
1C
ombi
nar
1.15
4.27
2,13
1,59
8C
ombi
nar
1.05
9.70
6,88
1,46
7C
ombi
nar
1.02
6.10
4,69
1,42
0C
ombi
nar
1.07
5.95
6,50
1,48
9C
ombi
nar
1.17
9.77
2,38
1,63
3C
ombi
nar
1.17
2.74
7,50
1,62
3C
ombi
nar
1.19
0.61
4,63
1,64
8C
ombi
nar
909.
311,
061,
259
Com
bina
r1.
086.
141,
381,
503
Com
bina
r1.
232.
673,
001,
706
Com
bina
r1.
171.
310,
501,
621
Com
bina
r1.
144.
097,
751,
584
Com
bina
r
977.
876,
631,
354
Com
bina
r1.
578.
430,
002,
185
Mut
ació
n1.
058.
109,
501,
465
Com
bina
r1.
678.
234,
632,
323
Mut
ació
n1.
443.
092,
381,
998
Mut
ació
n1.
618.
727,
752,
241
Mut
ació
n1.
038.
474,
381,
437
Com
bina
r
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
000
Repr
oduc
ir
1.09
5.85
1,63
1,51
7C
ombi
nar
953.
385,
131,
320
Com
bina
r1.
086.
756,
251,
504
Com
bina
r1.
317.
166,
501,
823
Com
bina
r1.
480.
824,
502,
050
Mut
ació
n1.
162.
900,
751,
610
Com
bina
r1.
063.
299,
631,
472
Com
bina
r
912.
971,
751,
264
Com
bina
r1.
107.
992,
881,
534
Com
bina
r93
3.16
6,69
1,29
2C
ombi
nar
1.11
6.32
3,00
1,54
5C
ombi
nar
957.
929,
941,
326
Com
bina
r1.
125.
120,
501,
557
Com
bina
r1.
809.
386,
382,
500
Mut
ació
n
1.13
7.32
0,13
1,57
4C
ombi
nar
1.37
4.28
8,13
1,90
2C
ombi
nar
1.11
2.91
1,38
1,54
1C
ombi
nar
1.79
6.26
3,38
2,48
6M
utac
ión
983.
729,
751,
362
Com
bina
r1.
288.
880,
631,
784
Com
bina
r1.
807.
598,
002,
500
Mut
ació
n
1.34
3.51
9,00
1,86
0C
ombi
nar
1.21
3.40
0,13
1,68
0C
ombi
nar
1.06
5.20
6,00
1,47
5C
ombi
nar
1.38
1.23
8,88
1,91
2C
ombi
nar
1.20
3.40
4,88
1,66
6C
ombi
nar
1.08
0.22
0,13
1,49
5C
ombi
nar
1.07
1.77
6,25
1,48
4C
ombi
nar
995.
431,
811,
378
Com
bina
r97
5.00
5,69
1,35
0C
ombi
nar
1.04
9.75
7,63
1,45
3C
ombi
nar
1.01
9.01
5,38
1,41
1C
ombi
nar
945.
836,
001,
309
Com
bina
r76
6.74
9,13
1,06
1Re
prod
ucir
1.02
1.67
3,94
1,41
4C
ombi
nar
1.33
8.56
0,50
1,85
3C
ombi
nar
1.36
0.93
7,88
1,88
4C
ombi
nar
1.55
8.94
2,75
2,15
8C
ombi
nar
1.40
3.83
6,38
1,94
3C
ombi
nar
1.26
4.43
6,88
1,75
0C
ombi
nar
1.15
0.15
4,75
1,59
2C
ombi
nar
1.77
5.93
8,38
2,45
8M
utac
ión
1.56
0.79
0,75
2,16
1C
ombi
nar
1.38
5.42
4,75
1,91
8C
ombi
nar
1.13
8.11
6,25
1,57
5C
ombi
nar
1.19
9.26
9,38
1,66
0C
ombi
nar
1.19
0.20
7,00
1,64
8C
ombi
nar
1.30
6.46
0,88
1,80
8C
ombi
nar
1.31
2.50
7,00
1,81
7C
ombi
nar
1.65
5.28
2,50
2,29
1M
utac
ión
1.21
7.14
0,88
1,68
5C
ombi
nar
1.15
1.29
5,88
1,59
4C
ombi
nar
1.18
4.08
5,88
1,63
9C
ombi
nar
1.34
3.83
4,38
1,86
0C
ombi
nar
1.44
9.89
9,75
2,00
7M
utac
ión
1.53
2.09
8,13
2,12
1C
ombi
nar
1.12
0.61
8,38
1,55
1C
ombi
nar
1.28
3.39
3,13
1,77
7C
ombi
nar
1.31
3.19
5,00
1,81
8C
ombi
nar
1.44
3.16
0,00
1,99
8C
ombi
nar
1.40
2.07
4,13
1,94
1M
utac
ión
1.13
7.58
6,25
1,57
5C
ombi
nar
1.37
7.03
9,00
1,90
6C
ombi
nar
1.34
2.73
9,00
1,85
9C
ombi
nar
1.30
7.63
1,75
1,81
0C
ombi
nar
1.63
6.85
7,00
2,26
6M
utac
ión
1.03
4.24
5,63
1,43
2C
ombi
nar
1.49
9.22
4,38
2,07
5M
utac
ión
1.24
0.66
9,50
1,71
7C
ombi
nar
1.03
9.23
1,63
1,43
9C
ombi
nar
1.29
5.32
1,75
1,79
3C
ombi
nar
1.15
9.59
9,38
1,60
5C
ombi
nar
1.13
6.67
6,63
1,57
3C
ombi
nar
896.
814,
501,
241
Com
bina
r1.
262.
580,
131,
748
Com
bina
r1.
488.
489,
632,
060
Mut
ació
n1.
595.
551,
132,
209
Com
bina
r
1.21
0.49
6,75
1,67
6C
ombi
nar
1.14
3.76
5,25
1,58
3C
ombi
nar
1.33
2.43
1,25
1,84
4C
ombi
nar
889.
711,
191,
232
Com
bina
r1.
261.
274,
501,
746
Com
bina
r91
7.30
8,69
1,27
0C
ombi
nar
898.
435,
001,
244
Com
bina
r
1.1
75
.61
6,5
31
.19
6.0
89
,36
1.1
55
.45
4,1
81
.22
6.9
80
,83
1.2
10
.67
3,8
91
.18
0.2
19
,31
1.2
63
.31
2,2
3
44
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
45CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
15
16
17
18
19
20
1.28
9.14
5,50
1,78
4C
ombi
nar
1.35
4.49
2,63
1,87
5C
ombi
nar
1.51
9.49
4,00
2,10
3M
utac
ión
1.15
2.75
8,88
1,59
6C
ombi
nar
1.18
3.73
2,88
1,67
9C
ombi
nar
1.43
3.47
7,13
2,23
8M
utac
ión
1.13
3.98
9,50
1,57
0C
ombi
nar
923.
577,
501,
278
Com
bina
r1.
457.
270,
502,
017
Mut
ació
n98
4.33
8,81
1,36
3C
ombi
nar
883.
819,
381,
254
Com
bina
r93
1.56
7,13
1,45
4C
ombi
nar
1.44
8.45
8,25
2,00
5C
ombi
nar
1.03
8.60
2,50
1,43
8C
ombi
nar
995.
567,
691,
378
Com
bina
r1.
195.
375,
631,
655
Com
bina
r1.
255.
113,
381,
780
Com
bina
r80
7.54
6,69
1,26
1C
ombi
nar
1.08
9.48
8,75
1,50
8C
ombi
nar
1.26
9.47
3,13
1,75
7C
ombi
nar
1.11
9.00
4,13
1,54
9C
ombi
nar
793.
532,
441,
098
Repr
oduc
ir1.
194.
334,
251,
694
Com
bina
r64
0.49
4,44
1,00
0Re
prod
ucir
1.04
4.96
9,06
1,44
6C
ombi
nar
1.19
4.81
3,63
1,65
4C
ombi
nar
1.11
5.86
0,63
1,54
5C
ombi
nar
1.26
4.91
4,75
1,75
1C
ombi
nar
1.47
8.14
9,50
2,09
7C
ombi
nar
1.33
3.88
7,88
2,08
3C
ombi
nar
1.26
3.05
0,25
1,74
8C
ombi
nar
1.15
8.08
6,25
1,60
3C
ombi
nar
1.15
4.40
1,25
1,59
8C
ombi
nar
1.02
1.92
5,88
1,41
5C
ombi
nar
1.17
1.77
0,38
1,66
2C
ombi
nar
863.
245,
631,
348
Com
bina
r
883.
337,
941,
223
Com
bina
r1.
250.
016,
251,
730
Com
bina
r1.
207.
412,
501,
671
Com
bina
r1.
268.
075,
751,
755
Com
bina
r1.
455.
980,
882,
065
Com
bina
r1.
781.
972,
882,
500
Mut
ació
n
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
000
Repr
oduc
ir72
2.41
8,25
1,00
0Re
prod
ucir
722.
418,
251,
025
Repr
oduc
ir97
2.05
2,00
1,51
8C
ombi
nar
1.65
5.61
1,63
2,29
2C
ombi
nar
1.02
5.49
9,06
1,42
0C
ombi
nar
1.34
3.52
3,63
1,86
0C
ombi
nar
1.16
1.65
3,38
1,60
8C
ombi
nar
1.15
5.63
3,00
1,63
9C
ombi
nar
1.22
4.39
9,50
1,91
2C
ombi
nar
1.37
6.51
4,88
1,90
5C
ombi
nar
1.13
9.97
9,25
1,57
8C
ombi
nar
1.30
7.11
7,38
1,80
9C
ombi
nar
1.10
0.03
9,75
1,52
3C
ombi
nar
956.
215,
251,
356
Com
bina
r97
4.50
4,81
1,52
1C
ombi
nar
919.
103,
501,
272
Com
bina
r1.
016.
783,
631,
407
Com
bina
r91
1.85
4,69
1,26
2C
ombi
nar
1.20
3.05
4,00
1,66
5C
ombi
nar
1.26
5.14
5,13
1,79
5C
ombi
nar
1.27
7.55
6,25
1,99
5C
ombi
nar
1.06
7.82
6,88
1,47
8C
ombi
nar
931.
737,
251,
290
Com
bina
r1.
248.
559,
631,
728
Com
bina
r1.
195.
881,
131,
655
Com
bina
r1.
273.
349,
631,
806
Com
bina
r93
6.78
7,25
1,46
3C
ombi
nar
959.
701,
381,
328
Com
bina
r90
8.36
7,31
1,25
7C
ombi
nar
1.08
2.91
4,88
1,49
9C
ombi
nar
5.30
5.51
3,00
2,50
0M
utac
ión
704.
976,
131,
000
Repr
oduc
ir70
4.97
6,13
1,10
1Re
prod
ucir
1.01
3.02
1,00
1,40
2C
ombi
nar
1.11
9.63
6,75
1,55
0C
ombi
nar
1.36
3.87
4,88
1,88
8C
ombi
nar
1.38
0.92
7,13
1,91
2C
ombi
nar
1.82
4.02
7,63
2,50
0M
utac
ión
1.46
7.14
0,63
2,29
1M
utac
ión
1.33
4.14
6,75
1,84
7C
ombi
nar
1.21
5.24
8,88
1,68
2C
ombi
nar
1.00
4.10
0,69
1,39
0C
ombi
nar
837.
066,
691,
159
Com
bina
r1.
054.
101,
631,
495
Com
bina
r1.
159.
699,
131,
811
Com
bina
r
1.93
9.66
5,00
2,50
0M
utac
ión
1.62
8.34
8,88
2,25
4M
utac
ión
1.58
7.99
6,38
2,19
8M
utac
ión
1.86
6.00
2,25
2,50
0M
utac
ión
2.15
4.41
7,75
2,50
0M
utac
ión
1.96
8.53
0,25
2,50
0M
utac
ión
1.04
4.76
5,25
1,44
6C
ombi
nar
1.21
8.43
1,00
1,68
7C
ombi
nar
1.17
5.98
0,63
1,62
8C
ombi
nar
1.19
4.72
4,38
1,65
4C
ombi
nar
1.18
9.89
5,88
1,68
8C
ombi
nar
1.05
0.41
3,38
1,64
0C
ombi
nar
995.
000,
501,
377
Com
bina
r1.
433.
092,
001,
984
Com
bina
r1.
160.
708,
001,
607
Com
bina
r1.
773.
716,
252,
455
Com
bina
r94
5.07
8,25
1,34
1C
ombi
nar
1.40
7.96
1,25
2,19
8C
ombi
nar
1.29
0.65
8,38
1,78
7C
ombi
nar
1.39
3.45
7,25
1,92
9C
ombi
nar
1.41
3.85
5,63
1,95
7C
ombi
nar
1.53
4.27
7,00
2,12
4C
ombi
nar
994.
434,
811,
411
Com
bina
r1.
093.
792,
881,
708
Com
bina
r
956.
052,
381,
323
Com
bina
r1.
137.
710,
501,
575
Com
bina
r1.
045.
612,
381,
447
Com
bina
r1.
480.
037,
382,
049
Com
bina
r1.
469.
415,
252,
084
Com
bina
r1.
125.
027,
131,
756
Com
bina
r
1.1
71
.34
6,2
51
.15
3.9
88
,59
1.1
96
.87
6,3
81
.42
1.8
11
,63
1.2
16
.60
0,4
61
.15
7.7
51
,62
21
22
23
24
25
1.32
1.20
0,50
2,06
3C
ombi
nar
1.24
0.01
6,00
1,93
6C
ombi
nar
1.17
5.37
4,25
1,83
5C
ombi
nar
1.77
8.85
4,75
2,50
0M
utac
ión
1.17
4.66
8,50
1,83
4C
ombi
nar
1.00
2.14
7,88
1,56
5C
ombi
nar
1.15
9.74
3,25
1,81
1C
ombi
nar
838.
597,
811,
309
Com
bina
r1.
354.
209,
632,
114
Com
bina
r1.
109.
202,
751,
732
Com
bina
r
878.
259,
881,
371
Com
bina
r85
1.97
8,44
1,33
0C
ombi
nar
1.15
2.44
4,25
1,79
9C
ombi
nar
1.39
6.01
4,88
2,18
0C
ombi
nar
923.
343,
881,
442
Com
bina
r
640.
494,
441,
000
Repr
oduc
ir64
0.49
4,44
1,00
0Re
prod
ucir
640.
494,
441,
000
Repr
oduc
ir64
0.49
4,44
1,00
0Re
prod
ucir
640.
494,
441,
000
Repr
oduc
ir
1.30
1.93
9,25
2,03
3C
ombi
nar
1.01
2.18
3,94
1,58
0C
ombi
nar
1.33
4.95
0,13
2,08
4C
ombi
nar
1.00
1.03
0,38
1,56
3C
ombi
nar
1.17
5.30
7,38
1,83
5C
ombi
nar
1.01
0.48
9,19
1,57
8C
ombi
nar
1.40
7.94
6,38
2,19
8C
ombi
nar
918.
007,
561,
433
Com
bina
r1.
472.
030,
632,
298
Com
bina
r1.
053.
953,
751,
646
Com
bina
r
900.
970,
311,
407
Com
bina
r1.
200.
845,
001,
875
Com
bina
r1.
583.
524,
132,
472
Mut
ació
n1.
551.
892,
502,
423
Mut
ació
n1.
317.
197,
002,
057
Com
bina
r
956.
748,
061,
494
Com
bina
r84
9.31
6,06
1,32
6C
ombi
nar
767.
582,
191,
198
Com
bina
r88
2.39
5,63
1,37
8C
ombi
nar
829.
141,
561,
295
Com
bina
r
1.12
4.26
2,25
1,75
5C
ombi
nar
900.
855,
191,
406
Com
bina
r92
4.25
6,19
1,44
3C
ombi
nar
1.10
1.03
3,38
1,71
9C
ombi
nar
940.
381,
501,
468
Com
bina
r
1.13
9.95
0,50
1,78
0C
ombi
nar
951.
735,
881,
486
Com
bina
r1.
170.
319,
001,
827
Com
bina
r1.
344.
508,
252,
099
Com
bina
r1.
104.
683,
131,
725
Com
bina
r
1.17
9.87
3,13
1,84
2C
ombi
nar
867.
697,
811,
355
Com
bina
r1.
512.
994,
502,
362
Com
bina
r1.
481.
038,
882,
312
Mut
ació
n1.
056.
914,
001,
650
Com
bina
r
917.
757,
131,
433
Com
bina
r1.
224.
047,
251,
911
Com
bina
r1.
074.
577,
631,
678
Com
bina
r1.
352.
198,
632,
111
Com
bina
r1.
200.
376,
501,
874
Com
bina
r
751.
310,
191,
173
Com
bina
r1.
103.
769,
381,
723
Com
bina
r85
1.99
8,50
1,33
0C
ombi
nar
1.20
7.95
0,13
1,88
6C
ombi
nar
966.
606,
631,
509
Com
bina
r
1.81
2.86
8,50
2,50
0M
utac
ión
1.50
2.57
2,88
2,34
6C
ombi
nar
1.30
8.14
3,50
2,04
2C
ombi
nar
1.77
3.48
7,63
2,50
0M
utac
ión
1.34
3.74
1,88
2,09
8C
ombi
nar
1.02
2.97
1,25
1,59
7C
ombi
nar
998.
740,
691,
559
Com
bina
r92
2.48
7,06
1,44
0C
ombi
nar
1.25
1.52
8,88
1,95
4C
ombi
nar
1.16
0.15
5,13
1,81
1C
ombi
nar
2.09
8.47
0,25
2,50
0M
utac
ión
2.05
5.39
6,63
2,50
0M
utac
ión
1.94
3.04
6,25
2,50
0M
utac
ión
1.85
3.22
8,63
2,50
0M
utac
ión
1.87
5.94
5,25
2,50
0M
utac
ión
1.59
1.35
0,50
2,48
5C
ombi
nar
1.28
2.94
4,75
2,00
3C
ombi
nar
1.06
7.62
1,00
1,66
7C
ombi
nar
1.00
5.46
2,69
1,57
0C
ombi
nar
1.42
7.54
3,88
2,22
9C
ombi
nar
1.14
3.23
4,38
1,78
5C
ombi
nar
1.37
9.42
7,63
2,15
4C
ombi
nar
1.62
7.81
5,63
2,50
0M
utac
ión
1.26
1.61
9,00
1,97
0C
ombi
nar
1.16
8.19
7,50
1,82
4C
ombi
nar
808.
949,
251,
263
Com
bina
r1.
400.
149,
252,
186
Com
bina
r1.
005.
395,
001,
570
Com
bina
r1.
086.
590,
001,
696
Com
bina
r1.
371.
484,
382,
141
Com
bina
r
1.82
4.03
2,13
2,50
0M
utac
ión
1.03
2.47
6,75
1,61
2C
ombi
nar
1.13
3.47
7,50
1,77
0C
ombi
nar
1.20
1.16
8,38
1,87
5C
ombi
nar
839.
891,
631,
311
Com
bina
r
1.1
71
.36
3,9
51
.15
3.1
16
,88
1.1
47
.65
5,3
31
.29
9.8
36
,86
1.1
33
.96
1,5
3
Porc
enta
jes
de s
elec
ción
Repr
oduc
ción
:
10%
Com
bina
ción
:
70%
Mut
ació
n :
20%
46
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
47CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
En la hoja llamada “Hoja1” se relaciona la siguiente información: tiempo de máquina, evolución de los in-dividuos durante las 20 generaciones, individuo rey en cada generación, los porcentajes de selección definidos por el usuario y finalmente los puntos piloto seleccio-nados por el programa. En la Tabla 2, se observan todos los datos anteriores a excepción de los puntos piloto, estos no se incluyeron debido a la magnitud de datos generados y a que estos no son un factor concluyente del método y se usan para revisar el comportamiento del programa.
En la hoja llamada “Cromosomas” se presentan las tablas de equivalencias generadas para los cromosomas binarios, estas tablas son utilizadas por el programa para codificar y descodificar los valores de las variables
que los utilizan (puntos finales de las permeabilidades relativas y los parámetros geoestadísticos).
Para validar el modelo obtenido se calcularon los percentiles 10, 20, 30, …, 80, 90 y 100 para todos los datos que componen cada propiedad petrofísica (porosidad, permeabilidad y saturación de agua) y se compararon con los percentiles 10, 20, 30, …, 80, 90 y 100 del individuo 2 de la población inicial, este in-dividuo fue seleccionado aleatoriamente entre los 20 individuos (realizaciones) generadas como población inicial. Los resultados obtenidos se pueden observar en la Tabla 3.
Adicional al cálculo de los percentiles, se hicieron gráficas donde se puede evidenciar el comportamiento
Figura 16. Representación variación permeabilidad
1 1 15
15
5 5
10
10
20
20
25
25
Permeabilidad_X0,000100 75 150 225 300
46
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
47CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Percentiles Porosidad
Individuo Inicial Individuo Final
10% 0,000 10% 0,000
20% 0,000 20% 0,000
30% 0,000 30% 0,000
40% 0,000 40% 0,000
50% 0,056 50% 0,010
60% 0,097 60% 0,086
70% 0,148 70% 0,124
80% 0,196 80% 0,174
90% 0,227 90% 0,210
100% 0,316 100% 0,320
Percentiles Permeabilidad
Individuo Inicial Individuo Final
10% 0,00 10% 0,00
20% 0,00 20% 0,00
30% 0,00 30% 0,00
40% 0,00 40% 0,00
50% 3,42 50% 1,38
60% 8,20 60% 6,70
70% 24,66 70% 15,77
80% 71,06 80% 48,87
90% 147,25 90% 102,82
100% 579,36 100% 994,79
Percentiles Sat. Agua
Individuo Inicial Individuo Final
10% 0,311 10% 0,315
20% 0,356 20% 0,351
30% 0,420 30% 0,409
40% 0,495 40% 0,479
50% 0,563 50% 0,543
60% 0,647 60% 0,607
70 0,933 70% 0,900
80% 1,000 80% 1,000
90% 1,000 90% 1,000
100% 1,000 100% 1,000
Tabla 3. Percentiles propiedades petrofísicas (modificadas)
de la propiedad, en el individuo inicial y en el individuo final. En este trabajo se muestra el comportamiento de la permeabilidad en seis capas del modelo de simula-ción (1, 5, 10, 15, 20 y 25). La Figura 16, muestra las capas que se relacionaron anteriormente, teniendo en cuenta que la capa 1 es la base (más profunda) y la capa 25 es el tope (superior) de la zona “C” después de ser escalada.
Se evidencia que las propiedades petrofísicas han sido modificadas dentro de ciertos rangos permitidos, logrando de esta forma mejorar el ajuste histórico (Fi-gura 14) y disminuyendo el valor de la función objetivo de 1’663,364 en la primera generación a 640,494 en la vigésima generación (Tabla 2).
CONCLUSIONES
• Se ha logrado aplicar una metodología algorítmica para obtener un modelo de distribución en 3D de las propiedades petrofísicas condicionado a los datos estáticos y dinámicos disponibles para una parte del área 7 del campo La Cira (yacimiento fluvial), con éxito.
• En este trabajo, se muestra que la computación evo-lutiva puede ser aplicada con éxito como herramienta para la caracterización de yacimientos de hidrocarbu-
48
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
49CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
ros condicionando el modelo geológico y estructural a los datos de producción.
• Tal como se evidencia en la Tabla 2 (Relación evolución individuos), la función objetivo se va re-duciendo, al igual que los valores promedio de las funciones objetivo en cada generación. Esta mejoría se puede ver en las Figuras 14 y 16.
• Para aplicar la metodología se elaboró un programa llamado INVERDINAMICA ©, desarrollado en pla-taforma PC y bajo ambiente Windows; el programa fue compilado, y presentado como un ejecutable.
• Esta metodología ayudará, a los ingenieros de yaci-mientos, a realizar la calibración (ajuste histórico) del modelo propuesto, disminuyendo las horas hombre requeridas para realizar este proceso.
• La combinación de los AGs (método de optimi-zación global), el método de los puntos piloto, la Simulación Gaussiana Secuencial (modelamiento geoestadístico) y la simulación streamline (simu-lador numérico) proporcionan una herramienta ro-busta para ser utilizada en la solución de problemas inversos aplicados a la industria del petróleo. Se genera de esta forma una solución muy económica en términos de costo computacional, eficiente y totalmente automática.
• La técnica definida en este trabajo puede ser aplica-da tanto a yacimiento en producción primaria como en yacimientos sometidos a la inyección de agua. Para yacimientos en recobro primario se haría ne-cesario el uso de un simulador convencional, como consecuencia, los tiempos de cómputo aumentarán considerablemente.
• Utilizar un modelo sin escalar, es posible, ya que el simulador de streamline permite su uso y, segura-mente, con el mejoramiento de las velocidades de cómputo en los equipos venideros, será posible obviar los errores que se generan al escalar los modelos.
• La metodología mostró ser una técnica apropiada para encontrar descripciones de yacimientos que hacen honor a los datos de producción, información de los pozos y a los rasgos geológicos; además es capaz de manejar gran cantidad de parámetros, lo cual es un problema cuando se trabaja en modelos de simulación de yacimientos.
• La metodología fue aplicada exitosamente en la parte superior de las arenas “C” de 3 patrones inyección/producción (12 pozos) del área 7 del Campo La Cira, este es un yacimiento heterogéneo de ambiente de depositación fluvial.
• El análisis de la conectividad de las arenas es un méto-do efectivo para clasificar las realizaciones litológicas y de esta forma se puede hacer una selección rápida de modelos ubicados en algún percentil en particular. En esta tesis el percentil utilizado fue el P-50.
• Para este estudio en particular, se encontró que la incertidumbre generada en la determinación de la Saturación de Agua, después de 30 años de produc-ción primaria, es mucho más significativa que la proveniente del cálculo del volumen de arcilla. Por esta razón, se cambió el volumen de arcilla por la saturación de agua.
• El concepto utilizado para diseñar la metodología de los AGs es tan universal que permite ser modificada fácilmente, y se adecua a las condiciones del proble-ma propuesto.
• Al comparar los resultados obtenidos por el simu-lador con los datos históricos, se observa que aún existen diferencias significativas en la producción al momento de la irrupción del agua, este problema se debe más a condiciones operativas que al método propuesto; se puede observar que el pozo La Cira 1210 es el que mejor ajusta tanto en la tasa de pro-ducción como en el corte de agua, esto se explica fácilmente al considerar su ubicación estratigráfica (parte baja de la estructura), lo que implica que la cantidad de gas libre al inicio de la inyección era
48
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
49CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
RECOMENDACIONES
• No se evaluaron los beneficios de aplicar un método de optimización local en el individuo rey al terminar la ejecución del Algoritmo Genético, por falta de tiempo. Se recomienda evaluar la implementación de un método de optimización local al finalizar la ejecución del AG.
• Aumentar las propiedades petrofísicas en conside-ración, incluyendo el volumen de arcilla como variable; incorporar parámetros adicionales tales como la presión capilar, propiedades PVT, características del acuífero, Skin, etc.
• La versión actual del software INVERDINAMI-CA© es la 2003.2, se recomienda continuar trabajando en el mejoramiento de este programa implementando más características que lo hagan más flexible. Durante el desarrollo de este trabajo de investigación se gene-raron varias aplicaciones programadas en VBA (Visual Basic para Aplicaciones) en Microsoft Excel, las cuales se deben continuar desarrollado y finalmente, deben ser integradas a INVERDINAMICA©, automatizando el proceso completo.
AGRADECIMIENTOS
Los autores quieren expresar su agradecimiento a Ecopetrol S.A. por la financiación de la totalidad de este proyecto y al Instituto Colombiano del Petróleo (ICP) por el apoyo logístico brindado. A los Ingenieros Jorge Mantilla V., Edgar Martínez H. y Edwin Pineda C. por su ayuda incondicional.
BIBLIOGRAFÍA
Bäck, T., 1996. “Evolutionary algorithms in theory and practice”. Oxford University Press, New York.
Baker, J. E., 1985. “Adaptive selection methods for genetic algorithms”. Grefenstette, J. J. ed., Proceedings of the First International Conference on Genetic Algorithms and their Applications, Lawrence Erlbaum Associates, Hillsdale, New Jersey.
Batycky, R. P., 1997. “A three-dimensional two-phase field scale streamline simulator”. PhD Thesis, Stanford Uni-versity, Dept. of Petroleum Engineering, Stanford, CA.
Batycky, R. P., Blunt, M. J., and Thiele, M. R., 1997. “A 3D Field Scale Streamline-Based Reservoir Simulator”, SPE Reservoir Engineering, 246-254.
Bernal, M. C., 1992. “Modelos de caracterización del ya-cimiento: geológico, de distribución de fluidos y petro-físicos”. Ecopetrol.
Bissell, R. C., Dubrule, O., Lamy, P., Swaby, P. and Lepine, O., 1987. “Combining geostatistical modelling with gradient information for history matching: the pilot point method”. Paper SPE 38730, presented at the 1997 SPE Annual Technical Conference and Exhibition, San Antonio, Texas, 5-8.
Bradley, H., 1987. “Petroleum engineering handbook”. So-ciety of Petroleum Engineers, Richardson, Texas.
Buckley, S. E. and Leverett, M. C. 1942. “Mechanisms of fluid displacement in sands”. Trans., AIME 146, 107-116.
Bush, M. D. and Carter, J. N., 1996. “Application of a mo-dified genetic algorithm to parameter estimation in the petroleum industry”. Dagli et al. eds., Intelligent Engi-neering Systems through Artificial Neural Networks 6, ASME Press, New York, 397-402.
Caers, J., Krishnan, S., Wang, Y. D., and Kovscek, A. R., 2001. “A geostatistical approach to streamline-based history matching”. In proceedings of the Stanford Center for Reservoir Forecasting, Stanford, CA.
Calvete, F., 2001. “Computación evolutiva: algoritmos ge-néticos y simulated annealing, técnicas meta-heurísticas de optimización en ingeniería de petróleos”. IX Congreso Colombiano del Petróleo, Asociación Colombiana de In-genieros de Petróleos (ACIPET), Bogotá.
menor que la de los otros dos pozos; el pozo con menor ajuste es La Cira 1214 donde se observa que el corte de agua del simulador aumenta rápidamen-te comparado con el corte de agua histórico. Puede decirse que este es el pozo estructuralmente más alto, y por lo tanto, el que tiene mayor efecto por presencia de gas. Esta condición no es compatible con el simulador streamline, ya que este no maneja fluidos compresibles.
50
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
51CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Carrillo, L. F., 1990. “Análisis petrofísicos especiales de corazones, pozo La Cira 1884”. Ecopetrol – ICP, Piedecuesta.
Chen, W. H., Gavalas, G. R., Seinfeld, J. H., And Wasserman, M. L., 1974. “A new algorithm for automatic history matching”. SPEJ, 593-608.
Corey, A. T., 1954. “The interrelation between gas and oil relative permeabilities”. Producers Monthly, 38-41.
Cuypers, M., Dubrule, O., Lamy, P. and Bissell, R., 1998. “Op-timal choice of inversion parameters for history-matching with the pilot point method”. Proceeding of the Sixth Euro-pean Conference on the Mathematics of Oil Recovery.
Dake, L. P., 1978. “Fundamentals Of Reservoir Enginee-ring”. Developments in Petroleum Science 8, Elsevier, Amsterdam.
Davis, L., 1991. “Handbook of genetic algorithms”. Van Nostrand Reinhold, New York.
De Jong, K. A., 1975. “An analysis of the behavior of a class of genetic adaptive systems”. Ph.D. dissertation, University of Michigan, Ann Arbor.
Dengo, C. A. and Covey, M. C., 1993. “Structure of the Eastern Cordillera of Colombia: Implications for Trap Styles and Regional Tectonics”. AAPG Bulletin, 77/8, 1315-1337.
Deschamps, T., Grussaute, T., Mayers, D. and Bissell, R., 1998. “The results of testing six different gradient opti-misers on two history matching problems”. Proceedings of the Sixth European Conference on the Mathematics of Oil Recovery.
Deutsch, C. V. and Journel, A. G. 1998. “GSLIB: Geostatis-tical software library and user’s guide”. second edition, Oxford University Press, New York.
Deutsch, C. V. 1998. “Fortran programs for calculating connectivity of three-dimensional numerical models and for ranking multiple realizations”. Computers & Geosciences, 241, 69.
Dykstra, H. and Parsons, H. L. 1950. “The prediction of oil recovery by waterflooding”. Secondary Recovery of Oil in the United States, 2nd ed., API, New York, 160-174.
Goldberg, D. E., 1989. “Genetic Algorithms in Search, Op-timization, and Machine Learning”. Addison-Wesley Publishing Company, Inc., Reading, Massachusetts.
Goldberg, D. E. and Deb, K., 1991. “A comparative analy-sis of selection schemes used in genetic algorithms”. Rawlins, G. J. E. (ed.), Foundations of Genetic Algori-thms, Morgan Kaufmann.
Gomez, S., Gosselin, O., and Barker, J. W., 1999. “Gradient-ba-sed history-matching with a global optimization method”. Paper SPE 56756 presented at the 1999 SPE Annual Tech-nical Conference and Exhibition, Houston, Texas, 3-6.
Holland, J. H., 1975. “Adaptation in natural and artificial systems: an introductory analysis with applications to biology, control, and artificial intelligence”. University of Michigan Press, Ann Arbor, Michigan.
Idrobo, E. A. 1999. “Characterization and ranking of reservo-ir models using geostatistics and streamline simulation”. PhD Dissertation, Texas A&M University.
Jensen, J., Lake, L., Corbett, P., and Goggin, D., 1997. “Sta-tistics for petroleum engineers and geoscientists”. Pren-tice Hall Petroleum Engineering Series, New Jersey.
Lobo, A., 1999. “Evaluación integrada de yacimientos, campo La Cira-Infantas”. Ecopetrol, Bogotá
Martínez, J., 1990. “Fundamentos de algoritmos genéticos”. II Congreso Colombiano de Neurocomputación, Univer-sidad Nacional de Colombia, Bogotá.
Mitchell, M., Forrest, S., and Holland, J. H., 1992. “The royal road for genetic algorithms: fitness landscapes and ga performance”. Varela, F. J. and Bourgine, P., (eds.), Towards a Practice of Autonomous Systems, Procee-dings of the First European Conference on Artificial Life, MIT Press.
Mitchell, M., 1996. “An introduction to genetic algorithms”. The MIT Press, Cambridge, Massachusetts.
Ouenes, A., 1992. “Application of simulated annealing to reservoir characterization and petrophysics inverse pro-blems”. PhD Thesis, New Mexico Institute of Mining and Technology.
Ramon, J.c., 1997. “Characterización and prediction of reser-voir architecture and petrophysical properties in fluvial channel sandstones, Middle Magdalena Basin, Colom-bia”. CT&F - Ciencia, Tecnología y Futuro, Piedecuesta, Santander: Ecopetrol, ICP , 1: 3.
Renard, Ph. and De Marsily, G., 1997. “Calculating equi-valent permeability: a review”. Advances in Water Re-sources, 20: 253-278.
50
SANTIAGO GONZÁLEZ et al.
CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
CARACTERIZACIÓN DINÁMICA DE YACIMIENTOS ESTRATIGRAFICAMENTE COMPLEJOS USANDO ALGORITMOS GENÉTICOS
51CT&F - Ciencia, Tecnología y Futuro - Vol. 2 Núm. 5 Dic. 2004
Romero, C.E., 2000. “A genetic algorithm for reservoir cha-racterisation using production data”. PhD Dissertation, University of London, Royal School of Mines.
Romero, C. E., Carter, J. N., Gringarten, A. C. and Zimmer-man, R. W., 2000. “A modified genetic algorithm for reservoir characterization”. SPE 64765, SPE Internatio-nal Oil and Gas Conference and Exhibition, Society of Petroleum Engineers, Beijing.
Romero, C. E., Carter, J. N., Gringarten, A. C. and Zimmer-man, R. W., 2000. “Improved reservoir characterization through evolutionary computation”. SPE 62942, SPE Annual Technical Conference and Exhibition, Society of Petroleum Engineers, Dallas.
Reyes, A., Fajardo, A. and Ballesteros, R., 1989. “Des-cripción macropetrográfica del pozo La Cira 1884”. Ecopetrol-ICP, Piedecuesta.
Satter, A. and Thakur, G. C., 1994. “Integrated petroleum reservoir management: a team approach”. Penn Well Publishing Company, Tulsa, Oklahoma.
Sen, M. K., Datta-gupta, A., Stoffa, P. L., Lake, L. W. and Pope, G. A., 1995. “Stochastic reservoir modeling using simulated annealing and genetic algorithms”. SPE For-mation Evaluation, 49-55.
Stoffa, P. L. and Sen, M. K., 1991. “Nonlinear multiparameter optimization using genetic algorithms: inversion of plane-wave seismograms”. Geophysics, 56, 1794-1810.
Surry, P. D. and Radcliffe, N. J., 1997. “Real representa-tions”. Belew, R.K. and Vose, M.D., eds.. Foundations of Genetic Algorithms 4, Morgan Kaufmann Publishers, San Francisco, California.
Thomas, L. K., Hellums, L. J. and Reheis, G. M., 1971. “A nonlinear automatic history matching technique for re-servoir simulation models”. Paper SPE 3475 presented at the SPE, 46th Annual Fall Meeting, New Orleans, Louisiana, 3-6.
Wang, Y. and Kovscek, A. R.,2000. “Streamline approach for history matching production data”. SPE Journal, (4): 353-362.
Whitley, D., 1999. “A free lunch proof for gray versus binary encodings”. Proceedings of the Genetic and Evolutiona-ry Computation Conference, Banzhat, W., eds., Morgan Kaufmann, San Francisco.
Wright, A. H., 1991. “Genetic algorithms for real parame-ter optimisation”. Rawlins, G. J. E. ed., Foundations of Genetic Algorithms, Morgan Kaufmann.
Wu, Z., Reynolds, A. C. and Oliver, D. S., 1998. “Condi-tioning geostatistical models to two-phase production data”. SPEJ, 4 (2).
Yang, P-h. and Watson, A. T., 1987. “Automatic history ma-tching with variable-metric methods”. Paper SPE 16977 presented at the SPE 62nd Annual Technical Conference and Exhibition, Dallas, Texas, 27-30.
Zalzala, A. M. and Fleming, P. J., 1997. “Genetic algorithms in engineering systems”. The Institution of Electrical En-gineers, London, U.K.