MX2014014770A - Analisis de datos sismicos utilizando coleccion de datos de nodos del fondo oceanico. - Google Patents

Abstract

La presente invención permite que se reduzca un error de tiempo de recorrido de RMS en una adquisición de datos sísmicos. Se recolectan primero las mediciones de campo de las coordenadas de fuente y receptor, velocidad del sonido en el agua como una función de la profundidad y el tiempo, sincronización del receptor y desfase del reloj. Se examinan entonces los datos sísmicos para medir el tiempo de recorrido de cada fuente hacia cada receptor. Puede entonces calcularse un tiempo de recorrido modelo en base a las mediciones de campo. Al perturbar de manera iterativa al menos uno de los datos medidos en campo usando una tabla de búsqueda y calcular el tiempo de recorrido después de cada perturbación hasta que se ha logrado un error aceptable de RMS, los datos sísmicos condicionados que toman en cuenta la naturaleza dinámica de la columna de agua proporcionarán la base para crear un mapa sísmico exacto que no se afecte por las condiciones cambiantes del agua.

Description

ANÁLISIS DE DATOS SÍSMICOS UTILIZANDO LA RECOLECCIÓN DE DATOS DE NODO DEL FONDO OCEANICO REFERENCIA CRUZADA CON LA SOLICITUD RELACIONADA La presente solicitud reivindica la prioridad de la Solicitud de Patente Provisional de los Estados Unidos Serie Número 61/655,325, presentada el 4 de junio de 2012, la totalidad de la cual se incorpora en la presente mediante la referencia.

ANTECEDENTES La presente invención se refiere en general a la adquisición y análisis de datos sísmicos a través de la conducción de un estudio sísmico y más particularmente a metodologías para mejorar la exactitud de los resultados en base al análisis de datos sísmicos.

La generación y registro de datos sísmicos involucra muchas diferentes configuraciones del receptor, incluyendo la colocación de geófonos o sismómetros en la superficie de la Tierra o fondo del mar, hidrófonos de remolque detrás de un recipiente sísmico marino, suspensión vertical de hidrófonos en el mar o colocación de geófonos en un sondeo (como en un perfil sísmico vertical) para registrar la señal sísmica. Una fuente, tal como una unidad vibradora, disparo de dinamita o una pistola de aire, genera vibraciones acústicas o elásticas que viajan hacia la Tierra, pasan a través de estratos con diferentes respuestas sísmicas y efectos de filtración, y regresan a la superficie para registrarse como datos sísmicos. La adquisición óptima varía de acuerdo con condiciones locales e involucra el empleo de la fuente adecuada (tanto el tipo como la intensidad), la configuración óptima de receptores, y la orientación de las líneas del receptor con respecto a las características geológicas. Esto asegura que pueda registrarse la máxima proporción de señal a ruido, que la resolución sea adecuada y que puedan reducirse o distinguirse efectos extraños tales como ondas de aire, ondas superficiales, múltiplos y difracciones, y retirarse a través del procesamiento.

La exploración sísmica submarina se utiliza ampliamente para localizar y/o analizar formaciones geológicas subterráneas respecto a la presencia de depósitos de hidrocarburos. Un tipo de estudio utiliza un bote que remolca una pluralidad de pistolas de aire y un conjunto de nodos del fondo oceánico (OBN) colocados en el fondo del océano. Los nodos se colocan en el fondo del océano por medio de un vehículo operado de manera remota (ROV) y cargador submarino; típicamente, el ROV y el cargador submarino se despliegan desde una lancha de equipo/ recuperación separada del remolque o pontón.

Para adquirir los datos, se liberan disparos de aire comprimido desde las pistolas de aire a intervalos periódicos conocidos y se registra la ubicación y sincronización de cada disparo. De igual modo, se registra la sincronización e intensidad de los nodos que detectan cada uno de los disparos de aire comprimido. Típicamente, los datos se recolectan durante al menos 30 días a fin de cubrir un ciclo completo de la marea, pero pueden recolectarse durante un mayor periodo de tiempo dependiendo de, por ejemplo, el tamaño del área por estudiarse.

La recolección de datos sísmicos de cuatro dimensiones generalmente comprende datos sísmicos tri dimensionales (3D) adquiridos en diferentes momentos sobre la misma área a fin de determinar cambios en la producción de un depósito de hidrocarburos a través del tiempo. Los cambios pueden observarse en la ubicación y saturación, presión y temperatura de un fluido. Los datos sísmicos en 4D son una de varias formas de datos sísmicos en un lapso de tiempo. Tales datos pueden adquirirse sobre la superficie o en un pozo de sondeo. Los datos sísmicos en un lapso de tiempo involucran la recolección de datos sísmicos a partir de la superficie o un pozo de sondeo, adquiridos en momentos diferentes, sobre la misma área, a fin de determinar cambios en la sub-superficie a través del tiempo, tal como el movimiento de fluido o efectos de recuperación secundaria. Los datos se examinan respecto a cambios en atributos relacionados a expresiones del contenido de fluido. Los datos sísmicos en un lapso de tiempo pueden repetir datos en 2D, 3D (lo cual se conoce como datos sísmicos en 4D)), datos Crosswell y VSP (perfil sísmico vertical).

El uso incrementado de adquisición OBN para monitoreo en un lapso de tiempo en el aguas profundas ha demostrado una necesidad de tener posiciones del receptor y de disparo muy exactas, así como también una buena comprensión de las variaciones de velocidad y altura de la columna de agua de frecuencia elevada, durante la adquisición, que afectan la sincronización del evento sísmico. Como se mencionó, los nodos se encuentran típicamente en el fondo del mar por más de treinta días y se observa así al menos un ciclo completo de las mareas y posiblemente grandes cambios en la velocidad de la columna de agua. Los ciclos de las mareas son problemáticos debido a que afectan la distancia vertical entre los nodos y disparos durante la adquisición. Deben comprenderse las variaciones de velocidad de frecuencia elevada con objeto de calcular las correcciones estadísticas de la columna de agua a fin de corregir los datos en una sola velocidad dependiente de la profundidad del agua para el procesamiento corriente abajo y la formación de imágenes de los datos OBN.

