Inversão de dados de sísmica de refração profunda a partir da curva tempo-distância

O trabalho em pauta tem como objetivo o modelamento da crosta, através da inversão de dados de refração sísmica profunda, segundo camadas planas horizontais lateralmente homogêneas, sobre um semi-espaço. O modelo direto é dado pela expressão analítica da curva tempo-distância como uma função que depende da distância fonte-estação e do vetor de parâmetros velocidades e espessuras de cada camada, calculado segundo as trajetórias do raio sísmico, regidas pela Lei de Snell. O cálculo dos tempos de chegada por este procedimento, exige a utilização de um modelo cujas velocidades sejam crescentes com a profundidade, de modo que a ocorrência das camadas de baixa velocidade (CBV) é contornada pela reparametrização do modelo, levando-se em conta o fato de que o topo da CBV funciona apenas como um refletor do raio sísmico, e não como refrator. A metodologia de inversão utilizada tem em vista não só a determinação das soluções possíveis, mas também a realização de uma análise sobre as causas responsáveis pela ambiguidade do problema. A região de pesquisa das prováveis soluções é vinculada segundo limites superiores e inferiores para cada parâmetro procurado, e pelo estabelecimento de limites superiores para os valores de distâncias críticas, calculadas a partir do vetor de parâmetros. O processo de inversão é feito utilizando-se uma técnica de otimização do ajuste de curvas através da busca direta no espaço dos parâmetros, denominado COMPLEX. Esta técnica apresenta a vantagem de poder ser utilizada com qualquer função objeto, e ser bastante prática na obtenção de múltiplas soluções do problema. Devido a curva tempo-distância corresponder ao caso de uma multi-função, o algoritmo foi adaptado de modo a minimizar simultaneamente várias funções objetos, com vínculos nos parâmetros. A inversão é feita de modo a se obter um conjunto de soluções representativas do universo existente. Por sua vez, a análise da ambiguidade é realizada pela análise fatorial modo-Q, através da qual é possível se caracterizar as propriedades comuns existentes no elenco das soluções analisadas. Os testes com dados sintéticos e reais foram feitos tendo como aproximação inicial ao processo de inversão, os valores de velocidades e espessuras calculados diretamente da interpretação visual do sismograma. Para a realização dos primeiros, utilizou-se sismogramas calculados pelo método da refletividade, segundo diferentes modelos. Por sua vez, os testes com dados reais foram realizados utilizando-se dados extraídos de um dos sismogramas coletados pelo projeto Lithospheric Seismic Profile in Britain (LISPB), na região norte da Grã-Bretanha. Em todos os testes foi verificado que a geometria do modelo possui um maior peso na ambiguidade do problema, enquanto os parâmetros físicos apresentam apenas suaves variações, no conjunto das soluções obtidas.

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.

Medidas de coerência para análise de velocidade na migração em tempo

Métodos iterativos para a análise de velocidade com migração dependem das funções objetivo medirem a horizontalização dos eventos de reflexão em famílias de ponto imagem comum (CIG). A migração em tempo é um método simples de imageamento para avaliar essas funções objetivo. Usando migração em tempo, estudamos a influência da função objetivo nos resultados da análise de velocidade com migração. Propomos duas novas funções objetivo para a análise velocidade com migração: Semblance Diferencial Estendido e o produto do Semblance Clássico com o Semblance Diferencial Estendido. Experimentos numéricos, utilizando os dados Marmousoft, mostram a eficácia das novas funções objetivo em estimar modelos velocidade produzindo eventos horizontalmente alinhados nas famílias de ponto imagem comum.

Migração 3D no tempo usando a teoria dos raios paraxiais

O trabalho aqui apresentado visa o estudo da migração 3D no tempo, em amplitudes verdadeiras, de dados da seção afastamento nulo ("zero-offset"), usando a aproximação de segunda ordem do tempo de trânsito de difração. O processo de migração no tempo usa a função difração determinada, corretamente, através de mesmos parâmetros que controlam a função tempo de trânsito de reflexão, além do auxílio dos raios paraxiais. Os parâmetros que controlam a função tempo de trânsito de reflexão são determinados através de pelo menos nove medidas de tempo de trânsito, obtidas a partir de dados pré-empilhados (tiro comum, ponto médio comum - CMP, etc.). Determinada a função difração, pode-se fazer a migração sem a necessidade de que se conheça o modelo de velocidades, pois a função tempo de trânsito de difração, como definido neste trabalho, não necessita do mesmo. O espalhamento geométrico é removido através do processo de migração e aplicação de uma função peso aos dados da seção afastamento nulo. Juntamente com a função peso, aplicou-se um filtro objetivando a recuperação da fase do sinal, ocorrida como consequência do processo de migração. Devido a limitação computacional, fez-se testes com dados sintéticos, em 2D, objetivando verificar a eficiência do programa. Os resultados obtidos foram satisfatórios, mostrando a eficiência e confiabilidade do processo.

