Mapeamento do relevo do embasamento de bacias sedimentares através da inversão gravimétrica vinculada

Apresentamos três novos métodos estáveis de inversão gravimétrica para estimar o relevo de uma interface arbitrária separando dois meios. Para a garantia da estabilidade da solução, introduzimos informações a priori sobre a interface a ser mapeada, através da minimização de um (ou mais) funcional estabilizante. Portanto, estes três métodos se diferenciam pelos tipos de informação físico-geológica incorporados. No primeiro método, denominado suavidade global, as profundidades da interface são estimadas em pontos discretos, presumindo-se o conhecimento a priori sobre o contraste de densidade entre os meios. Para a estabilização do problema inverso introduzimos dois vínculos: (a) proximidade entre as profundidades estimadas e verdadeiras da interface em alguns pontos fornecidas por furos de sondagem; e (b) proximidade entre as profundidades estimadas em pontos adjacentes. A combinação destes dois vínculos impõe uma suavidade uniforme a toda interface estimada, minimizando, simultaneamente em alguns pontos, os desajustes entre as profundidades conhecidas pelas sondagens e as estimadas nos mesmos pontos. O segundo método, denominado suavidade ponderada, estima as profundidades da interface em pontos discretos, admitindo o conhecimento a priori do contraste de densidade. Neste método, incorpora-se a informação geológica que a interface é suave, exceto em regiões de descontinuidades produzidas por falhas, ou seja, a interface é predominantemente suave porém localmente descontínua. Para a incorporação desta informação, desenvolvemos um processo iterativo em que três tipos de vínculos são impostos aos parâmetros: (a) ponderação da proximidade entre as profundidades estimadas em pontos adjacentes; (b) limites inferior e superior para as profundidades; e (c) proximidade entre todas as profundidades estimadas e um valor numérico conhecido. Inicializando com a solução estimada pelo método da suavidade global, este segundo método, iterativamente, acentua as feições geométricas presentes na solução inicial; ou seja, regiões suaves da interface tendem a tornar-se mais suaves e regiões abruptas tendem a tornar-se mais abruptas. Para tanto, este método atribui diferentes pesos ao vínculo de proximidade entre as profundidades adjacentes. Estes pesos são automaticamente atualizados de modo a acentuar as descontinuidades sutilmente detectadas pela solução da suavidade global. Os vínculos (b) e (c) são usados para compensar a perda da estabilidade, devida à introdução de pesos próximos a zero em alguns dos vínculos de proximidade entre parâmetros adjacentes, e incorporar a informação a priori que a região mais profunda da interface apresenta-se plana e horizontal. O vínculo (b) impõe, de modo estrito, que qualquer profundidade estimada é não negativa e menor que o valor de máxima profundidade da interface conhecido a priori; o vínculo (c) impõe que todas as profundidades estimadas são próximas a um valor que deliberadamente viola a profundidade máxima da interface. O compromisso entre os vínculos conflitantes (b) e (c) resulta na tendenciosidade da solução final em acentuar descontinuidades verticais e apresentar uma estimativa suave e achatada da região mais profunda. O terceiro método, denominado mínimo momento de inércia, estima os contrastes de densidade de uma região da subsuperfície discretizada em volumes elementares prismáticos. Este método incorpora a informação geológica que a interface a ser mapeada delimita uma fonte anômala que apresenta dimensões horizontais maiores que sua maior dimensão vertical, com bordas mergulhando verticalmente ou em direção ao centro de massa e que toda a massa (ou deficiência de massa) anômala está concentrada, de modo compacto, em torno de um nível de referência. Conceitualmente, estas informações são introduzidas pela minimização do momento de inércia das fontes em relação ao nível de referência conhecido a priori. Esta minimização é efetuada em um subespaço de parâmetros consistindo de fontes compactas e apresentando bordas mergulhando verticalmente ou em direção ao centro de massa. Efetivamente, estas informações são introduzidas através de um processo iterativo inicializando com uma solução cujo momento de inércia é próximo a zero, acrescentando, em cada iteração, uma contribuição com mínimo momento de inércia em relação ao nível de referência, de modo que a nova estimativa obedeça a limites mínimo e máximo do contraste de densidade, e minimize, simultaneamente, os desajustes entre os dados gravimétricos observados e ajustados. Adicionalmente, o processo iterativo tende a "congelar" as estimativas em um dos limites (mínimo ou máximo). O resultado final é uma fonte anômala compactada em torno do nível de referência cuja distribuição de constraste de densidade tende ao limite superior (em valor absoluto) estabelecido a priori. Estes três métodos foram aplicados a dados sintéticos e reais produzidos pelo relevo do embasamento de bacias sedimentares. A suavidade global produziu uma boa reconstrução do arcabouço de bacias que violam a condição de suavidade, tanto em dados sintéticos como em dados da Bacia do Recôncavo. Este método, apresenta a menor resolução quando comparado com os outros dois métodos. A suavidade ponderada produziu uma melhoria na resolução de relevos de embasamentos que apresentam falhamentos com grandes rejeitos e altos ângulos de mergulho, indicando uma grande potencialidade na interpretação do arcabouço de bacias extensionais, como mostramos em testes com dados sintéticos e dados do Steptoe Valley, Nevada, EUA, e da Bacia do Recôncavo. No método do mínimo momento de inércia, tomou-se como nível de referência o nível médio do terreno. As aplicações a dados sintéticos e às anomalias Bouguer do Graben de San Jacinto, California, EUA, e da Bacia do Recôncavo mostraram que, em comparação com os métodos da suavidade global e ponderada, este método estima com excelente resolução falhamentos com pequenos rejeitos sem impor a restrição da interface apresentar poucas descontinuidades locais, como no método da suavidade ponderada.

