RCAAP Repository

Com base em sísmica de alta resolução (perfilador de 7,0kHz e sonar de varredura lateral 100kHz), e amostras de sedimentos, foram investigados aspectos interessantes do fundo marinho na baía da Ilha Grande. Podem-se evidenciar efeitos combinados de ondas e correntes agindo sobre o fundo marinho, a partir dos índices de ripple (RI) e de paralelismo (PI2) de estruturas sedimentares (formas de fundo) observadas nos registros de sonar. Estruturas estas associadas a padrões sedimentares de fundo (ecofácies). Tais estruturas são predominantemente "mega ripples" localizadas em profundidades em torno de 16m, próximo à Enseada de Conceição de Jacareí onde ocorre uma transição textural de fundo de leste para oeste de um padrão de fundo grosso para um padrão fino. Suas cristas são orientadas a NE-SW e a ESE-WNW e bastante variáveis na forma (retas, onduladas, bifurcadas, descontínuas). Suas alturas variam de 0,2 a 0,55m, com valor médio em torno de 0.35 m. O comprimento de onda varia de 10 a 60m. Velocidades estimadas - unidirecionais para correntes de fundo e orbitais para ação de ondas - com base em parâmetros geológicos das estruturas observadas, sugerem uma possível gradação na magnitude das velocidades (35-80cm/s) na região dos padrões sedimentares de transição de fundo. Assim, as feições observadas (faciológicas e morfodinâmicas) do fundo marinho indicam que ondas e correntes ao interagirem com a geometria da baía podem provocar oscilações no nível de energia do meio, atuando no grau de selecionamento dos sedimentos e caracterizando diferentes tipos de fundo.

As imagens hiperespectrais são muito suscetíveis às interferências de ruídos, o que dificulta a identificação e a quantificação mineral. Portanto, torna-se fundamental o emprego de técnicas para minimizá-los. Um eficiente procedimento matemático para redução de ruídos é a Fração de Ruído Mínima (MNF). Esse método é derivado e análogo à transformação Análise de Componentes Principais, distinguindo-se por promover a maximização do ruído ao invés da variância dos dados, o que proporciona um ordenamento que reflete a qualidade das imagens. A transformação MNF pode ser subdividida em quatro etapas: a) obtenção de uma amostra do ruído; b) formulação de um índice de fracionamento do ruído; c) implementação de uma função de transformação linear no molde da APC, e d) inversão do MNF considerando apenas as informações de sinais. No presente trabalho é proposta uma metodologia de MNF seqüencial, utilizando referências internas e externas da imagem. Como referência externa foi utilizada a referência escura coletada durante o vôo simultaneamente à obtenção da imagem. A referência interna foi obtida por técnicas estatísticas para a segmentação da fração do ruído em relação ao sinal. Três tipos de ruídos podem ser definidos conforme a variância e correlação: a) ruídos não-correlacionados e com igual variância; b) ruídos altamente correlacionados; e c) ruídos com matriz de covariância desconhecida. Para as imagens AVIRIS pode-se esquematizar um procedimento seqüencial de MNF com duas etapas, sendo que, na primeira, são eliminados os ruídos de referência interna e na segunda, os de referência externa relativos ao ruído instrumental. Esse procedimento permite distinguir diferentes tipos de ruídos e eliminá-los seqüencialmente, de forma simples e eficiente.

