Alterações hidrogeomorfológicas devido à dinâmica de uso da terra na Bacia Hidrográfica do Ribeirão Santa Gertrudes (SP)

Agricultural and mining activities are directly related to changes in natural landscapes. From the perspective of anthropogenic geomorphology supported by general systems theory, this research has developed, in order to identify and analyze changes in the land use, mainly from agriculture and the exploitation of clay, and its implications in the hydro-geomorphological characteristics in the Santa Gertrudes Stream watershed (SP). This area is within the context of the Ceramic Pole Santa Gertrudes (SP), which besides its importance as a supplier of raw material, is characterized as the largest center of international reference in ceramic tiles on the American continent. For this purpose, we made land use and geomorphology maps of two scenarios, corresponding to the years 1962 and 2006, which allowed the identification of changes caused by human activities on the landforms of the area, such as the increase in area of parcels intended for mining activity, which went from 3.1% to occupy 19% of the catchment area of the respective period and that, among other changes, gave rise to new forms of relief as, for example, in levels of pit mining abrupt and smooth. The results indicate that the main features of representative of human changes in relief are represented for opening of large clay mining pits and agricultural activities, which intensified the denudation and sedimentation processes in the Santa Gertrudes Stream watershed.

Avaliação de riscos ambientais no Aeroporto Internacional de Viracopos, Campinas (SP)

Considering the current context of environmental concern, the main purpose of this research was to develop and apply indicators to evaluate the environmental risk for airports. As a result, 19 indicators related to storage and handling dangerous substances were created. These indicators represent five conditions varying from the worst (level 1) to the ideal condition (level 5). The indicators were applied, as a study of case, at Viracopos International Airport, Campinas (SP). The results showed a medium performance of 2.3, using a scale raging from 1.0 to 5.0, characterized as 'critical'. Then, the indicators were analyzed to determine the major weaknesses related to the environmental management at Viracopos International Airport, as well as to identify good practices and solutions. Among the suggested management strategies, we highlight the suitability of areas for storage of dangerous substances. In general, it can be inferred that the indicators proposed in this paper to evaluate the performance in the management of environmental risks at the Viracopos International Airport and propose actions to prevent, monitor, correct, and mitigate the environmental impacts of its operations.

Geotermometria das águas hipertermais do Sistema Aquífero Guarani no estado de São Paulo

Several geothermometers have been employed to determine the expected temperature in the groundwater reservoirs, most part of them based on the chemical equilibrium between thermal groundwaters and minerals which constitute the aquifer framework. Based on this principle, groundwater from Guarani Aquifer System, which presents temperatures upon to 38°C, classified as hyperthermal by the Brazilian Mineral Water Code, had their chemical composition evaluated with the main purpose to determine the expected temperatures in the reservoir based on the application of several geothermometers. These groundwaters can be classified as Na-HCO3 type evolving to Na-HCO3-Cl-SO4 type due the increase in sulphate and chloride concentrations. Measured groundwater temperatures show a good correlation with the depth of the aquifer, indicating an average geothermal gradient of 27.7°C/km. The lowest geothermal gradients, bellow 25°C/km, are observed in the central portion of the study area; meanwhile the highest gradients are observed in NE and W portions of the study area. Quartz equilibrium based geothermometers return temperature over the measured temperatures, meanwhile chalcedony geothermometer returns temperatures near the measured in the wells, except for the wells located in the area which presents the lowest geothermal gradients. The geothermometers based on the equilibrium between Na/K have returned temperatures diverse from the measured at the field, meanwhile the geothermometer based on the equilibrium between Na/K/Ca has presented good match with groundwater which present temperatures under 50°C, due the thermodynamic equilibrium between these waters and feldspars. Geothermometers show evidences about mixtures between groundwaters from GAS and underlying units, especially in wells located in the SW portion of the study area.

Compartimentação geomorfológica e processos deposicionais no megaleque fluvial do rio São Lourenço, Pantanal mato-grossense

