Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T12:07:56Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:42:09Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-10T15:30:42Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-18T11:32:11Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) / Made available in DSpace on 2014-09-18T11:32:11Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_AtenuacaoMultiplasCompressao.pdf: 24452115 bytes, checksum: 8631331284b78acc0b168236fbe6c6ea (MD5) Previous issue date: 2003 / O objetivo central deste trabalho é o estudo e a aplicação do método Kalman-Bucy no processo de deconvolução ao impulso e de deconvolução com predição, onde é considerado que os dados observados são classificados como não-estacionários. Os dados utilizados neste trabalho são sintéticos e, com isto, esta Tese tem características de um exercício numérico e investigativo. O operador de deconvolução ao impulso é obtido a partir da teoria de CRUMP (1974) fazendo uso das soluções das equações Wiener-Hopf apresentadas por KALMAN-BUCY (1961) nas formas contínuas e discretas considerando o processo como não estacionário. O operador de predição (KBCP) está baseado nas teorias de CRUMP (1974) e MENDEL ET AL (1979). Sua estrutura assemelha-se ao filtro Wiener-Hopf onde os coeficientes do operador (WHLP) são obtidos através da autocorrelação, e no caso (KBCP) são obtidos a partir da função bi(k). o problema é definido em duas etapas: a primeira consta da geração do sinal, e a segunda da sua avaliação. A deconvolução realizada aqui é classificada como estatística, e é um modelo fortemente baseado nas propriedades do sinal registrado e de sua representação. Os métodos foram aplicados apenas em dados sintéticos de seção fonte-comum obtida a partir dos modelos com interfaces contínuas e camadas homogêneas. / The main objective of this work is the study and the application of the Kalman-Bucy method in the processo f deconvolution to the impulse and deconvolution with prediction, considering the observed data as no stationary. The data used in this work are synthetic and, with this, this Thesis has characteristics of a numerical and search. The operator of deconvolution to the impulse is obtained from the Crump theory (1974), doing use of the solution of equation of Wiener-Holp presented by Kalman-Bucy in the continuoun and discrete forms considering the stacionary process. The prediction operator (KBCP) is based the Crump (1974) and Mendel et al (1979) theorics. Its structure resembles the Wiener-Hopf filter, where the coefficients of the operator are obtained through the autocorrelation, in the case (KBCP) are obtained from the function bi(k). The problem is defined in two steps: the first consists of the generation of the signal, and second of its evaluation. The deconvolution performed is classified as statistics, and is a model based in the properties of the registered signal and its representation. The method were applied only in synthetic data with common-shot section obtained from models with continuous interfaces and homogeneous layers.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-29T16:01:44Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-27T14:19:13Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) / Made available in DSpace on 2014-06-27T14:19:13Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_ProcessamentoInterpretativoDados.pdf: 4112570 bytes, checksum: 6957010292e9aff837e775e8d78ca0ed (MD5) Previous issue date: 2007 / ANP - Agência Nacional do Petróleo / Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área. / We present two potential field data interpretation methods applied to hydrocarbon prospecting. The first one uses aeromagnetic data to estimate the horizontal projection of the limit between the continental and oceanic crusts. This method is based on the existence of geological sources, which are magnetized and belongs exclusively to the continental crust. In this way, the estimates of the sources’ ends are used as estimates of the continental crust limits. The total-field anomaly measured above the continental shelf, continental slope and part of the continental rise is amplified using the downward continuation operator implemented in two different ways: the equivalent layer principle, and the Dirichlet boundary condition. Most of the computational load in computing the downward continued anomaly comes from the solution of a large-scale system of linear equations. This computational effort has been reduced not only by processing the whole area by moving windows of smaller dimensions, but also by the use of the conjugate gradient in the solution of the system of equations. Because the downward continuation operator is unstable, it was stabilized through the first-order Tikhonov stabilizing functional. Tests with noisecorrupted synthetic data have shown the efficiency of both implementations to enhance the termination of magnetic sources belonging to the continental crust, allowing, in this way, the estimation of the limit between the continental and oceanic crusts. Both implementations have been applied to two different areas offshore the Brazilian coast: Foz do Amazonas and Bacia do Jequitinhonha. The second method simultaneously delineates the basement topography and the geometry of salt structures occurring within the sedimentary rocks using gravity data. The interpretation models consist of a set of vertical, 2D juxtaposed prisms for the sedimentary pack, and of 2D horizontal prisms with polygonal cross sections for the salt structures. The solution has been stabilized by incorporating the geometric characteristics of the basement relief and of the salt structures, which are compatible with the a priori knowledge about the geological setting. To this end, we imposed inequality constraints to the parameters of the interpretation model, and used the stabilizing functionals known as global smoothness, weighted smoothness, and mass concentration along selected directions. We applied the method to synthetic data produced by simulated intracratonic and marginal basins, presenting density contrast with the basement varying with depth and containing salt structures. The results have shown that the method is potentially useful in simultaneously delineating the faulted basement relief and the salt structures geometries. We applied the method to real data along two gravity profiles across the Campos and Jequitinhonha Basins, and obtained interpretations in accordance with the known geology of the area.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-02T15:50:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-06-30T13:51:25Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) / Made available in DSpace on 2014-06-30T13:51:25Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_SimulacaoPerfisNucleares.pdf: 17341957 bytes, checksum: f62c8143abaa89b233e8512e81b09a51 (MD5) Previous issue date: 1993 / PETROBRAS - Petróleo Brasileiro S.A. / A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P1, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm3 e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas. / The identification and description of lithological parameters of a formation are essential in the evaluation of complex formations. Based on this, the combination of the nuclear tool response in uncased wells has been used systematically. The resultant logs can be considered as the interaction between two distinct phases: • The radiation transport phase from a source to one or more detectors through the formation. • The detection phase that consists of the collection of radiation, its transformation into current pulses, and the spectral distribuition of these pulses. As the presence of the detector does not strongly affect the radiation transport result, each phase can be simulated independent of the other, which allows us to introduce a new type of model in which the transport phase and the detection phase are uncoupled. In this work, the final response is simulated combining transport numerical solutions with a library of the detector responses to different incident energies and for each specific source - detector array. The radiation is calculated by the Finite Elements Method (FEM), as a 2½-D scalar flux derived from the numerical solution of the multigroup diffusion approximation of the Boltzmann transport equation in phase space. This is known as P1 approximation, where the variable direction is expanded in terms of the Legendre orthogonal polymonials, leading to the dimensionality reduction of the problem in such a way as to let it be more consistent with the FEM, where the flux depends only on the spatial variable and the physical properties of the formation. The NaI(Tl) response function is obtained separately by the Monte Carlo method (MC) where the life of a particle within the scintillator crystal is reconstructed simulating interaction by interaction the position, direction and energy of the different particles using a random number technique with associated appropriate probabilities laws. Each type of interaction (e.g., Rayleigh, Photo-electric, Compton and Pair production) is determined similarly and the simulation is concluded when the detector response functions are convolved with the scalar flux. The final response is the pulse-height spectrum of the simulated system. From this spectrum, a set of channels called detection windows are then selected. The count rates in each window show different dependencies on density and lithology. This fact allows one to use a combination of these windows in determining the density and photoelectric absorption factor of the formation. According to the method developed in this work, the logs in both thin and thick layers can be simulated. The performance of the method has been tested in complex formations, mainly where the presence of clay minerals, feldspars and micas have produced considerable effects sufficient to perturb the final response of the sonde. The results show that it is possible to identify physical and lithological parameters in formations having densities between 1.8 and 4.0 g/cm3 and photoelectric absorption factors in the interval of 1.5 to 5.0 barns/e-. The concentrations of Potassium, Uranium and Thorium can be obtained through the introduction of a new system of calibration which corrects the effects due to high variances and negative correlations observed on the mass concentration of Uranium and Potassium. In the simulation of the CNL response using the Tittle polynomial regression algorithm, it is verified that due to the limited vertical resolution of this sonde, the porosity value is poorly measured for most layers of thickness less than the source - far detector spacing, thus it has application only in thick layers. A new method was developed to solve this problem; the contribution of the relative area of each layer within the maximum information zone is determined. Thus, this neutron porosity makes possible an in-depth evaluation of expected CNL porosity-lithology response, convolving that area factor with the local formation porosity index, considering only thick layers. The presence of perturbating minerals is solved by considering the formation as formed by a predominant base matrix mineral, totally saturated by fresh water; the rest of the components are then considered as a perturbation of this base case. These results enable the calculation of synthetic well logs that can be used in inversion schemes in order to get a more detailed quantitative evaluation of complex formations.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T12:53:03Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-07-03T14:41:41Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) / Made available in DSpace on 2014-07-03T14:41:41Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AvaliacaoGeofisicaGeologica.pdf: 9834744 bytes, checksum: a2405c7a3ab5bd40721dc8d10090d688 (MD5) Previous issue date: 2013 / UFPA - Universidade Federal do Pará / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / INCT - Instituto Nacional de Ciência e Tecnologia / PETROBRAS - Petróleo Brasileiro S.A. / A utilização dos métodos de reflexão sísmica na exploração e desenvolvimento de reservatórios de hidrocarbonetos ocorre devido à sua vasta e densa amostragem, tanto em área quanto em profundidade, aliada ao refinamento de técnicas para o tratamento dos dados de reflexão sísmica, a partir destes dados, são geradas seções sísmicas, que após a aplicação de tratamento adequado, são utilizadas na interpretação dos estratos e/ou estruturas geológicas da subsuperfície. Neste trabalho é feita uma análise Geofísica Geológica de duas linhas sísmicas reais 2D marinhas da porção de quebra de talude da Bacia do Jequitinhonha. Para tanto, foi realizado um conjunto de processamento sísmico com objetivo de atenuar as reflexões múltiplas comuns em dados marinhos, além disso, foram estimados os modelos de velocidade em profundidade, utilizados para determinação das seções sísmicas migradas em profundidade. Nestas foram identificadas as superfícies refletoras. Através da análise dessas superfícies foram feitas as marcações de sismofácies, com base nos conceitos iniciais da sismoestratigrafia, com a finalidade de avaliar a qualidade do produto derivado do processamento sísmico, empregado neste estudo, para uma interpretação sismoestratigráfica, a qual está fundamentada na análise dos padrões de terminações dos refletores e padrão interno das sismofácies. / The seismic reflection methods use at the exploration and development hydrocarbon reservoirs is due to its vast and dense sampling in both area and depth, combined with techniques refinement to processing seismic reflection data, from these data, are generated seismic sections, which after applying appropriate treatment, are used in the interpretation of the strata and/or the subsurface geologic structures. At this paper have made one Geophysics Geological analysis at two real marine seismic lines 2D of break portion embankment Jequitinhonha Basin. Therefore, we performed a set of seismic processing in order to mitigate the multiple reflections common marine data, moreover, were estimated velocity models in depth, used to determine the depth migrated seismic sections. At these were identified reflective surfaces. Through analysis at these surfaces markings of sismofacies. have made based on the initial concepts seismic stratigraphy in order to evaluate the product quality derived from the seismic processing used at this study for seismic stratigraphic interpretation, which is based on the analysis termination patterns of reflectors and internal standard of sismofacies.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-20T14:53:12Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Ausência de palavras-chave on 2014-08-07T12:41:40Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-13T14:27:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T14:39:54Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) / Made available in DSpace on 2014-08-20T14:39:54Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_EmpilhamentoMetodoSuperficie.pdf: 6623577 bytes, checksum: 67b656d135334821f2f993c805d88214 (MD5) Previous issue date: 2003 / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / CTPETRO - Plano Nacional de Ciência e Tecnologia do Setor Petróleo e Gás Natural / FINEP - Financiadora de Estudos e Projetos / O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais. / The CRS stacking method simulates ZO seismic sections from multi-coverage data and does not dependente on a macro-velocity model. For 2-D medium the stacking traveltime depends on three parameters: the emergence angle of the normal ray (with respect to the measurement surface normal) and the wavefront curvatures of two hypothetical waves, called Normal-Incidence-Point (NIP) wave and Normal (N) wave. The CRS method consists of summing the amplitudes of the seismic traces in the multicoverage data along the surface defined by CRS stacking traveltime which that fits best the data set. The result of the CRS stack is assigned to points of a grid pre-defined in the ZO section. As the result obtain a simulated ZO section. This means that for each point of the ZO section must be estimated the three optimal parameters that yield the maximum coherence between the events of seismic reflection. In this Thesis I present formulae for the 2-D CRS method and for the NMO velocity that consider the topography of the measurement surface. The algorithm is based on the optimization strategy divided into three steps: 1) To search for the emergence angle and the curvature of the NIP wave, by applying a global optimization, 2) to search for the curvature of the N wave, by applying global optimization, and 3) to refine the initial parameters estimated in first two steps by applying local optimization. In the first two steps is used the Simulated Annealing (SA) algorithm and in the third step the Variable Metric (VM) algorithm. For the case of a measurement surface with smooth topography the curvature of this surface is included in the 2-D CRS stack formalism. This CRS algorithm implemented was applied to synthetic data set. The result is a simulated ZO section of high quality, with a high signal-to-noise ratio, and the estimative of the parameter triplet. It is performed a sensibility analysis for the new CRS stacking traveltime with respect to the curvature in several points of the curved measurement surface. This study showed that the CRS traveltime is more sensitive for fast midpoints of the central points and larger offsets. The expressions for the NMO velocities presented here is applied to estimate the interval velocities and the depth of the reflectors for 2-D model with a smooth topography. For the inversion of the velocities and the depth of the reflectors is considered the Dix-type inversion algorithm. The NMO velocity for a curved measurement surface deserves to best estimate the velocities and the depths of the reflectors than NMO velocities referred to planar surfaces. Also, I present an introduction to 3-D stack. In this case, the stacking traveltime depends on eight parameters. These parameters can be obtained by using some parameter-search strategies that I have showed in this Thesis. The combination of the strategy of the Traveltime Approximations and the strategy of Arbitrary Curvatures is used to apply 3-D CRS stack successful in synthetic and real data sets, respectively.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-09T12:43:26Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:13:45Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T16:47:31Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-20T17:00:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) / Made available in DSpace on 2014-08-20T17:00:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoMetodosEletrico.pdf: 10195935 bytes, checksum: a6776849fa107eb79f60ab784dd5f935 (MD5) Previous issue date: 1998 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / CNPq - Conselho Nacional de Desenvolvimento Científico e Tecnológico / DCT - Departamento de Cooperação Científica, Técnica e Tecnológica / Este trabalho trata da aplicação dos métodos geofísicos eletroresistivo e eletromagnético no município de Altamira, região central do Estado do Pará, com a finalidade de fazer o mapeamento da espessura das camadas aqüíferas, bem como a ocorrência de intrusões de diabásio na área do município, permitindo uma melhor locação de pontos para perfuração com fins de captação à água subterrânea. Na área estudada, verificam-se afloramentos de rochas cristalinas do embasamento coberto por sedimentos com idade variando desde o Paleozóico ao Quaternário, com esperança de se obter água nos sedimentos do Devoniano Inferior, mais precisamente na formação Maecuru. Foram executadas 34 sondagens elétricas verticais com arranjo Schlumberger com AB/2 variando de 100 a 500 metros. 108 sondagens eletromagnéticas transientes também foram executadas com arranjo "loop coincidente", com 50 metros de lado. Por problemas eletrônicos adicionados aos ruídos de fenômenos supraparamagnéticos, os resultados das sondagens eletromagnéticas apresentaram-se muito distorcidos, não tendo sido possível extrair deles informação geológica de relevância. Por isso, foram criados modelos eletromagnéticos sintéticos representando a situação geológica da área. Os dados obtidos no campo em forma de curvas de resistividade aparente, foram plotados em papel bilogarítmo e interpretados com a ajuda de algoritmos direto e inverso os quais considera modelos de camadas horizontalmente estratificadas, homogêneas e isotrópicas. Os resultados da interpretação permitiram estabelecer modelos geoelétricos aproximados refletindo a situação geológica na subsuperfície. Os resultados da interpretação das sondagens elétricas, foram confrontados com os resultados da interpretação qualitativa do mapa de isogamas do campo magnético total, permitindo delimitar com maior precisão a zona de incidência de soleiras de diabásio, situada a norte da área. Dos quatro horizontes geoelétricos distinguidos, o terceiro é essencialmente constituído por arenitos, apresentando, por isso, melhores condições para armazenamento de água subterrânea. / This work due application of the electroresistivity and electromagnetic geophysical methods in the municipal district of Altamira, in the central area of the State of Pará, with the purpose of mapping the thickness of the aquifer layers and occurrence of diabase intrusions for finding a groundwater wells. In this area are many outcrops of basement crystalline rocks covered with sediments of ages from the Paleozioc to Quaternary. The objective is to obtain water from the sediments of Inferior Devonian, more exactly in the Maecuru Formation. There were made thirty four vertical electric soundings using the Shlumberger arrangement, the AB/2 distance varying from 100 to 500 meters. Also there were made 108 transient electromagnetic soundings with "coincident loop" arrangement, using an square loop of 50 meters. Due to electronic problems and supraparamagnetic noises, the results of this electromagnetic soundings became very distorted, impossible to get from them important geologic information. Then, electromagnetic synthetic models were created to represent the geology of the area. The data were obtained from the field were presented as apparent resistivity curves. They were plotted using a log-log scale and interpreted using the EGSLIB/SEV1D first and then with the EGSLIB/SEV1DIN program. The models used were horizontal, isotropic and homogeneous layers. The interpretation results allowed to determine approximate geoelectrical models which reflects the underground geological situation. The interpretation results from the electric surveys were confronted with the results of the qualitative interpretation from the isogams map of the total magnetic field, allowing to delimit with more precision the zone of incidence of diabase sills, located at the north of the area. From the four distinguished geoeletric horizons, the third is essentially constituted by sands, therefore it is best for groundwater storage.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-10T12:28:13Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-07T16:13:21Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-14T14:23:37Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:04:51Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) / Made available in DSpace on 2014-08-22T17:04:51Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_AplicacaoTratamentoEstatistico.pdf: 16343162 bytes, checksum: 64dcc4f3d1c980c6fc446241fa786161 (MD5) Previous issue date: 1991 / Utilizando-se dados do campo de Camorim (Bacia de Sergipe-Alagoas), foi testado e aplicado um conjunto de técnicas estatísticas multivariantes (análises de agrupamentos, de componentes principais e discriminante) no intuito de identificar as fácies, previamente definidas em poços testemunhados, por meio dos perfis, viabilizando o reconhecimento das mesmas nos demais poços não testemunhados da área. A segunda etapa do processo de determinação das fácies consistiu no emprego de métodos auxiliares (análises composicional e de seqüência de fácies), que, combinados com as técnicas multivariantes, propiciaram melhores resultados na calibração rocha-perfil. A determinação das fácies, uma vez estabelecida, possibilitou o refinamento do processo de avaliação de formações ao viabilizar o exame de cada fácies-reservatório, isoladamente. Assim, esse procedimento tornou possível a escolha, para cada litologia, dos parâmetros utilizados na interpretação dos perfis ao mesmo tempo em que permitiu a totalização em separado dos valores de espessura, porosidade e saturação dos fluidos, bem como a adoção de diferentes valores de corte (cut-offs) para cada grupo considerado. Outras aplicações incluíram a melhoria na estimativa da porosidade e da permeabilidade, a adaptação de algoritmos para o cálculo preliminar de porosidade, a confecção de mapas de fácies e a geração automática de seções estratigráficas. Finalmente, foram destacadas a perspectiva de integração desse estudo com sistemas estatísticos de descrição de reservatórios, outras técnicas de determinação de fácies em desenvolvimento e a retomada da utilização de métodos estatísticos multivariantes em dados de perfis, como ferramenta de exploração. / A series of multivariate statistical techniques (cluster, principal component and discriminant analysis) was tested and applied to well log data from the Camorim field (offshore Sergipe State, Brazil) in order to identify facies previously defined through core description. The second step in the process of facies determination was supported by auxiliary methods (compositional and sequence facies analysis), which produced better results in the calibration of rock-log data, when combined with the multivariate techniques. The facies determination, once established, permits the refinement of the process of formation evaluation, enabling the examination of each facies separately. This procedure made it possible to choose, for each lithology, the parameters used in log interpretation. In parallel, this process allowed the summation of thickness, porosity, fluid saturation and the adoption of different cut-offs for each group, separately. Other applications included: enhancement in the estimation of porosity and permeability, the adaptation of algorithms designed for fast porosity estimation, the mapping of variables useful in the characterization of the vertical variability of the reservoir rocks and the automatic generation of stratigraphic sections. Finally, the possibility of integration of the work from this study with statistical systems of reservoir description, other facies determination techniques currently being developed and the utilization of multivariate statistical methods in well log data, as an auxiliary exploratory tool, were illustrated.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-11T16:47:14Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Definir palavras-chave on 2014-08-06T16:30:33Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-11T13:24:15Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-22T17:09:08Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) / Made available in DSpace on 2014-08-22T17:09:08Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_CaracterizacaoGeoeletricaArea.pdf: 3315484 bytes, checksum: 279133bbca3dc766b7e2c2c5bc9bb8b2 (MD5) Previous issue date: 2005 / CAPES - Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Foi desenvolvido um estudo da subsuperfície rasa para investigar a possibilidade de contaminação do subsolo por líquidos usados no tratamento de couros em área industrial. O estudo foi realizado na área da indústria Couro do Norte, Bairro de Maracacuera, no Distrito Industrial de Icoaraci (município de Belém-Pará). Durante o estudo foram feitas medidas para a determinação do fluxo de água do primeiro aqüífero e aplicada metodologia geofísica com a finalidade de detectar a possível contaminação do subsolo por poluentes liberados pelo curtume. As metodologias geofísicas utilizadas foram: eletrorresistividade (na forma de imageamento e través de sondagens verticais), potencial espontâneo, eletromagnético (Slingram) e radar de penetração no solo (GPR). A interpretação dos dados eletromagnéticos permitiu identificar a presença de zonas mais condutivas (provável pluma de contaminação) e zonas mais resistivas (blocos lateríticos). A distribuição dos potenciais medidos através do Método do Potencial Espontâneo indicou o sentido do fluxo subterrâneo local. Com as Sondagens Elétricas Verticais determinaram-se os valores de resistividade em subsuperfície para os diferentes horizontes correlacionando-os sempre com as informações obtidas no Caminhamento Elétrico. Através dos dados de GPR identificou-se a presença de blocos lateríticos, tubulações de metal e zonas com elevada atenuação (zona de sombra) que indicam a presença de contaminação. / A study covering the shallow subsurface was developed in order to investigate the possibility of underground contamination by liquids used in the treatment of animal skin. The study was carried out inside the area of the Couro do Norte industry, located in the Icoaraci Industrial District (Belém – Pará). During the study they were made measures for the determination of the flow of water of the first aquifer and applied geophysical methodology with the purpose of detecting the possible contamination of the underground for pollutant liberated by the tanning. The geophysical methodology applied was: resistivity (imaging and vertical electric soundings), spontaneous potential, electromagnetic (Slingram) and ground penetrating radar (GPR). Analyses of electromagnetic data lead to the identification of conductive zones (Contamination) and resistive zones (sandy material). Spontaneous potential gave information on the direction of local underground flow of water. Vertical electric soundings allowed to know the resistivity distribution in subsurface. The resistivity sounding models could be correlated to the electrical images obtained in electrical profiling. Through the GPR data was possible to identify the presence of lateritic blocks, metal piping and zones with high attenuation (shadow zones) that indicate the presence of contamination.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-06-17T11:50:38Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DesenvolvimentoConjuntoModelagem.pdf: 2394159 bytes, checksum: da7e66022f1b791ea4cec64ebf4c05a9 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T13:13:33Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-08-25T14:00:07Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DesenvolvimentoConjuntoModelagem.pdf: 2394159 bytes, checksum: da7e66022f1b791ea4cec64ebf4c05a9 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-08-29T13:34:37Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DesenvolvimentoConjuntoModelagem.pdf: 2394159 bytes, checksum: da7e66022f1b791ea4cec64ebf4c05a9 (MD5) / Made available in DSpace on 2014-08-29T13:34:37Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Dissertacao_DesenvolvimentoConjuntoModelagem.pdf: 2394159 bytes, checksum: da7e66022f1b791ea4cec64ebf4c05a9 (MD5) Previous issue date: 2004 / A Modelagem Eletromagnética Analógica tem sido uma poderosa ferramenta no estudo das respostas eletromagnéticas de uma grande variedade das técnicas de levantamentos e de modelos geológicos. Nas aplicações de métodos eletromagnéticos nas investigações de contaminação no subsolo, recentemente em grande demanda, são obtidas anomalias muito baixas (da ordem de 0,01% do campo primário). Com o objetivo de estudar em laboratório os problemas de levantamentos eletromagnéticos para alvos de baixo número de indução, um conjunto de modelagem em escala reduzida de alta sensibilidade, foi projetado e construído. Este conjunto opera na faixa de freqüência de 1 kHz até 1 MHz com excelente linearidade. O transdutor de transmissão gera suficiente densidade de fluxo magnético com corrente de até 3 App (pico a pico) circulando nele. O sistema de recepção, isto é, o conjunto de bobinas receptoras mais o pré-amplificador acoplado apresenta sensibilidade 3 V/A/m na freqüência de 10 kHz com inclinação de 20 dB/dec de freqüência. A alta sensibilidade do sistema de recepção combinada com a boa capacidade de geração do sistema transmissor, permite alcançar a resolução de 0,025% do campo primário. O conjunto dispõe ainda, de um sistema mecânico que permite posicionar o sistema de transdutores (transmissor e receptor) em posições arbitrárias com resolução de 1 mm e precisão de 0,05 mm. A precisão no deslocamento linear do sistema de sonda numa linha de 1715 mm de comprimento é alcançada por meio de um motor de passo controlado por um microcontrolador, que é rigidamente acoplado a um fuso que converte o movimento angular do motor em movimento linear. Com este conjunto, é possível simular uma grande variedade de técnicas de levantamentos eletromagnéticos no formato Slingram, isto é, dipolo-dipolo. Além da boa resolução mecânica e boa resolução de medida eletromagnética do conjunto, todo o processo de controle e aquisição de dados é operacionalizado por meio de um computador mestre. / Electromagnetic analog modelling has been a powerful tool to study the electromagnetic response of a great varieties of survey techniques of different geological models. In the applications of The electromagnetic methods in the investigations of subsurface contaminations of in great demand recently, very low anomalies (of the order of 0.01% Hp) are obtained. To study in the laboratory the problems of the electromagnetic surveys of low induction number targets, an electromagnet reduced scale model set-up is designed and construed. This set-up operates in the frequency range off 1 kHz to 1 MHz. The developed transmitting system generate a high magnetic flux density by passing a 3 AP-P current in it. The receiving system, consist of a set coils and amplifier coupled in it, has a sensibility of 3 V/A/m at 10 kHz with gradient of 20 dB/dec frequency. The overall resolution of the set-up is 0.025% HP, obtained due to the high sensitivity of the received system to generate primary field. The set-up also has a very precise mechanical system that permits locate the coils at any arbitrary position with a resolution of 1 mm and a precision 0.5 mm. The Tx-Rx system cam be moved linearly with great precision up to 1715 mm in length, by rotating a spindle employing a microprocessor controlled stepper motor attached to it. With set-up, it is possible to simulate a great variety of electromagnetic survey techniques of Slingram systems, that is dipole-dipole configurations. Besides good mechanical and electromagnetic measurements resolutions; all the process of data acquisition and other operations in the set-up cam be controlled through a computer.

Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-05-26T13:04:20Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Rejected by Irvana Coutinho (irvana@ufpa.br), reason: Indexar os assuntos on 2014-08-06T14:52:53Z (GMT) / Submitted by Cleide Dantas (cleidedantas@ufpa.br) on 2014-09-08T12:41:18Z No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Approved for entry into archive by Irvana Coutinho (irvana@ufpa.br) on 2014-09-08T14:44:47Z (GMT) No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) / Made available in DSpace on 2014-09-08T14:44:47Z (GMT). No. of bitstreams: 2 license_rdf: 23898 bytes, checksum: e363e809996cf46ada20da1accfcd9c7 (MD5) Tese_Migracao3DKirchhoff.pdf: 13210015 bytes, checksum: 55e96c2ad9e534d8dce3117b9d9f3623 (MD5) Previous issue date: 2013 / CENPES - Centro de Pesquisas Leopoldo Américo Miguez de Mello / PETROBRAS - Petróleo Brasileiro S.A. / O Feixe Gaussiano (FG) é uma solução assintótica da equação da elastodinâmica na vizinhança paraxial de um raio central, a qual se aproxima melhor do campo de ondas do que a aproximação de ordem zero da Teoria do Raio. A regularidade do FG na descrição do campo de ondas, assim como a sua elevada precisão em algumas regiões singulares do meio de propagação, proporciona uma forte alternativa na solução de problemas de modelagem e imageamento sísmicos. Nesta Tese, apresenta-se um novo procedimento de migração sísmica pré-empilhamento em profundidade com amplitudes verdadeiras, que combina a flexibilidade da migração tipo Kirchhoff e a robustez da migração baseada na utilização de Feixes Gaussianos para a representação do campo de ondas. O algoritmo de migração proposto é constituído por dois processos de empilhamento: o primeiro é o empilhamento de feixes (“beam stack”) aplicado a subconjuntos de dados sísmicos multiplicados por uma função peso definida de modo que o operador de empilhamento tenha a mesma forma da integral de superposição de Feixes Gaussianos; o segundo empilhamento corresponde à migração Kirchhoff tendo como entrada os dados resultantes do primeiro empilhamento. Pelo exposto justifica-se a denominação migração Kirchhoff-Gaussian-Beam (KGB). As principais características que diferenciam a migração KGB, durante a realização do primeiro empilhamento, de outros métodos de migração que também utilizam a teoria dos Feixes Gaussianos, são o uso da primeira zona de Fresnel projetada para limitar a largura do feixe e a utilização, no empilhamento do feixe, de uma aproximação de segunda ordem do tempo de trânsito de reflexão. Como exemplos são apresentadas aplicações a dados sintéticos para modelos bidimensionais (2-D) e tridimensionais (3-D), correspondentes aos modelos Marmousi e domo de sal da SEG/EAGE, respectivamente. / The Gaussian Beam (GB) is an asymptotic solution of the elastodynamic equation in the paraxial vicinity of a central ray, which approaches better the wave field than the standard zero-order ray theory. The GB regularity in the description of the wave field, as well as its high accuracy in some singular regions of the propagation medium, provide a strong alternative to solve seismic modeling and imaging problems. In this thesis, i presenty a new procedure for pre-stack depth migration with true-amplitude, combining the flexibility and robustness of Kirchhoff migration type using superposition of Gaussian beams to represent the wave field. The proposed migration algorithm comprises in two stacking process: the first is the beam stack is applied to subsets of seismic data multiplied by a weight function defined such that stack operator has the same formulation of the integral of the Gaussian beams superposition; the second is a weighted diffraction stack by means of the Kirchhoff type integral having as input the GB stacked data. For these reasons it is called Kirchhoff-Gaussian-Beam (KGB) migration. The main characteristics that distinguish the KGB migration, during the first stage stacking, with other migration methods that also use the theory of Gaussian beams, is the use of the first Fresnel zone projected to limit the width of the subset of seismic traces (beam) using a second-order approximation of the reflection travel time. Examples are shown for applications on two-dimensional (2-D) and three-dimensional (3-D) synthetic seismic data, respectively, to the models Marmousi and SEG/EAGE salt dome data sets.