A migração com amplitudes verdadeiras de dados de reflexão sísmica, em profundidade ou em tempo, possibilita que seja obtida uma medida dos coeficientes de reflexão dos chamados eventos de reflexão primária. Estes eventos são constituídos, por exemplo, pelas reflexões de ondas longitudinais P-P em refletores de curvaturas arbitrárias e suaves. Um dos métodos mais conhecido é o chamado migração de Kirchhoff, através do qual a imagem sísmica é produzida pela integração do campo de ondas sísmicas, utilizando-se superfícies de difrações, denominadas de Superfícies de Huygens. A fim de se obter uma estimativa dos coeficientes de reflexão durante a migração, isto é a correção do efeito do espalhamento geométrico, utiliza-se uma função peso no operador integral de migração. A obtenção desta função peso é feita pela solução assintótica da integral em pontos estacionários. Tanto no cálculo dos tempos de trânsito como na determinação da função peso, necessita-se do traçamento de raios, o que torna a migração em situações de forte heterogeneidade da propriedade física um processo com alto custo computacional. Neste trabalho é apresentado um algoritmo de migração em profundidade com amplitudes verdadeiras, para o caso em que se tem uma fonte sísmica pontual, sendo o modelo de velocidades em subsuperfície representado por uma função que varia em duas dimensões, e constante na terceira dimensão. Esta situação, conhecida como modelo dois-e-meio dimensional (2,5-D), possui características típicas de muitas situações de interesse na exploração do petróleo, como é o caso da aquisição de dados sísmicos 2-D com receptores ao longo de uma linha sísmica e fonte sísmica 3-D. Em particular, é dada ênfase ao caso em que a velocidade de propagação da onda sísmica varia linearmente com a profundidade. Outro tópico de grande importância abordado nesse trabalho diz respeito ao método de inversão sísmica denominado empilhamento duplo de difrações. Através do quociente de dois empilhamentos com pesos apropriados, pode-se determinar propriedades físicas e parâmetros geométricos relacionados com a trajetória do raio refletido, os quais podem ser utilizados a posteriori no processamento dos dados sísmicos, visando por exemplo, a análise de amplitudes.

O processamento convencional, sem um tratamento adequado de preservação de amplitude, mascara ou até inviabiliza o mapeamento de propriedades petrofísicas, prejudicando a correlação entre o dado sísmico e a informação de perfil. Uma forma de estimar corretamente as amplitudes e, por conseqüência, os coeficientes de reflexão nas interfaces sísmicas, é efetuar uma migração pré-empilhamento em verdadeira amplitude, onde a distorção das amplitudes devido ao espalhamento geométrico ao longo do raio de reflexão é compensado pela operação de migração. Entretanto, esse processo tem um custo elevado, sendo também dependente do modelo de velocidades utilizado. Uma rotina menos dispendiosa e mais estável em relação a erros no modelo de velocidade, consiste na transformação da seção sísmica obtida com a geometria de aquisição no afastamento comum em seções simuladas de afastamento nulo em amplitude verdadeira. A operação de transformação de uma seção sísmica de afastamento constante para a sua correspondente em afastamento nulo recebe o nome de migração para afastamento nulo em verdadeira amplitude (TA MZO). No caso de um meio com velocidade constante, a curva de empilhamento do MZO e a função peso se reduzem a fórmulas analíticas, minimizando o esforço computacional. O trabalho em questão tem dois objetivos principais: o primeiro é verificar a eficácia do algoritmo de TA MZO para velocidade constante em um modelo sintético com geologia complexa baseado em um reservatório turbidítico brasileiro de idade Neo-Albiana, onde a hipótese de velocidade constante é desrespeitada. O segundo objetivo é a realização de estudos quantitativos referentes à obtenção de amplitudes resultantes da aplicação da rotina descrita acima, bem como avaliar sua utilização para a construção de gráficos de amplitude com o afastamento (AVO) e com o ângulo (AVA) em interfaces relevantes para a caracterização de reservatórios.

Um levantamento experimental GPR (Ground Penetrating Radar) foi realizado no Rio Taquari, na Bacia do Pantanal Matrogrossense, Brasil, visando entender os processos de assoreamento na base do Rio e mapear os canais ao longo do Rio, que são essenciais para a navegação local. Muitos perfis GPR foram realizados no baixo Taquari, próximos à Fazenda Santa Maria e à Fazenda Porto Santa Luzia. Interpretações preliminares foram interessantes, permitindo localizar os canais do Rio e a ocorrência de grandes depósitos de areia, mostrando a viabilidade das aplicações do GPR em estudos de processos de sedimentação no Rio Taquari. Todos os dados adquiridos estão sendo analisados visando entender a dinâmica local nos processos de sedimentação, como parte de uma dissertação de mestrado em andamento por um estudante no IAG-USP. Os resultados deste levantamento geofísico servem com um guia para os estudos multidisciplinares do Programa Pantanal do governo do Estado do Mato Grosso do Sul.