The Pantanal wetland is an extensive depositional tract characterized by the presence of fluvial megafans, from which that of the Taquari River is the most extensive and well known. Located on the northern border of the Pantanal, the São Lourenço megafan is a large distributary fluvial system with an area of 16,000 km2, but poorly known regarding its geomorphological and geological features. The Sao Lourenço and the Taquari are coalescent megafans, having the Piquiri interfan meander river in between them. Based on the interpretation of satellite images and field validation, it was possible to identify channel patterns, to interpret depositional and erosional processes, to recognize the phenomena of river avulsion and to map depositional lobes. Three geomorphological zones were recognized on the São Lourenço megafan: 1) abandoned depositional lobes located in the upper/intermediate fan, composed of Pleistocene fluvial deposits and exhibiting distributary paleochannels on their surfaces, which are currently being dissected by tributary channels; 2) an active Holocene confined meander belt formed by fluvial aggradation in a N65E incised-valley; 3) active depositional lobes placed on the distal part of the system, which are the main site of sedimentation and are characterized by frequent avulsion events, channel bifurcation, distributary drainage pattern and the presence of depositional lobated landforms. The active depositional lobes were formed due to an important event of river avulsion that caused the lower portion of the meander belt to be abandoned.

Avaliação de áreas deposicionais e erosivas em cabos lamosos da zona costeira Amazônica através da análise multitemporal de imagens de sensores remotos

Os cabos lamosos Cassiporé e Orange no norte do Brasil constituem um ambiente dinâmico influenciado pelo rio Amazonas, onde as modificações na linha de costa estão sujeitas a severos processos de progradação e erosão. Imagens de sensores remotos ópticos e microondas foram coletadas de 1980 a 2003 e analisadas em um Sistema de Informação Geográfica (SIG), permitindo a identificação e quantificação da distribuição espacial das áreas de progradação e retrogradação ao longo da linha de costa. Durante este período, as maiores taxas de erosão ocorreram junto ao cabo Cassiporé, com recuo médio de 27,5 metros de distância linear e erosão de 1,37 km² de área de manguezal por ano. Por outro lado, os maiores índices de deposição de sedimentos ocorreram no cabo Orange, onde a planície costeira progradou 24,6 m ao ano, agregando 55,85 km² de manguezal à linha de costa nos últimos vinte e três anos. Os mecanismos de progradação determinaram um acréscimo na vegetação de manguezal de 50,8% ao longo das três últimas décadas. Um balanço sedimentar realizado na área pesquisada demonstrou que predominam os processos construtivos (61,3%) sob os processos erosivos (38,7%).

Estratigrafia da sucessão sedimentar Pós-Barreiras (Zona Bragantina, Pará) com base em radar de penetração no solo

Investigação pioneira aplicando Radar de Penetração no Solo (GPR) na área da Praia do Atalaia, norte do Brasil, nos permitiu caracterizar, pela primeira vez, fácies e estratigrafia dos depósitos conhecidos informalmente como Sedimentos Pós-Barreiras (Plioceno e mais recente). Esta sucessão sedimentar recobre discordantemente o embasamento miocênico, representado pelas formações Pirabas/Barreiras. Três unidades estratigráficas foram reconhecidas. A Unidade 1, inferior, consiste em um intervalo com 6 m de espessura média dominado por reflexões pobremente definidas e de baixa amplitude, e que intergradam com reflexões de média escala dos tipos tangencial, obliquo e hummocky. A Unidade 2, intermediária, possui cerca de 9 m de espessura e inclui, principalmente, reflexões oblíquas de larga escala, cujas configurações variam de paralela, tangencial, sigmoidal a sigmoidal-complexa. A Unidade 3, superior, corresponde a um intervalo entre 3,5 e 9 m de espessura, sendo dominada por reflexões hummocky, seguidas por reflexões de média escala dos tipos oblíquo, paralelo a sub-paralelo, e em corte-e-preenchimento. A análise da configuração das reflexões internas e geometria das reflexões nos leva a propor que a unidade correspondente aos Sedimentos Pós-Barreiras é mais variável faciologicamente que inicialmente imaginado, incluindo depósitos eólicos (dunas costeiras), bem como depósitos de cordão litorâneo, planície de maré, canal e mangue. Além disto, o mapeamento das três unidades descritas acima é importante para desvendar a complexidade de sedimentação versus erosão durante o Neógeno tardio no norte do Brasil.

2-D ZO CRS stack by considering an acquisition line with smooth topography

