Procedimento de otimização computacional na solução da inversão não linear com procedimentos de inversão linear

Apresentamos uma otimização ao procedimento de Inversão Gravimétrica regularizada proposta por RAMOS em 2013, para obtermos estimativas do relevo do embasamento de bacia sedimentares ainda mais eficientes computacionalmente. Este procedimento realiza a Inversão não linear utilizando formulação semelhante ao da inversão linear a cada iteração. Neste procedimento é utilizada uma matriz semelhante à matriz sensibilidade em nossa otimização computacional consiste em simplesmente não avaliar os elementos desta matriz que receberiam valores próximos de zero, atribuindo valor zero a estes elementos, dispensando assim dispendiosos cálculos. Mostramos em testes numéricos que a influência significativa da anomalia gravitacional somente é gerada pela massa de prova se esta se localiza a certa distância do ponto de observação. A partir dos valores da observação numérica estabelecemos quais elementos da matriz utilizada devem ser avaliados e quais simplesmente podem ser negligenciados recebendo valor zero. Com esta matriz construída, avaliamos em testes sintéticos a efetividade da estimação em recuperar o relevo verdadeiro. Comparamos o tempo computacional para a estimação do relevo do embasamento de modelos sintéticos entre o procedimento original proposto por RAMOS e o procedimento otimizado que propomos. Os resultados evidenciaram a economia de tempo computacional, embora sutil se o número de parâmetros for pequeno, demonstrando a eficiência da proposta.