Orietta Mata
Orietta Mata
Orietta Mata
INGENIERIA GEOFISICA
IMPLEMENTACIN DE ALGORITMOS PARA ESTUDIOS DE FACTIBILIDAD DE ATRIBUTOS SSMICOS UTILIZANDO DATOS SINTTICOS.
Por Br. Orietta Carolina Mata Pacheco. Realizado con la asesora del Profesor: Jos Regueiro.
Proyecto de Grado Presentado ante la Ilustre Universidad Simn Bolvar Como requisito Parcial para optar al Ttulo de Ingeniero Geofsico
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Este trabajo ha sido aprobado en nombre de la Universidad Simn Bolvar por el siguiente jurado calificador:
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 1: INTRODUCCIN 1.1 Antecedentes 1.2 Objetivo General 1.3 Software empleado
1 2 3 3
CAPTULO 2: MARCO TERICO 2.1 Aspectos Petrofsicos 2.1.1 Registros de pozo convencionales 2.1.1.1 Registro de gamma ray (GR) 2.1.1.2 Registro de densidad 2.1.1.3 Registro snico 2.1.1.4 Registro de resistividad 2.1.1.5 Registro de potencial espontneo (SP) 2.1.1.6 Registro calper 6 6 7 7 7 8 8 9
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
2.1.2 Registro de pozo no convencionales 2.2 Aspectos Ssmicos 2.2.1 Impedancias acsticas, coeficientes de reflexin y serie de reflectividad. 2.2.2 Sismograma Sinttico 2.2.3 Ondculas Fase Cero 2.2.4 Atributos Ssmicos 2.2.5 Clasificacin de los atributos ssmicos 2.2.5.1 Atributos pre apilamiento 2.2.5.2 Atributos post apilamiento 2.2.5.3 Anlisis de Traza Compleja 2.2.5.4 Atributos de Traza Compleja 2.2.5.5 Atributos Estadsticos de Amplitud 2.3 Aspectos Estadsticos 2.3.1 Estadstica. Estadstica descriptiva e inductiva 2.3.2 Variables. Tipos de variables. 2.3.3 Escalas de medidas 2.3.4 Medidas de tendencia central. Media, mediana y moda. 2.3.5 Rango percentil. Cuntiles, cuartiles, deciles y percentiles 2.3.6 Medidas de dispersin. Desviacin estndar y varianza
9 9 10
12 14 16 18 19 19 25 29 34 39 39 41 41 42 43 44
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 3: MARCO METODOLGICO 3.1 Rasgos Generales 3.2 Carga de los datos 3.3 Obtencin del Sismograma 3.3.1 Clculo de la escala en tiempo 3.3.2 Clculo de las impedancias acsticas, del coeficiente de reflexin y del sismograma sinttico 3.4 Aplicacin de los atributos ssmicos 3.4.1 Clculo de la traza compleja y de los atributos instantneos 3.4.2 Clculo de los atributos por ventana 3.5 Edicin de registros de pozo en tiempo 3.6 Herramientas de la interfase ASPASS / Atributos Instantneos 3.6.1 Herramienta Restar 3.7 Clculos estadsticos. 3.8 Herramientas de la interfase ASPASS / Estadsticas 52 52 52 53 55 55 56 56 46 47 49 49 51
CAPTULO 4: CARGA DE LOS DATOS DE POZO Y CLCULO DE SISMOGRAMAS SINTTICOS 4.1 Carga de los datos 4.2 Visualizacin de los registros y clculo del sismograma sinttico 4.3 Edicin de los registros en tiempo 4.4 Reclculo de sismogramas sintticos a partir de los registros 62 70 76 78
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
editados en tiempo
CAPTULO 5: CLCULO DE ATRIBUTOS 5.1 Clculo de Atributos Ssmicos Instantneos 5.2 Clculo de Atributos Ssmicos por Ventana. Atributos de Amplitud. CAPTULO 6: ESTADSTICAS 6.1 Inicio del mdulo ASPASS / Estadsticas 6.2 Comparacin de resultados 6.3 Herramientas del mdulo Estadsitcas 91 93 97 80 85
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
NDICE DE FIGURAS.
Registros de pozo convencionales. Proceso sintetizado de generacin del sismograma sinttico. Serie de ondculas pertenecientes a las diferentes familias. Clasificacin de los atributos ssmicos segn Alistair Brown. Anlisis de traza compleja La funcin amplitud instantnea calculada a partir de una traza ssmica compleja.
10 15 16 20 27 31
Figura 7.
32
Figura 8.
34
Figura 9
(A) Media, moda y mediana para distribuciones simtricas. (B) Media, mediana y moda para distribuciones no simtricas.
43
Figura 10.
Flujogramas del programa ASPASS. (A) Flujograma del mdulo ASPASS / Principal. (B) Flujograma del mdulo ASPASS / Atributos Instantneos. (C) Flujograma del mdulo ASPASS / Atributos de Amplitud. (D) Flujograma del mdulo ASPASS / Estadsticas
61
63 64
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 13. Figura 14. Figura 15. Figura 16. Figura 17. Figura 18. Figura 19. Figura 20. Figura 21. Figura 22. Figura 23. Figura 24. Figura 25. Figura 26.
Paso inicial para cargar un archivo. Seleccin del archivo a cargar. Proceso de seleccin de registros Interfase que permite organizar los registros e introducir los parmetros iniciales Mensaje de error que se despliega al introducir un registro en el campo errneo Parmetros iniciales de la ondcula. Diferentes maneras de activar el cuadro de Curvas. Cuadro de Curvas Pantalla principal mostrando los registros graficados Diferentes maneras de calcular el sismograma sinttico Sismograma calculado por las rutinas de ASPASS a partir de los registros de densidad y snico. Opciones relacionadas a la apariencia y forma del sismograma sinttico. Opcin Modificar Registros (A) Registros densidad y snico antes de la edicin de los valores en tiempo. (B) Registro densidad y snico despus de la edicin en tiempo. Cuadro de dilogo que permite escoger en que pista se encuentra cada registro. Sismograma sinttico calculado por ASPASS. Como inicializar el mdulo Atributos Instantneos. Pantalla principal del mdulo Atributos instantneos. Opcin Actualizar Sinttico
65 65 66 67 68 69 71 71 72 74 74 75 76 77
Figura 27. Figura 28. Figura 29. Figura 30. Figura 31.
78 79 80 81 82
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Opcin Calcular atributos Atributos instantneos calculados a partir del sismograma que se muestra en la pista 1. Barra de mens, opcin Herramientas de la pantalla de Atributos Instantneos. Salida de la herramienta Restar Sintticos, la cual permite restar dos sintticos e inclusive sus atributos para visualizar mejor las variaciones existentes Como inicializar el mdulo Atributos por Ventana, Atributos de Amplitud. Pantalla principal del mdulo Atributos por Ventana Opciones Actualizar Sinttico y Sinttico Original Opcin Definir Ventana Salida de los atributos calculados en la interfase de ASPASS / Atributos de Amplitud Lista de atributos de amplitud por ventana disponibles para cada pista. Basta con seleccionar un nombre de la lista para que l atributo sea calculado y mostrado por pantalla. Valores para cada atributo. Valores definidos del eje X para la pista 3. Diferentes maneras de iniciar el mdulo ASPASS / Estadsticas. Interfaz grfica del mdulo estadstico de ASPASS. Resultado de comparar el sismograma sinttico y los atributos seleccionados Lista de atributos disponibles para realizar la comparacin Campo de los valores de los atributos por ventana. Rojo denota original y azul modificado
82 83 83 83
Figura 36. Figura 37. Figura 38. Figura 39. Figura 40. Figura 41.
85 86 87 88 89 89
Figura 42. Figura 43. Figura 44. Figura 45. Figura 46. Figura 47. Figura 48.
90 90 92 93 94 95 96
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Herramienta de la interfaz ASPASS / Estadsticas que permite redimensionar el eje X de las pistas. Opcin Herramientas del mdulo ASPASS / Estadsticas. Lista de atributos a comparar dentro de la herramienta Variacin de los Atributos Variacin en los valores de los Atributos. Grfico de barras que muestra la variacin porcentual entre los valores obtenidos para original y modificado. Tabla representativa de los valores del grfico de barras.
96 97 98 99
Figura 53.
100
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
IMPLEMENTACIN DE ALGORITMOS PARA ESTUDIOS DE FACTIBILIDAD ALGORITMOS ESTUDIO DE ATRIBUTOS SSMICOS UTILIZANDO DATOS SINTTICOS.
Por Br. Orietta Carolina Mata P
RESUMEN. ESUMEN.
ASPASS (Atributos Ssmicos Post Apilamiento sobre Sismogramas Sintticos), fue diseado para realizar un estudio de factibilidad del uso de atributos ssmicos sobre sismogramas sintticos en la caracterizacin de yacimientos. Este software permite el clculo de sismogramas sintticos a partir de registros de pozo y la aplicacin de atributos ssmicos al mismo. Adems permite la edicin de registros en tiempo y el clculo del nuevo sismograma a partir de los registros modificados para despus comparar los cambios significativos y reflejar los resultados porcentualmente. Para demostrar la validez de las rutinas desarrolladas, se utilizaron un conjunto de datos de pozo, entre los cuales se encontraban registros snico y densidad. Los resultados obtenidos sirvieron como indicadores de la eficiencia y utilidad de las programas diseados, los cuales parecieron ser favorables.
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
DESTINY IS NOT A MATTER OF CHANCE, IT IS A MATTER OF CHOICE; IT IS NOT A THING TO BE WAITED FOR, IT IS A THING TO BE ACHIEVED..
WILLIAM JENNINGS BRYAN
Orietta C. Mata P.
ii
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
EL DESTINO NO ES CUESTIN DE OPORTUNIDADES SINO DE ELECCIONES, NO ES ALGO PARA SER ESPERADO, SINO PARA SER LOGRADO..
WILLIAM JENNINGS BRYAN
Orietta C. Mata P.
iii
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Dedicatoria.
Orietta C. Mata P.
iv
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Agradecimientos.
A lo largo de nuestras vidas, el camino que recorremos hasta llegar aqu es arduo y difcil. Sin embargo, como cada vez que alcanzamos la meta, las recompensas no se hacen esperar y mas pronto que tarde vemos el fruto del enorme sacrificio que empezamos unos cuantos aos atrs. Primero que nada quiero agradecerle a Dios, quien estuvo acompandome todo el trayecto y una que otra vez me llev en brazos hasta donde pens que no poda llegar. Por esas noches de iluminacin que me brindaste y por esos sabios consejos que al odo me susurraste. Por escucharme cuando me hizo falta y por secarme las lgrimas cuando las emociones se desbordabanGRACIAS! A mi familia, no caben en mi boca palabras para agradecerle toda su compaa y apoyo. A mi mami por orme, aconsejarme y regaarme tantas veces, para que no perdiera la perspectiva. Por todos los ratos desagradables que compartiste conmigo y por todas las lgrimas que derramaste en tu preocupacin de madre. GRACIAS! Sin eso no sera la mujer que soy ahora. Te amo. A mis hermanos, que de la mano me acompaaron por los caminos ms difciles que me tocaron vivirMuchachos, sin ustedes no hubiese podido hacerlo. A mis tos y mis abuelos, su ayuda y sus consejos no tienen precio. A mis primos, por su amistad sincera e invaluable. A todos de corazn, les agradezco el hecho de haber compartido conmigo tan grata experiencia. A mis amigos, los cuales por suerte no los puedo contar con los dedos de las manos. De ustedes aprend muchas cosas, incluso de las lgrimas que derram. A los que estuvieron para apoyarme cuando ms lo necesit y siguen estando. A los que compartieron conmigo momentos inolvidables y por razones de la vida ya no me acompaan. A todos aquellos que no se cansaron de mis llamadas de telfono, de mis malos humores y de mis historias que parecen de novela. Gracias a Mairim y a Ornella, son las mejores amigas que puedo desear. Gracias a Mayra, Gaby, Patricia,
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Grecia, Jenny y Ma. Antonieta , por compartir mis experiencias y haber salido airosas como para contarlasjajajaja!. Gracias a Oscar, Ernesto, Cho, Pibe, Ramn y Jacinto. Muchachos, gracias por brindarme una sonrisa cuando pens que no haba motivos para sonrer. A todos los adoro! Gracias a mis compaeros tesistas, los cuales no me dejaron desfallecer en esas noches de Messenger para no quedarme dormida, o simplemente por no dejar que me diera por vencida. Juancho y Cucha, especialmente gracias! Andrs Luke Landa, qu te puedo decir?... Te estar agradecida de por vida. Sin tus sabios consejos (y tu milagroso Petroluke, jajajaja!), no lo hubiese podido hacer. Literalmente, te has ganado mi respeto, mi cario y mi admiracin. Gracias! Gracias a mi tutor, el Prof. Jos Pepe Regueiro, ya que sin sus sabios consejos y su mucha pero mucha paciencia no habra podido terminar; y gracias al Prof. Mojtaba Taheri por sus horas de dedicacin en mis das difciles de geoestadstica. Y finalmente, pero no menos importante, a mi pap. Te fuiste cuando menos lo pensaba y me dejaste el corazn en pedacitos. Nos faltaron horas, das, meses y aos para compartir sonrisas, lgrimas y experiencias. Pero Dios sabe lo que hace y solo El entiende los motivos de sus acciones. Me enseaste un montn de cosas, las cuales no pude agradecer en su debido momento. Sin embargo, se que desde donde quieras que ests, estars viendo a la mujer en que me convert. Ojal y ests orgulloso. Se que no has dejado de acompaarme, pues en mis momentos de desesperacin escuch tus sabias palabras y tu carioso consuelo. Y en todos y cada uno de mis xitos te he sentido siempre junto a m. No pude decrtelo en su momento, pero quiero que sepas que te amo. Gracias por ayudar a mi mam a hacerme la mujer que soy. Te estar agradecida eternamente y ojal cuando volvamos a encontrarnos me recibas con un Bien hecho hija!. Espero que el da de mi graduacin ests junto a m, en primera fila, aplaudiendo como todos los dems, porque este logro es en tu honor. GRACIAS!
Orietta C. Mata P.
vi
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 1.
INTRODUCCIN.
En la caracterizacin de yacimientos nunca son suficientes las herramientas que se puedan utilizar, con tal de obtener una imagen lo ms parecida posible del subsuelo. Es as como encontramos mtodos que van desde la aplicacin de simples algoritmos para calcular atributos ssmicos hasta rutinas complicadas que involucran inversin ssmica.
Desde su aparicin en los aos 70, los atributos ssmicos han jugado un papel importante en el desarrollo de la caracterizacin, definicin y delimitacin de yacimientos. Se han encontrado numerosas aplicaciones a la determinacin de litologas, fluidos, porosidad, estratigrafa, efectos de entonacin etc. Sin embargo, es necesario el clculo de un conjunto de ellos para obtener la mayor informacin posible del subsuelo, lo que requiere tiempo y dedicacin. La aplicacin de estos atributos sobre sismogramas sintticos puede ayudar en gran medida a definir los parmetros esenciales para realizar un estudio sobre datos reales, delimitando la serie de estos que van a ser usados en el estudio. Este trabajo de grado presenta el desarrollo y aplicacin de un mtodo, que si quizs no es innovador, es al menos de gran utilidad para definir los cambios producidos en los atributos ssmicos a partir de variaciones en el sismograma sinttico. Estas variaciones a su vez, reflejan los cambios en los parmetros petrofsicos en presencia de algn rasgo especfico, como cambios de litologa por ejemplo.
Orietta C. Mata P.
CAPTULO 1: INTRODUCCIN
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
ASPASS (Atributos Ssmicos Post Apilamiento sobre Sismogramas Sintticos) permite el clculo de sismogramas sintticos a partir de registros de pozo y la posterior aplicacin de atributos ssmicos sobre los mismos. Adems realiza la edicin en tiempo de los datos de pozo para luego recalcular el sismograma y los atributos, reflejando porcentualmente las diferencias obtenidas despus de la modificacin. En el captulo 2, se explica detalladamente todas y cada una de las bases tericas que ataen a este estudio. En el captulo 3 se exponen todos los pasos que se siguieron para la realizacin del trabajo: desde la carga de los datos, su organizacin, su presentacin y su manipulacin hasta el posterior clculo de los sismogramas sintticos y de los atributos. Adicionalmente, la metodologa utilizada para el clculo estadstico permitir reflejar los resultados en forma porcentual. En los captulos 4, 5 y 6 se presentan el manual del usuario, donde se explica paso a paso como se usa el programa diseado, ASPASS, adems de mostrar los resultados obtenidos.
1.1.-Antecedentes. .1.-Antecedentes.
A pesar de que en la actualidad existen diversos paquetes que ofrecen dentro de sus funciones el clculo de atributos ssmicos a datos reales, no existe ninguno que especficamente estudie el resultados porcentualmente. Sin embargo, dentro de la Universidad Simn Bolvar se han desarrollado una serie de programas con el fin de automatizar y facilitar la enseanza de ciertos conceptos bsicos durante el desarrollo de la carrera de Ingeniera Geofsica. uso de atributos sobre datos sintticos y que refleje los
Orietta C. Mata P.
CAPTULO 1: INTRODUCCIN
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Dentro del rea de ssmica se encuentran: Implementacin de algoritmos para el clculo y aplicacin de atributos AVO instantneos en sismogramas sintticos, realizado por Carlos Cova en el ao 2003. Implementacin de algoritmos para el clculo de atributos de amplitud instantnea AVO en datos reales, realizado por Ramn Asuaje en el ao 2004. Y mas reciente, en el rea de petrofsica: Diseo de un algoritmo para la visualizacin y procesamiento automatizado de registros de pozo, por Andrs Landa en el ao 2004.
Orietta C. Mata P.
CAPTULO 1: INTRODUCCIN
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Es un lenguaje de alto desempeo diseado para realizar clculos tcnicos e integra el clculo, la visualizacin y la programacin en un ambiente fcil de utilizar donde los problemas y las soluciones se expresan en una notacin matemtica. Mat Lab es un sistema interactivo cuyo elemento bsico de datos es el arreglo que no requiere de un dimensionamiento previo. Esto permite resolver muchos problemas computacionales, especficamente aquellos que involucren vectores y matrices. (Atencia, 2001) Mat Lab se utiliza ampliamente en: Clculos numricos Desarrollo de algoritmos Modelado, simulacin y prueba de prototipos Anlisis de datos, exploracin y visualizacin Graficacin de datos con fines cientficos o de ingeniera Desarrollo de aplicaciones que requieran de una interfaz grfica de usuario (GUI, Graphical User Interface) Entre sus ms importantes caractersticas se encuentran: Uso de nmeros reales y complejos, adems de toda clase de
funciones matemticas. En relacin a esto provee funciones especficas (en forma de subrutinas) para realizar clculos complicados con matrices. Realizar grficos en 2D y 3D, con la posibilidad de ser modificados en
tiempo de ejecucin.
Orietta C. Mata P.
CAPTULO 1: INTRODUCCIN
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
clculos desde una pantalla de comandos, estos clculos pueden ser programados desde el editor de Mat Lab.
Orietta C. Mata P.
CAPTULO 1: INTRODUCCIN
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 2.
MARCO TERICO.
Debido a que el desarrollo del trabajo involucra varias disciplinas relacionadas todas a las geociencias, el marco terico se divide en varias categoras: aspectos petrofsicos, aspectos ssmicos y aspectos estadsticos.
2.1.- Aspectos petrofsicos. .1.Dentro de esta categora se hace referencia a todas aquellas definiciones y conceptos relacionados a registros de pozo.
2.1.1.- Registro convencionales de pozos. .1.Los registros convencionales de pozo son aquellos que provienen de mediciones directas hechas a lo largo de una formacin atravesada por un pozo, empleando para ello herramientas mecnicas cuyo principio de funcionamiento y calibracin depende de la propiedad fsica que se desea medir. Registro de pozo se define como un registro de una o mas medidas fsicas como funcin de la profundidad en un pozo. (Sheriff, 2001). Dentro de esta clase de registros se encuentran:
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
densidad. 2.1.1.2.- Registro de densidad. .1.2.El registro densidad refleja las mediciones de densidad de una formacin al medir la atenuacin de rayos gamma entre una fuente artificial y un detector. Los rayos gamma emitidos continuamente son canalizados dentro de la formacin, sufriendo mltiples colisiones con electrones, las cuales los hacen perder energa y formar una nube. El tamao de la nube depende principalmente de las propiedades difusivas de la formacin, es decir, de su densidad electrnica. Esta nube se encoge y se expande a medida que la densidad vara, encogindose cuando la densidad aumenta y viceversa. Las unidades son gramos/centmetros cbicos (gr/cm).
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
registro como tiempo de trnsito en la formacin, en microsegundos/pie. (s/pie). Los tiempos de viaje compresionales varan desde 40s/pie en formaciones duras hasta 150s/pie en rocas blandas. Las velocidades correspondientes, las cuales son el inverso de los tiempos de viaje, varan entre 25000 a 6600 pies/s.
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
principalmente para seleccionar zonas permeables, para estimar valores de Rw y, menos frecuente, para estimar arcillosidades.
diseada para medir el dimetro del hoyo construido durante la perforacin de un pozo, y el registro producido, dimetro vs. Profundidad, nos da una clara idea acerca del estado en el que se encuentra el hoyo al momento del perfilaje.
2.1.2.- Registro no convencionales de pozos. .1.2.Los registros no convencionales de pozo provienen de clculos realizados sobre los registros convencionales mediante frmulas o algoritmos matemticos que definen un parmetro elstico o una propiedad fsica y su variacin a lo largo del pozo. La serie de reflectividad a lo largo del pozo y el perfil de impedancias acsticas entran dentro de esta clasificacin.
2.2.- Aspectos Ssmicos. .2.En esta categora se exponen todos los conceptos y definiciones alusivas a ssmica de reflexin y a los atributos que se pueden obtener de los datos.
Orietta C. Mata P.
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 1: Registros de pozo convencionales. De izquierda a derecha: porosidad neutrn, gamma ray, densidad, velocidad de ondas de cizalla y calper. Tomado de Petroluke, 2004
(www.glossary.oilfield.slb.com/)
Z = .V
(1)
Orietta C. Mata P.
10
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
La impedancia acstica es importante ya que permite determinar los coeficientes de transmisin y reflexin en la interfaz de dos materiales que poseen dos impedancias acsticas diferentes. Debido a que no existe una continuidad en los valores de las impedancias de dos medios adyacentes se producen reflexiones de la onda. La intensidad de la fraccin de la onda incidente que se refleja puede ser obtenida a partir de la velocidad de la partcula y las presiones locales que son requeridas para tener un trayecto continuo a lo largo de la interfaz entre los materiales. El coeficiente de reflexin se define como la proporcin entre la amplitud de la onda reflejada y la onda incidente, es decir, cuanta energa es reflejada. Para ondas con incidencia normal, el coeficiente de reflexin se calcula de la siguiente manera:
R=
(2)
Donde R = coeficiente de reflexin, cuyos valores oscilan entre -1 y 1. 1 = densidad del medio 1 V1 = velocidad del medio 1. 2 = densidad del medio 2 V2 = velocidad del medio 2. Z1 = impedancia acstica del medio 1 Z2 = impedancia acstica del medio 2. Los valores tpicos de R son aproximadamente -1 si pasamos de agua a aire, significando que casi el 100% de la energa es reflejada y nada es transmitido; 0.5 de
Orietta C. Mata P.
11
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
agua a roca; y 0.2 de lutita a arena. Para incidencia no normal, es decir, con un ngulo de incidencia diferente de cero, el coeficiente de reflexin definido como la proporcin de las amplitudes depende de otros parmetros, como las velocidades de las ondas de cizalla, y son descritos como funcin del ngulo de incidencia por las ecuaciones de Zoeppritz. Una formacin sobre otra de mayor impedancia acstica, produce un coeficiente de reflexin positivo, mientras que una formacin sobre otra de menor impedancia acstica genera un coeficiente de reflexin negativo. Por lo tanto, una serie de capas pueden ser representadas por una serie de coeficientes de reflexin en sus bordes. Una onda ssmica generada en superficie se reflejar en las interfases con una amplitud proporcional al coeficiente de reflexin. (Requena, 2000). Los coeficientes de reflexin representan la ubicacin y magnitud de los cambios litolgicos del subsuelo por lo que constituyen el modelo geolgico del mismo La serie de reflectividad, no es ms que asociar los coeficientes de reflexin a un valor de tiempo, por lo tanto, al final obtendremos un perfil que asemeja bastante bien la distribucin de los estratos en el subsuelo.
Orietta C. Mata P.
12
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Otra definicin utilizada por los intrpretes ssmicos es que el sinttico, como comnmente se denomina, es un modelo unidimensional de energa acstica que se propaga a travs de los estratos de la Tierra. Se pueden comparar estratos definidos o realizar otra correlacin de algn punto escogido en un registro de pozo con las reflexiones importantes de una seccin ssmica. La calidad de la correspondencia en el sismograma sinttico depender de cun bueno es el registro de pozo, el procesamiento de los datos y la habilidad para extraer una ondcula representativa de la datos ssmicos, entre otros factores. (www.glossary.oilfield.slb.com/) Un sismograma sinttico completo, presentara los eventos primarios, primarios mas mltiples, mltiples solamente y primarios con prdidas de transmisin. Como el propsito es lograr la mayor similitud entre datos de perfiles y la traza ssmica de superficie, es necesario un buen control de los parmetros de la ondcula que lo define.(Requena, 2000). El sismograma sinttico puede ser graficado con ambas polaridades con la escala de tiempo que corresponda a la seccin de ssmica de superficie. Simultneamente con el sismograma sinttico se pueden graficar otros perfiles tales como snico ajustado, densidad de formacin, impedancia acstica y coeficientes de reflexin. Para efectos de correlacin e interpretacin se recomienda incluir perfiles tales como resistividad, gamma ray, calper, porosidad, conductividad o cualquier curva proveniente de una evaluacin petrofsica.
Orietta C. Mata P.
13
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
2.2.3.- Ondculas fase cero. .2.3.Una ondcula es un pulso unidimensional, usualmente como respuesta bsica de un reflector simple. Esta se origina como un paquete de energa desde la fuente, teniendo un origen especfico en tiempo, y retorna a los receptores como una serie de eventos distribuidos en tiempo y energa. La distribucin es una funcin de los cambios de velocidad y densidad en el subsuelo y la relativa posicin de estos con respecto a la fuente y el receptor. La energa que regresa no excede a la inicial, debido a que la energa en una ondcula decae a medida que se particiona a travs de las interfaces (Principio de particin de la energa). Las ondculas solo decaen debido a un proceso de prdida de energa como calor durante el proceso de propagacin. Esto se hace ms notable a altas frecuencias, ya que la ondcula tiende a mantener menores altas frecuencias en relacin a las bajas frecuencias. (http://www.glossary.oilfield.slb.com/). La ondcula elegida para el anlisis debera ser en principio, lo ms similar posible a las ondas generadas por la fuente ssmica usada.(Requena, 2000) Una ondcula fase cero es aquella construida simtricamente alrededor de cero y tiene un pico de mxima amplitud en este punto. (z Yilmaz, Seismic Data Processing). Ahora bien, existen diferentes ondculas fase cero, las cuales pertenecen a diferentes familias. Es as como encontramos: 1.- Ondcula Meyer. Esta familia de ondculas son simtricas, infinitamente derivables y poseen transformada discreta. 2.- Ondcula Gaussiana. Las ondculas pertenecientes a este grupo tambin son infinitamente derivables, pero slo las derivadas pares son simtricas. A esta familia
Orietta C. Mata P.
14
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 2: Proceso sintetizado de generacin del sismograma sinttico. Luego de obtener la escala vertical en tiempo a partir del registro de profundidad, con los registros de velocidad y densidad se genera el perfil de impedancias acsticas y la serie de reflectividad. Con la convolucin de una ondcula, en este caso fase cero, y los coeficientes de reflexin finalmente se obtiene la traza sinttica. Tomado de www.glossary.oilfield.slb.com/
pertenece la conocida ondcula Ricker, la cual es ampliamente utilizada, ya que sus propiedades espectrales son conocidas. La expresin matemtica correspondiente a la Ricker es:
w(tw) = 1 2( f p .t ) e
Donde:
( f p .t ) 2
(3)
fp = frecuencia pico para el espectro de amplitud de la ondcula. t = tiempo total ocupado por los n tiempos a igual intervalo de muestreo en los que se cumplen las muestras de la seal.
Orietta C. Mata P.
15
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
3.- Ondcula Mexican Hat. Provienen de la segunda derivada de la funcin de densidad de probabilidad de la Gaussiana. Son simtricas y su expresin matemtica es de la forma:
mexh ( x ) =
2 3
1/ 4
* e ( x
/ 2)
* (1 x 2 )
(4)
4.- Ondcula Morlet. Esta clase de ondculas son simtricas, y su frmula matemtica se expresa como:
morl ( x) = e ( x
/ 2)
* cos(5 x)
(5)
Figura 3: Serie de ondculas pertenecientes a las diferentes familias. (A) Ondcula Ricker, con 40 Hz de frecuencia, 40 muestras por ondcula y 1 ms de muestreo. (B) Ondcula Meyer. (C) Ondcula Morlet. (D) Ondcula Mexican Hat. Tomado de ASPASS, 2005.
Orietta C. Mata P.
16
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Los atributos ssmicos son medidas especficas de geometra, cinemtica o caractersticas estadsticas derivadas de la traza ssmica. (Valerio, 2001). Todos los parmetros manejados ssmicamente pueden ser denominados como atributos ssmicos. Ellos pueden ser velocidad, amplitud, frecuencia y la tasa de cambio de cualquiera de ellos con respecto al tiempo y al espacio. Los principales objetivos de estos atributos son proveer informacin detallada y certera al intrprete de los parmetros estructurales, estratigrficos y litolgicos del prospecto ssmico. (Turhan Taner, 2001). Estos atributos pueden ser medidos en un instante de tiempo o sobre una ventana; y pueden ser medidos tanto en una sola traza como en un conjunto de ellas. El anlisis de atributos incluye la evaluacin de varios parmetros de reservorios, incluyendo indicadores de hidrocarburos. (www.glossary.oilfield.slb.com) Un atributo es necesariamente una derivada de una medida ssmica bsica. Todos los atributos de horizontes y formaciones disponibles no son independientes de los otros, sino simplemente son diferentes maneras de presentar y estudiar una cantidad limitada de informacin bsica. Esta informacin bsica puede ser tiempo, amplitud, frecuencia y atenuacin y de esta forma establecemos la base para nuestra clasificacin de los atributos. . (Brown, 1996). Como una extensa generalizacin, los atributos derivados de tiempo proveen informacin estructural, los atributos derivados de amplitud proveen de informacin estratigrfica y de reservorio. Los atributos derivados de la frecuencia, an no son bien entendidos pero existe una amplia tendencia de que ellos proveen alguna informacin estratigrfica o de reservorio adicional. Los atributos relacionados a atenuacin, no se usan
Orietta C. Mata P.
17
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
hoy en da, pero existe la posibilidad de que el futuro ofrezcan informacin acerca de la permeabilidad. Ha sido notado que el contenido de amplitud de los datos ssmicos es el principal factor para la determinacin de parmetros fsicos, como impedancia acstica, coeficientes de reflexin, velocidades, absorcin, etc. El componente de fase es el principal factor para determinar la forma de los reflectores, su configuracin geomtrica, etc. (Brown, 1996).
Orietta C. Mata P.
18
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
2.2.5.2.- Atributos post apilamiento. 2.2.5.2.Los atributos post apilamiento pueden ser extrados a lo largo de un horizonte o resumidos sobre una ventana.
Orietta C. Mata P.
19
Orietta C. Mata P.
Figura 4: Clasificacin de los atributos ssmicos segn Alistair Brown. Modificado de Sheriff Dictionary, 2001.
AMPLITUD PRE
Intercepto AVO Gradiente AVO Interc vs. Gradiente Diferencia lejos cerca Factor Fluido
POST
POST
HORIZONTE Tiempo Isocrono Tendencia Residual Buzamiento Azimut Borde Iluminaci n Fase instant nea Coseno de la fase Curvatura
VENTANA
Coherencia Continuidad Semblanza Covariannza Dif. Valle Pico Corr. Max. Buzamiento Corr. Max. Azimut Seal a Ruido Indicador de estrato paralelo Indicador de estrato catico Diferencia de traza
HORIZONTE Amp . De reflexi n Amp . Compuesta Impedancia relativa Fuerza de reflexi n Proporci n de amp . Amp . Vs. background
VENTANA
HORIZONTE Frec . Instant nea Frec . Respuesta Envolvente Frec . Derivada del tiempo
VENTANA
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
20
CAPTULO 2: MARCO TERICO
HBRIDOS Forma de onda ESPESOR GROSOR SELECCI N Amp . M xima Mayor amp . Negativa Amp . Abs. M xima Diferencia pico valle DISTRIBUCI N
Energa tiempo medio Pendiente fierza de reflexin Pendiente de energ a media Proporcin positiva negativa
Mapa total abs. Energa total Promedio absoluto Energa promedio Fuerza de reflexi n promedio Amp. RMS Promedio de amp. Pico Varianza de amp.
ESPESOR GROSOR
Ancho de reflexi n Frec. Inst. promedio Frec. Inst. RMS No. De cruces por cero Frec. De pico espectral 1 frec. Esp. Dom . 2 frec. Esp. Dom 3 frec. Dom Ancho de banda espectral
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Las
mediciones
ssmicas
usualmente
involucran
una
apreciable
incertidumbre y no reflejan directamente una sola propiedad geolgica. Con algunas variables geolgicas, la correlacin con una propiedad en particular en una situacin es suficiente para no esperar una prxima situacin. Los atributos responden generalmente a una variedad de situaciones geolgicas y cambios geolgicos los cuales pueden significar un cambio en la correlacin. El problema es determinar la validez de una correlacin observada, especialmente cuando nosotros no entendemos la fsica que los rodea. Entre las maneras de calcular atributos estn el suavizado y el promediado por ventanas de varios tamaos, encontrando residuos, valores picos, midiendo la distribucin entre una ventana principal, sumas, dispersin, continuidad, bordes, linealidad de la curvatura, gradientes de las derivadas, valores absolutos, cambios en la polaridad, diferencias valles-picos, etc. Estas relaciones pueden ser medidas sobre ventanas, correlaciones, semblanza, covarianza, etc. Los atributos pueden ser medidos sobre una sola traza, un volumen de ellas o en diferentes maneras. Los primeros atributos identificados como tal fueron los atributos de traza compleja unidimensionales (envolvente de la amplitud, fase instantnea, frecuencia instantnea y polaridad aparente). Tambin pueden medirse a lo largo de una superficie seleccionada (atributos de horizonte), como extraccin de amplitud, magnitud de buzamiento, azimut del buzamiento, iluminacin artificial y coherencia. Las transformaciones de atributos reciben algunas veces nombres de propiedades fsicas como porosidad, saturacin de fluido, litologa, estratigrafa o discontinuidad estructural, etc. Usualmente basados en crossplots locales o
Orietta C. Mata P.
21
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
correlaciones locales con el registro de pozo u otras mediciones; ellas pueden tener una aproximacin razonable localmente pero pueden dar valores errneos bajo ciertas circunstancias. Antes de entrar de lleno en la clasificacin, definicin y expresiones matemticas de los atributos, revisemos rpidamente los siguientes conceptos, para entender ms claramente a qu tipos de cambios estn referidos.
(www.glossary.oilfield.slb.com) Amplitud: (1) La diferencia entre el mximo desplazamiento de la onda y el punto de no desplazamiento, o punto nulo. (2) La cantidad de desplazamiento de la ondcula ssmica medida de pico a valle. Frecuencia: La tasa de repeticin de longitudes de ondas completas de seales elctricas, luz, sonido y ondas ssmicas medidas en ciclos por segundos, o Hertz. Fase: Una descripcin del movimiento de, o significado de la comparacin de, ondas peridicas como las ondas ssmicas. Las ondas que tienen la misma forma, simetra o frecuencia y que alcanzan los valores mximos o mnimos simultneamente estn en fase. Ondas que no estn en fase son tpicamente descritas por diferencias angulares entre ellas, como por ejemplo, 180 grados fuera de fase. Las ondculas fase cero son simtricas en forma sobre el tiempo cero mientras que las ondculas no fase cero son asimtricas. Las ondculas no fase cero son convertidas a ondculas fase cero para lograr una mejor resolucin de la datos ssmicos. Conocidos los sintticos de pozos fase cero y los perfiles ssmicos verticales (vertical seismic profiles, VSP) pueden se comparados con ssmica de
Orietta C. Mata P.
22
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
superficie para determinar la fase relativa de la ondcula ssmica. Este conocimiento permite a los datos ssmicos de superficie ser corregidos a fase cero. Ahora bien, segn Turhan Taner, (2001) de acuerdo a sus caractersticas computacionales, los atributos post apilamiento se clasifican en: Atributos Instantneos: son calculados muestra a muestra y representan variaciones instantneas de varios parmetros. Valores instantneos de
atributos como envolvente de traza, sus derivadas, frecuencia y fase pueden ser determinados para la traza compleja. Atributos de ondcula: Esta clase comprende todos los atributos instantneos que son calculados en el pico de la envolvente de la traza y tienen diferentes relaciones directas con la Transformada de Fourier de la ondcula en la vecindad del pico de la envolvente. Por ejemplo, frecuencia instantnea en el pico de la envolvente es igual a la principal frecuencia del espectro de amplitud de la ondcula. Fase instantnea corresponde a interceptar fases de la ondcula. Estos atributos tambin son llamados atributos respuesta. A su vez, estos atributos pueden ser subdivididos en base a su relacin con la geologa: Atributos Fsicos: Relacionados a la cualidades y cantidades fsicas. La magnitud de la envolvente de la traza es proporcional al contraste de impedancia acstica; las frecuencias se relacionan al espesor del estrato, dispersin de onda y absorcin. Las velocidades instantneas y promedio se relacionan directamente con las propiedades fsicas de las rocas. Consecuentemente estos atributos son mayormente usados para clasificacin litolgica y caracterizacin de reservorios.
Orietta C. Mata P.
23
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Atributos Geomtricos: describen las relaciones espaciales y temporales de todos los otros atributos. La continuidad lateral medida por semblanza es un buen indicador de similaridad de los estratos as como de discontinuidad. El buzamiento de los estratos y la curvatura brindan informacin depositacional. Los atributos geomtricos son usados para interpretacin estratigrfica desde definir las caractersticas del evento y sus relaciones espaciales, y pueden ser usadas para cuantificar rasgos que asisten directamente en el reconocimiento de los patrones depositacionales y litologa relacionada. La mayora de los atributos, instantneos o de ondcula, son funciones de las
caractersticas de la ondcula ssmica reflejada. Sin embargo, nosotros podemos subdividir los atributos en dos categoras, como sigue a continuacin:
Atributos de Reflexin: corresponden a las caractersticas de la interfase. Todos los atributos instantneos o de ondcula pueden ser incluidos bajo esta denominacin.
Atributos de Transmisin: Relacionados a las caractersticas de un estrato entre dos interfases. Velocidades intervlicas, RMS y promedio, Q, absorcin y dispersin entran bajo esta categora. Debido a que los atributos individuales pueden ser representativos de varias
posibles condiciones, se puede minimizar esta incertidumbre inherente, o no unicidad, mediante la combinacin de diversos atributos en una manera lgica. Los atributos individuales que miden solo una cantidad son denominados Atributos Primitivos. Estos atributos primitivos pueden ser combinados lgicamente,
Orietta C. Mata P.
24
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
estadsticamente o matemticamente para formar Atributos Hbridos. (Brown, 1996) Un estudio ms detallado de las caractersticas, rasgos, comportamiento y definiciones de los atributos ssmicos post apilamiento ms usados se presenta a continuacin.
Orietta C. Mata P.
25
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
La tcnica de la Transformada de Hilbert es una tcnica para determinar la fase de una funcin de fase mnima a partir de su espectro de poder, implicando un clculo con operador de deconvolucin. Dado el espectro de poder P(f) y una ondcula de fase mnima, la representacin de la ondcula en el dominio de la frecuencia W(f) es:
W ( f ) = A( f )e j ( f ) = P( f )
1/ 2
e j ( f )
(6)
La amplitud A(f) es la raz cuadrada del espectro de poder. Tomando el logaritmo a ambos lados de la ecuacin dividimos la funcin en una parte imaginaria y una parte real:
ln[W ( f )] = (1 / 2) ln[ P( f )] + j ( f )
(7)
Para ser de fase mnima, la funcin debe ser analtica en el plano inferior. Entonces la Transformada de Hilbert puede ser usada para hallar la fase (f) a partir del ln P(f)/2:
( f ) = (1 / 2) ln[ P( f )] * 1 /(f )
(8)
Orietta C. Mata P.
26
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Tiempo
Tiempo
Tiempo
Tiempo
za Tra
po m Tie
Figura 5. Anlisis de traza compleja. Trazas real (a) y cuadratura (b) para un porcin de la traza ssmica. La envolvente es mostrada como una lnea punteada en (a) y (b). La fase instantnea es ploteada en (c), la frecuencia instantnea en (d), y la frecuencia promedio como curva punteada en (d). (e) es el diagrama isomtrico de la traza compleja. Modificado de Sheriff, 2001
Como las amplitudes y las fases son conocidas, la Transformada de Fourier puede ser calculada y la expresin en el dominio del tiempo para W(t) puede ser determinada. (Sheriff, 2001).
Orietta C. Mata P.
27
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Una frmula equivalente para f*(t) es dada por la Transformada de Hilbert (Rabiner y Gold, 1975; citado en Oria, 1997):
f * ( ) =
+
P.V .
f (t ) dt t
(9)
La tcnica de la Transformada de Hilbert permite generar la traza cuadratura a partir de la traza real observada. Esto permite que ambas porciones (la real y la imaginaria) estn disponibles para el anlisis. Sin embargo, la Transformada de Hilbert es muy sensible a pequeas variaciones de amplitud de la traza ssmica, consecuentemente su aplicacin irrestricta para el clculo de atributos puede generar a su vez ruido de alta frecuencia indeseable en las secciones de fase y frecuencia instantneas. (Sierra, 1990; citado en Oria, 1997). La Transformada de Hilbert puede ser expresada en forma digital como:
f * (t ) =
A
1 e jn f (t nt ) n
n0
(10)
y f * (t ) =
sen 2 (n / 2) f (t t ) n
f * (t ) =
f (t t )
sen 2 (n / 2) ,n 0 n
(11)
Orietta C. Mata P.
28
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
El operador de Hilbert normalizado en el tiempo (ecuacin A) es singular, desvanecindose en n, y decreciendo monotnicamente en magnitud mientras |n| crece para n singulares. (Oria; 1997). MatLab posee una funcin (y=hilb(x)) que calcula automticamente la Transformada de Hilbert, obteniendo de una manera rpida y sencilla la parte real de la traza y su cuadratura, evitando as esta serie clculos complejos, en los cuales se invierte mucho tiempo.
comportamiento estratigrfico a partir de los datos ssmicos, antes de la estimacin de los parmetros instantneos, el contenido de amplitud y frecuencia de la seal ssmica debe ser preservada en cada paso del procesamiento. Alguna variacin en la forma de la onda bsica que no es atribuible al subsuelo debe ser eliminada. Los
Orietta C. Mata P.
29
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
mltiples y todo tipo de ruido aleatorio limitan la integridad de los resultados.(Oria, 1997). El atributo Amplitud Instantnea es independiente de la fase. Provee informacin acerca de contrastes de impedancia acstica. Las reflexiones fuertes muchas veces estn asociadas a cambios litolgicos importantes entre estratos de roca adyacente, as como a travs de discontinuidades y lmites asociados con cambios agudos en el nivel relativo del mar o ambientes depositacionales. Cambios laterales en la amplitud instantnea (fuerza de reflexin) tambin pueden estar asociados con presencia de fluidos Variaciones laterales en el espesor de las capas cambian la interferencia de las reflexiones produciendo cambios graduales laterales en la potencia de reflexin. Cambios agudos locales pueden indicar fallamientos o acumulaciones de hidrocarburos donde las condiciones de entrampamiento son favorables.
Acumulaciones de hidrocarburos, especialmente gas, pueden mostrar reflexiones de gran amplitud que se conocen con el nombre de bright spots o puntos brillantes.
(Famiglietti, 1999).
El mdulo de la traza compleja F(t) se ha denominado Amplitud Instantnea y se expresa de la siguiente forma:
A(t ) = F (t ) = ( f 2 (t ) + h 2 (t ))1 / 2
(12)
La amplitud compuesta es la suma del valor absoluto de las amplitudes de las reflexiones identificadas en el tope y la base del reservorio u otro intervalo. La impedancia acstica derivada de la amplitud por inversin ssmica es otra manera
Orietta C. Mata P.
30
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
de combinar informacin del tope y la base del reservorio (con limitaciones de espesor). La amplitud derivada de la traza compleja es la fuerza de reflexin, amplitud instantnea, o envolvente de amplitud. Los atributos de amplitud calculados de alguna manera sobre una ventana de tiempo son muchos y variados. La mayora de ellos localizan propiedades netas del intervalo estudiado, algunos extraen una seleccin de las amplitudes dentro de una ventana, y algunos se avocan a
Figura 6. La funcin amplitud instantnea calculada a partir de una traza ssmica compleja. Tomado de Oria, 1997
Orietta C. Mata P.
31
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
(t ) = tg 1
h (t ) f (t )
(13)
La funcin fase incrementa montonamente con el tiempo de reflexin. Debido a que la fase es independiente de la fuerza de reflexin, con frecuencia los eventos dbilmente coherentes se clarifican. Sin embargo, en intervalos donde las amplitudes de reflexin son bajas, la determinacin de la fase mediante un examen visual de los datos es bastante difcil y en algunos casos, imposible. La fase instantnea tiende a ser lineal alrededor de la amplitud mxima.
Figura 7. La funcin fase instantnea calculada a partir de una traza ssmica compleja. Tomado de Oria, 1997
Orietta C. Mata P.
32
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
El atributo fase instantnea permite observar la presencia de estructuras en una seccin ssmica, como por ejemplo: fallas, discontinuidades, acuamientos y eventos con diferentes buzamientos. En ciertos yacimientos de gas, la fase instantnea puede ser usada para identificar y cartografiar reversiones de fase, las cuales podran indicar contenido de gas. (Famiglietti, 1999). Tambin se ha definido la Frecuencia Instantnea como:
(t ) =
(t ) (t )
(14)
Es decir, la derivada de la fase instantnea con respecto al tiempo. (Adriana Valerio, 2001) La frecuencia instantnea puede proveer informacin acerca de la frecuencia caracterstica de los eventos, efectos de absorcin y fracturamiento y espesores depositacionales. Cuando el comportamiento de la fase de una traza ssmica compleja puede ser calculado como una funcin continua en el tiempo, la rapidez del cambio de la fase puede entonces ser calculada, punto a punto y por consiguiente, el comportamiento de la frecuencia se conoce instantneamente. La frecuencia instantnea puede tener un grado de variacin que puede estar relacionado con la estratigrafa. Variaciones como las de los acuamientos y los bordes de interfases de fluidos (petrleo / agua), tienden a cambiar la frecuencia instantnea mas rpidamente. Un cambio hacia frecuencias bajas se observa muchas veces en
Orietta C. Mata P.
33
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
reflexiones de horizontes correspondientes a arenas con gas condensado y yacimientos de hidrocarburos. (Famiglietti, 1999).
Figura 8. La funcin frecuencia instantnea calculada a partir de una traza ssmica compleja. Tomado de Oria, 1997
Orietta C. Mata P.
34
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Espesor de litologa. Porosidad Arenas fluviales o deltaicas. Ciertos tipos de arrecifes. Discordancias. Efectos de entonacin. Cambios en la estratigrafa secuencial.
Cambios laterales en la amplitud han sido usados en estudios estratigrficos para separar reas de estratigrafa concordante de estratos caticos. En general, los estratos que son concordantes tienen amplitudes ms altas mientras que los estratos caticos tienen amplitudes ms bajas. Ventana de anlisis: Generalmente, se puede escoger una ventana que cubra una zona con pocos reflectores (entre 50 y 200ms), sobre los cuales se espera encontrar una anomala de amplitud o detectar variaciones laterales de amplitud. Para detectar cambios laterales de estratigrafa, se debe concentrar la ventana en un intervalo especfico de la secuencia de inters (tpicamente entre 100 y 200ms). Para detectar anomalas de amplitud, como puntos brillantes asociados a gas, se debe usar una ventana localizada, de aproximadamente entre 50 y 100ms. Para estudios especiales, se puede tomar una ventana extremadamente pequea. Concentrndose en una sola ondcula se puede proveer informacin del promedio de las caractersticas del reflector. Los atributos ms comunes dentro de esta categora son:
Orietta C. Mata P.
35
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
1 n 2 ai N i =1
Donde
1/ 2
(15)
n = nmero de muestras
1 n ai N i =1
Donde n = nmero de muestras
(16)
amplitud positiva de la traza se efecta un ajuste parablico. La salida del atributo es el mximo valor de la curva interpolada.
Orietta C. Mata P.
36
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
1 n ai ai > 0 N i =1
(17)
amplitud.
1 n ai ai < 0 N i =1
(18)
negativo de amplitud y la salida del atributo es el valor mximo negativo sobre la curva interpolada.
Amplitud Total Absoluta: para cada traza se suman el valor absoluto de las
(19)
Amplitud Total: para cada traza se calcula la amplitud total (suma aritmtica)
(20)
Orietta C. Mata P.
37
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
n 2 ai i =1 N
(21)
n 2 ai i =1
Amplitud Promedio:
(22)
Amp. principal =
a
N0
(23)
promedio de todas las muestras en la ventana. Luego sustrae el promedio a cada muestra y calcula el valor de la diferencia al cuadrado. La salida es la suma de todas estas variables divididas por el nmero de muestras.
V =
1 ( xi x)2 N
(24)
promedio de todas las muestras en la ventana. Luego se sustrae el promedio de cada ventana y posteriormente se resta el valor del promedio a cada muestra y calcula el valor de la diferencia al cubo. La salida es la suma de todos estos valores divididos por el nmero de muestras. (25)
Orietta C. Mata P.
38
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
S=
1 ( xi x)3 N
promedio de todas las muestras en la ventana. Luego se sustrae el promedio de cada ventana y despus sustrae el valor del promedio a cada muestra y calcula el valor de la diferencia a la cuarta potencia. La salida es la suma de todos estos valores divididos por el nmero de muestras.
K=
1 ( xi x)4 N
(26)
Orietta C. Mata P.
39
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Estadstica es una disciplina acadmica moderna, la cual provee de procedimientos cientficos para recopilar, organizar, resumir y analizar la informacin cuantitativa del tipo que comnmente se encuentra en las ciencias del comportamiento. Primero la estadstica permite el resumen y la presentacin de grandes cantidades de informacin en una manera que facilite su comunicacin e interpretacin. Segundo, la estadstica le permite al cientfico extender su investigacin mucho mas all de los parmetros restringidos, en la mayora de los cuales las investigaciones son desarrolladas. En la estadstica clsica se asume que todas las muestras estn tomadas aleatoriamente e independientemente de una sola distribucin. La posicin de las muestras o la relacin entre ellas no son consideradas.(Taheri y Viloria, 1998) Existen dos tipos principales de estadsticas: la estadstica descriptiva y la estadstica inductiva. La estadstica descriptiva es aquella que permite describir con precisin una serie de informacin cuantitativa; mientras que la estadstica inductiva es la que le permite al investigador establecer realizar inferencias y generalizaciones de pequeos grupos llamados muestra a un grupo mayor denominado poblacin, todo esto dentro de un muy bien definido grado de confianza. Una muestra es un pequeo grupo relativo de individuos escogidos en una manera cientfica para representar a un grupo mayor de los mismos individuos (este grupo mayor se define como poblacin), el cual el investigador est interesado en estudiar. En la estadstica inductiva, el trmino estadstica es usado de una manera especial. Un parmetro es alguna propiedad numrica de una poblacin. Entonces una estadstica es una propiedad numrica de una muestra, la cual es usada para estimar el valor del correspondiente parmetro de la poblacin. (Roscoe, 1974)
Orietta C. Mata P.
40
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Una poblacin puede ser finita o infinita. Por ejemplo, la poblacin consistente de todos los cerrojos producidos en una fbrica en un determinado da es finita, mientras que la poblacin formada por todos los posibles sucesos (caras, cruces) en tiradas de una moneda es infinita. Si una muestra es representativa de una poblacin, se pueden deducir importantes conclusiones acerca de sta, a partir del anlisis de la misma. (Spiegel, 1961).
Orietta C. Mata P.
41
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
42
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
A
Moda
Media Mediana
B
Moda Mediana Media
Figura 9. (A) Media, moda y mediana para distribuciones simtricas. (B) Media, mediana y moda para distribuciones no simtricas. Modificado de Roscoe, 1975
Si una serie de datos se colocan en orden de magnitud, el valor medio (o la media aritmtica de los dos valores medios) que divide al conjunto de datos en dos partes iguales es la mediana. Por extensin de esta idea se puede pensar en aquellos valores que dividen a los datos en cuatro partes iguales. Estos valores representados por Q1, Q2 y Q3 se llaman primero, segundo y tercer cuartil respectivamente; donde el valor de Q2 es igual a la mediana.
Orietta C. Mata P.
43
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Anlogamente los valores que dividen los datos en diez partes iguales se llaman deciles y se representan por D1, D2,,D9; mientras que los valores que dividen los datos en cien partes iguales se llaman percentiles y se representan por P1, P2, P3,P99. El quinto decil y el quincuagsimo percentil se corresponden con la mediana. Los percentiles P25 y P75 se corresponden con el primer y tercer cuartel, respectivamente. En conjunto, cuarteles, deciles, percentiles y otros valores obtenidos por subdivisiones anlogas de los datos se llaman cuntiles. (Siegel, 1961)
Orietta C. Mata P.
44
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Q=
Q3 Q1 ; 2
(27)
donde Q1 y Q3 son el primer y el tercer cuartel de los datos. El rango intercuartlico se emplea a veces, pero el rango semiintercuartlico es ms utilizado como medida de dispersin. El rango entre percentiles de un conjunto de datos viene definido por: Rango Percentil 10 90 =P90 P10, Donde P10 y P90 son los percentiles dcimo y nonagsimo de los datos. La desviacin media es una medida de la dispersin igual al promedio de los valores absolutos de la desviacin de los valores. Su expresin matemtica es la siguiente:
MD =
(28)
(29)
donde M es la media aritmtica de los nmeros y |Xi M| es el valor absoluto de las desviaciones de las diferentes Xi de M. La desviacin estndar o desviacin tpica es la raz cuadrada del promedio de las desviaciones media, y su expresin matemtica viene dada por:
1 n ( xi x) 2 n i =1
(30)
Finalmente, varianza se define como el cuadrado de la desviacin tpica y viene dada por:
2 =
1 ( xi x)2 n
(31)
Orietta C. Mata P.
45
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 3.
METODOLGICO. MARCO METODOLGICO.
A lo largo de este captulo se realiza una breve descripcin de la metodologa utilizada durante la realizacin de este trabajo, es decir, una explicacin detallada de cada una de las etapas que se emplearon para programar los diferentes algoritmos que permiten cumplir con los objetivos planteados.
Orietta C. Mata P.
46
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
47
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
nombres de los registros como los valores numricos para cada corrida contenidos en el archivo. ASPASS posee una interfase que permite establecer el X asociado al registro de profundidad y acta segn sea el caso. El usuario puede introducir registros obtenidos cada 0.25ft, cada 0.5 ft o cada pie y para cada uno de ellos actuar de manera similar para llegar al mismo resultado. Para el caso de X=0.25, se dise un algoritmo que aplica un filtro numrico a los datos, basado en el Mtodo de Diferencias Divididas de Newton, con el cual se reduce el tamao del archivo a la cuarta parte (tratando de conservar la mayor cantidad de informacin posible) haciendo el programa mas eficiente al minimizar la memoria de programa utilizada. Bsicamente el filtro coloca los valores de los registros a un pie de separacin, para lo cual aplica una interpolacin cada 4 muestras y le asigna el valor al punto medio. Para X=0.5 el algoritmo diseado es mas sencillo, ya que solo toma los valores
cada dos muestras, los promedia y asigna este nuevo valor calculado a la primera muestra. Igualmente, la salida esperada del algoritmo ser los registros correspondientes con muestras separadas cada pie. Si X=1, no se hace ninguna modificacin a los datos, ya que lo que se busca es precisamente eso, que los registros estn espaciados a un pie de separacin.
Orietta C. Mata P.
48
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
3.3.1.- Clculo de la escala en tiempo. .1.Despus de lograr los registros espaciados cada pie, se hace ms fcil el clculo de la funcin tiempo profundidad, a partir de la cual se obtienen los registros en escala de tiempo y finalmente el sismograma sinttico. Para realizar el clculo del sinttico es necesaria la utilizacin de los registros de profundidad y el registro snico. Gracias a que X=1, el registro snico est conformado por las velocidades intervlicas entre dos pies consecutivos, por lo que para su uso es preciso hacer una suma acumulativa a fin de obtener las velocidades para cada pie. Adems se necesita introducir una velocidad inicial (que para los efectos de diseo de la interfase se utiliz un valor correspondiente a 10.000 pie/s y el cual el usuario puede cambiar de acuerdo a su necesidad) que tambin es sumada a los datos del registro. Esta velocidad corresponde a la velocidad inicial por debajo de la primera velocidad registrada, es decir, para este caso corresponde a la velocidad intervlica desde profundidad 0 pie hasta la profundidad de 4300 pies; y puede ser definida en la misma interfaz grfica donde se define el valor de X.
Orietta C. Mata P.
49
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
A partir de la ecuacin
t=
correspondiente a los datos de profundidad del registro, consiguiendo los tiempos asociados a cada pie. El siguiente paso es muestrear la escala de tiempo a 1milisegundo (t igual a 1), por lo que se utiliza una funcin que redondea los valores al entero ms cercano. De esta manera se consiguen varios valores de profundidad para un mismo tiempo, siendo ahora necesario tener un solo valor de profundidad por milisegundo. Para lograr esto, se emplea un algoritmo que promedia acumulativamente los valores de profundidad para cada tiempo, es decir, el algoritmo busca los dos primeros valores asociados a un mismo tiempo, los promedia y verifica si el valor que sigue en el registro tambin posee el mismo tiempo asociado. De ser as, entonces el algoritmo promedia el valor obtenido anteriormente con el nuevo valor hallado y as sucesivamente hasta obtener un valor de profundidad por segundo. Esta operacin aplicada al registro de profundidad se aplica simultneamente a todos los registros cargados, teniendo finalmente un valor de cada registro asociado a cada t. Si el usuario lo define as en los parmetros iniciales, luego de obtener los registros en tiempo cada milisegundo, se remuestrean los registros a 2ms (esto es t igual a 2). Empleando una rutina que lee los valores de los registros cada 2 muestras y crea una columna aparte con estos datos, se logran solo los valores asociados a cada 2ms. Finalmente, el tiempo obtenido se multiplica por dos ya que se est considerando tiempo doble de viaje.
Orietta C. Mata P.
50
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
acsticas, 3.3.2.- Clculo de las impedancias acsticas, del coeficiente .3. de reflexin y del sismograma sinttico.
Para el clculo de las impedancias se obtienen primero las velocidades instantneas para cada milisegundo, las cuales se tienen al invertir el registro acumulado del snico. A partir de la frmula, V = velocidad y =densidad. El coeficiente de reflexin se calcula utilizando las impedancias obtenidas anteriormente a partir de la relacin R1 =
2 1 . Para obtener la serie de reflectividad, 2 + 1
Z = .V
basta con asociar los coeficientes obtenidos a los valores de tiempo. La convolucin de la serie de reflectividad con una ondcula, finalmente da como resultado el sismograma. En este aspecto debemos hacer dos acotaciones importantes: la primera, la ondcula utilizada en este caso es la ondcula Ricker, la cual es una ondcula fase cero (a igual intervalo de muestreo de la serie de reflectividad) calculada a partir de una rutina programada en MatLab por CREWES Project; y la segunda, , la funcin utilizada para realizar la convolucin es una rutina provista por MatLab denominada convz.m, la cual est especialmente diseada para realizar la convolucin de una traza con una ondcula fase cero, evitando el retraso (delay) que se produce al convolver la ondcula y la traza mediante la funcin de convolucin convencional.
Orietta C. Mata P.
51
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
52
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
para lo que se dise un algoritmo especfico que permite tomar los valores de tiempo del sismograma y los asocia a los valores de amplitud. Estos valores de amplitud seleccionados son introducidos en las respectivas frmulas de los atributos y la salida del algoritmo es un valor promedio para cada ventana de anlisis, mostrado relacionado al valor de tiempo promedio de la ventana, es decir, se muestra asociado al valor medio de la ventana en el eje tiempo. En la interfase diseada, adems, se muestran los valores calculados para cada atributo y se puede modificar la escala de cada pista, para apreciar mejor la diferencia que existe en un mismo atributo calculado para la traza original y para la traza modificada. Los valores tope y base de la ventana de estudio tambin se presentan en la interfase, a manera solo de referencia, para ofrecerle un usuario una informacin mas detallada del intervalo en tiempo que se est trabajando. Entre los atributos de amplitud programados para este mdulo encontramos: amplitud RMS, amplitud absoluta promedio, amplitud pico promedio, amplitud valle promedio, energa total, amplitud total absoluta, amplitud total, energa promedio, amplitud promedio, varianza en amplitud, simetra en amplitud y kurtosis en amplitud.
Orietta C. Mata P.
53
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
que se realizaran las modificaciones. La rutina asocia los valores de tiempo escogidos a los valores del registro mostrado, definiendo el tope y la base de la ventana. Estos dos valores de promedian para obtener un valor ponderado al cual se le suma y se le resta una pequea diferencia para as obtener dos valores relacionados al promedio inicial. Con una funcin aleatoria se producen el mismo nmero de valores que posee la ventana, los cuales oscilan entre los dos promedios ponderados calculados anteriormente. Estos valores vendrn a sustituir los valores del registro original, produciendo las modificaciones deseadas por el usuario del programa. Esta edicin en tiempo podr introducir un pequeo margen de error, el cual no es significante en comparacin con la longitud total del registro. Los nuevos registros obtenidos a partir de la edicin, se almacenan en un nuevo vector, (denominados copia1, copia2 y copia3 dependiendo de en que pista se ubica el registro), lo que implica que los valores originales siempre estn disponibles si se quiere volver al sismograma sinttico original. Posteriormente, se puede recalcular el sinttico tomando como referencia los valores de copia1, copia2 y/o copia3, calculando un sismograma igual al original a excepcin de la ventana que fue modificada. Este nuevo sismograma se almacena de igual manera en otro vector para que los valores del sismograma original puedan ser reutilizados, bien sea para desplegarlos por pantalla o para realizar una nueva edicin. Con este sismograma producto de la edicin en tiempo, es posible recalcular todos los atributos ssmicos disponibles, con lo que ahora se establecen puntos de comparacin y se da paso al ltimo mdulo de ASPASS.
Orietta C. Mata P.
54
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Para facilitar la delimitacin grfica del intervalo modificado, los nuevos registros que aparecen por pantalla, ahora presentaran dos lneas rojas horizontales que marcan el alcance en tiempo del intervalo modificado, y todos los clculos derivados realizados con estos registros, arrastraran esta delimitacin en tiempo.
3.6.-
Herramientas
de
la
interfase
ASPASS
Atributos
Instantneos. Instantneos.
Dentro de la interfase de atributos instantneos se program una herramienta til que permite observar mejor los cambios ocurridos dentro de los valores que se presentan en la interfase.
Orietta C. Mata P.
55
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
56
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
En estas rutinas se toman en cuenta los valores calculados durante los otros tres mdulos y por una serie de operaciones aritmticas sencillas (como sumas, restas y reglas de tres) se calculan los porcentajes de los valores modificados con respecto a los originales. Es decir, si el valor original es considerado 100%, entonces el porcentaje del modificado es obtenido por una regla de tres. La diferencia de 100 menos el porcentaje del modificado, es justamente la diferencia o variacin porcentual de los datos. Para los atributos de amplitud (atributos de ventana), fue relativamente fcil implementar este procedimiento, ya que se dispone de dos puntos para cada atributo, original y modificado respectivamente. Ahora bien, con los atributos instantneos es necesario implementar otra metodologa que permita, antes del clculo de los porcentajes, obtener un valor puntual por traza para la ventana de anlisis. Para ello, se hizo uso del concepto de varianza, con el cual se resume la informacin de toda una ventana en un punto representativo del resto de los valores. Calculando un valor de varianza para la traza original y luego un valor de varianza para la traza modificada en el intervalo, ya se disponen de los dos puntos necesarios para realizar los clculos porcentuales. Los resultados obtenidos se resumen en un grfico de barras, Atributo vs. Porcentaje, donde se identifican rpidamente aquel (aquellos) atributo(s) mas sensibles a las variaciones realizadas en los registros de pozo durante la edicin en tiempo. Adicionalmente se program una tabla, como parte de esta herramienta, donde se visualizan las cifras correspondientes a los valores de cada atributo, permitiendo no solo identificar las variaciones visualmente sino adems matemticamente.
Orietta C. Mata P.
57
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
ASPASS / Principal
INICIO
Introduccin de los parmetros bsicos del muestreo de los datos, velocidad inicial y los parmetros de la ondcula.
No
Si
Modificar registros de pozo en un intervalo de inters Grficos de los registros en las pistas
Orietta C. Mata P.
58
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
INICIO
No
Sinttico Original
Si
Presentacin de los atributos calculados por pantalla Inicio del mdulo atributos por ventana Inicio del mdulo estadstico
Orietta C. Mata P.
59
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
INICIO
Presentacin del sismograma por pantalla Definir ventana sobre el intervalo de inters
No
Sinttico Original
Si
Presentacin de los atributos calculados por pantalla Inicio del mdulo atributos por ventana Inicio del mdulo estadstico
Orietta C. Mata P.
60
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
ASPASS / Estadsticas
INICIO
No Mostrar Tabla
Si
FIN
FIN
Figura 10. Flujogramas del programa ASPASS. (A) Flujograma del mdulo ASPASS / Principal. (B) Flujograma del mdulo ASPASS / Atributos Instantneos. (C) Flujograma del mdulo ASPASS / Atributos de Amplitud. (D) Flujograma del mdulo ASPASS / Estadsticas
Orietta C. Mata P.
61
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 4.
C CARGA DE LOS DATOS DE POZO Y CLCULO DE SISMOGRAMAS SINTTICOS.
ASPASS posee interfaces interactivas para cada una de sus funciones. Desde cargar los datos hasta desplegar el sinttico y calcular los respectivos atributos, ASPASS ofrece un modo amigable y sencillo de interactuar con el usuario. Es importante recordar, que antes de realizar cualquier accin, se debe chequear que el directorio actual de MatLab sea el adecuado, es decir, colocar en el directorio actual la ubicacin exacta desde donde estn almacenados todos los algoritmos. Para una mayor comodidad, se recomienda tambin que los datos a emplear se encuentren dentro de este mismo directorio, para as evitar mayores complicaciones.
Orietta C. Mata P.
62
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Inmediatamente se despliega la pantalla del men principal, desde la cual se realizan todas las operaciones relacionadas al clculo del sismograma sinttico. Antes de iniciar el proceso de carga de cualquier archivo es importante sealar que los nombres de los registros de pozo deben poseer ciertas siglas dentro del nombre, lo que hace la carga de los datos eficiente. El registro de profundidad debe tener las siglas DEPTH, el snico DT, el registro de densidad RHOB y el gamma ray GR. Esto con la finalidad de que los algoritmos asociados a la verificacin de los registros necesarios para el clculo del sismograma sinttico, funcionen correctamente y no den ningn error. As mismo, en caso de querer cargar un archivo ASCII, se debe tener en cuenta una
Orietta C. Mata P.
63
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
consideracin especial: se debe eliminar del archivo el encabezado con la informacin de pozo y en la lnea que contiene los nombres de los registros se debe colocar al inicio los caracteres ~A, con la finalidad de que el algoritmo diseado para la lectura de datos no emita ningn error. A partir de esta pantalla se realiza el proceso de carga de los datos. Para ello, se presiona en la barra de men la opcin Archivo y luego la opcin Cargar. Deben aparecer dos opciones: Archivo LAS y Archivo ASCII. Se selecciona la opcin deseada y seguidamente se abrir un cuadro como el que se muestra a continuacin para seleccionar el archivo que contiene los datos de pozo.
Orietta C. Mata P.
64
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
65
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Luego de seleccionar el archivo a abrir, se despliega una lista que contiene los nombres de los registros que se encuentran dentro de este. El usuario podr escoger cules registros cargar y continuar con el proceso. A continuacin el programa preguntar si se desea cargar otro archivo, lo cual es muy til en el supuesto caso de que la informacin requerida se encuentre en archivos separados.
El siguiente paso consiste en organizar los registros e introducir los parmetros iniciales de muestreo de los datos. Para esto, la interfase muestra listas que contienen los
Orietta C. Mata P.
66
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
nombres de los registros seleccionados, los cuales deben ser colocados en el campo que lleve su nombre. Adems de ser una manera eficiente de organizar los datos, es una manera efectiva de comprobar si los registros cargados realmente son los necesarios para realizar el clculo del sismograma sinttico y la conversin tiempo profundidad. En caso de no ser as, se desplegar un cuadro de advertencia para notificarle al usuario que ha ocurrido un error. Por otra parte, los parmetros iniciales que deben introducirse corresponden a la distancia X en el registro de profundidad, el intervalo t al que se van a muestrear los datos; y la velocidad inicial, que corresponder a la velocidad intervlica desde profundidad cero pies hasta el valor de la primera medicin.
Figura 16. Interfase que permite organizar los registros e introducir los parmetros iniciales
Orietta C. Mata P.
67
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 17. Mensaje de error que se despliega al introducir un registro en el campo errneo
Al presionar el botn Cargar, se mostrar la interfase que permite al usuario introducir los parmetros de la ondcula usada en la generacin del sinttico. Para la ondcula Ricker, es necesario introducir la frecuencia necesaria (en Hz), el intervalo de muestreo, el cual se recomienda que sea el mismo intervalo de muestreo que se selecciono para los datos y finalmente la cantidad de muestras por ondcula. En caso de que la ondcula no sea Ricker, se ofrece una lista con diferentes opciones de ondculas fase cero que el usuario puede usar a manera de comparacin cualitativa. Estas ondculas son: meyer, morlet y mexican hat y como parmetro slo se debe incluir el nivel de suavizado
Orietta C. Mata P.
68
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
(refinement) de la funcin. Presionando el botn Mostrar se puede observar la ondcula diseada mediante los parmetros introducidos.
Cuando presionemos el botn Ok automticamente se inicia el proceso de carga de los datos, el cual tarda unos pocos segundos. Cuando este proceso finaliza, se retorna a la pantalla de men principal de ASPASS.
Orietta C. Mata P.
69
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Ahora bien, si se desea abrir una sesin de trabajo previa los pasos a seguir son ms sencillos. Basta con presionar la opcin de la barra de mens Archivo y luego Abrir. Inmediatamente aparecer un cuadro de dilogo donde se puede seleccionar el archivo a abrir y posteriormente los registros en las pistas correspondientes.
Visualizacin c 4.2.- Visualizacin de los registros y clculo del sismograma sinttico. sinttico.
El primer paso a seguir a la hora de calcular el sismograma sinttico es visualizar los registros introducidos en la pantalla del men principal de ASPASS. Para esto, solo se debe presionar la opcin Edicin de la barra de mens, y seguidamente la opcin Visualizar Registros o tambin se puede presionar directamente el botn Visualizar Registros que se encuentra en la pantalla principal. Inmediatamente se abrir un cuadro como el que se muestra a continuacin:
Orietta C. Mata P.
70
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 20. Cuadro de Curvas. Interfase que permite al usuario organizar los registros en las pistas y colocar la apariencia de estos.
Orietta C. Mata P.
71
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Esta interfase permite al usuario ordenar los registros de pozo en las pistas y adems colocarle la apariencia deseada. Se pueden colocar varios registros en una misma pista, o si se prefiere uno en cada una de ellas. Tambin se puede seleccionar el color, el estilo y el grosor de los registros. Eventualmente, pero no para esta parte del proyecto, se habilitaran las opciones de Log,Grid e Inv, con las cuales se podrn personalizar mas an los registros graficados en las pistas. Luego de escoger los registros y asignarles una determinada apariencia, se mostraran en la pantalla principal de ASPASS de la siguiente manera:
Orietta C. Mata P.
72
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Ahora bien, en caso de que se desee reorganizar los registros en las pistas o cambiar las caractersticas de uno de ellos, se usa la opcin Reordenar Registros, de la barra de mens, Edicin. Esta borrar automticamente los registros que se tenan graficados y volver a graficar los registros seleccionados con las modificaciones correspondientes. Los registros que no se hayan reordenado, aparecern de nuevo en la pantalla en las mismas condiciones iniciales. Sin embargo, en caso de que se haya editado algn registro en tiempo, al aplicar la opcin Reordenas Registros se perder cualquier cambio que se haya realizado. La interfase de Reordenar Registros es exactamente igual a la Visualizar Registros, con la excepcin que cuando se abra no aparecer en blanco sino que mostrar la informacin previamente cargada. Para calcular el sismograma sinttico se pueden utilizar dos maneras: seleccionando la opcin Sintetico de la barra de mens y luego la opcin Calcular Sinttico o presionando directamente el botn Calcular Sinttico ubicado en la pantalla principal de ASPASS.
Orietta C. Mata P.
73
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Inmediatamente aparecer en la pantalla del men principal el sismograma sinttico calculado a partir de los registros densidad y snico.
Figura 23. Sismograma calculado por las rutinas de ASPASS a partir de los registros de densidad y snico.
Orietta C. Mata P.
74
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Existen otras opciones relacionadas al sinttico, mediante las cuales se puede modificar su forma y su apariencia. Una de ellas permite cambiar la ondcula fase cero utilizada para generar el sinttico. Presionando la opcin Modificar ondcula desde la pantalla o utilizando la barra de mens a travs de Sintetico, Parmetros de la ondcula, se desplegar el la interfase Parmetros de la ondcula en la cual se puede seleccionar una de las opciones que all aparecen. La otra opcin permite cambiar el color, el grosor o el estilo de lnea bajo el cual aparecer la sismograma en la pantalla. Se puede acceder al cuadro Preferencias del sinttico presionando el botn Preferencias ubicado en la interfase del men principal o en la barra de mens Sintetico, Preferencias.
Figuras 24. Opciones relacionadas a la apariencia y forma del sismograma sinttico. Modificar ondcula y Parmentros de ondcula permiten modificar el tipo de ondcula para la convolucin y Preferencias permite modificar el color, el grosor y el estilo de lnea del sismograma
Orietta C. Mata P.
75
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figuras 25. Opcin Modificar Registros, con la cual se puede realizar la edicin en tiempo.
Al activar esta herramienta, aparecer un cursor que le permite al usuario escoger los valores tope y base del intervalo a modificar. La salida de los algoritmos se mostrar por pantalla de manera similar a la original, y solo se observar variacin en la ventana sealada.
Orietta C. Mata P.
76
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
B
Figura 26. (A) Registros densidad y snico antes de la edicin de los valores en tiempo. (B) Registro densidad y snico despus de la edicin en tiempo.
Orietta C. Mata P.
77
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Rec 4.4.- Reclculo sismogramas sintticos a partir de los registros editados editados en tiempo.
Posterior a la edicin en tiempo, se puede recalcular el sismograma a partir de los registros editados. Mientras que no haya ocurrido una modificacin en los datos de pozo, al presionar la opcin Calcular Sinttico el sismograma que aparecer es el original. Sin embargo, luego de la edicin, al activar Calcular Sinttico se abrir un cuadro de dilogo como se muestra a continuacin:
Figura 27. Cuadro de dilogo que permite escoger en que pista se encuentra cada registro.
Este cuadro permite organizar en que pista se encuentra cada registro, asignando un nmero de pista al registro de densidad y un nmero de pista al registro snico, para facilitar as los clculos del sismograma sinttico, el cual aparecer por pantalla una vez que se halla cerrado el cuadro de dilogo.
Orietta C. Mata P.
78
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
B
Figura 28. Sismograma sinttico calculado por ASPASS. (A) Calculado antes de la edicin de los registros en tiempo. (B) Calculado a partir de los registros de pozo modificados. Notese el leve cambio que existe entre uno y otro en la zona sealada a raz de la modificacin realizada.
Mediante el botn Sinttico Original, ubicado en la pantalla del men principal, se puede volver a visualizar rpidamente el sismograma original, es decir, el sismograma calculado antes de realizar la edicin de los registros de pozo.
Orietta C. Mata P.
79
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 5.
CLCULO DE ATRIBUTOS.
Para poder realizar el clculo de atributos ssmicos es necesario haber cargado los datos de pozo y calculado con anterioridad el sismograma sinttico.
La interfase diseada para el clculo de los atributos instantneos, es ligeramente parecida a la interfase del men principal. En esta pantalla, adems de visualizar los diferentes atributos de traza compleja, se puede aplicar una herramienta que permite
Orietta C. Mata P.
80
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
caracterizar mejor cualitativamente las variaciones observadas en los atributos producidos por diferentes sismogramas sintticos con alguna caracterstica similar.
El primer paso para calcular los atributos es visualizar en la pantalla el sismograma sinttico a partir de cual sern calculados. Presionando el botn Actualizar Sinttico ubicado en la pantalla principal, se mostrar en la pista correspondiente a sismograma sinttico el ltimo sismograma que se visualiz en la pantalla del men principal de ASPASS. Es decir, el botn Actualizar Sinttico bsicamente lo que hace es actualizar el sismograma que se estaba visualizando en la pantalla anterior.
Orietta C. Mata P.
81
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 31. Opcin Actualizar Sinttico la cual permite visualizar el sismograma en la pantalla de Atributos Instantneos
La opcin Sinttico Original, la cual se encuentra justo debajo de Actualizar Sinttico en la pantalla principal de Atributos Instantneos, permite visualizar el sinttico generado a partir de los registros sin edicin. Es una herramienta muy til para comparar los atributos producidos por el sismograma original y el obtenido a partir de los registros editados. Calcular atributos permite calcular los atributos ssmicos de traza compleja y visualizarlos por pantalla, en las pistas correspondientes a cada uno de ellos.
Figura 32. Opcin Calcular atributos, con la cual se visualizan los atributos de traza compleja en la pantalla de Atributos Instantneos
Orietta C. Mata P.
82
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 33. Atributos instantneos calculados a partir del sismograma que se muestra en la pista 1.
Presionado el botn Preferencias Atributos se puede cambiar, de manera similar a las preferencias del sinttico, la apariencia en pantalla de los atributos, modificando el color, grosor o el estilo de la lnea que representa a los atributos. Ahora bien, la barra de mens de esta interfase presenta una opcin de Herramientas la cual puede ser til desde el punto de vista de un anlisis cualitativo.
Orietta C. Mata P.
83
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Restar Sintticos permite la resta de dos sismogramas e inclusive de sus atributos de traza compleja, para visualizar mejor las variaciones que se producen entre unos y otros.
Figura 35. Salida de la herramienta Restar Sintticos, la cual permite restar dos sintticos e inclusive sus atributos para visualizar mejor las variaciones existentes entre uno y otro.
Para poder hacer uso de esta herramienta adecuadamente, se deben haber calculado los atributos asociados a la traza original y a la traza modificada. De lo contrario, la herramienta no podr ser implementada debido a la ausencia de uno de los valores necesarios para el clculo.
Orietta C. Mata P.
84
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 36. Como inicializar el mdulo Atributos por Ventana, Atributos de Amplitud.
Orietta C. Mata P.
85
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Seguidamente aparecer la pantalla de ASPASS / Atributos de Amplitud, la cual bsicamente mantiene el diseo de los otros dos mdulos. Adems de las cuatro pistas caractersticas, (la correspondiente al sismograma y las tres correspondientes a los atributos), esta interfase posee una serie de botones y listas desplegables que interactan con el usuario para hacer ms fcil el clculo de los valores de los atributos.
Para iniciar el clculo de los atributos por ventana, el primer paso consiste en visualizar el sismograma sinttico (original o modificado, segn sea el caso) en la pantalla. Los botones del panel Actualizar Sinttico y Sinttico Original permiten colocar el sismograma en la primera pista, la cual es la diseada para este fin.
Orietta C. Mata P.
86
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 38. Opciones Actualizar Sinttico y Sinttico Original, a partir de las cuales se puede visualizar el sismograma en la interfase.
La opcin Actualizar Sinttico muestra el ltimo sismograma que se visualiz en la pantalla del men principal, mientras que la opcin Sinttico Original muestra el sinttico calculado a partir de los datos originales (sin edicin en tiempo). Eventualmente, si el ltimo sismograma visualizado en la interfase ASPASS / Principal es el sinttico original, entonces al presionar cualquiera de estos dos botones se obtendr el mismo resultado. Luego de visualizar el sismograma sinttico en la interfase de ASPASS / Atributos de Amplitud, el siguiente paso es definir la ventana de anlisis en tiempo. Para ello es necesario presionar el botn Definir Ventana, y posteriormente escoger sobre el sinttico los puntos base y tope en tiempo del intervalo. Obviamente, si queremos calcular los atributos de ventana para el intervalo editado en tiempo, se recomienda escoger los valores inferior y superior de este, sealados por la lnea roja que delimita la ventana. Sin embargo, no existe ninguna restriccin dentro del programa en relacin a este aspecto y se puede escoger cualquier intervalo deseado siguiendo el procedimiento anteriormente descrito.
Orietta C. Mata P.
87
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Base
Tope
Figura 39. Opcin Definir Ventana, con la cual se delimita el intervalo de anlisis.
Inmediatamente despus de seleccionar los valores de la ventana aparecern los valores de los atributos seleccionados para cada pista. Los valores de tope y base de la ventana tambin se mostrarn en los cuadros destinados para este fin en la interfase., a manera de informacin para el usuario. El atributo preestablecido para cada pista es el de amplitud RMS, pero basta con desplegar la lista de atributos para seleccionar un nuevo el cual se calcular automticamente para los valores del intervalo del sismograma que se muestra en la pista 1.
Orietta C. Mata P.
88
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 40. Salida de los atributos calculados en la interfase de ASPASS / Atributos de Amplitud. (En este caso se observa para toadas las pistas el atributo preestablecido: Amplitud RMS)
Figura 41. Lista de atributos de amplitud por ventana disponibles para cada pista. Basta con seleccionar un nombre de la lista para que l atributo sea calculado y mostrado por pantalla.
Orietta C. Mata P.
89
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
De igual manera, los valores numricos de los atributos correspondientes a cada pista son mostrados en la parte inferior derecha de la pantalla, en los campos asignados para este fin. Esto le permite al usuario visualizar mejor los resultados obtenidos.
Figura 42. Valores para cada atributo. (En este caso, como en las tres pistas se calculo el mismo atributo preestablecido, Amplitud RMS, el valor es el mismo)
En la parte superior derecha de la interfase se encuentra una herramienta til que permite modificar la escala del eje de las abcisas para cada pista si el usuario as lo dispusiese. Basta con introducir los lmites del intervalo y seleccionar la pista donde se requiere el cambio.
Orietta C. Mata P.
90
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CAPTULO 6.
ESTADSTICAS.
El mdulo ASPASS / Estadsticas, es el ltimo mdulo desarrollado dentro de ASPASS y es el que permite comparar estadsticamente los resultados obtenidos a lo largo de los otros tres mdulos. Es decir, para que sea posible a utilizacin del estudio estadstico se debe disponer previamente de algn resultado que nos permita establecer una comparacin entre original y modificado.
Orietta C. Mata P.
91
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 44. Diferentes maneras de iniciar el mdulo ASPASS / Estadsticas. (A) Opcin Estadsticas del men principal. (B) Opcin Estadsticas de la pantalla de atributos instantneos y (C) Opcin Estadsiticas del men de la pantalla de atributos de amplitud.
La interfase correspondiente a ASPASS / Estadsticas mantiene el diseo de las pantallas anteriores, presentando cuatro pistas y la serie de botones del panel, los cuales permitirn realizar las operaciones deseadas por el usuario. Como en el caso de la interfase de atributos instantneos y la interfase de atributos por ventana, la primera pista corresponde al campo sismograma sinttico y las otras tres se usan para mostrar la informacin de los atributos calculados.
Orietta C. Mata P.
92
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Orietta C. Mata P.
93
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
atributos instantneos, la traza modificada se encuentra justo encima de la traza original a excepcin de la ventana donde se produjo la modificacin, haciendo mas evidente la variacin producida.
Figura 46. Resultado de comparar el sismograma sinttico y los atributos seleccionados La lnea roja denota la traza original y sus derivados mientras que la lnea azul representa la traza modificada y sus derivados.
En la parte superior de las pistas 2, 3 y 4 se encuentran las listas desplegables con los nombres de los atributos que ASPASS dispone. Seleccionado el nombre del atributo y luego presionando el botn Comparar es posible visualizar en la misma pista el valor del
Orietta C. Mata P.
94
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
atributo para la traza original y para la traza modificada. El atributo preestablecido es amplitud instantnea, sin embargo la lista posee todos los atributos instantneos adems de los de ventana.
Figura 47. Lista de atributos disponibles para realizar la comparacin. Se incluyen todos los atributos instantneos ms los atributos de amplitud
Para aquellos atributos no instantneos, se dispone de un campo en la parte inferior derecha de la pantalla donde se muestran los valores de los atributos para la traza original y para la traza modificada. Esto a manera de herramienta, para que el usuario pueda no solo ver la diferencia grficamente sino que adems cuantitativamente.
Orietta C. Mata P.
95
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Figura 48. Campo de los valores de los atributos por ventana. Rojo denota original y azul modificado
En la esquina superior derecha se encuentra una herramienta que permite modificar la escala del eje x de las pistas. Basta con colocar los lmites y luego seleccionar la pista, para redimensionar la pista. Esta herramienta es particularmente til para observar mejor cuan lejanos o cercanos estn los valores entre ellos.
Figura 49. Herramienta de la interfaz ASPASS / Estadsticas que permite redimensionar el eje X de las pistas.
Orietta C. Mata P.
96
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Al presionar esta opcin y luego la pestaa de Variacin de los Atributos aparecer una lista que contiene todos los atributos disponibles en el paquete. Se pueden seleccionar tantos atributos como el usuario desee y para realizar una seleccin mltiple se debe presionar al mismo tiempo la tecla Control y hacer click con el mouse sobre el nombre del atributo que se quiere comparar. Esta herramienta permite el clculo de un grfico de barras, en el cual se muestran comparativamente los porcentajes de los atributos modificados con respecto a los porcentajes de los atributos originales. Adems se presenta la variacin porcentual entre los valores originales y los modificados, pudiendo establecer claramente aquel (aquellos) atributo (s) que son ms sensibles a los cambios realizados en la edicin de registros.
Orietta C. Mata P.
97
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
La lnea roja que se muestra en el grfico de barras denota el valor de 100%, representando todos los valores originales, valor sobre el cual se calcularon los porcentajes para cada atributo. Las barras azules representan el porcentaje de los valores modificados con respecto a los originales obtenidos por una simple regla de tres y finalmente las barras verdes representan la diferencia calculada al restarle a 100% el porcentaje del modificado, constituyendo una manera fcil y sencilla de cuantificar la variacin producida. Si el valor
Orietta C. Mata P.
98
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
modificado excede al valor original, la resta de 100 menos el valor dar como resultado un nmero negativo, por lo que se tomar en cuenta el valor absoluto de la operacin.
Figura 52. Variacin en los valores de los Atributos. Grfico de barras que muestra la variacin porcentual entre los valores obtenidos para original y modificado.
De esta pequea interfase es posible desplegar una tabla que contiene todos los valores del grfico de barras, permitindole al usuario saber cuantitativamente cuanto han variado los resultados. La tabla contiene los atributos calculados, el valor original, el valor modificado, la diferencia entre estos dos valores, el valor porcentual del modificado y la
Orietta C. Mata P.
99
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
diferencia porcentual entre 100% y el valor del modificado. Vale la pena recordar que para los atributos de ventana se dispone de un punto representativo del intervalo para cada traza y que para los atributos instantneos fue implementado el concepto estadstico de varianza para obtener de igual manera un punto por traza para la ventana. Los valores presentados en azul corresponden a las barras azules de la interfase Variaciones en los valores de los Atributos, mientras que los valores en verde representan las barras de color verde.
Orietta C. Mata P.
100
CAPTULO 6: ESTADSTICAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
CONCLUSIONES. CONCLUSIONES.
El paquete ASPASS fue desarrollado para ser usado por cualquier usuario y con los datos de cualquier pozo, siempre y cuando se encuentren en el formato establecido (*.LAS o *.TXT). Se disearon cuatro mdulos en total, cada uno con una funcin diferente pero a la vez capaz de complementar las salidas de los otros mdulos. El primer mdulo denominado ASPASS / Principal es el responsable de la carga de los datos de pozo, de su organizacin y su posterior presentacin grfica. Adems, es posible realizar el clculo del sismograma sinttico y la edicin de registros en tiempo para generar nuevos sismogramas con caractersticas diferentes en un intervalo dado. Durante el segundo el mdulo se pueden calcular los atributos instantneos de los sismogramas sintticos originales y modificados. Tambin se puede implementar una herramienta til desde el punto de vista cualitativo, la cual permite restar los valores de las trazas original y modificada, visualizando claramente las variaciones producidas por las modificaciones realizadas. Los atributos por ventana, que por ahora solo estn constituidos por aquellos relacionados a amplitud, se pueden obtener haciendo uso del tercer mdulo. En esta interfase se encuentran unos campos que le permiten al usuario visualizar el valor del atributo seleccionado para la ventana de anlisis, la cual constituye una herramienta cuantitativa muy til para cuantificar variaciones. Y finalmente, en el cuarto mdulo se realiza una comparacin estadstica para los resultados obtenidos a lo largo de los otros mdulos. A partir de este, no solo es posible
Orietta C. Mata P.
101
CONCLUSIONES
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
establecer una comparacin sino que adems definir aquel (aquellos) atributo(s) ms sensibles a las modificaciones implementadas en la traza original. Cada uno de los algoritmos diseados realiza una funcin especfica, la cual desempea satisfactoriamente. Los resultados obtenidos utilizando el pozo de prueba son los esperados, lo que confirma la validez del paquete. Luego de cargar los registros de pozo, calcular el sismograma sinttico, editar registros en tiempo y calcular atributos para ambas trazas (original y modificada) es posible establecer, basndonos en la comparacin de las variaciones porcentuales de los valores, los atributos ms sensibles a los cambios realizados. Para probar el programa, los datos originales de pozo fueron sometidos a dos tipos de cambios fundamentales: (1) variacin en espesores y (2) variacin en litologa. Para ello fue necesario contar con un registro que pudiera aportar este tipo de informacin, como lo es gamma ray, el cual se us estrictamente como apoyo para definir el yacimiento de inters. ASPASS respondi satisfactoriamente y arroj resultados que se ajustan a la naturaleza de los cambios. Se usaron tres pozos diferentes para probar los algoritmos de carga y lectura de datos. Para cada uno de los casos, de igual manera los resultados obtenidos fueron exitosos.
Orietta C. Mata P.
102
CONCLUSIONES
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
RECOMENDACIONES. RECOMENDACIONES.
ASPASS permite realizar un modelado unidimensional, debido a que usa una traza sinttica extrada de informacin de pozos para realizar todo el anlisis. Por esta razn un estudio geoestadstico que integrara la informacin de varios pozos a lo largo de un yacimiento no fue posible. En lo sucesivo, ASPASS debera expandir su alcance para que le permita al usuario utilizar mltiples pozos y mediante la implementacin de algunos conceptos claves de geoestadstica, establecer simulaciones y estimaciones de ciertas propiedades fundamentales del yacimiento. Eventualmente se recomienda implementar una metodologa que permita trabajar y manipular datos ssmicos reales, donde se puedan emplear los resultados obtenidos para datos sintticos verificando as su efectividad. Adicionalmente, en el mdulo de atributos por ventana deberan agregarse ms atributos, para que el usuario no solo caracterice el intervalo en base a amplitud sino que adems incorpore una caracterizacin con base en atributos de frecuencia y de fase. La implementacin de toda una metodologa en programacin que incluya interfases no es sencilla. Se necesitan al menos conocimientos bsicos de los comandos y funciones necesarias para establecer las bases del estudio. Conseguir bibliografa que permita adiestrarse un poco sera lo ms imperativo, por lo que sera de mucha utilidad de que al menos la biblioteca de geofsica incorporara dentro de sus libros algunos de programacin en Mat Lab, sobretodo si se piensa seguir la misma lnea de desarrollo de
Orietta C. Mata P.
103
RECOMENDACIONES
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
tesis que se ha tenido hasta el momento, donde el uso de un lenguaje de programacin es necesario y fundamental.
Orietta C. Mata P.
104
RECOMENDACIONES
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
REFERENCIAS.
Fuentes bibliogrficas: Asuaje, Ramn. Implementacin de algoritmos para el clculo de atributos de amplitud instantnea AVO en datos reales. Proyecto de Grado, USB. 2004. Brown, Alistair. Interpretation of Three Dimensional Seismic Data (AAPG Memoir 42). 4th edition. American Asociation of Petroleum Geologist.1996. Brown, A. Seismic attributes and their classification. The Leading Edge, Volumen 15, Nmero 10, p. 1090. 1996. Brown, Alistair. Picking philosophy for 3D stratigraphic interpretation. The Leading Edge. Septiembre, 1998. p. 1198. Chen, Q y S. Sidney. Seismic attribute technology for reservoir forecasting and monitoring. The Leading Edge, Volumen 16, Nmero 5, p. 445. 1997. Cova, Carlos. Implementacin de algoritmos para el clculo y aplicacin de atributos AVO instantneos en sismogramas sintticos. Proyecto de Grado, USB. 2003. Famiglietti, Nathaly. Estudio de atributos ssmicos post apilamient en un volumen de datos 3D del rea de Acema. Proyecto de Grado. USB. 1999. Garca, Eduardo. Curso bsico de perfiles de pozo en hoyo desnudo. Gua de curso de petrofsica en USB. Caracas, 2001.
Orietta C. Mata P.
105
REFERENCIAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Gonzlez, Mayra. Interpretacin estructural del Campo Corozo, norte de Monagas, utilizando ssmica 3D migrada en profundidad. Proyecto de Grado, USB. 2004 Kallweit, R. S y L.C. Word. The Limits of resolution of zero phase wavelets. Geophysics. Vol. 47, 1035 1046.1982. Landa, Andrs. Diseo de un algoritmo para la visualizacin y procesamiento automatizado de registros de pozo. Proyecto de Grado, USB. 2004. Lewis, Catherine. Seismic attributes for reservoir monitoring: A feasibility study using forward modeling. The Leading Edge, Vol. 16. N 5. Mendoza, Jorge. Apuntes del curso de mtodos numricos aplicados a geofsica. Universidad Simn Bolvar, Septiembre Diciembre 2003. Oria, Dorian. Determinacin de atributos ssmicos no convencionales como posibles indicadores de hidrocarburos. Proyecto de Grado, UCV. 2004. zdogan, Yilmaz. Seismic Data Processing. Society of Exploration Geophysicist. Volumen 2. EEUU, 2000. Requena, Neptal. Anlisis de velocidad de ssmica de superficie y de pozos para la migracin en profundidad, utilizando datos 3D de los campos el Furrial y Carito. Proyecto de Grado, UCV.2000. Roscoe, John. Fundamental Research Statistics for the Behavioral Sciences. 2nd edition. Editorial HRW. EEUU, 1975. Sheriff, R. E. Encyclopedic dictionary of exploration geophysics: Society of Exploration Geophysicist, 1991 Spiegel, Murray. Estadstica. Editorial Mc. Graw Hill, Mxico. 1961.
Orietta C. Mata P.
106
REFERENCIAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Taheri, M. y R. Viloria. Caracterizacin Geoestadstica de los Yacimientos. CIED PDVSA. Caracas, 1998. Taner, M. y R. Sheriff. Seismic Stratigraphy Applications to Hydrocarbon Exploration. American Association of Petroleum Geologists Memoir 26, 301 327. EEUU, 1977. Taner, M. Attributes Revisited. Rock Solid Images, Houston, Texas.1992 Valerio, Adriana. Interpretacin Estructural en tiempo y estinacin de anomalas de yacimiento con atributos ssmicos post apilamientos en el Campo Oritpano 600 700. Proyecto de Grado, USB. 2000. p.
@ Fuentes no bibliogrficas:
www.glossary.oilfield.slb.com/ www.spwla.org www.edge-online.com www.rocksolidimages.com www.lgc.com Atencia, Javier y Ral Nestar. Aprenda MatLab 6.0 como si estuviera en primero. Universidad de Navarra. Espaa. Junio, 2001. Barnes, Arthur. Instantaneous frequency and amplitude at the envelope peak of a constant phase wavelet. Geophysics. Vol. 56 N 7, Julio 1991. SEG. pp 1058 1060.
Orietta C. Mata P.
107
REFERENCIAS
Implementacin de algoritmos para estudios de factibilidad en atributos ssmicos utilizando datos sintticos.
Esqueda, Jos. MatLab e Interfases Grficas. Universidad Autnoma de Baja California. Noviembre, 2002. Landmark. PostStack Family Reference Manual. Septiembre, 2001. Robertson, J y Nogami, H. Complex seismic trace analysis of thin beds. Geophysics. SEG. Vol. 49, N 4. Abril, 1984. pp 344 352. Taner, T; Koehler, F y Sheriff, E. Complex seismic trace analysis. Geophysics. SEG. Vol. 44, N 6. 1979. pp. 1041 1063.
Orietta C. Mata P.
108
REFERENCIAS