Investigações geofísicas foram realizadas sobre a Suíte Granítica de Itu, localizada próxima à cidade de Itu, Estado de São Paulo-Brasil, com o objetivo de localizar zonas de fraturas nas rochas graníticas e determinar o topo da rocha sã, através da integração dos métodos GPR-Ground Penetrating Radar e Resistividade Elétrica. Os perfis geofísicos foram realizados sobre uma mesma linha para estudos comparativos entre estas duas metodologias, com a filosofia de interpretação integrada. Um perfil GPR de 80m foi adquirido com as antenas de 50, 100 e 200 MHz, além de dois perfis de caminhamento elétrico com dipólos de 2 e 10m. As análises dos perfis geofísicos permitiram identificar dois fortes refletores inclinados: um em torno de 10m de profundidade, interpretado como fraturas preenchidas com água e o outro entre 12 e 17m de profundidade, que corresponde a uma região resistiva mergulhando depois 15m de profundidade, interpretada como o topo da rocha granítica sã. Além disso, duas regiões anômalas foram identificadas: uma em torno de 50m e a outra em torno de 80m. Nos perfis GPR, essas regiões são caracterizadas por uma região sem refletores, devido à elevada atenuação das ondas eletromagnéticas. Nos perfis elétricos, essas regiões correspondem às regiões condutivas que podem estar relacionadas com a presença de uma fratura subvertical.

This paper illustrates a new reservoir characterization approach using seismic pattern recognition methodology based on principal component analysis, trace shape and 3D multi-attributes classification on a set of 3D seismic data volumes. We develop the reservoir interpretation workflow (<A HREF="#fig01">figure 1</A>) as follows: 1) Seismic facies maps obtained by Kohonen's Self Organizing Map (SOM) Neural Network method to seismic facies volume from a set of various seismic attributes, obtained by Hierarchical classification method; and 2) Statistical analysis of available seismic attributes by application of Principal Component Analysis (PCA) statistical method before classification. How to discriminate the amount of information that each geophysicist is dealing with, in a reservoir characterization process and what is the suitable procedure to help to discriminate between data and information? We present some results of its application on the Paleocene/Eocene reservoirs of a Campos Basin field.

A análise da viabilidade técnica para monitoramento do fluxo de fluidos em reservatórios, através de levantamentos sísmicos do tipo time-lapse, ou 4D, deve envolver não apenas a avaliação laboratorial das modificações provocadas no reservatório pela substituição de fluidos com a conseqüente variação dos parâmetros elásticos mas, principalmente, uma fase prévia de modelagem sísmica cuja importância pode ser, dentre outras, quantificar o impacto do efeito tuning que pode afetar as esperadas variações de amplitude após o time-lapse. Nesse caso, tal efeito precisa ser investigado com profundidade, pois devido à substituição de fluídos ao nível dos reservatórios, a variação do tuning observado em camadas delgadas pode introduzir um fator-modificação das amplitudes sísmicas que inviabilize a correta interpretação da sísmica 4D, especialmente quando se procura fazer uma interpretação mais quantitativa com a análise do comportamento de AVO e outras técnicas que exigem a investigação da variação de amplitudes.

An Eulerian model and a Lagrangian particle model are used to study the dispersion of a contaminant released from a low height source in the Stable Boundary Layer (SBL). The Eulerian model is based on the solution of the advection-diffusion equation by the Laplace transform technique. The Lagrangian model is based on a generalized form of the Langevin equation. Turbulence inputs are parameterised according to two procedures capable of generating continuous values in all stability conditions and in all heights of Planetary Boundary Layer (PBL). Statistical indices were calculated to compare the predicted and observed values of ground-level concentration. According to the statistical analysis, predicted concentration values agree well with the observed ones.

Accounting for the current knowledge of the Convective Boundary Layer (CBL) structure and characteristics, a new formulation for eddy diffusivities has been derived, to be used in atmospheric dispersion models. Expressions are proposed for eddy diffusivities depending on source distance for inhomogeneous turbulence calculated directly with the Batchelor theory. The classical statistical diffusion theory, the observed spectral properties and the observed characteristics of energy-containing eddies are used to estimate these parameters. In addition, a vertical eddy diffusivity was introduced in an air pollution model and validated with the data of Copenhagen experiments and compared with the results of the obtained from classic formulation of the literature. A statistical evaluation show that, the results of the dispersion model with the K-parameterisation included, produces a good fit of the measured ground-level concentration data in the experimental conditions considered. Furthermore, the current work suggests that wind and eddy diffusivities profiles physically more realistic are important in regions near of an elevated source, and improves the description of the transport process of atmospheric pollutants.