Controle de vibração em reatores elétricos através de absorvedores dinâmicos viscoelásticos

Reatores elétricos utilizados em subestações de energia têm funções específicas na operacionalização e manutenção da funcionalidade destas instalações, sendo equipamentos aplicados ao controle de reativos, de correntes de curto circuito e surtos de manobra em sistemas de transmissão de energia elétrica. Detectou-se que alguns reatores em funcionamento apresentaram elevados níveis de vibração estrutural. Dois destes reatores operam na subestação da cidade de Rurópolis, situada a oeste do estado do Pará. Outro problema decorrente do elevado nível de vibração destas estruturas é a geração de pontos quentes no núcleo eletromagnético do reator, podendo gerar bolhas de gases altamente inflamáveis. A existência destes problemas operacionais pode ocasionar ao Sistema de Transmissão do Norte do Brasil interrupções na transmissão de energia elétrica. Isto geraria na região Norte e em outras interligadas (Nordeste, Sudeste e Centro-Oeste) elevados prejuízos sociais e financeiros. A aplicação de ferramentas para a minimização do problema vibro-acústico (redução dos níveis vibração e, por conseguinte, ruído) pode evitar danos estruturais e operacionais que resultem na sua indisponibilidade operativa reduzindo, conseqüentemente, perdas e a necessidade do desenvolvimento de equipamentos de custo elevado com baixos índices de vibração. No intuito de se reduzir tais níveis, propõe-se à aplicação de absorvedores dinâmicos na estrutura externa dos reatores. Uma vez dimensionados em quantidade, posição e parâmetros definidores (massa, rigidez e amortecimento) espera-se criar um sistema, que sintonizado com a freqüência da excitação de natureza eletromagnética, absorverá parte da energia vibratória do reator, minimizando-se, os elevados índices detectados e, por conseqüência, todos os problemas subseqüentes. São indicados compostos viscoelásticos nos absorvedores dinâmicos (com propriedades dinâmicas dependentes tanto da temperatura como da freqüência de trabalho), em virtude de estes materiais aumentarem o campo de ação dos absorvedores, adequando-se à elevada densidade modal verificada nas chaparias metálicas do reator.

Análise da potencialidade a condução do BDT através da densidade eletrônica de estados via tight binding e desenvolvimento de software (B3J)

