Modelagem e simulação de um separador bifásico com inversão de fases induzida

The present work has the main goal to study the modeling and simulation of a biphasic separator with induced phase inversion, the MDIF, with the utilization of the finite differences method for the resolution of the partial differencial equations which describe the transport of contaminant s mass fraction inside the equipment s settling chamber. With this aim, was developed the deterministic differential model AMADDA, wich was admensionalizated and then semidiscretizated with the method of lines. The integration of the resultant system of ordinary differential equations was realized by means of a modified algorithm of the Adam-Bashfort- Moulton method, and the sthocastic optimization routine of Basin-Hopping was used in the model s parameter estimation procedure . With the aim to establish a comparative referential for the results obtained with the model AMADDA, were used experimental data presented in previous works of the MDIF s research group. The experimental data and those obtained with the model was assessed regarding its normality by means of the Shapiro-Wilk s test, and validated against the experimental results with the Student s t test and the Kruskal-Wallis s test, depending on the result. The results showed satisfactory performance of the model AMADDA in the evaluation of the MDIF s separation efficiency, being possible to determinate that at 1% significance level the calculated results are equivalent to those determinated experimentally in the reference works

Projeto e Operação em Escala Semi-Industrial de um Equipamento para Tratamento de Águas Produzidas na Indústria do Petróleo Utilizando Nova Tecnologia: misturador-decantador à inversão de fases - MDIF

The aim of this work is to use a new technology in the treatment of produced wastewaters from oil industry. An unit for treat produced waters called UTMDIF, was designed, installed and operated in an industrial plant for treatment of effluents from oil industry. This unit operates by means of the method of separation of phase inversion and can become a promising alternative to solve the problem of oil/water separation. This method constitutes the basis of the working of a new design of mixersettler of vertical configuration which occupies small surface area. The last characteristic becomes specially important when there is limitation on the lay-out of the plant, for example, over maritime platforms to explore oil. This equipment in a semi-industrial scale treats produced wastewaters contaminated with oil at low concentrations (ranging from 30 to 150 mg/L) and throughputs of 320 m3/d (47,4 m3 m-2 h-1). Good results were obtained in oil/water separation which leads to the necessary specification to discharge those wastewaters. Besides, the non dependence of the efficiency of separation in spite of the salinity of the medium becomes the equipment an attractive new technology to treat wastewaters containing oil at low concentrations

Desempenho de escolares de 1ª série na bateria de identificação de erros de reversão e inversão na escrita: estudo preliminar

OBJETIVO: caracterizar o desempenho de escolares da 1ª série na bateria de identificação de erros de reversão e inversão na escrita. MÉTODOS: participaram deste estudo 30 escolares de 1ª série de ensino público do município da cidade de Marília-SP, de ambos os gêneros, na faixa etária de 7 anos a 7 anos e 11 meses de idade. Como procedimento foi realizada a adaptação brasileira da bateria de identificação dos erros de reversão e inversão na escrita. Esta bateria é composta por 3 testes: teste de desempenho contínuo, teste de memória de curto prazo e teste de controle de escrita. Cada teste é composto por subtestes que visam identificar a capacidade dos escolares em identificar as inversões e reversões de letras e números isolados e letras e números em sequências. RESULTADOS: os resultados revelaram que os escolares da 1ª série apresentaram desempenho superior para identificação de reversão e inversão de letras dentro de palavras do que isoladamente, identificação da reversão para letras e números com desempenho superior se comparado a inversão de letras na palavra. Além disso, os escolares deste estudo apresentaram menor tempo para identificação de reversão e inversão de letras e números isolados do que em sequência. CONCLUSÃO: a Bateria adaptada para este estudo permitiu conhecer o perfil dos escolares de 1ª série quanto a identificação dos erros de inversão e reversão, entretanto, é necessário novas aplicações para verificar a ocorrência de identificação destes erros em populações de diferentes fase de alfabetização.

Arcabouço geofísico, isostasia e causas do magmatismo cenozóico da província Borborema e de sua margem continental (Nordeste do Brasil)