Modelagem do mCSEM no domínio do tempo usando transformada discreta de Fourier

A modelagem do mCSEM é feita normalmente no domínio da frequência, desde sua formulação teórica até a análise dos resultados, devido às simplificações nas equações de Maxwell, possibilitadas quando trabalhamos em um regime de baixa frequência. No entanto, a abordagem através do domínio do tempo pode em princípio fornecer informação equivalente sobre a geofísica da subsuperfície aos dados no domínio da frequência. Neste trabalho, modelamos o mCSEM no domínio da frequência em modelos unidimensionais, e usamos a transformada discreta de Fourier para obter os dados no domínio do tempo. Simulamos ambientes geológicos marinhos com e sem uma camada resistiva, que representa um reservatório de hidrocarbonetos. Verificamos que os dados no domínio do tempo apresentam diferenças quando calculados para os modelos com e sem hidrocarbonetos em praticamente todas as configurações de modelo. Calculamos os resultados considerando variações na profundidade do mar, na posição dos receptores e na resistividade da camada de hidrocarbonetos. Observamos a influência da airwave, presente mesmo em profundidades oceânicas com mais de 1000m, e apesar de não ser possível uma simples separação dessa influência nos dados, o domínio do tempo nos permitiu fazer uma análise de seus efeitos sobre o levantamento. Como parte da preparação para a modelagem em ambientes 2D e 3D, fazemos também um estudo sobre o ganho de desempenho pelo uso do paralelismo computacional em nossa tarefa.

Modelagem de perfis de indução no domínio do tempo por diferenças finitas

Apresentamos aqui uma metodologia alternativa para modelagem de ferramentas de indução diretamente no domínio do tempo. Este trabalho consiste na solução da equação de difusão do campo eletromagnético através do método de diferenças finitas. O nosso modelo consiste de um meio estratificado horizontalmente, através do qual simulamos um deslocamento da ferramenta na direção perpendicular às interfaces. A fonte consiste de uma bobina excitada por uma função degrau de corrente e o registro do campo induzido no meio é feito através de uma bobina receptora localizada acima da bobina transmissora. Na solução da equação de difusão determinamos o campo primário e o campo secundário separadamente. O campo primário é obtido analiticamente e o campo secundário é determinado utilizando-se o método de Direção Alternada Implícita, resultando num sistema tri-diagonal que é resolvido através do método recursivo proposto por Claerbout. Finalmente, determina-se o valor máximo do campo elétrico secundário em cada posição da ferramenta ao longo da formação, obtendo-se assim uma perfilagem no domínio do tempo. Os resultados obtidos mostram que este método é bastante eficiente na determinação do contato entre camadas, inclusive para camadas de pequena espessura.

Modelamento e correção de descentralização das imagens de tempo de trânsito

As imagens da parede do poço obtidas através da ferramenta UBI (Ultrasonic Borehole Imager) são amplamente utilizadas por geólogos e engenheiros de petróleo para identificar eventos geológicos nas paredes de poços abertos e na inspeção dos tubos de revestimento, uma vez que este perfil praticamente fornece uma fotografia da parede do poço. As ferramentas de imageamento acústico produzem imagens do tempo de trânsito e da amplitude do pulso acústico gerado pela ferramenta e refletido na parede do poço. Entretanto, estas imagens podem ter uma interpretação não realista, uma vez que elas podem ter seu aspecto alterado em razão do movimento da ferramenta no interior do poço. Este trabalho apresenta o modelamento das imagens de tempo de trânsito a partir da aplicação do critério de Coulomb para a ruptura da parede do poço submetida a um estado plano de tensões, a qual fornecerá a seção do poço, que é a forma geométrica que será traçada pela ferramenta de imageamento acústico do poço. O deslocamento ascensional da ferramenta e as imperfeições da parede do poço, normalmente, são os responsáveis pelo deslocamento da posição do transdutor em relação ao eixo do poço. Este efeito pode ter grande responsabilidade nas imperfeições das imagens de tempo de trânsito. Assim, a correção dessas imagens, chamada de correção da descentralização, busca o reposicionamento do transdutor para a posição do eixo do poço. Apresenta-se, também, um método de correção do efeito da descentralização da ferramenta baseado neste modelamento. O método é proposto com base na geometria analítica plana e no método do raio para a definição do tempo de trânsito do pulso acústico, com o objetivo de reconstruir as imagens de tempo de trânsito obtidas com a ferramenta descentralizada, ou seja, corrigir estas imagens tornando-as como se fossem obtidas com a ferramenta centralizada em relação ao eixo do poço.

