Tempos de transito multiparametricos : estimação e inversão

Abstract

Orientadores : Martin Tygel, Lucio Tunes dos SantosTese (doutorado) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação CientificaResumo: Nesta tese, desenvolvemos um método, baseado no Algoritmo de Dix, para a estimação de um modelo de velocidade em profundidade. O método de Dix é capaz de construir um modelo com várias camadas homogêneas, separadas por interfaces curvas. Nosso método também gera modelos com essa estrutura. Porém, a velocidade em cada camada é uma função afim da profundidade. Desta forma, permitimos a presença de heterogeneidades, possibilitando um ajuste muito melhor dos dados. Além disso, utilizamos, como dados de entrada, os vários parâmetros cinemáticos, obtidos de maneira automática, pela técnica conhecida como Common Reflection Surface (CRS). Essa técnica é capaz de lidar diretamente com dados de multicobertura, fazendo um uso consistente de toda a informação disponível no levantamento sísmico. As estratégias convencionais usam quantidades cuja obtenção já se mostra um problema em si. Dado o maior grau de liberdade dos dados usados pelo método CRS, o tempo de trânsito depende de vários parâmetros. Estudamos também o problema de como estimar os parâmetros cinemáticos do tempo de trânsito que seriam usados para posterior inversão. Para tanto, foi aplicado um recente algoritmo de otimização, conhecido com Gradiente Espectral Projetado. Resultados promissores indicam o potencial de tal abordagem. Além disso, apresentamos uma aplicação para o modelo de velocidade construído. Desenvolvemos uma técnica para a obtenção de curvas de amplitude por afastamento (AVO) e amplitude por ângulo (AVA). Estas curvas são usadas pela indústria do petróleo para a caracterização de reservatórios. Sua obtenção implica um elevado custo computacional. O método proposto consiste em estimar o fator de espalhamento geométrico através de traçamento de raios no modelo invertido. Com isso, obtém-se as curvas AVO e A VA, para pontos de interesse em profundidade, a baixos custos computacionais. Os resultados para modelos sintéticos mostraram-se muito promissoresAbstract: For a fixed, central ray in an isotropic elastic or acoustic media, traveltime moveouts of rays in its vicinity can be described in terms of a certain number of parameters that refer to the central ray only. The determination of these parameters out of multicoverage data leads to very powerful algorithms that can be used for several imaging and inversion processes. Assuming two-dimensional propagation, the travei time expressions depend on three parameters directly related to the geometry of the unknown model in the vicinity of the central ray. We present a new method to extract these parameters out of coherency analysis applied directly to the data. Application of the method on a synthetic example shows an excellent performance of the algorithm both in accuracy and efficiency. The results obtained so far indicate that the algorithm may be a feasible option to solve the corresponding, harder, full three-dimensional problem, in which eight parameters, instead of three, are required. In conventional processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macro-velocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with NMO velocities derived from a previous velocity analysis conducted on CMP data. We presents a modified version of the Hubral and Krey algorithm that is adapted to advantageously use the above described parameters (the CRS attributes) as its input. Moreover, each velocity layer is no longer restricted to be a constant, being now a affine function on depth. Finally, we present a method to obtain a true-amplitude migration and amplitude-versusangle (AVA) at selected points using the attributes generated by the CRS attributes. Our approach combines the CRS stackjinversion process applied to multicoverage data, together with the use of a kinematic Kirchhoff migration, to achieve true-amplitudes (TA) at assigned depth points of the migrated images. Our method is designed to aggregate amplitude information on selected points of a reflector, after a purely kinematic image (migration) has been obtained. The method is tested on a synthetic inhomogeneous layered model with good resultsDoutoradoDoutor em Matemática Aplicad

    Similar works