Los nodos del fondo oceánico tienen todos relojes internos que son independientes al reloj principal (típicamente hora GPS), lo cual sirve como la referencia para todas las observaciones en un estudio dado. Cada reloj interno experimenta cierta cantidad de derivación del reloj, debido a factores tales como envejecimiento del oscilador de cristal, error humano en la calibración, y simplemente activación retrasada en comparación con el reloj GPS principal, durante el despliegue del nodo. Para cada nodo, el valor de derivación del reloj se mide cuando el nodo se sincroniza con el reloj principal tras la recuperación del nodo a partir del fondo del mar. Sin embargo, existen casos donde esta medición no puede realizarse o se realiza de manera no confiable, lo cual conduce a incertidumbre en la medición de derivación del reloj. Para determinar derivación del reloj en nodos con mediciones de campo no confiables, se comparan los rastros del disparo espacialmente circundante con un tiempo de adquisición significativamente diferente y se selecciona una derivación del reloj que retira discontinuidades en el arribo directo. Sin embargo, si existe un cambio correspondiente en la velocidad del agua, debido a cambios en el tiempo de adquisición, entonces un análisis de discontinuidad puede engañar o conducir a ajustes de derivación incorrectos.

Un proceso que corrige solo uno de colocación de disparo, colocación de receptor, sincronización del receptor o velocidad de columna de agua, será sujeto a fuga de errores en las otras dimensiones. Debe utilizarse una téenica que resuelva simultáneamente la posición de la fuente, la posición del receptor, la derivación del reloj receptor y la velocidad del agua como una función del tiempo de adquisición.

Es un objeto primario y ventaja de la presente invención el proporcionar un método, sistema y producto que resuelva el modelo de posiciones de disparo y de receptor, sincronización del receptor y la velocidad de la columna de agua, de manera simultánea.

Otros objetos y ventajas de la presente invención serán en parte aparentes a aquellos expertos en la materia y en parte lo serán en lo sucesivo.

BREVE DESCRIPCIÓN DE LA INVENCIÓN El tiempo que lleva tomar energía para recorrer desde la fuente hasta el nodo se determina por la posición relativa del disparo (I) y el receptor (J) y la velocidad del sonido, o lentitud s, a través del agua (Ecuación 1) (Mitchell et al . ) . El tiempo de arribo directo, registrado por el nodo, es el tiempo de recorrido de la energía más el error de desfase del reloj al momento del disparo (Ecuación 2). Un proceso de inversión de tiempo de recorrido reduce el error entre los tiempos de arribo directo, recolectados, reales, y los tiempos de arribo directo, calculados por el modelo (por la Ecuación 1 la Ecuación 2) para producir una solución a las posiciones de disparo y de receptor, derivación del reloj receptor y velocidad de la columna de agua que mejor se ajusta a los datos.

Para acondicionar los datos a la recolección de arribo directo, se lleva a cabo una designación angular que corrige el arribo directo a la fase cero. A mayores desplazamientos de linea transversal, la separación lateral de cada secuencia fuente provoca que el arribo directo parezca un triplete en lugar de un solo tren de ondas. Los diseños direccionales estabilizan la recolección de los tiempos de arribo directos a mayores desplazamientos de la linea transversal, mediante colapso del triplete en un solo tren de ondas. La eliminación de espectros secundarios dependientes del ángulo del arribo directo permite que la recolección de realice en el pico de un tren de ondas de fase cero. Típicamente, el tiempo de arribo directo se elige como el menor de 3500m o tres veces la profundidad del receptor, con objeto de evitar la recolección de tiempos de arribo directos cuando la energía reflejada interfiere con el arribo directo. La des-marcación del arribo directo mueve eficazmente la fuente hacia la superficie marítima, elevación de O. Estas elevaciones se establecen entonces como la altura de la marea modelo, con relación al nivel del mar promedio, al momento del disparo. Las profundidades de receptor iniciales también se corrigen respecto a la altura de la marea en base al tiempo de despliegue de nodo, cuando la medición de profundidad inicial se realizó por el ROV (Vehículo Operado de Manera Remota).

El proceso de inversión de tiempo de recorrido requiere de una solución inicial para las posiciones del receptor y disparo, la sincronización del receptor y la velocidad de la columna de agua. Las posiciones y sincronización se obtienen generalmente mediante el uso de métodos no iterativos que utilizan una suposición de velocidad de agua constante. Las soluciones iniciales de sincronización y posición pueden ser una solución los "suficientemente aproximada" para comenzar el proceso de inversión de tiempo de recorrido. El modelo de velocidad inicial se deriva de las mediciones de presión, conductividad y temperatura registradas durante los tránsitos de columna de agua del equipo de adquisición mientras se despliegan y recuperan los nodos.

Cada rastro, dentro del rango de desplazamiento seleccionado, se correlaciona con el desplazamiento más próximo del nodo dado. Cuando el coeficiente de correlación es inaceptablemente bajo, la recolección puntual se rechaza del proceso de inversión.

Por consiguiente, en un aspecto de la invención, se proporciona un método para reducir errores de RMS en datos recolectados durante un estudio sísmico sobre un área predeterminada, comprendiendo las etapas de a) obtener mediciones de campo de i) posiciones x, y, y z de fuente, ii) posiciones x, y, y z de receptor, iii) velocidad del sonido en el agua como una función de tiempo y profundidad, iv) sincronización del receptor y v) desfase del reloj; b) medir el tiempo de recorrido de cada fuente hacia cada receptor; c) determinar un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y d) perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

En otro aspecto de la invención, un programa de computadora almacenado en un medio de almacenamiento no transitorio y adaptado para ejecutarse en un procesador de computadora con propósitos de reducir el error de RMS en datos recolectados durante un estudio sísmico sobre un área predeterminada, comprendiendo un código de programa para obtener mediciones de campo de las posiciones de fuente, posiciones del receptor, velocidad del sonido en el agua como una función de tiempo y profundidad, sincronización del receptor y desfase del reloj; medir el tiempo de recorrido de cada fuente hacia cada receptor; determinar un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

En otro aspecto de la presente invención, se proporciona un sistema de computadora que incluye un programa de computadora almacenado en un medio de almacenamiento no transitorio y adaptado para ejecutarse en un procesador, reduciendo el programa de computadora, de manera iterativa, el error de RMS en datos recolectados durante un estudio sísmico sobre un área predeterminada, comprendiendo un código de computadora para: obtener mediciones de campo de las posiciones de fuente, posiciones del receptor, velocidad del sonido en el agua como una función de tiempo y profundidad, sincronización del receptor y desfase del reloj; medición del tiempo de recorrido de cada fuente hacia cada receptor; determinación de un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

