INVERSÃO SÍSMICA TOMOGRÁFICA USANDO NORMA DE INTEGRAL DE FUNÇÃO

Abstract

O alto consumo de tempo em processamento computacional é um grande problema geralmente associado aos métodos de tomografia de raios sísmicos. Isto ocorre porque, em cada passo do processo iterativo da inversão definida pelo método tomográfico, o problema do "two-point ray-tracing" deve ser resolvido para cada par de fonte-receptor. A fim de resolver este tipo de problema, a norma vetorial, comumente usada nas funções de erro a ser minimizado no processo de inversão, é substituída por uma norma da integral de função, que nos permite estimar parâmetros do modelo através da minimização da área entre as curvas de tempos observados e calculados que são interpoladas (ou ajustadas) aos pontos dos dados. Desenvolvimentos matemáticos e experimentos numéricos relativamente simples, com modelos de campos de velocidade sísmica compressional (polinomialmente parametrizados), mostram que a norma da integral de função nos permite poupar uma enorme quantidade de tempo de processamento sem uma importante perda de precisão. Em alguns casos, os parâmetros do modelo podem ser melhor estimados com a norma de função do que com a norma vetorial, tradicionalmente utilizada na inversão de tomografia sísmica.