Modelo fractal para resistividade complexa de rochas: interpretação petrofísica e aplicação à exploração geoelétrica

Rochas contendo metálicos disseminados ou partículas de argila em ambiente natural onde soluções eletrolíticas normalmente preenchem os poros das rochas, exibem um tipo de polarização em baixas freqüências conhecido como polarização induzida. Nesta tese foi desenvolvido um novo modelo para descrever o fenômeno de polarização das rochas, não apenas em baixas freqüências, mas compreendendo todo o espectro eletromagnético, possível de utilização na prospecção geoelétrica. Este novo modelo engloba a maioria dos modelos utilizados até o momento como casos especiais, além de superar as limitações dos mesmos. Seu circuito analógico inclui uma impedância não linear do tipo r (iwtf)^-n que simula o efeito das superfícies rugosas das interfaces entre os grãos bloqueadores (partículas metálicas e/ou de argilas) e o eletrólito. A impedância de Warburg generalizada está em série com a resistência dos grãos bloqueadores da passagem de corrente e em paralelo com a impedância da dupla camada associada a essas interfaces. Esta combinação está em série com a resistência do eletrólito nas passagens dos poros bloqueados. Os canais não bloqueados são representados por uma resistência que corresponde à resistividade normal CC da rocha. A combinação desta resistência com a capacitância "global" da rocha é finalmente conectada em paralelo ao resto do circuito mencionado acima. Os parâmetros deste modelo incluem a resistividade CC (p₀), a cargueabilidade (m), três tempos de relaxação (t, Tf and T2), um fator de resistividade de grãos (δ_r), e o expoente de freqüência (η). O tempo de relaxação fractal (Tf), e o expoente de frequencia (η) estão relacionados à geometria fractal das interfaces rugosas entre os minerais condutivos (grãos metálicos e/ou partículas de argila bloqueando os canais dos poros) e o eletrólito. O tempo de relaxação (T) é um resultado da relaxação em baixa freqüência das duplas camadas elétricas formadas nas interfaces eletrólito-cristais, enquanto (T₀) é o tempo de relaxação macroscópico da amostra como um todo. O fator de resistividade dos grãos (δ_r) relaciona a resistividade dos grãos condutivos com o valor de resistividade CC da rocha. A resistividade CC da rocha (p₀), e δ_r estão relacionados à porosidade, à condutividade do eletrólito e às relações mineralógicas entre a matriz e os grãos condutivos. O modelo foi testado sobre um intervalo largo de freqüências contra dados experimentais de amplitude e fase da resistividade bem como para dados de constante dielétrica complexa. Os dados utilizados neste trabalho foram obtidos a partir da digitalização de dados experimentais publicados, obtidos por diversos autores e englobando amostras de rochas sedimentares, ígneas e metam6rficas. É mostrado neste trabalho que os parâmetros deste modelo permitem identificar diferenças texturais e mineralógicas nas rochas. Bote modelo foi introduzido, primeiramente, como propriedade intrínseca de um semiespaço homogêneo sendo demonstrado, neste trabalho, que a resposta observada em superfície reflete as propriedades intrínsecas do meio polarizável, sendo o acoplamento eletromagnético desprezível em freqüências menores que 104 Hz. Em seguida, o meio polarizável foi embebido em um pacote de N camadas sendo demonstrado que os parâmetros fractais do meio polarizável podem ser obtidos do levantamento em superfície para diferentes espessuras dessa camada. Isto justifica a utilização pura e simples de modelos de polarização desenvolvidos para amostras em laboratório para ajustar dados de campo, o que vem sendo feito sem uma justificativa bem fundamentada. Estes resultados demonstram a importância para a prospecção geolétrica do modelo proposto nesta tese.

Efeitos do strike das estruturas geológicas de duas dimensões nas pseudos-seções de IP resistividade