Os dados sísmicos terrestres são afetados pela existência de irregularidades na superfície de medição, e.g. a topografia. Neste sentido, para obter uma imagem sísmica de alta resolução, faz-se necessário corrigir estas irregularidades usando técnicas de processamento sísmico, e.g. correições estáticas residuais e de campo. O método de empilhamento Superfície de Reflexão Comum, CRS ("Common-Reflection-Surface", em inglês) é uma nova técnica de processamento para simular seções sísmicas com afastamento-nulo, ZO ("Zero-Offset", em inglês) a partir de dados sísmicos de cobertura múltipla. Este método baseia-se na aproximação hiperbólica de tempos de trânsito paraxiais de segunda ordem referido ao raio (central) normal. O operador de empilhamento CRS para uma superfície de medição planar depende de três parâmetros, denominados o ângulo de emergência do raio normal, a curvatura da onda Ponto de Incidência Normal, NIP ("Normal Incidence Point", em inglês) e a curvatura da onda Normal, N. Neste artigo o método de empilhamento CRS ZO 2-D é modificado com a finalidade de considerar uma superfície de medição com topografia suave também dependente desses parâmetros. Com este novo formalismo CRS, obtemos uma seção sísmica ZO de alta resolução, sem aplicar as correições estáticas, onde em cada ponto desta seção são estimados os três parâmetros relevantes do processo de empilhamento CRS.

A quick review of 2D topographic traveltimes

O método de empilhamento Superfície de Reflexão Comum (SRC) foi originalmente introduzido como um método data-driven para simular seções afastamento-nulo a partir de dados sísmicos de reflexão pré-empilhados 2-D adquiridos ao longo de uma linha de aquisição reta. Este método está baseado em uma aproximação de tempos de trânsito hiperbólica de segunda ordem parametrizada com três atributos cinemáticos do campo de onda. Em dados terrestres, os efeitos topográficos desempenham um papel importante no processamento e imageamento de dados sísmicos. Assim, esta característica tem sido considerada recentemente pelo método SRC. Neste trabalho apresentamos uma revisão das aproximações de tempos de trânsito SRC que consideram topografia suave e rugosa. Adicionalmente, nós revemos também a aproximação de tempos de trânsito Multifoco para o caso da topografia rugosa. Por meio de um exemplo sintético simples, nós fornecemos finalmente as primeiras comparações entre as diferentes expressões de tempos de trânsito.

Migração (2,5-D) com amplitudes verdadeiras em meios com gradiente constante de velocidade

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.

3D raytracing through homogeneous anisotropic media with smooth interfaces

O raio conectando dois pontos em um meio anisotrópico, homogêneo por partes e com variação lateral, é calculado utilizando-se técnicas de continuação em 3D. Se combinado com algoritmos para solução do problema de valor inicial, o método pode ser estendido para o cálculo de eventos qS₁ e qS₂. O algoritmo apresenta a mesma eficiência e robustez que implementações de técnicas de continuação em meios isotrópicos. Rotinas baseadas neste algoritmo têm várias aplicações de interesse. Primeiramente, na modelagem e inversão de parâmetros elásticos na presença de anisotropia. Em segundo lugar, as iterações de Newton-Raphson produzem atributos da frente de onda como vetor vagarosidade e a matrix hessiana do tempo de trânsito, quantidades que permitem determinar o espalhamento geométrico e aproximações de segunda ordem para o tempo de trânsito. Estes atributos permitem calcular as amplitudes ao longo do raio e investigar os efeitos da anisotropia no empilhamento CRS em modelos de velocidade simples.

Estimation of fractures orientation from qP reflectivity using multiazimuthal AVO analysis

Investigamos a estimativa da orientação de fratura; mergulho e direção, através da análise de AVO multiazimutal de ondas qP e de suas convertidas qS₁ e qS₂. Assumimos fraco contraste de impedância, fraca anisotropia e ainda que, o meio fraturado comporta-se efetivamente como um meio transversalmente isotrópico (TI). Dentro destas hipóteses, a estimativa da orientação de fratura reduz-se a estimativa da orientação do eixo de simetria a partir de dados de refletividade de ondas qP. Aproximações lineares da refletividade da onda qP são utilizadas na inversão.

Wave equation depth migration using complex Padé approximation

Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domínio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido às ondas evanescentes e produz imagens com menos ruídos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos íngremes. Os resultados foram comparados àqueles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD).

