CRS seismic data imaging system: a case study for basib reevaluation

Este trabalho apresenta resultados práticos de uma atenção sistemática dada ao processamento e à interpretação sísmica de algumas linhas terrestres do conjunto de dados do gráben do Tacutu (Brasil), sobre os quais foram aplicadas etapas fundamentais do sistema WIT de imageamento do empilhamento CRS (Superfície de Reflexão Comum) vinculado a dados. Como resultado, esperamos estabelecer um fluxograma para a reavaliação sísmica de bacias sedimentares. Fundamentado nos atributos de frente de onda resultantes do empilhamento CRS, um macro-modelo suave de velocidades foi obtido através de inversão tomográfica. Usando este macro-modelo, foi realizado uma migração à profundidade pré- e pós-empilhamento. Além disso, outras técnicas baseadas no empilhamento CRS foram realizadas em paralelo como correção estática residual e migração de abertura-limitada baseada na zona de Fresnel projetada. Uma interpretação geológica sobre as seções empilhadas e migradas foi esboçada. A partir dos detalhes visuais dos painéis é possível interpretar desconformidades, afinamentos, um anticlinal principal falhado com conjuntos de horstes e grábens. Também, uma parte da linha selecionada precisa de processamento mais detalhado para evidenciar melhor qualquer estrutura presente na subsuperfície.

Migração em profundidade pré-empilhamento utilizando os atributos cinemáticos do empilhamento por superfície de reflexão comum

O empilhamento por superfície de reflexão comum (ou empilhamento SRC), conhecido como empilhamento CRS, do inglês Commom reflection surface, constitui-se em um novo método para o processamento sísmico na simulação de seções afastamento nulo (AN) e afastamento comum (AC). Este método é baseado em uma aproximação paraxial hiperbólica de segunda ordem dos tempos de trânsito de reflexão na vizinhança de um raio central. Para a simulação de seção AN, o raio central é um raio normal, enquanto que para a simulação de uma seção AC o raio central é um raio de afastamento finito. Em adição à seção AN, o método de empilhamento SRC também fornece estimativas dos atributos cinemáticos do campo de onda, sendo aplicados, por exemplo, na determinação (por um processo de inversão) da velocidade intervalar, no cálculo do espalhamento geométrico, na estimativa da zona de Fresnel, e também na simulação de eventos de tempos de difrações, este último tendo uma grande importância para a migração pré-empilhamento. Neste trabalho é proposta uma nova estratégia para fazer uma migração em profundidade pré-empilhamento, que usa os atributos cinemáticos do campo de onda derivados do empilhamento SRC, conhecido por método CRS-PSDM, do inglês CRS based pre-stack depth migration. O método CRS-PSDM usa os resultados obtidos do método SRC, isto é, as seções dos atributos cinemáticos do campo de onda, para construir uma superfície de tempos de trânsito de empilhamento, ao longo da qual as amplitudes do dado sísmico de múltipla cobertura são somadas, sendo o resultado da soma atribuído a um dado ponto em profundidade, na zona alvo de migração que é definida por uma malha regular. Similarmente ao método convencional de migração tipo Kirchhoff (K-PSDM), o método CRS-PSDM precisa de um modelo de velocidade de migração. Contrário ao método K-PSDM, o método CRS-PSDM necessita apenas computar os tempos de trânsito afastamento nulo, ao seja, ao longo de um único raio ligando o ponto considerado em profundidade a uma dada posição de fonte e receptor coincidentes na superfície. O resultado final deste procedimento é uma imagem sísmica em profundidade dos refletores a partir do dado de múltipla cobertura.

Migração 3-D Kirchhoff-Gaussian-Beam (KGB) pré-empilhamento no domínio da profundidade

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.

Migração Kirchhoff pré-empilhamento em profundidade modificada usando o operador de feixes gaussianos