Os levantamentos de IP-resistividade efetuados na Serra dos Carajás não foram executados ortogonalmente às estruturas geológicas, pois utilizaram linhas anteriormente abertas pelas equipes de geoquímica. Este fato motivou este estudo teórico da influência da direção das linhas de medidas de IP-resistividade em relação ao "strike" da estrutura. Foi usado o programa de elementos finitos de Rijo (1977), desenvolvido para levantamentos perpendiculares às estruturas com as adaptações necessárias. A modificação principal foi na rotina de transformação inversa de Fourier. Para o caso simples dos levantamentos perpendiculares, a transformada inversa é uma integral discreta com apenas sete pontos. No entanto, para as medidas obliquas, o integrando é oscilatório, e portanto, a integral a ser calculada é mais complexa. Foi adaptado um método apresentado por Ting e Luke (1981), usando dezoito pontos em cada integração. Foi constatado que o efeito da direção da linha em relação ao "strike" é desprezível para ângulos maiores que 60 graus. Para ângulos menores, o efeito consiste no alongamento da anomalia, com pequenas alterações em seu centro. Não há uma maneira simples de compensar este efeito com mudanças nos parâmetros do modelo.

Estimativas da condutividade térmica dos minerais e rochas e influência de parâmetros térmicos e petrofísicos na resistividade aparente da formação

O presente estudo realiza estimativas da condutividade térmica dos principais minerais formadores de rochas, bem como estimativas da condutividade média da fase sólida de cinco litologias básicas (arenitos, calcários, dolomitos, anidritas e litologias argilosas). Alguns modelos térmicos foram comparados entre si, possibilitando a verificação daquele mais apropriado para representar o agregado de minerais e fluidos que compõem as rochas. Os resultados obtidos podem ser aplicados a modelamentos térmicos os mais variados. A metodologia empregada baseia-se em um algoritmo de regressão não-linear denominado de Busca Aleatória Controlada. O comportamento do algoritmo é avaliado para dados sintéticos antes de ser usado em dados reais. O modelo usado na regressão para obter a condutividade térmica dos minerais é o modelo geométrico médio. O método de regressão, usado em cada subconjunto litológico, forneceu os seguintes valores para a condutividade térmica média da fase sólida: arenitos 5,9 ± 1,33 W/mK, calcários 3.1 ± 0.12 W/mK, dolomitos 4.7 ± 0.56 W/mK, anidritas 6.3 ± 0.27 W/mK e para litologias argilosas 3.4 ± 0.48 W/mK. Na sequência, são fornecidas as bases para o estudo da difusão do calor em coordenadas cilíndricas, considerando o efeito de invasão do filtrado da lama na formação, através de uma adaptação da simulação de injeção de poços proveniente das teorias relativas à engenharia de reservatório. Com isto, estimam-se os erros relativos sobre a resistividade aparente assumindo como referência a temperatura original da formação. Nesta etapa do trabalho, faz-se uso do método de diferenças finitas para avaliar a distribuição de temperatura poço-formação. A simulação da invasão é realizada, em coordenadas cilíndricas, através da adaptação da equação de Buckley-Leverett em coordenadas cartesianas. Efeitos como o aparecimento do reboco de lama na parede do poço, gravidade e pressão capilar não são levados em consideração. A partir das distribuições de saturação e temperatura, obtém-se a distribuição radial de resistividade, a qual é convolvida com a resposta radial da ferramenta de indução (transmissor-receptor) resultando na resistividade aparente da formação. Admitindo como referência a temperatura original da formação, são obtidos os erros relativos da resistividade aparente. Através da variação de alguns parâmetros, verifica-se que a porosidade e a saturação original da formação podem ser responsáveis por enormes erros na obtenção da resistividade, principalmente se tais "leituras" forem realizadas logo após a perfuração (MWD). A diferença de temperatura entre poço e formação é a principal causadora de tais erros, indicando que em situações onde esta diferença de temperatura seja grande, perfilagens com ferramentas de indução devam ser realizadas de um a dois dias após a perfuração do poço.