The Borborema Province (BP) is a geologic domain located in Northeastern Brazil. The BP is limited at the south by the São Francisco craton, at the west by the Parnaíba basin, and both at the north and east by coastal sedimentary basins. Nonetheless the BP surface geology is well known, several key aspects of its evolution are still open, notably: i)its tectonic compartmentalization established after the Brasiliano orogenesis, ii) the architecture of its cretaceous continental margin, iii) the elastic properties of its lithosphere, and iv) the causes of magmatism and uplifting which occurred in the Cenozoic. In this thesis, a regional coverage of geophysical data (elevation, gravity, magnetic, geoid height, and surface wave global tomography) were integrated with surface geologic information aiming to attain a better understanding of the above questions. In the Riacho do Pontal belt and in the western sector of the Sergipano belt, the neoproterozoic suture of the collision of the Sul domain of the BP with the Sanfranciscana plate (SFP) is correlated with an expressive dipolar gravity anomaly. The positive lobule of this anomaly is due to the BP lower continental crust uplifting whilst the negative lobule is due to the supracrustal nappes overthrusting the SFP. In the eastern sector of the Sergipano belt, this dipolar gravity anomaly does not exist. However the suture still can be identified at the southern sector of the Marancó complex arc, alongside of the Porto da Folha shear zone, where the SFP N-S geophysical alignments are truncated. The boundary associated to the collision of the Ceará domain of the BP with the West African craton is also correlated with a dipolar gravity anomaly. The positive lobule of this anomaly coincides with the Sobral-Pedro II shear zone whilst the negative lobule is associated with the Santa Quitéria magmatic arc. Judging by their geophysical signatures, the major BP internal boundaries are: i)the western sector of the Pernambuco shear zone and the eastern continuation of this shear zone as the Congo shear zone, ii) the Patos shear zone, and iii) the Jaguaribe shear zone and its southwestern continuation as the Tatajuba shear zone. These boundaries divide the BP in five tectonic domains in the geophysical criteria: Sul, Transversal, Rio Grande do Norte, Ceará, and Médio Coreaú. The Sul domain is characterized by geophysical signatures associated with the BP and SFP collision. The fact that Congo shear zone is now proposed as part of the Transversal domain boundary implies an important change in the original definition of this domain. The Rio Grande do Norte domain presents a highly magnetized crust resulted from the superposition of precambrian and phanerozoic events. The Ceará domain is divided by the Senador Pompeu shear zone in two subdomains: the eastern one corresponds to the Orós-Jaguaribe belt and the western one to the Ceará-Central subdomain. The latter subdomain exhibits a positive ENE-W SW gravity anomaly which was associated to a crustal discontinuity. This discontinuity would have acted as a rampart against to the N-S Brasiliano orogenic nappes. The Médio Coreaú domain also presents a dipolar gravity anomaly. Its positive lobule is due to granulitic rocks whereas the negative one is caused by supracrustal rocks. The boundary between Médio Coreaú and Ceará domains can be traced below the Parnaíba basin sediments by its geophysical signature. The joint analysis of free air anomalies, free air admittances, and effective elastic thickness estimates (Te) revealed that the Brazilian East and Equatorial continental margins have quite different elastic properties. In the first one 10 km < Te < 20 km whereas in the second one Te ≤ 10 km. The weakness of the Equatorial margin lithosphere was caused by the cenozoic magmatism. The BP continental margin presents segmentations; some of them have inheritance from precambrian structures and domains. The segmentations conform markedly with some sedimentary basin features which are below described from south to north. The limit between Sergipe and Alagoas subbasins coincides with the suture between BP and SFP. Te estimates indicates concordantly that in Sergipe subbasin Te is around 20 km while Alagoas subbasin has Te around 10 km, thus revealing that the lithosphere in the Sergipe subbasin has a greater rigidity than the lithosphere in the Alagoas subbasin. Additionally inside the crust beneath Sergipe subbasin occurs a very dense body (underplating or crustal heritage?) which is not present in the crust beneath Alagoas subbasin. The continental margin of the Pernambuco basin (15 < Te < 25 km) presents a very distinct free air edge effect displaying two anomalies. This fact indicates the existence in the Pernambuco plateau of a relatively thick crust. In the Paraíba basin the free air edge effect is quite uniform, Te ≈ 15 km, and the lower crust is abnormally dense probably due to its alteration by a magmatic underplating in the Cenozoic. The Potiguar basin segmentation in three parts was corroborated by the Te estimates: in the Potiguar rift Te ≅ 5 km, in the Aracati platform Te ≅ 25 km, and in the Touros platform Te ≅ 10 km. The observed weakness of the lithosphere in the Potiguar rift segment is due to the high heat flux while the relatively high strength of the lithosphere in the Touros platform may be due to the existence of an archaean crust. The Ceará basin, in the region of Mundaú and Icaraí subbasins, presents a quite uniform free air edge effect and Te ranges from 10 to 15 km. The analysis of the Bouguer admittance revealed that isostasy in BP can be explained with an isostatic model where combined surface and buried loadings are present. The estimated ratio of the buried loading relative to the surface loading is equal to 15. In addition, the lower crust in BP is abnormally dense. These affirmations are particularly adequate to the northern portion of BP where adherence of the observed data to the isostatic model is quite good. Using the same above described isostatic model to calculate the coherence function, it was obtained that a single Te estimate for the entire BP must be lower than 60 km; in addition, the BP north portion has Te around 20 km. Using the conventional elastic flexural model to isostasy, an inversion of crust thickness was performed. It was identified two regions in BP where the crust is thickened: one below the Borborema plateau (associated to an uplifting in the Cenozoic) and the other one in the Ceará domain beneath the Santa Quitéria magmatic arc (a residue associated to the Brasiliano orogenesis). On the other hand, along the Cariri-Potiguar trend, the crust is thinned due to an aborted rifting in the Cretaceous. Based on the interpretation of free air anomalies, it was inferred the existence of a large magmatism in the oceanic crust surrounding the BP, in contrast with the incipient magmatism in the continent as shown by surface geology. In BP a quite important positive geoid anomaly exists. This anomaly is spatially correlated with the Borborema plateau and the Macaú-Queimadas volcanic lineament. The integrated interpretation of geoid height anomaly data, global shear velocity model, and geologic data allow to propose that and Edge Driven Convection (EDC) may have caused the Cenozoic magmatism. The EDC is an instability that presumably occurs at the boundary between thick stable lithosphere and oceanic thin lithosphere. In the BP lithosphere, the EDC mechanism would have dragged the cold lithospheric mantle into the hot asthenospheric mantle thus causing a positive density contrast that would have generated the main component of the geoid height anomaly. In addition, the compatibility of the gravity data with the isostatic model, where combined surface and buried loadings are present, together with the temporal correlation between the Cenozoic magmatism and the Borborema plateau uplifting allow to propose that this uplifting would have been caused by the buoyancy effect of a crustal root generated by a magmatic underplating in the Cenozoic