A teoria dos feixes gaussianos foi introduzida na literatura sísmica no início dos anos 80 por pesquisadores russos e tchecos, e foi originalmente utilizada no cálculo do campo de ondas eletromagnéticas, baseado na teoria escalar da difração. Na teoria dos feixes gaussianos, o campo de ondas sísmicas é obtido por uma integral, cujo o integrando é constituído de duas partes, a saber: (1) as amplitudes dos campos das ondas na vizinhança do ponto de observação e (2) a função fase de cada um desses campos de ondas, que neste caso é representada por um tempo de trânsito paraxial complexo. Como ferramenta de imageamento, mais precisamente como operador de migração, os primeiros trabalhos usando feixes gaussianos datam do final da década de 80 e início dos anos 90. A regularidade dos campos de ondas descritos pelos feixes gaussianos, além de sua alta precisão em regiões singulares do modelo de velocidades, tornaram o uso de feixes gaussianos como uma alternativa híbrida viável para a migração. Nesse trabalho, unimos a flexibilidade da migração tipo Kirchhoff em profundidade em verdadeira amplitude com a regularidade da descrição do campo de ondas, representado pela sobreposição de feixes gaussianos. Como forma de controlar de forma estável quantidades usadas na construção de feixes gaussianos, utilizamos informações advindas do volume de Fresnel, mais precisamente a zona de Fresnel ao redor do ponto de reflexão e a zona de Fresnel projetada, localizada ao redor do ponto de registro do sismograma e cuja a informação se encontra nas curvas de reflexão de dados sísmico. Nosso processo de migração pode ser chamado como uma migração Kirchhoff em verdadeira amplitude usando um operador de feixes gaussianos.

Aproximações hiperbólicas dos tempos de trânsito de raios paraxiais a um raio central refletido e difratado

A simulação de uma seção sísmica de afastamento nulo (AN) a partir de dados de cobertura múltipla para um meio 2-D, através do empilhamento, é um método de imageamento de reflexão sísmica muito utilizado, que permite reduzir a quantidade de dados e melhor a razão sinal/ruído. Baseado na aproximação hiperbólica dos tempos de trânsito dependente de três parâmetros ou atributos cinemáticos de frentes de onda, recentemente, vem desenvolvendo-se um novo método para simular seções (AN) chamado método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC). Também, seguindo este novo conceito de imageamento sísmico está surgindo um método para simular seções com afastamento comum (AC) a partir de dados de cobertura múltipla usando aproximações dos tempos de trânsito paraxiais na vizinhança de um raio central com afastamento finito. Esta nova aproximação dos tempos de trânsito depende de cinco atributos cinemáticos. Neste trabalho, a partir da aproximação dos tempos de trânsito paraxiais em relação a um raio central com afastamento finito, derivamos uma nova equação do tempo de trânsito usando a condição de um ponto difrator em profundidade, reduzindo a equação original para quatro parâmetros. Para ambas aproximações (reflexão e difração), mostramos a superfície de empilhamento SRC com afastamento finito (SRC-AF). Considerando um modelo sintético, realizamos um estudo comparativo das aproximações dos tempos de trânsito para as quatro configurações sísmicas (fonte comum (FC), receptor comum (RC), ponto-médio-comum (PMC) e afastamento comum (AC)). Para analisar o comportamento do operador SRC-AF, quando este é perturbado, discutimos sua sensibilidade em relação a cada um dos cinco parâmetros (K₁, K₂, K₃, β_S e β_G). Esta análise de sensibilidade é realizada em duas formas: Sensibilidade através da primeira derivada e Sensibilidade no Empilhamento SRC-AF. Após realizar a análise de sensibilidade utilizamos uma nova condição, K₂ = 0 e assim, obtemos uma nova aproximação, agora dependente de três parâmetros. Usando essas aproximações hiperbólicas (em função de cinco, quatro e três parâmetros), propomos um algoritmo para a simulação de seções AC a partir de dados de cobertura múltipla. Finalmente, é apresentado um estudo da zona de Fresnel, com o objetivo de determinar a delimitação da abertura da superfície de empilhamento SRC-AF.

Ecomorphological patterns of the fishes inhabiting the tide pools of the Amazonian Coastal Zone, Brazil