Resistividade e polarização induzida (IP): modelagem analógica

Foi desenvolvida neste trabalho uma nova técnica de modelagem analógica para IP-resistividade em laboratório. Para tanto, foram utilizados modelos geométricos cilíndricos e tabulares. Para a construção destes modelos usou-se cimento, areia quartzítica e grafita. Como a grafita apresenta um intenso fenômeno IP, variou-se o teor de grafita contido nos modelos. Assim, para diferentes formas e medidas dos modelos, experimentaram-se diferentes conteúdos de partículas polarizáveis. As medidas foram feitas no sentido de verificar o comportamento da resposta IP e resistividade dos modelos quanto à forma, conteúdo de grafita, profundidade e orientação do modelo estudado.

Modelagem unidimensional de dados do método MCSEM com variação na resistividade da água do mar e na orientação dos sensores

O Marine Controlled Source ElectroMagnetic (MCSEM) é um novo método geofísico que tem sido empregado para a detecção de camadas resistivas no sub-solo oceânico e tem se mostrado uma ferramenta promissora para a indústria do petróleo. Neste trabalho, é modelado levantamentos do método MCSEM incluindo variações nos valores de resistividade da água do mar. Tais variações podem resultar da influência de correntes marinhas, gradientes de temperatura ou de qualquer fator que possa afetar a salinidade da água. Observamos que os dados do método são afetados de maneira significativa por tais variações. Tanto a variação na resistividade quanto na espessura da camada de água pode afetar os resultados do MCSEM de maneira semelhante. É modelado também o posicionamento aleatório dos receptores, mostrando um decaimento significativo da amplitude do sinal quando os receptores se posicionam a uma certa distância da linha de medidas.

DesG: uma nova metodologia para intervenção e coleta de dados magnetotelúricos

Esta tese apresenta duas contribuições para a interpretação geofísica, voltadas ao Método Magnetotelúrico. A primeira trata de uma nova abordagem para a interpretação MT, denominada MÉTODO DESCRITIVO-GEOLÓGICO (DesG), em alusão à incorporação explícita de informação a priori de correlação fácil com a descrição geológica tradicional. O intérprete define por meio de elementos geométricos (pontos e linhas) o arcabouço de feições geológicas bem como fornece valores de resistividade aos corpos geológicos presumidos. O método estima a distribuição de resistividade subsuperficial em termos de fontes anômalas próximas aos elementos geométricos, ajustando as respostas produzidas por estes corpos às medidas de campo. A solução obtida fornece então informações que podem auxiliar na modificação de algumas informações a priori imprecisas, permitindo que sucessivas inversões sejam realizadas até que a solução ajuste os dados e faça sentido geológico. Entre as características relevantes do método destacam-se: (i) os corpos podem apresentar resistividade maior ou menor do que a resistividade do meio encaixante, (ii) vários meios encaixantes contendo ou não corpos anômalos podem ser cortados pelo perfil e (iii) o contraste de resistividade entre corpo e encaixante pode ser abrupto ou gradativo. A aplicação do método a dados sintéticos evidencia, entre outras vantagens, a sua potencialidade para estimar o mergulho de falhas com inclinação variável, que merece especial interesse em Tectônica, e delinear soleiras de diabásio em bacias sedimentares, um sério problema para a prospecção de petróleo. O método permite ainda a interpretação conjunta das causas do efeito estático e das fontes de interesse. A aplicação a dados reais é ilustrada tomando-se como exemplo dados do COPROD2, cuja inversão produziu soluções compatíveis com o conhecimento geológico sobre a área. A segunda contribuição refere-se a desenho de experimento geofísico. Por meio de indicadores diversos, em especial a matriz densidade de informação, é mostrado que a resolução teórica dos dados pode ser estudada, que guia o planejamento da prospecção. A otimização do levantamento permite determinar os períodos e as posições das estações de medida mais adequados ao delineamento mais preciso de corpos cujas localizações são conhecidas aproximadamente.