Neste trabalho nos propomos a fazer um estudo acerca da potencialidade de condução eletrônica no polímero BDT (1,3-benzoditiol 4H-ciclopenta[2,1-b:3,4b’]). O estudo usual de polímeros conjugados é feito de modo a obter sua densidade de estados com diversos tipos e níveis de dopagem. O método de Huckel é o mais utilizado e se baseia na separabilidade das ligações sigma e pi que é possível quando a molécula estudada é plana. Os polímeros conjugados são em sua maioria planos e estão inseridos nesta aproximação. O monômero do BDT apresenta sua geometria fora do plano por apresentar ligações com orbitais sp3. Para contornar esse problema foi desenvolvido o programa B3J, que considera todos os orbitais de valencia (s, px, py e pz). O programa B3J calcula a densidade de estados de sistemas poliméricos. O estudo das bandas do BDT foi feito com este software. Calculamos a densidade de estados do sistema neutro e com diversos níveis de dopagem, com distribuição aleatória e ordenada dos defeitos, dopagem do tipo n e do tipo p. O comportamento do quadrado do coeficiente da expansão da função de onda foi obtido para polímeros de até 20 monômeros. Estes cálculos foram obtidos com geometrias dos métodos AM1 e PM3. Obtivemos os espectros de absorção de oligômeros a fim de inferir seu comportamento para um polímero. Foram utilizados cálculos de otimização de geometria através dos métodos semi-empíricos AM1 e PM3 e ZINDO/S e o método DFT. Em outro objetivo desta monografia há o estudo do aproveitamento de tetrâmeros de BDT como dispositivos eletrônicos. Tais oligômeros foram otimizados em diversos valores de potencial elétrico, com a inserção em suas cadeias de moléculas doadoras e aceitadoras para induzir um aumento no momento de dipolo da mesma.

Estudo de solvatocromismo em líquidos moleculares orgânicos via método seqüencial Monte Carlo/Mecânica Quântica

Utilizou-se o método seqüencial Monte Carlo / Mecânica Quântica para obterem-se os desvios de solvatocromismo e os momentos de dipolo dos sistemas de moléculas orgânicas: Uracil em meio aquoso, -Caroteno em Ácido Oléico, Ácido Ricinoléico em metanol e em Etanol e Ácido Oléico em metanol e em Etanol. As otimizações das geometrias e as distribuições de cargas foram obtidas através da Teoria do Funcional Densidade com o funcional B3LYP e os conjuntos de funções de base 6-31G(d) para todas as moléculas exceto para a água e Uracil, as quais, foram utilizadas o conjunto de funções de base 6-311++G(d,p). No tratamento clássico, Monte Carlo, aplicou-se o algoritmo Metropólis através do programa DICE. A separação de configurações estatisticamente relevantes para os cálculos das propriedades médias foi implementada com a utilização da função de auto-correlação calculada para cada sistema. A função de distribuição radial dos líquidos moleculares foi utilizada para a separação da primeira camada de solvatação, a qual, estabelece a principal interação entre soluto-solvente. As configurações relevantes da primeira camada de solvatação de cada sistema foram submetidas a cálculos quânticos a nível semi-empírico com o método ZINDO/S-CI. Os espectros de absorção foram obtidos para os solutos em fase gasosa e para os sistemas de líquidos moleculares comentados. Os momentos de dipolo elétrico dos mesmos também foram obtidos. Todas as bandas dos espectros de absorção dos sistemas tiveram um desvio para o azul, exceto a segunda banda do sistema de Beta-Caroteno em Ácido Oléico que apresentou um desvio para o vermelho. Os resultados encontrados apresentam-se em excelente concordância com os valores experimentais encontrados na literatura. Todos os sistemas tiveram aumento no momento de dipolo elétrico devido às moléculas dos solventes serem moléculas polares. Os sistemas de ácidos graxos em álcoois apresentaram resultados muito semelhantes, ou seja, os ácidos graxos mencionados possuem comportamentos espectroscópicos semelhantes submetidos aos mesmos solventes. As simulações através do método seqüencial Monte Carlo / Mecânica Quântica estudadas demonstraram que a metodologia é eficaz para a obtenção das propriedades espectroscópicas dos líquidos moleculares analisados.