Identificamos padrões ecomorfológicos que refletem a ecologia de espécies encontradas em poças de maré na Zona Costeira Amazônica (ZCA). Indivíduos de 19 espécies foram coletados no estado do Pará durante duas expedições em 2011. Foram estabelecidas dominância, grau de residência, guildas tróficas e tomadas medidas morfométricas de até 10 indivíduos de cada espécie. Calculou-se 23 atributos ecomorfológicos relacionados à locomoção, posição e forrageio, utilizados para o cálculo da distância ecomorfológica. Foram utilizadas Análises de Componentes Principais (PCA) para avaliar que atributos ecomorfológicos explicaram a variação entre as espécies. O teste de Mantel foi utilizado para testar a correlação da distância taxonômica com a morfologia das espécies e um teste de Mantel parcial para avaliar a correlação das guildas tróficas com os padrões ecomorfológicos, controlando-se o efeito da distância taxonômica entre as espécies. Nas análises formaram-se dois eixos principais para variação em relação aos padrões de locomoção, correlacionados à largura do pedúnculo caudal e formato da nadadeira anal, ocorrendo influência da distância taxonômica entre as espécies nos padrões ecomorfológicos. Espécies dominantes e residentes apresentaram menor capacidade de natação contínua. Quanto à posição na coluna d'água, formaram-se dois eixos principais da variação, correlacionados à posição do olho, área da nadadeira pélvica e formato do corpo, ocorrendo influência da distância taxônomica entre as espécies nas dissimilaridades morfológicas. A PCA agrupou espécies de hábito pelágico com espécies de hábito bentônico. Em relação ao forrageio, formaram-se dois eixos principais da variação, correlacionados ao tamanho da boca, tamanho do olho e comprimento do trato digestório. Espécies de diferentes guildas permaneceram agrupadas, sugerindo fraca relação da morfologia com o forrageio e não houve influência da distância taxonômica nas dissimilaridades nas guildas tróficas. Espécies residentes e dominantes em poças de maré na ZCA apresentam hábito sedentário, ocorrendo pouca influência da distância taxonômica nos padrões ecomorfológicos que se referem à posição na coluna d'água e locomoção, demonstrando que espécies distantes filogeneticamente podem apresentar padrões ecomorfológicos similares, e a morfologia demonstrou-se como fraca preditora das táticas de forrageio.

Edge effects and the impact of wildfires on populations of small non-volant mammals in the forest-savanna transition zone in Southern Amazonia

O efeito de borda e do fogo sobre a comunidade de pequenos mamíferos não-voadores foi investigada em uma área de transição entre Cerrado e Floresta em uma área de matriz de soja na Amazônia Oriental. Os indivíduos foram coletados em 24 transectos, dos quais 16 foram distribuídos em área sem efeito do fogo e oito distribuídos com efeito do fogo. Um total de 11 espécies foi registrado, incluindo seis roedores e cinco marsupiais. A espécie Hylaeamys megacephalus foi a mais abundante em áreas sem efeito do fogo. A abundância e riqueza de pequenos mamíferos não-voadores apresentaram uma diminuição em áreas queimadas, entretanto o efeito do fogo parece mascarar o efeito de borda nestas mesmas áreas. Em relação ao efeito de borda, sem nenhum efeito de fogo, a relação entre a abundância de pequenos mamíferos não-voadores com a distância da borda foi positiva. A relação entre a diversidade de pequenos mamíferos e efeito de borda pode ser ligada à vegetação matriz e características ecológicas de cada espécie.

Estudo da estrutura geo-elétrica da Região do Juruá, AM, pelo método magnetotelúrico

Utilizando-se dados magnetotelúricos (MT), foi obtida uma imagem geo-elétrica nítida da região do Juruá, Bacia do Solimões, na forma de seções geo-elétricas. Os dados de campo foram registrados ao longo de três linhas de 15 km, espaçadas de 3.5 km, recobrindo uma área de 100 km². O espaçamento entre as 35 estações é irregular, variando de 400 m a 3500 m. A faixa de freqüências utilizada cobriu de 0.001 Hz até 300 Hz, o que permitiu investigar de 100 m até 60 km de profundidade. Os dados apresentam-se afetados pelo efeito de distorção estática. Para corrigir este efeito foi utilizada a mediana da resistividade do primeiro condutor, correspondente à Formação Solimões. Foi utilizado o invariante do tensor MT para interpretar a estrutura geo-elétrica do Juruá. As seções geo-elétricas foram obtidas a partir do agrupamento dos dados resultantes da transformação de Bostick e da inversão 1D de Occam, para cada estação. Foi identificada uma seqüência de camadas condutivas e resistivas, correspondentes ao pacote sedimentar, uma zona de falhas e o topo do embasamento geo-elétrico, caracterizando a Bacia do Solimões. Abaixo do embasamento geo-elétrico foram também identificados uma zona condutora, seguida por uma camada de baixa condutividade, a profundidades iguais ou superiores a 20 km. Esta camada é interpretada como sendo de composição de gabro, estando associada a processos de acreção vertical, intimamente ligados à estabilização crustal e espessamento da litosfera. Os resultados apresentam uma boa concordância com os perfis de resistividade de poços e dados sísmicos de superfície.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