Ecobiometria e Triplex Doppler para a avaliação da vitalidade fetal em Cebus apella

Os objetos do trabalho foram: 1) realizar a ecobiometria do concepto para acompanhar o crescimento fetal e determinar a idade gestacional em Cebus apella; 2) descrever o momento em que os órgãos do feto são observados; 3) realizar a sexagem fetal; 4) avaliar o fluxo sanguíneo das artérias uterinas (AU) e umbilical (Aum), mensurando os índices de resistividade (IR) e pulsatilidade (IP); 5) observar a presença ou ausência da incisura protodiastólica (IPD) e do componente diastólico nas ondas de fluxo sanguíneo das AU e Aum, respectivamente, durante a gestação em C. apella. Seis fêmeas adultas da espécie C. apella foram selecionadas e, posteriormente, condicionadas para os procedimentos de colpocitologia ou contenção química. Para o monitoramento do ciclo reprodutivo das fêmeas e crescimento folicular, foram realizados exames colpocitológicos e ultrassonografia (US), respectivamente, para a escolha do melhor dia da monta natural ou inseminação artificial. Vinte dias após a cópula ou inseminação, era realizado o diagnóstico da gestação por meio da US. O dia de cada exame em relação ao parto foi contado em retrospecto (nascimento = dia 0). Os exames ultrassonográficos foram feitos nos dias -133, -113, -83, -53, -21 e -1 antes do parto. A US bidimensional em modo B foi utilizada para mensurar o saco gestacional do dia - 133 até -113; os diâmetros biparietal, fronto-occiptal, circunferência da cabeça e circunferência abdominal do dia -113 a -1; e comprimento do fêmur do dia -53 até -23, mostrando o aumento desses parâmetros com o avanço da gestação. O coração e o estômago começaram a ser visualizados no dia -113 e, os demais órgãos e a sexagem no dia -83. O Triplex Doppler foi empregado para avaliar o fluxo sanguíneo durante o período gestacional, mostrando uma diminuição nos IR e IP das AU e Aum, do dia -133 a -1, bem como observar o desaparecimento da IPD (dia -1) e aparecimento do componente diastólico (dia -53) na onda de fluxo sanguíneo das AU e Aum, respectivamente. Entre os dias -113 a -1, a média da freqüência cardíaca fetal, obtida pelo Triplex Doppler, foi de 189 ± 2,43 bpm. O presente trabalho permitiu determinar a idade gestacional, avaliar o crescimento anatômico do feto, descrever o momento em que os órgãos são visualizados e realizar a sexagem em C. apella. Constatou-se, também, que o fluxo sanguíneo das AU e Aum são importantes parâmetros para avaliar a vitalidade fetal em C. apella.

Caracterização da liga 6101 refinada com a adição de diferentes teores de cobre e solidificadas em molde unidirecional horizontal e em molde "U"

A liga 6101 (série 6xxx) foi modificada através da adição de cobre nos percentuais de 0,05%Cu e 0,3%Cu, e também de refinadores de grão à base de uma pré-liga Al-Ti-B. As ligas foram caracterizadas termicamente (coeficiente de transferência de calor metal/molde, Velocidade de deslocamento da isoterma liquidus e taxa de resfriamento), mecanicamente (limite de resistência à tração e módulo de tenacidade) e eletricamente (resistência, resistividade e condutividade elétrica) utilizando duas formas de vazamentos: uma sendo através de solidificação em molde unidirecional horizontal e a outra através do molde em “U”. Para a solidificação em molde unidirecional horizontal as caracterizações mostraram-se bastante influenciadas pelos refinadores de grão e pela viscosidade que o cobre exerce na eficiência do contato metal/molde das ligas, alterando significativamente suas propriedades térmicas, mecânicas e elétricas. Por outro lado, na solidificação em molde “U” as ligas foram avaliadas apenas através de caracterização mecânica e elétrica, tendo em vista a preocupação sobre as variáveis do processo de conformação da liga metálica, tais como, a avaliação da taxa de deformação, do limite de resistência a tração, da tenacidade e da resistência, resistividade e condutividade elétrica, do perfil metálico produzido, para corpos de prova de diferentes diâmetros. Para o molde em “U” houve ainda a adição e avaliação de mais duas ligas sendo a Al-0,6Mg-0,4Si-0,1Cu e Al-0,6Mg-0,4Si-0,2Cu. As propriedades mecânicas e elétricas sofreram forte influência do encruamento, do teor de cobre e das microcavidades existentes nas estruturas das ligas.

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

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

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

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