Estudo comparativo de abordagens de inversão do fator de qualidade Q

Seismic wave dispersion and attenuation studies have become an important tool for lithology and fluid discrimination in hydrocarbon reservoirs. The processes associated to attenuation are complex and are encapsulated in a single quantitative description called quality factor (Q). The present dissertation has the objective of comparing different approaches of Q determination and is divided in two parts. Firstly, we made performance and robustness tests of three different approaches for Q determination in the frequency domain. They are: peak shift, centroid shift and spectral ratio. All these tests were performed in a three-layered model. In the suite of tests performed here, we varied the thickness, Q and inclination of the layers for propagation pulses with central frequency of 30, 40 and 60 Hz. We found that the centroid shift method is produces robust results for the entire suíte of tests. Secondly, we inverted for Q values using the peak and centroid shift methods using an sequential grid search algorithm. In this case, centroid shift method also produced more robust results than the peak shift method, despite being of slower convergence

Evidências da herança geotectônica pré-cambriana na geração da barcia potiguar: um estudo geofísico multdisciplinar

The 3D gravity modeling of the Potiguar rift basin consisted of a digital processing of gravity and aeromagnetic data, subsidized by the results of Euler deconvolution of gravity and magnetic data and the interpretation of seismic lines and wells descriptions. The gravity database is a compilation of independent geophysical surveys conducted by several universities, research institutions and governmental agencies. The aeromagnetic data are from the Bacia Potiguar and Plataforma Continental do Nordeste projects, obtained from the Brazilian Petroleum Agency (ANP). The solutions of the Euler Deconvolution allowed the analysis of the behavior of the rift main limits. While the integrated interpretation of seismic lines provided the delimitating horizons of the sedimentary formations and the basement top. The integration of these data allowed a 3D gravity modeling of basement topography, allowing the identification of a series of internal structures of the Potiguar rift, as well intra-basement structures without the gravity effect of the rift. The proposed inversion procedure of the gravity data allowed to identify the main structural features of the Potiguar rift, elongated in the NE-SW direction, and its southern and eastern faulted edges, where the sedimentary infill reachs thicknesses up to 5500 m. The southern boundary is marked by the Apodi and Baixa Grande faults. These faults seem to be a single NW-SE oriented fault with a strong bend to NE-SW direction. In addition, the eastern boundary of the rift is conditioned by the NE-SW trending Carnaubais fault system. It was also observed NW-SE oriented faults, which acted as transfer faults to the extensional efforts during the basin formation. In the central part of the residual anomaly map without the gravity effect of the rift stands out a NW-SE trending gravity high, corresponding to the Orós-Jaguaribe belt lithotypes. We also observe a gravity maximum parallel to the Carnaubais fault system. This anomaly is aligned to the eastern limit of the rift and reflects the contact of different crustal blocks, limited by the eastern ward counterpart of the Portalegre Shear Zone