Modelagem do tempo de execução de obras civis: estudo de caso na Universidade Federal do Pará

Este artigo propõe a adoção de um modelo numérico prognóstico para estimar a variável "tempo de execução" para empreendimentos públicos de forma objetiva. O trabalho de campo consistiu na aplicação de métodos estatísticos para analisar dados de obras licitadas e executadas durante o período de 2006 a 2009 na Universidade Federal do Pará (UFPA). A análise de dados envolveu cálculos de regressões lineares e transformadas das funções. Após estratificação e tratamento inicial dos dados, os elementos adotados para construção do modelo final se restringiram a 102 obras de um total de 225 originariamente pesquisadas, resultando nos seguintes parâmetros estatísticos: coeficiente de correlação (R) de 0,899; coeficiente de determinação (R²) de 0,808; coeficiente de determinação ajustado (R² ajustado) de 0,796 e erro padrão (Se) de 0,41. Estes parâmetros indicam forte correlação linear entre as variáveis, indicando que 79,60% da variação do tempo para executar uma obra pública podem ser causadas pela variação, em conjunto, das variáveis área construída, custo orçado, capacidade técnica operacional do contratante, capacidade operacional da empresa, tipologia de serviço e estação do ano.

Uso de árvore de decisão para avaliação da segurança estática em tempo real de sistemas elétricos de potência

As técnicas utilizadas para avaliação da segurança estática em sistemas elétricos de potência dependem da execução de grande número de casos de fluxo de carga para diversas topologias e condições operacionais do sistema. Em ambientes de operação de tempo real, esta prática é de difícil realização, principalmente em sistemas de grande porte onde a execução de todos os casos de fluxo de carga que são necessários, exige elevado tempo e esforço computacional mesmo para os recursos atuais disponíveis. Técnicas de mineração de dados como árvore de decisão estão sendo utilizadas nos últimos anos e tem alcançado bons resultados nas aplicações de avaliação da segurança estática e dinâmica de sistemas elétricos de potência. Este trabalho apresenta uma metodologia para avaliação da segurança estática em tempo real de sistemas elétricos de potência utilizando árvore de decisão, onde a partir de simulações off-line de fluxo de carga, executadas via software Anarede (CEPEL), foi gerada uma extensa base de dados rotulada relacionada ao estado do sistema, para diversas condições operacionais. Esta base de dados foi utilizada para indução das árvores de decisão, fornecendo um modelo de predição rápida e precisa que classifica o estado do sistema (seguro ou inseguro) para aplicação em tempo real. Esta metodologia reduz o uso de computadores no ambiente on-line, uma vez que o processamento das árvores de decisão exigem apenas a verificação de algumas instruções lógicas do tipo if-then, de um número reduzido de testes numéricos nos nós binários para definição do valor do atributo que satisfaz as regras, pois estes testes são realizados em quantidade igual ao número de níveis hierárquicos da árvore de decisão, o que normalmente é reduzido. Com este processamento computacional simples, a tarefa de avaliação da segurança estática poderá ser executada em uma fração do tempo necessário para a realização pelos métodos tradicionais mais rápidos. Para validação da metodologia, foi realizado um estudo de caso baseado em um sistema elétrico real, onde para cada contingência classificada como inseguro, uma ação de controle corretivo é executada, a partir da informação da árvore de decisão sobre o atributo crítico que mais afeta a segurança. Os resultados mostraram ser a metodologia uma importante ferramenta para avaliação da segurança estática em tempo real para uso em um centro de operação do sistema.