Modelagem 2,5D dos campos usados no Método Eletromagnético a Multi-Frequência - EMMF

Esta tese mostra a modelagem 2,5D de dados sintéticos do Método Eletromagnético a Multi-frequência (EMMF). O trabalho é apresentado em duas partes: a primeira apresenta os detalhes dos métodos usados nos cálculos dos campos gerados por uma bobina horizontal de corrente colocada sobre a superfície de modelos bidimensionais; e a segunda, usa os resultados obtidos para simular os dados medidos no método EMMF, que são as partes real e imaginária da componente radial do campo magnético gerado pela bobina. Nesta segunda parte, observamos o comportamento do campo calculado em diversos modelos, incluindo variações nas propriedades físicas e na geometria dos mesmos, com o intuito de verificar a sensibilidade do campo observado com relação às estruturas presentes em uma bacia sedimentar. Com esta modelagem, podemos observar as características dos dados e como as duas partes, real e imaginária, contribuem com informações distintas e complementares. Os resultados mostram que os dados da componente radial do campo magnético apresentam muito boa resolução lateral, mesmo estando a fonte fixa em uma única posição. A capacidade desses dados em distinguir e resolver estruturas alvo será fundamental para o trabalho futuro de inversão, bem como para a construção de seções de resistividade aparente.

Modelagem numérica da influência do eletrojato equatorial em dados magnetotelúricos produzidos por estruturas tridimensionais