BREVE DESCRIPCIÓN DE LOS DIBUJOS La presente invención se entenderá de manera más completo y se apreciará mediante lectura de la siguiente Descripción Detallada en conjunto con los dibujos acompañantes, en los cuales: La Figura 1 es una vista esquemática de un sistema de estudio sísmico; La Figura 2 es un esquema que muestra la geometría en un sistema de estudio sísmico que se basa en la presente invención; Las Figuras 3a-3c son diagramas estructurales de tablas de búsqueda usadas con la presente invención; La Figura 4 es un diagrama de flujo de alto nivel; Las Figuras 5a-5b son diagramas de flujo que ilustran un aspecto de la presente invención; La Figura 6 es una serie de gráficas que representan una vista en mapa de un resultado disparo a disparo para un receptor dado; y Las Figuras 7a-7d son vistas gráficas de un mapa de receptor/nodo que muestra la diferencia de RMS.

DESCRIPCIÓN DETALLADA Haciendo ahora referencia a los dibujos, en donde los números de referencia similares se refieren a partes similares a través de los mismos, se observa en la Figura 1 una vista ilustrativa del paradigma usado para recolectar datos sísmicos usando nodos del fondo oceánico (OBNs). En la vista de la Figura 1, una embarcación 10, jala a un conjunto de fuentes de energía (por ejemplo, pistolas de aire comprimido) 12 en una forma serpentina a través de una rejilla 14 de área predeterminada sobre una estructura de tiempo predeterminada (por ejemplo, 60-90 dias). Los OBNs 16 se colocan en el fondo del océano 18 por medio de un cargador 20 que los transporta desde una embarcación 17 (podría ser también la embarcación 10) hasta el fondo del océano, y un ROV 22 que los lleva desde el cargador 20 y los coloca en una posición predeterminada en el fondo del océano. Una vez que se ha colocado un cierto número de los OBNs 16, la embarcación 10 inicia su movimiento en serpentina alrededor de la rejilla 14. Las fuentes de energía 12 expulsan una fuente de energía hacia el fondo del océano 18 a intervalos predeterminados, y la energía se recibe por los OBNs 16 que se encuentran en la trayectoria de la fuente de energía. Una vez que se han colocado todos los OBNs 16 (también referidos como "nodos" o "receptores") en el fondo del océano, el cargador 20 y el ROV 22 regresan a la superficie donde pueden recuperarse los datos de medición de la velocidad de la columna de agua y ciertos datos de posición del receptor. En la Figura 2, se ilustra la geometría de la fuente/disparos 12 y el nodo/receptor/OBN 16.

Es bien sabido que cada uno de los OBNs tiene un reloj atómico así como también los sensores de energía (por ejemplo, hidrófonos y geófonos) para detectar la energía en las ondas de recorrido descendente (energía directa), así como también las ondas de recorrido ascendente (energía reflejada). Después del despliegue y recuperación hacia la embarcación, cada reloj atómico se sincroniza con un reloj principal en base a GPS. La sincronización requiere de la determinación de la frecuencia de inicio del reloj atómico de tal manera que la sincronización concuerde con la del reloj principal en base a GPS. Según se conoce en la materia, los relojes atómicos experimentarán cierta derivación mesurable; la derivación se debe (1) a un desplazamiento lineal provocado por falla al determinar perfectamente la frecuencia de inicio durante la calibración, (2) derivación cuadrática (o derivación que se representa por una relación cuadrática) debido al decaimiento del cristal del reloj atómico a través del tiempo, y (3) una forma de derivación volumétrica que se provoca simplemente por error humano por tiempos de inicio inferiores a la perfección. En general, el desfase del reloj puede medirse en el campo cuando se recupera cada nodo al final de la recolección de datos y el tiempo transcurrido para cada reloj atómico se compara a través del tiempo transcurrido en el reloj principal de GPS.

Otros puntos de datos medidos en campo incluyen las coordenadas para cada OBN (receptor) 16. Las posiciones La x y la y para los OBNs 16 se determinan usando herramientas de hallazgo de rango utilizadas en la superficie, pero se sujetan a errores relacionados a la distancia siendo mayores los errores que se producen a mayor profundidad del agua. La posición Z de cada OBN 16 se determina en base a una lectura de presión al momento de la colocación. Las coordenadas de cada fuente de energía cada vez que se libera un "disparo" (de energía) se registran también usando teenología GPS, pero se sujetan a ligeras variaciones debido a ondulaciones provocadas por las ondas superficiales en el océano. La velocidad de la columna de agua o la velocidad (o la lentitud) del sonido en el agua como una función de tiempo y profundidad, es de naturaleza dinámica debido a las mareas, corrientes, cambios de salinidad del agua, cambios de temperatura del agua y otros factores conocidos, pero sin importar que se mida en campo al momento de despliegue del OBN mediante el uso de tecnología conocida, presente en el cargador 20 , esos factores en presión, conductividad y temperatura y esta suposición de velocidad constante del agua se utilizan en el cálculo de un tiempo de recorrido modelo.

Cada una de estas mediciones de campo permite que se calcule un tiempo de recorrido modelo mediante el uso de ecuaciones conocidos (ecuaciones (1) y (2) anteriores y también reproducidas en la Figura 2). Una vez que se recupera un OBN 16 y los datos adquiridos en cada uno se cargan en una computadora para análisis, el tiempo de recorrido modelo puede compararse a través del tiempo de recorrido real o "recolectado" que se calcula en base a los datos reales presentes en cada OBN 16. Con objeto de modelar de manera más precisa los datos y de cuantificar las variaciones/fluctuaciones debidas a la naturaleza dinámica de la columna de agua y otros errores en las mediciones de campo, se lleva a cabo sobre los datos una gradiente conjugada (o puede utilizarse otro modelo matemático) utilizando una serie de tablas de búsqueda, según se representa de manera estructural en las Figuras 3a-3c, y una computadora que se programa para ejecutar iteraciones de las ecuaciones de inversión (1) y (2) anteriores al perturbar o dar libertad primero a los datos de posición La x y la y de cada receptor, perturbar o dar libertad, en segundo lugar, a derivación cuadrática de reloj, y al perturbar, en tercer lugar, los datos de posición X, Y, y Z de las fuentes, según se describirá de manera más completa a continuación. Los datos pueden además uniformarse y reducirse el error de RMS entre el tiempo de arribo modelo y el tiempo de arribo recolectado al ejecutar perturbaciones de los datos de posición Z de cada receptor, derivación lineal de reloj y derivación volumétrica de reloj. Después de que se ejecuta cada perturbación, puede producirse una representación gráfica del error de RMS mediante el programa de computadora y observarse para análisis, un ejemplo de lo cual se proporciona en la Figura 7d. Además, la creciente mejora en la reducción de error de RMS puede observarse gráficamente en una serie de 10 gráficas circulares, según se representa en la Figura 6, en donde, el error de RMS para cada receptor puede observarse gráficamente después de perturbar cada una de las variables medidas en campo a través del uso de las tablas de búsqueda y aplicación de la gradiente conjugada.