Avaliação radiográfica da altura mandibular anterior após vestibuloplastia por inversão de retalhos: estudo em humanos

A proposta do presente estudo foi de analisar, através de análise cefalométrica por telerradiografias de perfil, a altura óssea na região mandibular anterior após a realização de vestibuloplastia pela técnica de inversão de retalhos. Foram avaliados 13 pacientes (8 mulheres e 5 homens), com média igual a 54,3 anos que apresentavam, em média, 25,01 mm de altura mandibular anterior. As medidas foram obtidas utilizando-se traçado cefalométrico e comparando-se o período pré-operatório com os períodos pós-operatórios, que consistiram de 30, 90 e 180 dias. Durante o período analisado, houve uma redução da altura mandibular anterior, estatisticamente não significante, de 0,6 mm.

Modelagem gravimétrica da Bacia de Taubaté-vale do Rio Paraíba do Sul, leste do estado de São Paulo

This paper presents the gravity survey refinement of the Paraíba do Sul river valley, east region of São Paulo state, by acquisition of 646 gravity stations. Further than the traditional Free-Air and Bouguer corrections, it was necessary to apply the Terrain Correction due the proximity of the Mantiqueira and Mar mountain ranges. An enhanced understanding of Taubaté Basin structural framework was the main objetive of this work. As result, we present the Bouguer Gravity Map of the east region of São Paulo State and the Residual Bouguer Gravity Map generated by Upward Continuation and High-pass filtering. The 2D Gravity Modeling was applied in 4 sections, using bi-dimensional models for the bodies that produced the anomalies (Talwani Method). After integrating the gravity data with the interpretation of 11 seismic sections, it was possible to delineate the Taubate Basin structural framework, which shows an utmost sedimentary thickness of 800 m in the depocenters. It was also interpreted an ultramaphic alkaline plug in the subsurface of Caçapava city.

Avaliação da eficiência de corte de fresas Carbide em função dos métodos de esterilização físicos e químico. Análise gravimétrica

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Da invenção à inversão do autor: copyleft, all rights reversed

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Informação a priori na inversão de dados magnetotelúricos