A América do Sul apresenta várias peculiaridades geomagnéticas, uma delas, é a presença do Eletrojato Equatorial, o qual se estende de leste para oeste no Brasil ao longo de aproximadamente 3500 km. Considerando-se o fato de que a influência do Eletrojato Equatorial pode ser detectada a grandes distâncias do seu centro, isto suscita o interesse em se estudar os seus efeitos na exploração magnetotelúrica no Brasil. A influência do eletrojato equatorial na prospecção magnetotelúrica tem sido modelada para meios geológicos uni e bidimensionais valendo-se para isto de soluções analíticas fechadas e de técnicas numéricas tais como elementos finitos e diferenças finitas. Em relação aos meios geológicos tridimensionais, eles tem sido modelados na forma de "camadas finas", usando o algoritmo "thin sheet". As fontes indutoras utilizadas para simular o eletrojato equatorial nestes trabalhos, tem sido linhas de corrente, eletrojatos gaussianos e eletrojatos ondulantes. Por outro lado, o objetivo principal da nossa tese foi o modelamento dos efeitos que o eletrojato equatorial provoca em estruturas tridimensionais próprias da geofísica da prospecção. Com tal finalidade, utilizamos o esquema numérico da equação integral, com as fontes indutoras antes mencionadas. De maneira similar aos trabalhos anteriores, os nossos resultados mostram que a influência do eletrojato equatorial somente acontece em frequências menores que 10^-1 Hz. Este efeito decresce com a distância, mantendo-se até uns 3000 km do centro do eletrojato. Assim sendo, a presença de grandes picos nos perfis da resistividade aparente de um semi-espaço homogêneo, indica que a influência do eletrojato é notável neste tipo de meio. Estes picos se mostram com diferente magnitude para cada eletrojato simulado, sendo que a sua localização também muda de um eletrojato para outro. Entretanto, quando se utilizam modelos geo-elétricos unidimensionais mais de acordo com a realidade, tais como os meios estratificados, percebe-se que a resposta dos eletrojatos se amortece significativamente e não mostra muitas diferenças entre os diferentes tipos de eletrojato. Isto acontece por causa da dissipação da energia eletromagnética devido à presença da estratificação e de camadas condutivas. Dentro do intervalo de 3000 km, a resposta eletromagnética tridimensional pode ser deslocada para cima ou para baixo da resposta da onda plana, dependendo da localização do corpo, da frequência, do tipo de eletrojato e do meio geológico. Quando a resposta aparece deslocada para cima, existe um afastamento entre as sondagens uni e tridimensionais devidas ao eletrojato, assim como um alargamento da anomalia dos perfis que registra a presença da heterogeneidade tridimensional. Quando a resposta aparece deslocada para baixo, no entanto, há uma aproximação entre estes dois tipos de sondagens e um estreitamento da anomalia dos perfis. Por outro lado, a fase se mostra geralmente, de uma forma invertida em relação à resistividade aparente. Isto significa que quando uma sobe a outra desce, e vice-versa. Da mesma forma, comumente nas altas frequências as respostas uni e tridimensionais aparecem deslocadas, enquanto que nas baixas frequências se mostram com os mesmos valores, com exceção dos eletrojatos ondulantes com parâmetros de ondulação α = —2 e —3. Nossos resultados também mostram que características geométricas próprias das estruturas tridimensionais, tais como sua orientação em relação à direção do eletrojato e a dimensão da sua direção principal, afetam a resposta devido ao eletrojato em comparação com os resultados da onda plana. Desta forma, quando a estrutura tridimensional é rotacionada de 90°, em relação à direção do eletrojato e em torno do eixo z, existe uma troca de polarizações nas resistividades dos resultados, mas não existem mudanças nos valores da resistividade aparente no centro da estrutura. Ao redor da mesma, porém, se percebe facilmente alterações nos contornos dos mapas de resistividade aparente, ao serem comparadas com os mapas da estrutura na sua posição original. Isto se deve à persistência dos efeitos galvânicos no centro da estrutura e à presença de efeitos indutivos ao redor do corpo tridimensional. Ao alongar a direção principal da estrutura tridimensional, as sondagens magnetotelúricas vão se aproximando das sondagens das estruturas bidimensionais, principalmente na polarização XY. Mesmo assim, as respostas dos modelos testados estão muito longe de se considerar próximas das respostas de estruturas quase-bidimensionais. Porém, os efeitos do eletrojato em estruturas com direção principal alongada, são muito parecidos com aqueles presentes nas estruturas menores, considerando-se as diferenças entre as sondagens de ambos tipos de estruturas. Por outro lado, os mapas de resistividade aparente deste tipo de estrutura alongada, revelam um grande aumento nos extremos da estrutura, tanto para a onda plana como para o eletrojato. Este efeito é causado pelo acanalamento das correntes ao longo da direção principal da estrutura. O modelamento de estruturas geológicas da Bacia de Marajó confirma que os efeitos do eletrojato podem ser detetados em estruturas pequenas do tipo "horst" ou "graben", a grandes distâncias do centro do mesmo. Assim, os efeitos do eletrojato podem ser percebidos tanto nos meios estratificados como tridimensionais, em duas faixas de freqüência (nas proximidades de 10^-1 Hz e para freqüências menores que 10^-3 Hz), possivelmente influenciados pela presença do embasamento cristalino e a crosta inferior, respectivamente. Desta maneira, os resultados utilizando o eletrojato como fonte indutora, mostram que nas baixas freqüências as sondagens magnetotelúricas podem ser fortemente distorcidas, tanto pelos efeitos galvânicos da estrutura tridimensional como pela presença da influência do eletrojato. Conseqüêntemente, interpretações errôneas dos dados de campo podem ser cometidas, se não se corrigirem os efeitos do eletrojato equatorial ou, da mesma forma, não se utilisarem algoritmos tridimensionais para interpretar os dados, no lugar do usual modelo unidimensional de Tikhonov - Cagniard.