Para entender aún más la naturaleza de la invención, puede ser útil una descripción del problema a resolverse.

El tiempo de arribo ty de la onda directa de la Fuente i medida por el Nodo j se escribe donde ti es el tiempo de la Fuente i, XSÍ = (XSÍ, ysi, ZSÍ) es la ubicación de la Fuente i, XNj = (XNj, yNj, ZNj) es la ubicación en el lecho marino del Nodo j, y s es lentitud en la columna de agua. La ecuación es una familiar excepto tal vez por la cantidad dy que es el desfase del reloj del Nodo j en el momento ti. El desfase del reloj se provoca por variaciones en la frecuencia de un oscilador de cristal de un nodo. Para este estudio tomamos que la variación sea la suma de un error de desplazamiento de frecuencia fijo y el asi llamado término de envejecimiento: donde fj es la frecuencia de oscilador del Nodo j en Hertz al momento t, fj es la frecuencia designada del oscilador, aj caracteriza el error de frecuencia fija y bj la velocidad de envejecimiento del oscilador por día (división entre 86,400 da la velocidad de envejecimiento por segundo). Finalmente, Tj es el momento en el cual se sincroniza el reloj. Es directo para mostrar que el desfase del reloj resultante, en segundos, al momento de la Fuente i se da por Claramente, un error de desplazamiento de frecuencia conduce a la derivación lineal a través del tiempo: un valor de aj igual a 10, por ejemplo, contribuiría 86.4 segundos al día a la derivación total. El envejecimiento da como resultado la variación cuadrática a través del tiempo: un valor de bj igual a 2, por ejemplo, agregaría 21.6 milisegundos de derivación después de 50 días. (Ver Olofsson y Woje (2010) para una discusión comprensiva de la sincronización de reloj, incluyendo comportamiento del oscilador más complejo que la ecuación (2).

Observe que la lentitud s, en la ecuación (1) depende del tiempo de la fuente así como también la trayectoria de recorrido entre la fuente y el nodo.

Motivados por la disponibilidad de las mediciones de velocidad del sonido del ROV en estudios de nodo (ver, por ejemplo, Hays et al., 2008, Figura 3) parametrizamos la lentitud en la columna de agua como En esta ecuación, so es una lentitud de fondo fija, obtenida, por ejemplo, mediante interpolación de datos de ROV registrados en momentos de tránsito de ROV. La cantidad s es una corrección constante (aplicada en toda profundidad) a la lentitud de fondo en el momento t. La aproximación en la ecuación (4) indica que las variaciones laterales en la columna de agua en cualquier momento dado son insignificantes. (Esta suposición puede violarse, por ejemplo, a medida que las corrientes oceánicas empujan agua nueva hacia el área de estudio).

Las cantidades desconocidas que típicamente buscamos obtener a través de la inversión de los tiempos de arribo directos, se mencionan a continuación: en donde ns es el número de fuentes, nN es el número de nodos, y la duración de tiempo del estudio se ha dividió a intervalos n0-l, con el intervalo de tiempo kw unido por las correcciones de lentitud constante Ok y s¾+i. Las variables conforman los elementos de un vector modelo m. Dada la primera suposición, mp, en el modelo linealizamos la ecuación (1) obteniendo T (m)«T(raa)-A(L,)ám . (6) Aquí, T es un vector de una primera recolección de tiempo de arribo. A es una matriz de derivados parciales y Am es un vector de perturbaciones a la primera suposición. (Las recolecciones son relativas a los momentos de fuente; en otras palabras, ti se toma hacia el lado izquierdo de la ecuación (1) con objeto de formar T). La ecuación (6) se resuelve para Am. Para problemas grandes, típicamente se utiliza un método iterativo tal como gradientes conjugadas. Para problemas pequeños, la descomposición de valor singular o svd es factible y reveladora. Aplicamos este procedimiento después de un pequeño estudio sintético.

Ejemplo Sinté tico La Figura la muestra el diseño "planeado" del estudio. Se colocan 25 nodos separados a 1000 m y 225 fuentes separadas a 500 m. El estudio se adquiere desde el fondo hacia la parte superior y de izquierda a derecha, un disparo cada seis horas durante 56 dias. Los relojes se sincronizan un día antes de que se inicie el disparo. Las variables en el problema son aquellas mencionadas en (5) con anterioridad. La primera suposición, mo, utiliza las posiciones x-y mostradas en la figura. La primera suposición en la profundidad del nodo se toma a 1000 m para todos los nodos y la primera suposición en los coeficientes de derivación es cero para todos los nodos. La velocidad de fondo se asume constante a 1510 rns-1 y una corrección constante al fondo se considera a intervalos de un día. El número total de variables en el problema es por lo tanto igual a 2x225 + 5x25 + 57 =632. Finalmente, el permitir desplazamientos de hasta 3000 m produce un total de 2442 recolecciones. Los valores de parámetro verdaderos, aquellos usados para calcular los verdaderos tiempos de arribo, se describen posteriormente.

Svd revela cuatro vectores singulares en el espacio nulo de A en este problema. Tres de los vectores singulares describen la traslación y rotación del estudio en el plano x-y. La Figura Ib ilustra uno de estos. Para obtener esta figura, simplemente representamos a escala el vector singular y le agregamos mo. Se seleccionó una escala grande para propósitos de ilustración. Esta parte del espacio nulo se entiende fácilmente, ya que la traslación y/o rotación de todos los nodos y fuentes tiene efecto nulo sobre el tiempo de recorrido. De hecho, existe el mismo espacio nulo cuando se resuelve solo para x-y de la fuente y receptor.

El cuarto vector nulo se ilustra en la Figura le. En la inspección detallada de esta figura puede observarse que se ha expandido el estudio. También, los nodos son más profundos y la velocidad del agua es mayor. Para entender esta parte del espacio nulo, se considera extender el estudio lateralmente a fin de que cualquiera de los dos puntos inicialmente separados por la distancia lateral R se separen por la distancia aR. Enseguida, al mantener fijas las profundidades de disparo, se tensa verticalmente por el mismo factor a. Finalmente, se multiplica la velocidad por a. El siguiente efecto es que permanecen sin cambio los tiempos de recorrido. La Figura le tiene a=1.04.

No existe componente diferente de cero de los coeficientes de derivación en el espacio nulo arriba descrito. Sin embargo, la derivación se caracteriza prominentemente en un número de vectores singulares con valores singulares muy pequeños (no mostrados). Estos vectores singulares indican que los coeficientes de derivación a y b pueden relacionarse, pero existen contribuciones diferentes de cero provenientes también de las otras variables. Ilustramos al resolver el problema sintético, primero libre de ruido y después con datos ruidosos.

La fila superior de la Figura 2 resume las diferencias entre el modelo verdadero y el modelo de inicio usado en la inversión. El modelo verdadero tiene nodos y fuentes desplazados de manera aleatoria de los lugares planeados: los nodos se desplazan de manera lateral cuando mucho 10 m, mientras que las fuentes se desplazan hasta 5 m. La profundidad del nodo se encuentra dentro de 3 m por encima o por debajo de 1000 m. Los coeficientes de desfase del reloj se seleccionan de manera aleatoria en el rango de ±10 para a, y ±2 para b. La velocidad del agua es constante a través de todo el estudio a 1500 ms_1.

La fila media de la Figura 2 muestra las diferencias entre el modelo verdadero y el resultado de una inversión cuando los datos son primeros tiempos de arribo libres de ruido. Hemos utilizado nuevamente svd (se excluyen los cuatro vectores nulos). La inversión es exacta en este caso. La fila inferior de la Figura 2 muestra el resultado de una inversión con datos ruidosos. Aqui los primeros tiempos de arribo se contaminaron con ruido aleatorio de manera uniforme en el rango de ±1 ms. La figura indica que las ubicaciones del nodo se determinan de manera exacta dentro de 2 m de manera lateral en este ejemplo. Lo mismo es verdadero para la mayoría de las fuentes. (Los mayores errores ocurren cerca de los bordes del estudio en donde son pocas las recolecciones). Las pequeñas discrepancias en la profundidad del nodo (principalmente demasiado superficiales) parecen compensarse por pequeños errores en la velocidad del agua (principalmente demasiado lenta). Claramente, los coeficientes de derivación son sensibles al ruido; sin embargo, cuando se calcula la derivación total al final del estudio se encuentra que es exacta en un margen de unos cuantos milisegundos - los coeficientes de derivación se relacionan.

- : · · - = -; ! , ~ . . . -- · - · · : = . : - . ; ; :: - Figura 1 : (a) patrón del estudio de nodo sintético planeado, Los nodos se muestran en amarillo: las fuentes en rojo. Vw es la velocidad del agua, ZN es la profundidad del nodo y a y b son coeficientes de derivación de nodo. Esta es la primera suposición m0 en la inversión. (b) m0 + 2000Vi, donde Vi es uno de los tres vectores nulos que describen la traslación y rotación de la geometría del estudio. Permanecen sin cambio la velocidad, profundidades de nodo y coeficientes de derivación de nodo. (c) m0 + 2000Vi, donde Vi es un vector nulo que provoca expansión o contracción del estudio, combinado con cambio en la velocidad. Permanecen sin cambio los coeficientes de derivación de nodo.

NotfcXv NodoZ OasptmHoib- Oeqptarieito-li RaritXT Veteüa) i l ! | l Figura 2. Histogramas que muestran los resultados de inversiones de datos libres de ruido y ruidoso. De izquierda a derecha, las cantidades mostradas son: colocación de nodo lateral, profundidad de nodo, coeficiente de derivación de nodo a, coeficiente de derivación de nodo b, colocación lateral de fuente y velocidad en el agua. Fila superior: la diferencia entre el modelo verdadero y el modelo de inicio, es decir, m - mo. Fila media: la diferencia entre el modelo verdadero y el modelo invertido para recolecciones libres de ruido. Fila inferior: la diferencia entre el modelo verdadero y el modelo invertido cuando *1 ms de ruido se agrega a las recolecciones.