Esta Tese apresenta dois estudos aplicados à inversão de dados magnetotelúricos. No primeiro deles, os parâmetros obtidos na inversão são as dimensões da malha de parametrização da subsuperfície, sendo conhecida, a priori, a resistividade de uma heterogeneidade e a do seu meio envolvente; no outro estudo, é abordado o uso de operadores de derivadas de ordem maior do que um com a finalidade de estabilizar o problema inverso. No primeiro estudo, os resultados podem ser considerados satisfatórios somente se a informação sobre as resistividades tem erro menor do que 20%. No segundo estudo, os resultados demonstram que o uso de operadores de ordem maior do que um podem ser mais eficazes do que o uso convencional do operador de primeira derivada, pois além de estabilizarem o problema inverso, esses operadores contribuem para melhorar a resolução das heterogeneidades de resistividade da subsuperfície. Ambos os estudos são inéditos, pois a prática de inversão de dados magnetotelúricos consiste de obter como resultado do problema inverso a resistividade dos prismas de uma malha de parametrização de dimensões fixas, usando como estabilizador o operador de primeira derivada. Os modelos usados nos estudos são bidimensionais e representam uma subsuperfície com uma e duas heterogeneidades de forma prismática envolvidas por ambiente homogêneo. O desempenho das técnicas foi testado com dados sintéticos com e sem ruído gaussiano, bem como dados reais do perfil COPROD2. Durante o trabalho, são, ainda, descritas as técnicas de inversão denominadas creeping e jumping e feita uma comparação e avaliação sobre elas. Mostra-se aqui que, ao contrário do que afirmam muitos pesquisadores, a inclusão de informação a priori sore os parâmetros pode ser feita na técnica do creeping com a mesma facilidade com que é feita na técnica do jumping.

Inversão de dados eletromagnéticos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada

Neste trabalho apresentamos um estudo da aplicação do regularizador “Variação Total” (VT) na inversão de dados geofísicos eletromagnéticos. O regularizador VT reforça a proximidade entre os parâmetros adjacentes, mas, quando a influência de uma descontinuidade é sentida nos dados, este permite mudanças abruptas sobre os parâmetros. Isso faz com que o método seja uma alternativa válida, quando os dados observados usados na inversão provém de um ambiente geológico com uma distribuição suave de condutividade, mas que pode apresentar descontinuidades em lugares como as interfaces entre as camadas geoelétricas, como na margem de uma zona de óleo ou de um corpo de sal, que podem ser zonas muito resistivas no interior de sedimentos condutivos. Quando, devido a baixa resolução nos dados, o método não tem informações o suficiente para identificar a interface, o regularizador variação total reforça a proximidade entre os parâmetros adjacentes fazendo um transição suave entre as condutividades camadas, da mesma forma que é apresentado pela suavidade global. O método de Variação Total permite que modelos menos suaves sejam alcançados porque na norma L₁ a medida de desajuste entre os pares de parâmetros adjacentes, dará o mesmo valor se a variação dos parâmetros é suave ou se a variação é abrupta, o que não é o caso se o mesmo desajuste é medido na norma L₂, pois em uma distribuição suave a medida do desajuste é menor, sendo assim favorecida pela minimização desta norma. O uso deste regularizador permite uma melhor estimativa do tamanho de um corpo, seja ele resistivo ou condutivo. O trabalho está apresentado na forma de três artigos, cada um descrevendo uma etapa no desenvolvimento do problema da inversão, seguindo uma sequência de complexidade crescente no problema direto. O primeiro artigo neste trabalho é intitulado “Inversão de dados do CSEM marinho 1D de meio estratificado anisotrópico com o regularizador Variação Total”. Este descreve o passo inicial no desenvolvimento do problema: a inversão de dados do CSEM marinho de modelos estratificados 1D com anisotropia na condutividade das camadas. Este problema se presta bem para este desenvolvimento, porque tem solução computacional muito mais rápida do que o 2D, e nele já estão presentes as características principais dos dados do método CSEM marinho, como a largura muito grande da faixa de amplitudes medidas em um levantamento, e a baixa resolução, inerente às baixas frequências empregadas. A anisotropia acrescenta uma dificuldade a mais no problema, por aumentar o nível de ambiguidade nos dados e demandar ainda mais informação do que no caso puramente isotrópico. Os resultados mostram que a aplicação dos vínculos de igualdade do método VT permite a melhor identificação de uma camada alvo resistiva do que a simples aplicação dos vínculos tradicionais de suavidade. Até onde podemos aferir, esta solução se mostra superior a qualquer outra já publicada para este problema. Além de ter sido muito importante para o desenvolvimento de códigos em paralelo. O segundo artigo apresentado aqui, “Inversão de dados Magnetotelúricos com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, trata da inversão de dados do método Magnetotelúrico em ambientes 2D. Este problema demanda um esforço computacional muito maior do que o primeiro. Nele, estudamos a aplicação do método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão. A construção da matriz de sensibilidade é a etapa que demanda mais tempo no processo de inversão, e o uso do método de estados adjuntos foi capaz de reduzir muito este tempo, gerando derivadas com um bom nível de aproximação. Esta etapa da pesquisa foi fundamental pelo problema direto ser matematicamente e computacionalmente muito mais simples do que o do CSEM marinho 2D. Novamente em comparação com a aplicação do regularizador de suavidade global, o regularizador de Variação Total permitiu, neste problema, uma melhor delimitação das bordas de heterogeneidades bidimensionais. A terceira parte deste trabalho, apresentada no artigo “Inversão de dados do CSEM marinho 2.5D com o regularizador Variação Total e o uso da matriz de sensibilidade aproximada”, apresenta a apliação do método de Variação Total ao problema da inversão de dados CSEM marinho 2.5D. Usamos o método dos estados adjuntos para gerar uma boa aproximação para as derivadas necessárias para a construção da matriz de sensibilidade usada na inversão, acelerando assim o processo de inversão. Para deixar o processo de inversão ainda mais rápido, lançamos mão da programação em paralelo com o uso de topologia. A comparação entre a aplicação do regularizador de suavidade global, e o regularizador de Variação Total permitiu, assim como nos casos anteriores, uma melhor delimitação das bordas de heterogeneidades bidimensionais.