Análise da resolução do modelo de Hughes em sondagens magnetotelúricas

Na maioria dos métodos de exploração geofísica, a interpretação é feita assumindo-se um modelo da Terra uniformemente estratificado. Todos os métodos de inversão, inclusive o de dados eletromagnéticos, exigem técnica de modelamento teórico de modo a auxiliar a interpretação. Na literatura os dados são geralmente interpretados em termos de uma estrutura condutiva unidimensional; comumente a Terra é assumida ser horizontalmente uniforme de modo que a condutividade é função somente da profundidade. Neste trabalho uma técnica semi-analítica de modelagem desenvolvida por Hughes (1973) foi usada para modelar a resposta magnética de duas camadas na qual a interface separando as camadas pode ser representada por uma expansão em série de Fourier. A técnica envolve um método de perturbação para encontrar o efeito de um contorno senoidal com pequenas ondulações. Como a perturbação é de primeira ordem a solução obtida é linear, podemos então usar o princípio da superposição e combinar soluções para várias senoides de forma a obter a solução para qualquer dupla camada expandida em série de Fourier. Da comparação com a técnica de elementos finitos, as seguintes conclusões podem ser tiradas: • Para um modelo de dupla camada da Terra, as camadas separadas por uma interface cuja profundidade varia senoidalmente em uma direção, as respostas eletromagnética são muito mais fortes quando a espessura da primeira camada é da ordem do skin depth da onda eletromagnética no meio, e será tanto maior quanto maior for o contraste de condutividade entre as camadas; • Por outro lado, a resistividade aparente para este modelo não é afetada pela mudança na frequência espacial (v) do contorno; • Em caso do uso da solução geral para qualquer dupla camada na Terra cuja interface possa ser desenvolvida em série de Fourier, esta técnica produziu bons resultados quando comparado com a técnica de elementos finitos. A linerização restringe a aplicação da técnica para pequenas estruturas, apesar disso, uma grande quantidade de estruturas pode ser modelada de modo simples e com tempo computacional bastante rápido; • Quando a dimensão da primeira camada possui a mesma ordem de grandeza da estrutura, esta técnica não é recomendada, porque para algumas posições de sondagem, as curvas de resistividade aparente obtidas mostram um pequeno deslocamento quando comparadas com as curvas obtidas por elementos finitos.

Análise teórica do problema de Weaver da falha infinita, modo TE magnetotelúrico

É apresentada uma solução totalmente analítica do modelo da falha infinita para o modo TE magnetotelúrico, levando em conta a presença do ar, utilizando como base o trabalho de Sampaio apresentado em 1985, que apresenta uma solução parcialmente analítica e parcialmente numérica – solução híbrida. Naquela solução foram aplicadas oito condições de contorno, sendo que em quatro delas foram encontradas inconsistências matemáticas que foram dirimidas com alterações adequadas nas soluções propostas por Sampaio. Tais alterações propiciaram a chegarse à solução totalmente analítica aqui apresentada. A solução obtida foi comparada com a solução de Weaver, com a de Sampaio e com o resultado do método numérico dos elementos finitos para contrastes de resistividade iguais a 2, 10 e 50. A comparação da solução analítica, para o campo elétrico normalizado, com a solução de elementos finitos mostra que a solução analítica proporcionou resultados mais próximos, em comparação aos fornecidos por Weaver e por Sampaio. Este é um problema muito difícil, aberto para uma solução analítica definitiva. A solução apresentada aqui é, nesta direção, um grande passo.

Aplicação de análise de séries temporais em perfis de poço