Os métodos GPR-Ground Penetrating Radar e Resistividade Elétrica foram utilizados com o objetivo de localizar zonas de fraturas e o topo da rocha sã, numa pedreira de granito ornamental, localizada próxima à cidade de Capão Bonito, Estado de São Paulo-Brasil. Oito perfis de reflexão GPR e 19 sondagens de velocidade WARR-Wide Angle Reflection and Refraction foram adquiridas utilizando-se as antenas de 25, 50, 100 e 200 MHz. Também foram adquiridas duas SEV's-Sondagens Elétricas Verticais. As análises dos perfis GPR permitiram localizar fortes refletores inclinados e subhorizontais, provavelmente relacionados com zonas de fraturas preenchidas com água. Além disso, foram identificadas duas zonas sem reflexões rasas que podem estar relacionadas com intrusões de diques de diabásio, e um refletor hiperbólico relacionado com um matacão. A rocha sã foi identificada a 23 m de profundidade. Os resultados obtidos irão auxiliar o projeto da engenharia de minas orientando o avanço na frente de trabalho.

Atualmente existem várias técnicas geofísicas de mapeamento de fraturas como GPR, e caminhamentos eletromagnético e geoelétrico dipolo-dipolo. Esses métodos não fornecem informações acerca do conteúdo da zona de fratura, que pode estar preenchida de argila. Usando métodos nucleares, tais como, espectrometria gama e emanomeria de radônio, torna-se possível identificar a natureza da composição mineral do material da fratura e estudar a circulação de água usando o transporte de radônio pela água de percolação. Neste trabalho, a integração das técnicas nuclear, GPR e geoelétrica nas rochas intrusivas graníticas de Itu mostrou-se muito interessante porque torna possível encontrar e identificar fraturas em rocha granítica preenchidas com argila ou contendo água doce.

Os equipamentos foram desenvolvidos para efetuarem medidas de potenciais espontâneo e induzido na superfície. Em ambos os casos foi usado um cabo multifio com entradas múltiplas de tensão ao longo do comprimento, que tornou possível conectá-lo a 6 eletrodos de aço inox cravados no solo. Um sistema digital de seqüenciamento automático controlado por mini interruptor de pressão permite acoplar os eletrodos de potencial seqüencialmente para um voltímetro eletrônico DC com circuito de autozeramento baseado em circuito de servocontrole PID. Um filtro ativo tipo rejeição de 60 Hz e um tipo passa-baixa com constante de tempo de 1 segundo são usados para minimizarem os ruídos de linha e os esféricos. O sistema de controle eletrônico torna também possível acoplar o cabo multifio com o receptor datalogger do SYSCAL R-2 no arranjo para caminhamento dipolo- dipolo.