Inversão de dados de ondas de superfície Rayleigh

Neste trabalho compilamos informações sobre um grande número de medidas de velocidade de grupo para ondas Rayleigh do modo fundamental, com período até 100 segundos. Tais dados consistiram de informações retiradas da literatura geofísica e cobriram toda a Terra. Parte dos dados foi organizada em trabalhos anteriores e uma segunda parte foi apresentada aqui de forma inédita. Para a América do Sul, selecionamos os principais conjuntos de dados de tais ondas e elaboramos diversos perfis onde a distribuição de velocidade de ondas cisalhantes foi obtida a partir da inversão das curvas de dispersão de velocidade de grupo. Tais perfis serviram para termos uma ideia inicial da estrutura interna da Terra em nosso continente. Com o conjunto global de dados de velocidade de grupo foi possível obtermos os mapas de distribuição lateral de valores de velocidade para cada período referencial entre 20 e 100 segundos. Tais mapas foram produzidos da mesma forma que os mapas de velocidade de fase de ROSA (1986), onde a amostragem for para realizada para blocas medindo 10x10 graus, englobando toda a Terra, em projeção mercator. O valor de velocidade de grupo em cada bloco, para cada período, foi obtido a partir da inversão estocástica dos dados de anomalia de velocidade em relação aos modelos regionalizados de JORDAN (1981) com os valores de velocidade de grupo de ROSA et al. (1992). Os mapas de velocidade de grupo obtidos aqui foram então empregados, na América do Sul, com os valores de velocidade de fase dos mapas obtidos por ROSA (1986). Assim, foi possível determinarmos, em profundidade, os mapas de variação de velocidade de onda cisalhante e os mapas de distribuição de valores de densidade. Com isto, pudemos construir o primeiro mapa de profundidade do Moho (todo do Manto Superior) da América do Sul.

Inversão vinculada de dados de EMAP e MT-1D anisotrópico