A un nivel elevado, según se muestra en el diagrama de flujo de la Figura 3, un aspecto de la presente invención proporciona un método que involucra tomar mediciones de campo de coordenadas de fuente y de receptor, velocidad del sonido en el agua como una función de profundidad y tiempo, sincronización del receptor, y desfase del reloj se recolectan primero. Los datos sísmicos se examinan entonces para medir el tiempo de recorrido de cada fuente hacia cada receptor. Puede entonces calcularse un tiempo de recorrido modelo en base a las mediciones de campo. Al perturbar de manera iterativa al menos uno de los datos medidos en campo usando una tabla de búsqueda y calcular el tiempo de recorrido después de cada perturbación hasta que se ha logrado un error de RMS aceptable, los datos sísmicos acondicionados que toman en cuenta la naturaleza dinámica de la columna de agua proporcionarán la base para crear un mapa sísmico exacto que no se afecta por las condiciones cambiantes del agua.

Haciendo referencia a los diagramas de flujo de las Figuras 4a y 4b, el proceso de gradiente, conjugada, más detallado, mediante el cual se reduce el error de RMS entre el tiempo de arribo modelo y el tiempo de arribo recolectado es el siguiente: Para comenzar, en la etapa 100, se determinan los coeficientes de correlación para cada receptor a fin de realizar una determinación cualitativa de recolección, y si existe cualquier par de fuente-receptor con coeficientes de correlación por debajo de un cierto umbral (por ejemplo, 95%), se retiran de la consideración aquellos pares de fuente-receptor y no se analizan. Sin embargo, en el aguas profundas existe típicamente muy poco ruido y los coeficientes de correlación generalmente son muy elevados. Enseguida, en la etapa 102, se determina un tiempo de recorrido inicial mediante el uso de una suposición de velocidad de agua constante (por ejemplo, suponiendo que la velocidad de la columna de agua no cambia a través del tiempo).