Esta tese descreve a aplicação de análise de séries temporais em perfis de poço. Através desta técnica é possível avaliar-se a repetição e a resolução vertical de perfis, e determinar-se o intervalo de amostragem e a velocidade de perfilagem ideais para diferentes perfis. A comparação entre três poços é também feita, baseada num mesmo tipo de perfil. Para tanto, na seqüência utilizada, procurou-se manter num mesmo domínio os dados cuja quantidade total de amostras (N) por perfil não ultrapassou 2048. Desses dados, foram inicialmente retirados o valor médio das amostras e o alinhamento polinomial algébrico eventualmente nelas embutido. Em seguida, foi efetuada a aplicação do ponderador cossenoidal, do filtro passa-alta, da janela Hanning, do cálculo da função coerência, do espectro de fase, da razão sinal-ruído e dos espectros de potência do sinal e do ruído, nesta ordem. Para a função coerência, fez-se necessário o cálculo dos níveis de confiança de 50%, 90% e 95%. O cálculo do primeiro nível teve por base a necessidade de se determinar a resolução vertical de alguns perfis, e dos demais, a fim de que fosse obtida uma informação referente à localização daqueles níveis para a coerência calculada. Em relação ao espectro de fase, seu cálculo surgiu da necessidade de se obter uma informação adicional a respeito dos perfis manipulados, ou seja, o conhecimento da ocorrência ou não de deslocamento relativo de profundidade entre a seção principal e a seção repetida. A razão sinal-ruído foi calculada no sentido de possibilitar a comparação, como elemento avaliador dos diversos tipos de perfis, com a coerência e o cálculo dos espectros de potência. Os espectros de potência do sinal e do ruído foram calculados para se ter mais um parâmetro de avaliação da seção repetida, já que em tese, os espectros de potência do sinal e do ruído da seção repetida devem ser iguais aos respectivos espectros da seção principal. Os dados utilizados na aplicação da metodologia proposta foram fornecidos pela PETROBRÁS e oriundos de quatro poços da Bacia Potiguar emersa. Por questões de sigilo empresarial, os poços foram identificados como poços A, B, C e D. A avaliação da repetição entre diferentes tipos de perfis indica que, para o poço A, o perfil micro-esférico (MSFL) tem repetição melhor do que o perfil de porosidade neutrônica (CNL), o qual tem, por sua vez, repetição melhor do que o perfil de raios gama normal (GR). Para os perfis do poço D, uma diminuição da velocidade de perfilagem de 550 m/h para 275 m/h é vantajosa apenas para o perfil de porosidade neutrônica. Já a velocidade de perfilagem de 920 m/h, utilizada, na obtenção dos perfis do poço C, é totalmente inadequada para os perfis de resistividade (MSFL, ILD e ILM). A diminuição do intervalo de amostragem de 0,20 m para 0,0508 m, nos perfis de raios gama e de porosidade neutrônica, e 0,0254 m para o perfil de densidade, apresenta bons resultados quando aplicada no poço D. O cálculo da resolução vertical indica, para o perfil de porosidade neutrônica, uma superioridade qualitativa em relação ao perfil de raios gama normal, ambos pertencentes ao poço A. Para o poço C, o perfil micro-esférico apresenta uma resolução vertical na mesma ordem de grandeza da resolução do perfil de raios gama do poço B, o que evidencia ainda mais a inconveniência da velocidade de perfilagem utilizada no poço C. Já para o poço D, o cálculo da resolução vertical indica uma superioridade qualitativa do perfil de densidade de alta resolução em relação ao perfil de raios gama de alta resolução. A comparação entre os poços A, B e D, levada a efeito através dos respectivos perfis de porosidade neutrônica normais, comprova que a presença de ruído aleatório, em geral, está diretamente ligada à porosidade da formação - uma maior porosidade indica uma presença maior de ruído e, por conseguinte, uma queda qualitativa no perfil obtido. A análise do espectro de fase de cada perfil indica um deslocamento em profundidade, existente entre as seções principal e repetida de todos os perfis do poço C. E isto pôde ser confirmado com a posterior superposição das seções.