Organofilização de uma Mg-bentonita da bacia do Parnaíba-Sul do Maranhão e sua utilização em poli (metacrilato de metila)

Bentonitas são argilas que tem como seu principal constituinte argilominerais do grupo da esmectita, predominantemente montmorillonita. De acordo com o cátion predominante no espaço intercamada da esmectita, a bentonita pode ser classificada como sódica, cálcica ou magnesiana. Essas argilas possuem vasta aplicação industrial, como fluidos de perfuração, pelotização, moldes de fundição, dentre outros. Para algumas aplicações mais específicas e que agregam maior valor ao produto final, como na síntese de nanocompósitos polímero/argila, faz-se necessário à intercalação de íons orgânicos na intercamada do argilomineral. No Brasil, a produção industrial de argilas organofílicas é pequena e voltada para os mercados de tintas, graxas e resinas de poliéster. Empresas do setor de bentonitas, que ainda não estão produzindo esse tipo de material, vêm mostrando crescente interesse nesta aplicação. Dentro desse contexto, este trabalho buscou avaliar o potencial da Bentonita Formosa, uma Mg-bentonita recentemente descrita e relativamente abundante no nordeste do Brasil, na produção de argilas organofílicas e sua aplicação em síntese de nanocompósitos polímero/argila. Para isso, foram realizadas sínteses variando a concentração dos íons surfactantes hexadeciltrimetilamônio (HDTMA⁺) e dodeciltrimetilamônio (DTMA⁺) em 0,7, 1,0 e 1,5 vezes o valor de CEC, com tempo de reação de 12 horas e variação de temperatura de 25 ºC e 80 ºC. A Mg-Bentonita in natura e ativada com carbonato de sódio foi utilizada como material de partida. Tanto o material de partida como as argilas organofílicas obtidas foram caracterizadas por DRX, DTA/TG e IV. As argilas que apresentaram melhores resultados de intercalação foram utilizadas nas proporções de 1%, 3% e 10% para a síntese de nanocompósitos poli(metacrilato de metila) (PMMA)/argila. As análises de DRX confirmaram a intercalação dos íons orgânicos no espaço intercamada da Mg-esmectita com e sem ativação. Com os resultados de IV foi possível observar que a razão de confôrmeros gauche/trans diminui com o aumento do espaçamento basal. Os resultados de DTA/TG confirmaram a estabilidade térmica das argilas organofílicas à temperatura máxima de 200 °C, o que possibilita a utilização desse material em síntese de nanocompósitos polímero/argila obtidos por processo de fusão. A análise de DRX confirmou a intercalação do PMMA no espaço intercamada da Mg-esmectita em todos os nanocompósitos produzidos. Com as análises de DSC foi possível observar o aumento da temperatura de transição vítrea para todos os nanocompósitos, quando comparados com PMMA puro. Com isso, é possível concluir que a Mg-Bentonita pode ser intercalada com íons alquilamônio, sem a necessidade prévia de ativação sódica, formando argilas organofílicas, assim como sua utilização em síntese de nanocompósitos. Essa possibilidade de utilização da Mg-bentonita in natura pode representar uma importante diferença em termos de custos de processo, na comparação com as bentonitas cálcicas existentes no Brasil, ou mesmo as importadas, que precisam ser ativadas durante o beneficiamento. Finalmente, acredita-se que a pesquisa deve avançar com a avaliação das propriedades mecânicas dos nanocompósitos produzidos neste trabalho, visando as diferentes possibilidades de aplicações desses materiais.

Perfil genotípico de resistência do vírus Influenza A (H1N1) pandêmico aos inibidores da neuraminidase em pacientes procedentes da mesorregião de Belém no período de maio de 2009 a maio de 2012