La siguiente etapa 104 utiliza las tablas de búsqueda para ejecutar perturbaciones mediante actualización o dando libertad a los datos de posición de receptor la x y la y. Para ciertos receptores que tienen mayor error de RMS, de manera simultánea a la actualización de datos de posición de receptor la x y la y, también se actualizará el desfase del reloj (lineal, cuadrática y volumétrica). Las perturbaciones se ejecutarán de manera continua hasta que se haya logrado un error de RMS aceptable. Si después de que se ha ejecutado un umbral predeterminado de perturbaciones aún existen ciertas recuperaciones con error de RMS inaceptablemente elevado, aquellos receptores simplemente se retiran de la tabla de búsqueda y no vuelven a introducirse en el análisis hasta después de que se han ejecutado por completo todas las perturbaciones, como se explicará en lo sucesivo.

Enseguida, en la etapa 106, los datos de posición x, y, y z, para cada fuente, se permiten la libertad de ser actualizados y se ejecutan perturbaciones usando datos provenientes de la tabla de búsqueda. Estas perturbaciones se ejecutan continuamente hasta que se ha logrado un error de RMS aceptable, según se ilustra en la etapa 106.

Para refinar aún más el error de RMS, los datos de posición Z de los receptores se permiten entonces la libertad de ser actualizados usando datos atraídos desde las tablas de búsqueda, según se ilustra en la etapa 108. Después de cada actualización, se ejecuta la perturbación y se determina el error de RMS en la etapa 108. Una vez que el error de RMS es aceptablemente bajo en los datos de posición de receptor z, se utilizan entonces las tablas de búsqueda para permitir libertad a el desfase del reloj cuadrática, en la etapa 110. Después de que se han ejecutado suficientes perturbaciones con libertad dada a el desfase del reloj cuadrática, se ejecuta el mismo modelo de perturbación mientras se da libertad a el desfase del reloj lineal en la etapa 112.

Una vez que el error de RMS es lo suficientemente bajo después de haber dado libertad a el desfase del reloj lineal, los receptores previamente retirados (de la etapa 104) se agregan nuevamente a las tablas de búsqueda en la etapa 114, y las perturbaciones se ejecutan de manera continua solo sobre aquellos receptores hasta que el error de RMS para cada uno se aproxima al error de RMS para los receptores circundantes.

Con relación a las Figuras 7a-7d, se muestran gráficas que pueden producirse mediante el uso de la presente invención, que reflejan los limites de sincronización de nodo/receptor. En la Figura 8 se muestra que un pequeño número de receptores/nodos reciben la habilidad de actualizar sus valores de derivación cuadráticos. Los circuios negros representan ubicaciones de nodo donde fue correcta la derivación medida en campo y no se modificó en la gradiente conjugada/de inversión.

En la Figura 7b, el error de tiempo de RMS inicial para cada nodo muestra que necesitan hacerse mejoras a las posiciones de disparos y receptores, sincronización de los nodos y el modelo de velocidad. La mayoría de los nodos bajo las condiciones iniciales tuvo un error de RMS de entre 0.8 y 1.3 mseg. En la Figura 7c, la gráfica representa los nodos después de que se ha ejecutado la inversión de la línea de tiempo y muestra un error de RMS muy pequeño. Finalmente, la Figura 7d muestra la gráfica después de que se ha ejecutado la gradiente conjugada, con errores de RMS de menos de 0.08 mseg para cada nodo.

Con respecto a la gráfica de diez círculos de la Figura 6, cada uno de los diez círculos muestra gráficamente el error de RMS del nodo para el cual se creó la gráfica en diferentes etapas de acondicionamiento de datos. El primer círculo 200 muestra gráficamente los coeficientes de correlación para los disparos recibidos por el nodo, y en el ejemplo de esta Figura, la correlación es muy elevada para todos los pares de disparo. El siguiente círculo 202 muestra gráficamente el error de RMS inicial con la suposición de velocidad constante que se había hecho, siendo los errores demasiado elevados para una representación gráfica, sísmica, confiable. El siguiente círculo 204 muestra gráficamente el error de RMS después de que se ha ajustado el nodo a la sincronización volumétrica, pero manteniendo la suposición de velocidad constante; produce un índice de error de RMS también que es inaceptablemente elevado para una representación gráfica, sísmica, confiable. El siguiente círculo 206 muestra gráficamente el error de RMS después de que se ha aplicado el modelo inicial de velocidad dinámica en base a los datos de medición de campo de PCT; el error de RMS después del modelo inicial de velocidad dinámica se ha aplicado en base a los datos de medición de campo de PCT; el índice de error de RMS permanece inaceptablemente elevado para representación gráfica sísmica en 4D confiable. El siguiente círculo 208 muestra gráficamente el modelo de velocidad dinámica del círculo 206, pero con la actualización de lentitud aplicada; estos mejora el error de RMS, pero el error aún es inaceptablemente elevado para representación gráfica sísmica en 4D confiable. El siguiente círculo 208 muestra gráficamente el error de RMS después de que se han actualizado las posiciones del receptor la x y la y; nuevamente, el error de RMS se hace más pequeño, pero aún requiere de mejora para tener una representación gráfica sísmica en 4D confiable. El siguiente círculo 210 muestra gráficamente el error de RMS después de que se han actualizado las posiciones de disparo x, y, y z; el error de RMS ha disminuido nuevamente y se está haciendo más uniforme, pero podría mejorarse después a través de perturbaciones adicionales. El siguiente circulo 212 muestra gráficamente el error de RMS después de que se ha actualizado la posición de receptor z; nuevamente mejora del error de RMS. El siguiente circulo 214 muestra gráficamente el error de RMS después de que se ha actualizado el desfase del reloj lineal, y entonces el circulo final 216 muestra gráficamente el error de RMS después de que se ha actualizado la derivación cuadrática de reloj. Existe poca diferencia entre los últimos tres circuios, ya que la derivación lineal y cuadrática no se tomaron en cuenta a si mismas para demasiados errores de sincronización, pero después de que la gradiente conjugada se ha ejecutado por completo el error de RMS alcanzó un nivel aceptablemente bajo para proporcionar capacidades exactas y confiables de representación gráfica sísmica en 4D.