Uma técnica para a inversão de dados magnetotelúricos é apresentada neste trabalho. Dois tipos de dados são tratados aqui, dados gerados por modelos unidimensionais com anisotropia na condutividade das camadas e dados bi-dimensionais de levantamentos do método EMAP (ElectroMagnetic Array Profiling). Em ambos os casos fazemos a inversão usando vínculos aproximados de igualdade para estabilizar as soluções. Mostramos as vantagens e as limitações do uso destes vínculos nos processos de inversão. Mesmo vinculada a inversão ainda pode se tornar instável. Para inverter os dados 2-D do EMAP, apresentamos um processo que consiste de três partes: 1 – A construção de um modelo interpretativo e da aproximação inicial para a inversão a partir dos dados de seções de resistividade aparente filtradas pelo processo de filtragem do EMAP; 2 – a inclusão de uma camada de corpos pequenos aflorantes, chamada de camada destatic shift, aos modelos interpretativos para resolver as fontes de distorções estáticas que contaminam os dados; 3 – o uso dos vínculos aproximados de igualdade absoluta para estabilizar as soluções. Os dois primeiros passos nos permitem extrair o máximo de informação possível dos dados, enquanto que o uso dos vínculos de igualdade nos permite incluir informação a priori que possua significado físico e geológico. Com estes passos, obtemos uma solução estável e significativa. Estudaremos o método em dados sintéticos de modelos bi-dimensionais e em dados reais de uma linha de EMAP feita na Bacia do Paraná.

Modelagem sísmica e inversão na presença de anisotropia

A necessidade da adoção de modelos elásticos anisotrópicos, no contexto da sísmica de exploração, vem crescendo com o advento de novas técnicas de aquisição de dados como VSP, walkway VSP, tomografia poço a poço e levantamentos sísmicos com grande afastamento. Meios anisotrópicos, no contexto da sísmica de exploração, são modelos efetivos para explicar a propagação de ondas através de meios que apresentam padrões de heterogeneidade em escala muito menor que o comprimento de onda das ondas sísmicas. Particularmente, estes modelos são muito úteis para explicar o dado sísmico mais robusto que são as medidas de tempo de trânsito. Neste trabalho, são investigados aspectos da propagação de ondas, traçado de raios e inversão de tempos de trânsito em meios anisotrópicos. É estudada a propagação de ondas SH em meios anisotrópicos estratificados na situação mais geral onde estas ondas podem ocorrer, ou seja, em meios monoclínicos com um plano vertical de simetria especular. É mostrado que o campo de ondas SH refletido a partir de um semi-espaço estratificado, não apresenta qualquer informação sobre a possível presença de anisotropia em subsuperfície. São apresentados métodos simples e eficientes para o traçado de raios em 3D através de meios anisotrópicos estratificados, baseados no princípio de Fermat. Estes métodos constituem o primeiro passo para o desenvolvimento de algoritmos de inversão de tempos de trânsito para meios anisotrópicos em 3D, a partir de dados de VSP e walkaway VSP. Esta abordagem é promissora para determinação de modelos de velocidade, que são necessários para migração de dados sísmicos 3D na presença de anisotropia. É efetuada a análise da inversão tomográfica não linear, para meios estratificados transversalmente isotrópicos com um eixo de simetria vertical(TIV). As limitações dos dados de tempo de trânsito de eventos qP para determinação das constantes elásticas, são estabelecidas e caracterizados os efeitos da falta de cobertura angular completa na inversão tomográfica. Um algoritmo de inversão foi desenvolvido e avaliado em dados sintéticos. A aplicação do algoritmo a dados reais demonstra a consistência de meios TIV. Esta abordagem é útil para casos onde há informação a priori sobre a estratificação quase plana das formações e onde os próprios dados do levantamento poço a poço apresentam um alto grau de simetria especular em relação a um plano vertical. Também pode ser útil em interpretações preliminares, onde a estimativa de um meio estratificado, serve como modelo de fundo para se efetuar análises mais detalhadas, por exemplo, como um modelo de velocidades anisotrópico para migração, ou como um modelo de calibração para análises de AVO.