Linearização dos coeficientes de reflexão de ondas qP em meios anisotrópicos

Os coeficientes de reflexão de onda qP em uma interface plana que separa dois meios anisotrópicos dependem dos parâmetros elásticos dos meios envolvidos de maneira não linear. Aproximações lineares para a refletividade de ondas qP pressupondo fraco contraste de impedância e fraca anisotropia levam a uma forma simples para se fazer à análise de AVO/AVD. Neste trabalho a solução das equações de Zoeppritz é rescrita explicitamente em função de matrizes de impedância e polarização. Além disso, utilizando-se esta abordagem, é apresentada uma metodologia geral mais simples para se obter as formas linearizadas. Estas formas linearizadas para os coeficientes de reflexão da onda qP e das ondas convertidas para uma onda qP incidente apresentam resultados muito próximos da formulação exata para ângulos de incidência menores que 30° considerando-se um contraste de impedância moderado e anisotropia em limites geologicamente aceitáveis.

Deconvolução de Euler: passado, presente e futuro – um tutorial

Neste tutorial apresentamos uma revisão da deconvolução de Euler que consiste de três partes. Na primeira parte, recordamos o papel da clássica formulação da deconvolução de Euler 2D e 3D como um método para localizar automaticamente fontes de campos potenciais anômalas e apontamos as dificuldades desta formulação: a presença de uma indesejável nuvem de soluções, o critério empírico usado para determinar o índice estrutural (um parâmetro relacionado com a natureza da fonte anômala), a exeqüibilidade da aplicação da deconvolução de Euler a levantamentos magnéticos terrestres, e a determinação do mergulho e do contraste de susceptibilidade magnética de contatos geológicos (ou o produto do contraste de susceptibilidade e a espessura quando aplicado a dique fino). Na segunda parte, apresentamos as recentes melhorias objetivando minimizar algumas dificuldades apresentadas na primeira parte deste tutorial. Entre estas melhorias incluem-se: i) a seleção das soluções essencialmente associadas com observações apresentando alta razão sinal-ruído; ii) o uso da correlação entre a estimativa do nível de base da anomalia e a própria anomalia observada ou a combinação da deconvolução de Euler com o sinal analítico para determinação do índice estrutural; iii) a combinação dos resultados de (i) e (ii), permitindo estimar o índice estrutural independentemente do número de soluções; desta forma, um menor número de observações (tal como em levantamentos terrestres) pode ser usado; iv) a introdução de equações adicionais independentes da equação de Euler que permitem estimar o mergulho e o contraste de susceptibilidade das fontes magnéticas 2D. Na terceira parte apresentaremos um prognóstico sobre futuros desenvolvimentos a curto e médio prazo envolvendo a deconvolução de Euler. As principais perspectivas são: i) novos ataques aos problemas selecionados na segunda parte deste tutorial; ii) desenvolvimento de métodos que permitam considerar interferências de fontes localizadas ao lado ou acima da fonte principal, e iii) uso das estimativas de localização da fonte anômala produzidas pela deconvolução de Euler como vínculos em métodos de inversão para obter a delineação das fontes em um ambiente computacional amigável.

Metodologia geofísica aplicada ao estudo arqueológico dos sítios Bittencourt e Jambuaçu, Estado do Pará

Medidas magnéticas, de radiação gama e eletromagnéticas com o Radar de Penetração no Solo (GPR) foram realizadas em dois sítios arqueológicos localizados nas proximidades de uma área destinada a construção de um mineroduto. O objetivo das medidas foi identificar locais para escavação arqueológica visando a descoberta de remanescentes de habitantes que povoaram a região em tempos pretéritos, auxiliando nos estudos para recompor a sua história. As medidas magnéticas evidenciaram anomalias que levaram a descoberta de uma grande quantidade de fragmentos cerâmicos e de outros artefatos. As medidas de radiação gama permitiram a discriminação entre solo arenoso e argiloso, o que representa um resultado importante, pois o solo arenoso geralmente é correlacionado a camada de ocupação arqueológica onde são encontrados os remanescentes. O Radar de Penetração no Solo confirmou as anomalias magnéticas e mostrou outras feições anômalas não relacionadas aos remanescentes arqueológicos.