Aunque el proceso se ha descrito como una operación secuencial, el procesamiento a través de dimensiones puede ocurrir, y preferentemente lo hace, de manera simultánea, con objeto de reducir el tiempo de procesamiento. La gráfica de diez círculos ejemplificada en la Figura 6 se utiliza para permitir un disparo difícil y deducir las libertades que se le dan, pero el círculo más importante es el último, 218. Para llegar al último círculo, el proceso de inversión preferentemente resuelve todos los términos a los que se dio libertad de manera simultánea. Además, aunque se ha descrito una gradiente conjugada como una manera preferida de resolver la ecuación lineal, también podrían utilizarse otras fórmulas matemáticas sin apartarse del alcance y espíritu de la presente invención, según se define por las reivindicaciones anexas.

Una vez que se ha reducido el error de RMS a través de todas dimensiones, los datos pueden utilizarse para generar mapas sísmicos exactos. Además, debido al bajo error de tiempo de recorrido de RMS asociado con los datos acondicionados, los posteriores estudios realizados en la misma rejilla pueden compararse de manera más exacta con los estudios anteriores a fin de permitir una determinación más exacta de los cambios que ocurren en los depósitos de hidrocarburo que se revelaron en los mapas sísmicos. Sin la reducción del error de RMS del tiempo de recorrido, los cambios en la velocidad de la columna de agua entre los estudios podrían conducir a determinaciones inexactas en los cambios que han ocurrido en el depósito de hidrocarburos. Por lo tanto, el uso del modelo de perturbación iterativa de la presente invención proporciona una herramienta útil para conducir estudios sísmicos en 4D.

Las tablas de búsqueda y procesamiento de las perturbaciones se realizan usando una computadora o un sistema computarizado, en donde una computadora almacena las tablas de búsqueda en una base de datos y se ha programado para leer la base de datos y procesar los datos de acuerdo con el proceso descrito en la presente. Esto podría realizarse en una estación de trabajo independiente o a través de múltiples computadoras que se encuentran en comunicación entre sí.

Ejemplo de Funcionamiento Un estudio de nodo de agua profunda, con profundidades de agua que varían desde 1600m hasta 2300m, se utiliza para demostrar el proceso de inversión del tiempo de recorrido. La esquina noroeste del estudio ejemplar se afectó por un fondo de agua inclinado, en lugares mayores de 8 grados. El fragmento de observación en el estudio se adquirió en aproximadamente 64 días.

Para el estudio ejemplar, las recolecciones con coeficientes de correlación menores de 0.92 se rechazaron de la inversión del tiempo de recorrido. La causa típica para coeficientes de correlación menores fue la presencia de arribos de refracción interferente. De otro modo, el OBN que se registra en ambientes de agua profunda es típicamente muy quieto, con pocos rastros de "ruido" eliminados por esta etapa.

Con relación a la Figura 6, esta gráfica muestra la velocidad de intervalo medida a 1350m de profundidad de cada tránsito de columna de agua de los ROVs y cargadores submarinos. El eje vertical es la velocidad y el eje horizontal es el día Julián de la medición. Se utilizaron dos ROVs, los rombos azul obscuro y circuios azul claro, y un cargador submarino sobre el cual el dispositivo de medición se cambió al estudio medio, cuadrados rosas y triángulos amarillos.

El parámetro inicial que tiene el mayor efecto sobre los resultados de la inversión del tiempo de recorrido es el modelo de velocidad de agua inicial. Las iteraciones tempranas de la inversión del tiempo de recorrido usaron todas las mediciones de velocidad de agua realizadas durante la adquisición. Durante cada tránsito de columna de agua de los ROVs (usados para desplegar y recuperar los nodos) y del cargador submarino (un dispositivo que eleva y disminuye muchos nodos hacia el fondo del mar para un despliegue y recuperación más eficiente de los nodos por ROV) se realiza una medición de las velocidades de la columna de agua. Estos tránsitos dieron como resultado 79 perfiles registrados de velocidad de la columna de agua, realizados a partir de cuatro dispositivos diferentes. La velocidad de intervalo a 1350m de profundidad de agua puede observarse en la Figura 1 para todas las mediciones realizadas a partir de cada dispositivo de medición. Después de la investigación, se determinó utilizar solo las mediciones del cargador submarino, los cuadrados rosa y triángulos amarillos, debido a que eran más estables y había mediciones realizadas con aquellos dispositivos. La diferencia de velocidad medida entre el extremo derecho de los cuadrados rosa y el extremo izquierdo de los triángulos amarillos se supuso inicialmente era una diferencia de medición volumétrica entre los dos dispositivos, pero el análisis de los tiempos de recorrido sísmico registrados mostró un cambio de velocidad significativo ocurrido entre los días Julián 147 y 152. Las velocidades medidas se asumieron espacialmente constantes para cualquier momento dado, variando así solo por profundidad de agua y tiempo de medición.

Cuando es posible, los valores medidos en campo para el desfase del reloj de receptor se reforzaron a través de todo el proceso de inversión. Para este estudio, las derivaciones de reloj medidas en campo se utilizaron para aproximadamente 85% de los nodos. Cuando se limita la inversión para permitir actualizaciones hasta el 15% de nodos con mediciones de derivación no confiables observadas, los resultados para las posiciones de disparo y de receptor, sincronización del receptor y velocidades de la columna de agua a través del estudio se hacen más estables.

La Figura 8 muestra la diferencia de RMS final de los tiempos de arribo directo modelados y los tiempos de arribo directo recolectados reales para el estudio ejemplar. En esta figura, el mapa de nodos muestra la diferencia de RMS de los tiempos de arribo directos recolectados y los tiempos de arribo modelados en base a los resultados finales de la inversión de tiempo de recorrido. La mayoría de los errores de RMS es menor a 0.1 mseg. Los nodos en el área superior izquierda del mapa se influyen por una pendiente inclinada en el fondo del agua y de este modo se refleja más energía a desplazamientos más breves. Las barras de color varían desde 0.0 mseg hasta 0.5 mseg de error de RMS.

