PREDICCIÓN GENERALIZADA DE MÚLTIPLOS DE SUPERFICIE 3-D
REFERENCIA CRUZADA CON SOLICITUDES RELACIONADAS Esta solicitud reclama el beneficio de la Solicitud provisional de Patente Norteamericana Serie Número 60/560, 129, presentada en Abril 7, 2004, la cual está incorporada a la presente descripción como referencia. Campo de la Invención Las modalidades de la presente invención se refieren generalmente a la investigación sísmica marina y más particularmente, a un método para atenuar el efecto de los múltiplos de superficie en señales sísmicas. Antecedentes de la Invención La investigación sísmica es un método para determinar la estructura de formaciones subterráneas en la tierra. La investigación sísmica generalmente utiliza fuentes de energía sísmica, las cuales generan ondas sísmicas y receptores sísmicos los cuales detectan las ondas sísmicas. Las ondas sísmicas se propagan en las formaciones de la tierra en donde una porción de la sonda se refleja de las interfases entre las formaciones subterráneas. La amplitud de polaridad de las ondas reflejadas es determinada por la diferencia y la ¡mpedancia acústica entre las capas de rocas que comprenden las formaciones subterráneas. La impedancia acústica de una capa de rocas es el producto de la velocidad de la propagación acústica dentro de la capa y la densidad de la capa. Los receptores sísmicos detectan las ondas sísmicas reflejadas y convierten las ondas reflejadas en señales eléctricas representativas. Las señales generalmente son transmitidas por medios eléctricos, ópticos, de radio y otros medios, a aparatos los cuales graban las señales. A través del análisis de las señales grabadas (o trazos) puede ser determinada la forma, posición y composición de las formaciones subterráneas. La investigación sísmica marina es un método para determinar la estructura de formaciones subterráneas subyacentes a los cuerpos de agua. La investigación sísmica marina generalmente utiliza fuentes de energía sísmica y receptores sísmicos localizados en el agua los cuales son, ya sea remolcado detrás de una embarcación o colocados en el fondo del agua desde una embarcación. La fuente de energía generalmente es un aparato explosivo o un sistema de aire comprimido el cual genera energía sísmica, la cual entonces se propaga en la forma de ondas sísmicas a través del cuerpo de agua y dentro de las formaciones de la tierra debajo del fondo del agua Conforme las ondas sísmicas golpean las interfases entre las formaciones subterráneas, una porción de las ondas sísmicas se vuelve a reflejar a través de la tierra y el agua a los receptores sísmicos para ser detectadas, transmitidas y grabadas. Los receptores sísmicos generalmente utilizados en la investigación sísmica marina son sensores de presión, tales como hidrófonos. Adícionalmente, aunque se pueden utilizar sensores de movimiento, tales como los acelerómetros. Tanto las fuentes sísmicas como los receptores pueden ser colocados nuevamente de manera estratégica para cubrir el área de investigación. Sin embargo, las ondas sísmicas se reflejan de interfases diferentes a aquellas que se encuentran entre las formaciones subterráneas, como podría ser deseado. Las ondas sísmicas también se reflejan del fondo del agua y de la superficie del agua, y las ondas reflejadas resultantes mismas continúan reflejándose. Las ondas las cuales reflejan tiempos múltiples reflejados son denominadas "múltiplos". Las ondas que reflejan tiempos múltiples reflejados en la capa de agua entre la superficie del agua arriba y el fondo del agua debajo, son denominadas "múltiplos del fondo del agua". Los múltiplos del fondo del agua han sido reconocidos desde hace tiempo como un problema en el procesamiento sísmico marino y la interpretación, de modo que se han desarrollado múltiples métodos de atenuación basados en la ecuación de la onda para manipular los múltiplos del fondo del agua. Sin embargo, puede ser definido un conjunto más grande de múltiplos que contienen múltiplos del fondo de agua como un subconjunto. El conjunto más grande incluye múltiplos con reflejos ascendente de las interfases entre las formaciones subterráneas además de reflejos ascendente del fondo del agua. Los múltiplos en el conjunto más grande tienen en común sus reflejos descendentes en la superficie del agua y por lo tanto, son denominados "múltiplos de superficie". La figura uno, que se explica más adelante, proporciona ejemplos de los diferentes tipos de reflejos.
La figura 1 muestra una vista en forma de diagrama de la investigación sísmica. El procedimiento es designado generalmente con el número 100. Las formaciones subterráneas que van a ser exploradas tales como 102 y 104, yacen debajo de un cuerpo de agua 106. Las fuentes de energía sísmica 108 y los receptores sísmicos 1 10 son colocados en el cuerpo de agua 106, generalmente por una o más embarcaciones sísmicas (no mostradas). Una fuente sísmica 108, tal como una pistola de aire, crea ondas sísmicas en el cuerpo de agua 106 y una porción de las ondas sísmicas viaja descendente a través del agua hacia las formaciones subterráneas 102 y 104 más allá del cuerpo de agua 106. Cuando las ondas sísmicas alcanzan un reflector sísmico, una porción de las ondas sísmicas se refleja ascendente y una porción de las ondas sísmicas continúa descendente. El reflector sísmico puede ser el fondo del agua 1 12 o una de las interfases entre dos formaciones subterráneas, tales como la interfase 1 14 entre las formaciones 1 12 y 104. Cuando las ondas reflejadas que viajan ascendente alcanzan la interfase de agua/aire en la superficie del agua 1 16, la mayor parte de la porción de las ondas se refleja descendente nuevamente. Continuando de este modo, las ondas sísmicas pueden reflejar tiempos múltiples entre reflectores ascendente, tales como el fondo del agua 1 12 o las interfases de la formación más abajo, y el reflector descendente en la superficie del agua 1 16 arriba, como se describe más delante de una manera más completa. Cada vez que las ondas reflejadas se propagan pasan la posición de un receptor sísmico 1 10, el receptor 1 10 percibe las ondas reflejadas y genera señales representativas. Los reflejos primarios son aquellas ondas sísmicas las cuales se han reflejado solamente una vez del fondo del agua 1 12 y una interfase entre las formaciones subterráneas, antes de ser detectadas por un receptor sísmico 1 10. Un ejemplo de un reflejo primario se muestra en la figura 1 , por medio de las trayectorias del rayo 120 y 122. Los reflejos primarios contienen la información deseada acerca de las formaciones subterráneas, la cual es la meta de la investigación sísmica marina. Los múltiplos de superficie son aquellas ondas que se han reflejado múltiples veces entre las superficies del agua 1 16 y cualesquiera reflectores ascendente, tales como el fondo del agua 1 12 y o las interfases de la formación antes de ser percibidas por un receptor 1 10. Un ejemplo de un múltiplo de superficie el cual es específicamente un múltiplo del fondo del agua, se muestra por medio de las trayectorias del rayo 130, 132, 133 y 136. Al punto en la superficie del agua 1 16 en el cual es reflejada la onda descendente es al que generalmente nos referimos, como un punto de reflejo descendente 133. Los múltiplos de superficie que inician en la trayectoria del rayo 130 son el múltiplo del orden 1 , ya que el múltiplo contiene un reflejo de la superficie del agua 1 16. Dos ejemplos generales de múltiplos de superficie con reflejos ascendente desde el fondo del agua 1 12 y las interfases de la formación se muestran por medio de las trayectorias del rayo 140, 142, 144, 146, 148 y 150 y por las trayectorias del rayo 160, 162, 164, 166, 168 y 170. Ambos de estos últimos dos ejemplos de múltiplos de superficie son múltiplos del orden dos, ya que los múltiplos contienen dos reflejos desde la superficie del agua 1 16. En general un múltiplo de superficie del orden i si el múltiplo contiene i reflejos desde la superficie del agua 1 16. Los múltiplos de superficie son ruidos extraños los cuales obscurecen ia señal deseada de reflejo primario. La atenuación del múltiplo de superficie es una inversión previa al apilado del campo de onda grabado, la cual elimina todos los órdenes de todos los múltiplos de superficie presentes dentro de la señal sísmica marina. A diferencia de algunos algoritmos de atenuación de múltiplos basados en la ecuación de onda, la atenuación múltiple no requiere diseño alguno de o supuestos con respecto a las posiciones, formas y coeficientes de reflejo de los reflectores que causan los múltiplos. En vez de ello, la atenuación del múltiplo de superficie depende de la consistencia física interna entre múltiples eventos primarios que deben existir en cualquier conjunto de datos marinos grabado correctamente. La información necesaria para el proceso de atenuación del múltiplo de superficie, ya está contenido dentro de los datos sísmicos. Varios métodos de la técnica anterior han tratado de eliminar los múltiplos de superficie de los trazos grabados. Por ejemplo, se ha observado que el tiempo de viaje para un múltiplo del fondo del agua es una función de la "compensación", la distancia entre la fuente y el receptor y el número de veces que se reflejan múltiplos desde la superficie. Por ejemplo, si el múltiplo se refleja desde la superficie una vez antes de ser recibido por el micrófono y la compensación es 0, el tiempo de viaje del múltiplo es exactamente dos veces el de las ondas principales. Este hecho ha sido utilizado en varios esquemas para eliminar los múltiplos. Otros métodos comprenden esquemas de trazo de rayos complejos los cuales generan una onda de múltiplos sintéticos y la restan de la onda real para obtener una grabación supuestamente libre de múltiplos. Sin embargo, estos datos son problemáticos, debido a que requieren un conocimiento importante de la estructura submarina, así como la configuración del fondo del océano antes de que pueda ser generada la onda sintética. Se pueden generar múltiplos sintéticos similares utilizando métodos más exactos que no comprenden directamente el trazado del rayo, por ejemplo, técnicas de propagación en el campo, pero nuevamente estos métodos requieren un conocimiento detallado de por lo menos el fondo del océano, así como la forma de las interfases del fondo del mar y por lo tanto, no son tan prácticos como se habría deseado. Los algoritmos actuales de predicción de múltiplos de superficie requieren circunvuluciones de pares de trazos para las cuales la localización del receptor para un trazo del par coincide con la localizacíón del disparo para el otro trazo. Debido a que las localizaciones de la fuente y el receptor en los datos registrados, rara vez coinciden de manera precisa, es una práctica común regularizar los conjuntos de datos a una geometría nominal, de modo que sea lograda esta coincidencia de las ubicaciones de la fuente y el receptor. Los múltiplos entonces son pronosticados por medio de esta geometría regular y luego des-regularizados a la geometría original antes de la resta. Desafortunadamente, la regularizacíón y (especialmente) los procesos de des-regularización, con frecuencia son inexactos, los cuales con frecuencia conducen a errores importantes en los múltiplos pronosticados. Por lo tanto, existe una necesidad en la técnica de un método mejorado para eliminar la grabación de los eventos de reflejos de múltiples superficies de las grabaciones sísmicas para propósitos de procesamiento de datos sísmicos. Breve Descripción de la Invención Una o más modalidades de la presente invención tienen como objetivo minimizar los errores ocasionados por los procesos de regularización y des-regularización minimizando el involucramiento de los procesos de regularización y evitando completamente el proceso de des-regularización. En una modalidad, la presente invención funciona con datos de muchas de las subsuperficies a la vez para minimizar los procesos de regularización. Como tal, la presente invención puede comprender circunvoluciones de trazos de líneas diferentes de subsuperficie. En otra modalidad, con el objeto de mantener la eficiencia asociada con el procesamiento basado en la línea de la subsuperficie, el algoritmo de la presente invención está formulado para operar en pares de líneas de subsuperficie en un momento. Todavía en otra modalidad, la presente invención se refiere a la implementación generalizada de un algoritmo de predicción de múltiplos relacionados con la superficie tridimensional (SMP), el cual no hace supuestos acerca de la regularidad o distribución de los trazos en los conjuntos de datos grabados, es decir, no existe el concepto de una geometría nominal. Como tales, los conjuntos de datos grabados son simplemente tratados como una colección de trazos definidos por las ubicaciones de su fuente y de su receptor. De esta manera, una o más modalidades de la presente invención están configuradas para pronosticar múltiplos para cualquier conjunto de trazos definidos de manera similar por sus localizaciones de la fuente y receptor. En particular, la presente invención puede pronosticar múltiplos en la localización correcta, la compensación y el azimuth. También pueden ser diseñadas una o más modalidades de la presente invención que toman en cuenta todas las irregularidades en la geometría de adquisición, las cuales pueden ser particularmente benéficas para resolver la pluma del cable. La exactitud de los múltiplos pronosticados puede depender de la distribución de los trazos en el conjunto de datos grabado. Por consiguiente, una o más modalidades de la presente invención se refieren a un método para pronosticar una pluralidad de múltiplos de superficie para una pluralidad de datos objetivo en una grabación de datos sísmicos. En una modalidad, el método incluye crear un archivo que contiene información con respecto a la pluralidad de pares de trazos grabados. Cada par de trazos grabados es substancialmente más cercano a un trazo deseado del lado del disparo y a un trazo deseado del lado del receptor. El método incluye además la circunvolución de pares de trazos grabados para generar una pluralidad de circunvoluciones y apilar las circunvoluciones de cada trazo objetivo. En otra modalidad, el método incluye (a) seleccionar un trazo objetivo; (b) seleccionar un punto de reflejo descendente potencial para el trazo objetivo seleccionado; (c) calcular por lo menos uno de un punto medio deseado del lado del disparo, la compensación, el azimuth, y por lo menos uno de un punto medio deseado del lado del receptor, la compensación y el azimuth utilizando el punto de reflejo ascendente potencial seleccionado y el trazo objetivo seleccionado. El punto medio deseado del lado del disparo, la compensación y el azimuth definen un trazo deseado del lado del disparo. El punto medio deseado del lado del receptor, la compensación y el azimuth definen un trazo deseado del lado del receptor. El método incluye además (d) determinar un par de trazos grabados substancialmente más cercanos al trazo deseado del lado del disparo y el trazo deseado del lado del receptor; y (e) circunvolucionar el par de trazos grabados para generar una circunvolución. Breve Descripción de los Dibujos De modo que la manera en la cual pueden ser entendidas de manera detallada las características anteriormente mencionadas de la presente invención, una descripción más particular de la presente invención, brevemente resumida anteriormente, pueden tener como referencia las modalidades, algunas de las cuales están ilustradas en los dibujos adjuntos. Sin embargo, deberá observarse que los dibujos adjuntos ilustran solamente modalidades físicas de la presente invención y por lo tanto, no deberán considerarse como limitantes de su alcance, ya que la invención puede admitir otras modalidades igualmente efectivas. La figura 1 ilustra una vista diagramática de una investigación sísmica marina. La figura 2 ilustra un diagrama de flujo de la primera etapa en un método para realizar una predicción tridimensional de múltiplos de superficie de acuerdo con una o más modalidades de la presente invención . Las figuras 3A y 3B ilustran un diagrama de flujo de la segunda etapa en el método para realizar la predicción tridimensional de múltiplos de superficie de acuerdo con una o más modalidades de la presente invención. La figura 4 ¡lustra una vista de planta de una geometría de adquisición de acuerdo con una o más modalidades de la presente invención. La figura 5 ilustra una red de cómputo dentro de la cual pueden ser implementadas varias modalidades de la presente invención. Descripción Detallada de la Invención La figura 2 ¡lustra un diagrama de flujo de la primera etapa 200 en un método para realizar la predicción tridimensional de múltiplos de superficie de acuerdo con una o más modalidades de la presente invención. En el paso 210, es seleccionado el trazo objetivo. Un ejemplo de un trazo objetivo seleccionado se ilustra en la figura 4 como el trazo (S, R) en los trazos objetivo definen las ubicaciones en las cuales deberán ser pronosticados los múltiplos. En el paso 220, se determina o define una apertura 305 para el trazo objetivo seleccionado. La apertura 305 puede ser un área rectangular y centrada en la localización del punto medio M del trazo objetivo. Otras formas geométricas para la apertura 305 están contempladas por las modalidades de la presente invención. La apertura 305 es definida como que incluye substancialmente todos los puntos de reflejo descendente potenciales (DRPs) de los múltiplos de superficie para el trazo objetivo. Como un ejemplo, un punto de reflejo descendente potencial X se ilustra en la figura 4. En el paso 230, la apertura 305 está colocada dentro de una pluralidad de celdas. En una modalidad, el punto medio del trazo objetivo está localizado en uno de los nodos de rejilla (centros de celdas). La separación de las rejillas puede ser arbitraria. Los nodos de rejilla definen los DRPs potenciales para el trazo objetivo. En el paso 240, se selecciona un DRP potencial, tal como un primer DRP para el trazo objetivo seleccionado. En el paso 250, son calculados el punto medio deseado del lado del disparo Ms, la compensación Xs y el azimuth ?s y el punto medio deseado del lado receptor MR, la compensación XR y el azimuth TR. El Ms es la localizacíón del punto medio entre la fuente y el DRP potencial seleccionado. La compensación Xs es la distancia horizontal entre el DRP potencial seleccionado y la fuente S. El azimuth TR es definido como el ángulo entre la línea que conecta la fuente S y el DRP potencial seleccionado y alguna dirección fija, la cual es generalmente la dirección en línea. El MR es la ubicación del punto medio entre el receptor R y el DRP potencial seleccionado. La compensación XR es la distancia horizontal entre el DRP potencial seleccionado y el receptor R. El azimuth TR es definido como el ángulo entre la línea que conecta el receptor R y el DRP potencial seleccionado y alguna dirección fija, la cual es generalmente la dirección en línea. En una modalidad, el punto medio deseado del lado del disparo Ms, la compensación Xs y el azimuth ?s y el punto medio deseado del lado del receptor MR, la compensación XR y el azimuth TR son calculados basados en el trazo objetivo seleccionado y el DRP potencial seleccionado. Los puntos medios, compensaciones y azimuths juntos definen el trazo deseado del lado del disparo (S,X) y el trazo deseado del lado del receptor (X R). En el paso 260, son determinados el trazo de entrada más cercano al trazo deseado del lado del disparo y al trazo de entrada más cercano al trazo deseado del lado del receptor. En una modalidad, se determinan las entradas de los trazos más cercanos minimizando una función objetivo, la cual define la cercanía de dos trazos basada en sus puntos medios, compensaciones y azimuths. Un ejemplo de una función de objetivo es: D2 = |?m|2 + wx|?x|2 + w? |?T|2 en donde D mide la cercanía entre los trazos, ?m, ?x y ?? son las diferencias en el punto medio, la compensación y el azimuth, respectivamente, y wx y w? son los pesos que definen la importancia relativa de los errores en las compensaciones y azimuths comparadas con el error en los puntos medios. De manera notable, wx no tiene dimensión, mientras que w? tiene dimensiones de L2. En una modalidad w? es ajustado en cero debido a una cobertura deficiente del azimuth en el conjunto de datos de entrada. En otra modalidad, puede existir un valor mínimo para la función objetivo minimizada arriba de la cual se considera que no existen trazos que coincidan. En el paso 270, la información con respecto a los trazos de entrada más cercanos es almacenada en un archivo, al cual nos podemos referir como un archivo de índice de circunvolución (CIF). Por ejemplo, dicha información puede incluir identificadores para los trazos de entrada más cercanos, sus líneas de subsuperficie asociadas, el punto de reflejo descendente potencial seleccionado X, el punto medio deseado del lado del disparo Ms, la compensación Xs y el azimuth ?s, el punto medio deseado del lado del receptor MR, la compensación XR y el azimuth ?R y se podrá seleccionar el trazo objetivo seleccionado. En el paso 280, se toma una determinación con respecto a si la apertura incluye otro DRP potencial para el trazo objetivo seleccionado. Si la respuesta es en el sentido afirmativo, entonces el procesamiento regresa al paso 240 en el cual es seleccionado otro DRP potencial. Si la respuesta es negativa, entonces el procesamiento continúa al paso 285, en el cual se toma una determinación con respecto a si existe otro trazo objetivo. Si la respuesta es afirmativa, entonces el procesamiento regresa al paso 210, en el cual es seleccionado otro trazo objetivo. Si la respuesta es negativa, entonces el procesamiento continúa al paso 290, en el cual el CIF es dividido en uno o más subarchivos de acuerdo con los pares de líneas de subsuperficíe que contienen los trazos de entrada más cercanos. De esta manera, cada subarchivo contiene información dirigida solamente a un par de líneas de subsuperficie, en donde cada línea de subsuperficie contiene un trazo de entrada más cercano a cualquiera del trazo deseado del lado del disparo o el trazo deseado del lado del receptor. El orden de las líneas de subsuperficie en el par no es crítico. Las figuras 3A y 3B ilustran un diagrama de flujo de la segunda etapa 300 en un método para realizar una predicción tridimensional de múltiplos de superficie de acuerdo con una o más modalidades de la presente invención. En el paso 310, se selecciona el primer subarchivo. En el paso 320, la información con respecto a un par de trazos de entrada más cercanas al trazo deseado del lado del disparo y al trazo deseado del lado del receptor para un trazo objetivo seleccionado se lee del subarchivo seleccionado. En el paso 330, un par de trazos de entrada correspondientes a la información con respecto al par de trazos de entrada más cercanos es extraída de un conjunto de datos sísmicos grabado. El conjunto de datos sísmicos grabado puede ser almacenado en cualquier archivo o almacenamiento de datos generalmente conocidos por las personas expertas en la técnica. El conjunto de datos sísmicos grabado puede ser extrapolado a la compensación cero. El conjunto de datos sísmicos grabado puede ser una recolección de trazos previos al apilado definidos por el punto medio, la compensación y el azimuth. Cada trazo y el conjunto de datos sísmicos grabados pueden tener un identificador de línea de subsuperficie y un identificador único del trazo que puede ser utilizado para identificar el trazo de entrada dentro del conjunto de datos sísmicos grabados. El conjunto de datos sísmicos grabados puede ser organizado en líneas de subsuperficie o cualesquiera de otras subdivisiones, tales como líneas de navegación. En el paso 340, una corrección del diferencial del movimiento de salida es aplicada al par de trazos registrados extraídos para corregir las compensaciones de los trazos registrados extraídos para la compensación del lado del disparo deseado, en la compensación del lado del receptor deseado. En el paso 350, el par de trazos registrados corregidos y extraídos son círcunvolucionados. En el paso 355, la circunvolución es almacenada. En el paso 360, se toma una determinación con respecto a si el subarchivo seleccionado contiene otro par de trazos de entrada que van a ser circunvolucíonados. Si ia respuesta es afirmativa, entonces el procesamiento regresa al paso 320. Si la respuesta es negativa, entonces el procesamiento continúa al paso 365, en el cual son clasificadas las circunvoluciones de acuerdo con los trazos objetivo. En el paso 370, todas las circunvoluciones por cada tres objetivos son apiladas juntas para tener una sola circunvolución apilada por trazo objetivo para el subarchívo seleccionado. En el paso 375, se toma una determinación con respecto a si existe otro subarchivo de CI F. Si la respuesta es afirmativa, entonces es seleccionado su archivo (paso 378) y el procesamiento regresa al paso 320. Si la respuesta es negativa, entonces el procesamiento continúa al paso 380 en el cual todas las circunvoluciones apiladas por cada subarchivo en el CIF son clasificadas de acuerdo con su trazo objetivo. En el paso 385, todas las circunvoluciones apiladas por cada subarchivo en el CIF son apiladas por cada trazo objetivo para obtener una sola circunvolución apilada por trazo objetivo para todos los subarchivos. En el paso 390, la firma de la fuente es descircunvolucíonada de acuerdo con técnicas generalmente conocidas por aquellos expertos en el arte. En el paso 395, puede ser aplicado un filtro-p tridimensional para corregir el efecto del apilado en la onda de acuerdo con las técnicas generalmente conocidas por los expertos en el arte. La figura 5 ilustra una red de cómputo 500 en la cual pueden ser implementadas las modalidades de la presente invención. La red de cómputo 500 incluye un sistema de computadora 530, el cual puede ser implementado como cualquier computadora personal convencional o estación de trabajo, tal como la estación de trabajo basada en UNIX. El sistema de computo 530 está en comunicación con los dispositivos de aimacenamíento de disco 529, 531 y 533 los cuales pueden ser externos a los dispositivos de almacenamiento de disco duro. Está contemplado que los dispositivos de almacenamiento de disco 529, 531 y 533 son unidades de disco duro convencionales, y como tales, serán implementadas por medio de una red de área local o por acceso remoto. Por supuesto, aunque los dispositivos de almacenamiento de disco 529, 531 y 533 están ¡lustrados como dispositivos separados, puede ser utilizado un solo dispositivo de almacenamiento de disco para almacenar cualquiera y todas las instrucciones del programa, datos de medición y resultados que se desean. En una modalidad, los datos sísmicos de los hidrófonos son almacenados en el dispositivo de almacenamiento de disco 531. El sistema de cómputo 530 puede recuperar los datos apropiados del dispositivo de almacenamiento de disco 531 para realizar la predicción de múltiplos de superficie 3-D de acuerdo con instrucciones de programa que corresponden a los métodos aquí descritos. Las instrucciones de programa pueden ser escritas en un lenguaje de programación de computación, tal como C++, Java y similares. Las instrucciones del programa pueden ser almacenadas en una memoria legible por computadora, tal como un dispositivo de almacenamiento de disco del programa 533. Por supuesto, el medio de memoria que almacena las instrucciones del programa puede ser cualquier tipo convencional o utilizado para el almacenamiento de programas de cómputo, incluyendo, unidades de disco duro, diskettes, CD-ROMs y otros medios ópticos, cintas magnéticas y similares. De acuerdo con la modalidad preferida de la presente invención, el sistema de cómputo 530 presenta la salida principalmente en un despliegue de gráficos 527 o alternativamente, por medio de la impresora 528. El sistema de computo 530 puede almacenar los resultados de los métodos descritos anteriormente en el almacenamiento de disco 529 para su uso posterior y el análisis adicional. El teclado 526 y el dispositivo de señalamiento (por ejemplo, un ratón, bola de rastreo o similar) 525 puede ser provisto con el sistema de cómputo 530 para hacer posible la operación interactiva. El sistema de cómputo 530 puede estar localizado en un centro de datos remoto a la región de la investigación. La computadora del sistema 530 está en comunicación con los hidrófonos (ya sea directamente o por medio de la unidad de grabación, no mostrada) para recibir señales indicadoras de la energía sísmica reflejada. Estas señales, después del formateo convencional y otro procesamiento inicial, son almacenadas por la computadora del sistema 530, como datos digitales en el almacenamiento de disco 531 para la recuperación posterior y procesamiento de la manera descrita anteriormente. Aunque la figura 5 ¡lustra el almacenamiento de disco 531 como conectado directamente a la computadora del sistema 530, también está contemplado que el dispositivo de almacenamiento de disco 531 puede ser accesible a través de la red de área local o por acceso remoto. Además, aunque los dispositivos de almacenamiento de disco 529, 531 están ilustrados como dispositivos separados para almacenar los datos sísmicos de entrada y los resultados del análisis, los dispositivos de almacenamiento de disco 529, 531 pueden ser implementados dentro de una sola unidad de disco (junto con o separados del dispositivo de almacenamiento de disco del programa 533) o de cualquier otra manera convencional como seria perfectamente entendido por los expertos en la técnica refiriéndose a esta especificación . Aunque lo anterior está relacionado con las modalidades de la presente invención, se pueden prever otras modalidades adicionales de la presente invención, sin salirse del alcance básico de la misma, y el alcance de la misma es determinado por las reivindicaciones siguientes.