O vírus Influenza é o responsável pela gripe, uma doença que ocasiona milhões de mortes e hospitalizações todos os anos. Nas infecções severas, especialmente em pessoas com risco para complicações, os antivirais tornam-se os principais meios para o manejo clínico, merecendo especial destaque os inibidores da neuraminidase (INAs). De fato, na pandemia de 2009 a Organização Mundial da Saúde (OMS) recomendou o uso do oseltamivir para o tratamento dos doentes. Porém, devido à evolução genética viral, surgiram cepas com mutações no gene codificador da neuraminidase (NA) responsáveis por substituições aminoacídicas que levam à resistência aos fármacos INAs. Assim, a OMS passou a recomendar a vigilância de resistência genotípica para os vírus Influenza. Este trabalho teve como objetivos verificar a ocorrência de mutações no gene codificador da NA dos vírus Influenza A (H1N1) pandêmico que possam estar relacionadas à resistência aos INAs em cepas circulantes na mesorregião metropolitana de Belém no período de maio de 2009 a maio de 2012 e analisar, através da modelagem de proteínas, as substituições aminoacídicas da NA que possam estar influenciando na conformação protéica. Durante o período de estudo, foram recebidas no Laboratório de Vírus Respiratórios 2619 amostras clínicas de pacientes que apresentavam sinais e sintomas de infecção respiratória aguda com até cinco dias de evolução. Para a detecção do genoma viral foi feita a extração do RNA viral, seguida de RT-PCR em tempo real utilizando marcadores específicos para Influenza A H1N1pdm, resultando em 744 (28,4%) positivas. Parte das amostras positivas foram então inoculadas em células MDCK. Para as amostras isoladas em cultura de células, foi feita uma nova extração do RNA viral seguida de uma RT-PCR e semi-nested (PCR) utilizando iniciadores específicos para o gene NA, e posterior análise em sequenciador automático ABI Prism 3130xl (Applied Biosystems). A modelagem molecular da NA foi realizada através dos softwares SWISS-MODEL, MODELLER 9.10, PROCHECK, VERIFY3D e PYMOL. A análise parcial das sequências da neuraminidase nas amostras sequenciadas mostrou que não houve a circulação de cepas de vírus H1N1pdm com a mutação H275Y, a principal envolvida na resistência ao oseltamivir. Porém, em duas amostras foi identificada a substituição D199N que já foi relatada em vários estudos mostrando uma possível associação com o aumento da resistência ao oseltamivir. As amostras de 2012 apresentaram duas substituições (V241I e N369K) que estão relacionadas com um possível papel na compensação dos efeitos negativos causados pela mutação H275Y. A modelagem molecular mostrou que na mutação D199N houve uma alteração na estrutura da proteína NA próxima ao sítio de ligação ao antiviral. A análise filogenética revelou que as amostras de 2012 formaram um cluster isolado, demonstrando uma variação muito mais temporal do que geográfica. Este representa o primeiro estudo de resistência dos vírus Influenza H1N1pdm na mesorregião metropolitana de Belém, representando um importante instrumento para que os profissionais de saúde adotem estratégias mais eficazes no manejo da doença e no desenvolvimento de novos fármacos anti-influenza.

Memórias de trabalho nas paisagens missioneiras do "antes-tempo"

O artigo em questão estabelece reflexões em torno do tema do trabalho na região missioneira gaúcha, tomando como referência o legado das missões jesuítico-guaranis e seus desdobramentos no contemporâneo, a partir das memórias das pessoas que vivem naquela porção do Estado. Para tanto, o artigo se baseia numa pesquisa de campo de cerca de um ano, que teve como um dos interesses etnográficos a tentativa de compreender os processos de transformação das paisagens na região através das ações técnico-culturais relacionadas ao universo do trabalho.

Partial time use of anterior repositioning splints in the management of TMJ pain and dysfunction: a one-year controlled study

This study aimed at evaluating the effectiveness of partial use of anterior repositioning appliances in the management of TMJ pain and dysfunction when compared to stabilization splints and a control group in a one-year follow-up. Sample was initially constituted by 60 patients, randomly divided into three groups: I- stabilization splints, II- repositioning splints and III- no treatment. The whole sample was evaluated by means of TMJ and muscle palpation, mandibular AROM, analysis of occlusal contacts, joint sounds inspection and Visual Analogue Scale (VAS) for one year; 52 patients composed the final sample. A significant (after 15 days) improvement in pain report (VAS) and palpation index was found for group II (p≤0.01). The occurrence of occlusal alterations as posterior open bite or gross interferences after the splint therapy and increased muscle tenderness were not problems in this study. Similar results in joint noises reduction were observed for the entire sample. It was concluded that controlled partial use of repositioning splints is a beneficial tool in the management of intra-articular pain and dysfunction, with no risks of irreversible occlusal changes.