Conclusiones en base al Ejemplo La inversión del tiempo de recorrido es un proceso vital para identificar y corregir los errores de colocación de disparo y receptores en adquisición de OBN. Antes de la inversión del tiempo de recorrido, los errores de RMS de los tiempos de arribo directos, recolectados versus modelados, se encontraron entre 1 mseg y 2 mseg, lo cual podría ser un nivel aceptable para formación de imagen en 3D pero no para análisis en 4D. Los resultados de la inversión de tiempo de recorrido tienen errores de RMS para la mayoría de los nodos menores de 0.1 mseg. El proceso construye un modelo de velocidad de alta frecuencia que varía por tiempo de disparo y profundidad del agua.

El modelo de velocidad de la columna de agua de alta frecuencia se utilizará para corregir los datos para estadísticas de la columna de agua, lo cual es vital para el procesamiento de los datos de nodo a través de regularización, des-marcación y formación de imagen.

El proceso de inversión del tiempo de recorrido determina de manera exacta las correcciones de posición, correcciones de sincronización y actualiza la velocidad de la columna de agua de alta frecuencia en profundidades de agua que exceden aproximadamente 400m. Cuando son superficiales las profundidades de agua, el proceso tiene cierta dificultad en la determinación de las actualizaciones de sincronización de el desfase del reloj, correctas. Esto se debe principalmente al arribo directo en los nodos que se influyen por otros eventos a desplazamientos mucho más breves, conduciendo a una ventana de tiempo de adquisición mucho menor que se observa por los rastros utilizados para la inversión del tiempo de recorrido. Las actualizaciones de velocidad del agua derivadas por el proceso en estudios superficiales son correctas y pueden utilizarse para procesamiento ulterior corriente abajo.

Claims (15)

REIVINDICACIONES

1. Un método para reducir error de RMS en datos recolectados durante un estudio sísmico sobre un área predeterminada, que comprende las etapas de: a) obtener mediciones de campo de: i) posiciones de fuente x, y, y z; ii) posiciones del receptor x, y, y z; iii) velocidad del sonido en el agua como una función de tiempo y profundidad; iv) sincronización del receptor; y v) desfase del reloj; b) medir el tiempo de recorrido de cada fuente hacia cada receptor; c) determinar un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y d) perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

2. El método de acuerdo con la reivindicación 1, que comprende la etapa adicional de perturbar al menos una segunda de tales mediciones de campo, diferente de la primera, y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

3. El método de acuerdo con la reivindicación 2, que comprende la etapa adicional de perturbar una tercera de tales mediciones de campo, diferente de la primera y la segunda, y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

4. El método de acuerdo con la reivindicación 3, en donde tal primera medición de campo son las posiciones la x y la y de los receptores, la segunda de tales mediciones de campo son las posiciones x, y, y z de las fuentes, y la tercer medición de campo es el desfase del reloj.

5. El método de acuerdo con la reivindicación 1, que comprende además la etapa de: después de que cada perturbación se lleva a cabo para tal al menos una primera de tales mediciones de campo, desplegar gráficamente el error de RMS para cada posición de receptor.

6. Un programa de computadora almacenado en un medio de almacenamiento no transitorio y adaptado para ejecutarse en un procesador de computadora para propósitos de reducir el error de RMS en datos recolectados durante un estudio sísmico sobre un área predeterminada, comprendiendo un código de programa para: a) obtener mediciones de campo de: i) posiciones de fuente; ii) posiciones del receptor; iii) velocidad del sonido en el agua como una función de tiempo y profundidad; iv) sincronización del receptor; y v) desfase del reloj; b) medir el tiempo de recorrido de cada fuente hacia cada receptor; c) determinar un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y d) perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

7. El programa de computadora de la reivindicación 6, en donde el programa de computadora se programa además para perturbar al menos una segunda de tales mediciones de campo, diferente de la primera, y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

8. El programa de computadora de acuerdo con la reivindicación 7, que comprende además el código para perturbar una tercera de tales mediciones de campo, diferente de la primera y segunda mediciones de campo, y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

9. El programa de computadora de acuerdo con la reivindicación 8, en donde tal primera medición son las posiciones la x y la y de los receptores, la segunda de tales mediciones de campo son las posiciones x, y, y z de las fuentes y la tercer medición de campo es el desfase del reloj.

10. El programa de computadora de acuerdo con la reivindicación 6, que comprende además un código de computadora que provoca error de RMS para que cada posición de receptor se despliegue gráficamente después de que se lleva a cabo cada perturbación para tal al menos una primera de tales mediciones de campo.

11. Un sistema de computadora que incluye un programa de computadora almacenado en un medio de almacenamiento no transitorio y adaptado para ejecutarse en un procesador, reduciendo el programa de computadora de manera iterativa el error de RMS en los datos recolectados durante un estudio sísmico sobre un área predeterminada, comprendiendo el código de computadora para: a) obtener mediciones de campo de: i) posiciones de fuente; ii) posiciones del receptor; iii) velocidad del sonido en el agua como una función de tiempo y profundidad; iv) sincronización del receptor; y v) desfase del reloj; b) medir el tiempo de recorrido de cada fuente hacia cada receptor; c) determinar un tiempo de recorrido modelo en base a cada una de tales mediciones de campo; y d) perturbar al menos una primera de tales mediciones de campo y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

12. El sistema de computadora de acuerdo con la reivindicación 11, en donde el programa de computadora se programa además para perturbar por lo menos una segunda de tales mediciones de campo, diferente de la primera, y re-calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

13. El sistema de computadora de acuerdo con la reivindicación 6, que comprende además un código de computadora para perturbar una tercera de tales mediciones de campo, diferente de la primera y segunda mediciones, y re calcular el tiempo de recorrido hasta que se ha calculado un error de RMS aceptable.

14. El sistema de computadora de acuerdo con la reivindicación 13, en donde tal primera medición son las posiciones del receptor la x y la y de los receptores, la segunda de tales mediciones de campo son las posiciones x, y, y z de las fuentes, y la tercer medición de campo es el desfase del reloj.

15. El sistema de computadora de acuerdo con la reivindicación 11, que comprende además un código de computadora que provoca que el error RMS para cada posición de receptor se despliegue gráficamente después de que se lleva a cabo cada perturbación para tal al menos una primera de tales mediciones de campo.