Aspectos texturais do magmatismo e tramas da tectônica impostas ao Granito Chaval na Zona de Cisalhamento Santa Rosa, extremo Noroeste da Província Borborema

A região Noroeste da Província Borborema apresenta uma diversidade de corpos graníticos de natureza e evolução tectônica diversificadas, do Paleoproterozoico ao Paleozoico, com maior incidência relacionada ao Neoproterozoico e alojamento em diferentes fases da orogenia Brasiliana. Um desses exemplos é o Granito Chaval, que representa um batólito aflorante próximo à costa Atlântica do Ceará e Piauí, intrusivo em ortognaisses do Complexo Granja e supracrustais do Grupo Martinópole. Ele é, em parte, coberto por depósitos cenozoicos costeiros e rochas sedimentares paleozoicas da Bacia do Parnaíba. O Granito Chaval tem como característica marcante a textura porfirítica, destacando-se megacristais de microclina, em sienogranitos e monzogranitos, e outras feições texturais/estruturais de origem magmática, Essas permitiram interpretar sua evolução como de alojamento relativamente raso do plúton, conduzido por processos de cristalização fracionada, mistura de magmas com fluxo magmático e ação gravitacional em função da diferença de densidade do magma, levando à flutuação e ascensão de megacristais de microclina no magma residual, com alojamento de leucogranitos e pegmatitos nos estágios finais da evolução deste plutonismo. Por outro lado, em toda a metade Leste do plúton, encontra-se um rico acervo de estruturas tectógenas de cisalhamento, relacionada à implantação da Zona de Cisalhamento Transcorrente Santa Rosa, que levou a transformações tectonometamórficas superpostas às feições magmáticas, as quais atingiram condições metamórficas máximas na fácies anfibolito baixo. Cartograficamente, foram individualizados três domínios estruturais em que estão presentes uma gama de variações petroestruturais do Granito Chaval, sejam feições texturais/estruturais ígneas e tectônicas. As rochas plutônicas foram deformadas e modificadas progressivamente à medida que se dirige para Leste, no qual as rochas mudam-se para tonalidades mais escuras do cinza e os processos de cominuição e recristalização dinâmica reduzem, progressivamente, a granulação grossa desses granitos bem como o tamanho dos fenocristais para dimensões mais finas, mantendo-se suas características porfiroides. Desse modo, a trama milonítica se torna evidente, acentuando-se ao atingir a porção principal da Zona de Cisalhamento Transcorrente Santa Rosa. Como principais feições estruturais, destacam-se extinção ondulante forte; encurvamento e segmentação de cristais; geminação de deformação; rotação de cristais; microbudinagem; foliação anastomosada, inclusive S-C; lineação de estiramento; formas amendoadas de porfiroclastos, fitas e folhas de quartzo e recristalização. Os produtos desses processos de cisalhamento resultam na formação de protomilonitos, milonitos e ultramilonitos. Essas faixas miloníticas representam os locais de maior concentração da deformação, por isso é possível acompanhar progressivamente suas modificações texturais e mineralógicas, configurando uma sequência clássica de deformação progressiva heterogênea, por cisalhamento simples, em condições frágil-dúctil e dúctil. O alojamento do Granito Chaval aconteceu no final do Criogeniano (aproximadamente 630 Ma) e pode ser interpretado como magmatismo sin a tardi-tectônico em relação ao evento Brasiliano. O processo de cisalhamento que gerou a Zona de Cisalhamento Transcorrente Santa Rosa se formou nos incrementos finais da deformação de uma colisão continental em um sistema de cavalgamento oblíquo, em que se edificou o Cinturão de Cisalhamento Noroeste do Ceará, devido ao extravasamento lateral de massas crustais em fluxo dúctil acontecido no final da orogenia Brasiliana no Noroeste da Província Borborema.