O presente trabalho tem por objetivo a utilização da classificação não-supervisionada de dados geofísicos aéreos (magnéticos e de espectrometria de raios gama) e gravimétricos em Cuba Centro-Oriental, com o intuito de gerar subsídios para interpretação geológica regional e tectono-estrutural na área de estudo. A região centro-oriental de Cuba abrange amplos setores das principais unidades geológicas documentadas no arquipélago, incluindo o Arco Vulcânico Cretácico (AVC), os ofiolitos Mesozóicos, a margem continental da plataforma das Bahamas (Jurássico Superior-Cretáceo Superior) e a cobertura sedimentar Cenozóica. As características geológicas da área a tornam um interessante `laboratório' para utilização de dados geofísicos e para avaliação de métodos estatísticos multivariados. O estudo revelou que as classes geofísicas, derivadas da aplicação desses métodos aos dados, possuem uma boa correlação com as formações geológicas e ambientes tectono-estruturais da área.

We investigate the estimation of fractures orientation, strike and dip, through multiazimuthal AVO analysis of qP and its converted waves qS1 and qS2. We assume weak impedance contrast, weak anisotropy and that the fractured medium behaves as an effective transversally isotropic (TI) medium. Under these assumptions, the estimation of fractures orientation is reduced to the estimation of the orientation of an axis of symmetry from qP reflectivity data. Linearized approximations of qP reflectivity are used for inversion.

Two-point raytracing problem is solved for events in a piecewise homogeneous and laterally varying 3D anisotropic media by continuation techniques. In conjunction with the shooting method the algorithm can be used for computation of qP, qS1, and qS2 events. The algorithm has the same performance and robustness as previous implementations of the continuation method for tracing rays in isotropic models. Routines based on our algorithm have several useful applications. First, an efficient forward problem solver for traveltime inversion of elastic parameters in the presence of anisotropy. Second, Newton-Raphson iterations during two-point raytracing produce wavefront attributes, slowness and wavefront curvature. These attributes allows the computation of geometrical spreading and second order approximations for traveltimes. Therefore it can be used to investigate the effects of anisotropy on CRS, in simple velocity models.

Uma das maiores preocupações em projetos geoambientais é o estudo da contaminação dos solos e da água subterrânea causada pela disposição de lixo urbano. Os métodos geofísicos aplicados são um bom caminho para pegar algumas informações requeridas num programa de caracterização. Este trabalho mostra alguns resultados obtidos através de levantamentos de resistividade, eletromagnético indutivo e GPR (Ground Penetrating Radar), na busca de informações sobre o impacto causado pela área de disposição do lixo em Cuiabá, Estado do Mato Grosso, Brasil.

Este artigo mostra resultados parciais de investigações geofísicas realizadas na borda da Bacia Sedimentar de São Paulo, dentro do campus da Universidade de São Paulo. O principal objetivo desta pesquisa é a caracterização dos sedimentos e do topo das rochas do embasamento da bacia, utilizando os métodos GPR-Ground Penetrating Radar e Resistividade Elétrica. Os perfis geofísicos foram realizados ao longo de uma mesma linha, visando comparar os resultados com as duas metodologias. Um perfil de 216m foi adquirido com as antenas de 50 MHz e um perfil de resistividade com dipólos de 20m. Além disso, duas Sondagens Elétricas Verticais - SEV's foram feitas sobre a mesma linha. A análise do perfil GPR permitiu identificar o limite entre a base do solo e uma camada de areia. Os resultados estão de acordo com as informações de poços. O perfil elétrico identificou uma zona anômala de condutividade entre 80 e 110 m, relacionada provavelmente com uma zona de cizalhamento.

Este trabalho visa o estudo de contaminação causada por um lixão, através dos métodos geofísicos da resistividade, eletromagnético indutivo. Os trabalhos foram realizados em Ribeirão Preto, São Paulo, Brasil, em um lixão desativado. Os dados obtidos foram tratados com softwares apropriados que possibilitaram o realce das anomalias, as quais auxiliaram na caracterização do problema ambiental causado.

Os métodos numéricos são comumente utilizados na simulação de casos difíceis da geofísica de exploração, e, para esta razão, é sempre necessário desenvolver algoritmos computacionais rápidos e eficientes para atingir esses objetivos. Neste artigo, trabalhamos num algoritmo computacional unidimensional com a finalidade de simular perfis condutivos em poço, que foram coletados em ambientes geológicos complexos. O programa foi usado para analisar efeitos nos perfis de indução, escolhendo como modelo diversas seqüências intercaladas de camadas paralelas horizontais finas (folhelho-areia), com alto contraste de condutividades, o que caracteriza uma situação geológica muito comum nos reservatórios, e, sempre, de acordo com a indústria de petróleo, apresenta dificuldades de interpretação. Dessa forma, para realizar nossa modelagem, utilizamos as ferramentas ILD, 2C40 e 6FF40, desenvolvidas pela SCHLUMBERGER e muito usadas na perfilagem de poços. Nossos resultados mostram que as ferramentas focalizadas do tipo 6FF40 possuem melhores respostas nos casos analisados, mas, na maioria das sondas, os valores de condutividades e espessuras das camadas ficam muito longe dos modelos teóricos.