Métodos geofísicos aplicados a prospecção de água subterrânea no município de Ponta de Pedras - Ilha de Marajó - Pará

Este trabalho foi realizado numa área de 25 km² em Ponta de Pedras, Ilha de Marajó. Tem por objetivo o dimensionamento geométrico e seleção de zonas favoráveis à água subterrânea. Foram utilizados dois métodos geofísicos: sísmica de refração e eletroresistividade. O método sísmico foi empregado em caráter experimental, visando verificar a viabilidade de sua aplicação na área. Os resultados da interpretação confirmaram a suspeita inicial, da existência de horizontes geológicos pouco profundos, não possíveis de serem distinguidos pelo método sísmico. Com o método de eletroresistividade foram feitas sondagens verticais e perfis horizontais. Os resultados da interpretação indicaram áreas promissoras na parte central e numa faixa ao norte da área prospectada. A espessura do aquífero raso varia entre 0 e 50 m.

Environmental sensitivity index (ESI) mapping of oil spill in the amazon coastal zone: the PIATAM Mar project

A importância do monitoramento ambiental é medida pelos vários casos de derramamentos de óleo ocorridos no mundo durante as últimas três décadas. Isto tem incentivado as empresas e órgãos do governo envolvidos na prevenção e combate a estes acidentes a aperfeiçoarem cada vez mais os métodos, tanto preventivos como corretivos, para a minimização dos danos gerados por acidentes com derramamento de óleo. Este trabalho objetiva contextualizar de forma histórica como os acidentes com derramamento de óleo propiciaram o desenvolvimento de pesquisa tecnológica a partir de parcerias entre empresas de petróleo, agências de governo, universidades e institutos de pesquisa no Brasil, em especial na zona costeira Amazônica. Como resultado, índices de sensibilidade ambiental ao derramamento de óleo (ISA) foram definidos especialmente para a Amazônia costeira, onde processos fluviais e marinhos se encontram na foz do maior rio do mundo, o rio Amazonas. Perspectivas de pesquisa e respostas de emergência a acidentes são apresentadas, a fim de se conservar a diversidade socioambiental da mais importante região tropical do planeta.

Spatio-temporal distribution of fish larvae in relation to ontogeny and water quality in the oligohaline zone of a North Brazilian estuary

O padrão de distribuição das larvas de peixes na Baía do Guajará foi estudado através de coletas trimestrais e discutido em relação aos estágios de desenvolvimento dos indivíduos e a exposição dos mesmos aos contaminantes despejados pela Capital do Estado do Pará, Belém. A densidade e a diversidade larval foram baixas com forte dominância dos clupeídeos, engraulídeos e, em menor grau, os cianídeos. O principal período de reprodução foi definido no início do período chuvoso. Larvas de Clupeiformes em pré-flexão e flexão foram encontradas nos pontos de coletas mais afastados da cidade, enquanto aquelas em pós-flexão prevaleceram nas margens da cidade. Por outro lado, os cianídeos em pré-flexão e flexão foram capturados perto dos centros de atividade urbana, enquando aqueles em pós-flexão foram pouco abundantes. É sugerido que a baía se encontra na rota migratória dos cianídeos entre a área de desova e os berçários mais distantes. Apesar das águas no entorno da cidade de Belém mostraram sinais de contaminação, a qualidade ambiental na Baía do Guajará no momento do estudo estava apropriada para a vida das larvas de peixes. Nitrato com pH foram as variáveis que melhor explicaram a distribuição larval.