Morfotectônica da região centro-leste do Estado de São Paulo e áreas adjacentes de Minas Gerais: termocronologia e paleotensões

Pós-graduação em Geologia Regional - IGCE

Neotectônica do baixo vale do rio Jaguaribe - Ceará

Pós-graduação em Geologia Regional - IGCE

Avaliação estocástica de afundamentos de tensão através de análise de curtos-circuitos deslizantes e banco de dados relacional

Atualmente, é grande a preocupação com relação aos distúrbios que afetam a qualidade da energia elétrica, em especial os afundamentos de tensão, devido aos prejuízos associados aos usuários. Isto decorre da crescente presença de cargas sensíveis a estes distúrbios, principalmente nos processos industriais. Para se avaliar corretamente o problema, é fundamental dispor de estimativas confiáveis do número e das características dos afundamentos de tensão que podem afetar determinada carga ou processo. Dentro deste contexto, este trabalho propõe a implementação de uma metodologia de estimação dos afundamentos de tensão. Esta metodologia calcula a amplitude, duração e freqüência de ocorrência do fenômeno, tendo como base um estudo estatístico estocástico. Neste estudo a posição da falta ao longo das linhas é sorteada, a fim de melhor emular a aleatoriedade da ocorrência das faltas no sistema. Para avaliar a aplicabilidade da metodologia implementada, realizaram-se estudos de casos utilizando um sistema elétrico real pertencente a uma concessionária da região norte (CELPA). A análise permite afirmar que a metodologia implementada pode ser eficazmente utilizada para auxiliar engenheiros no planejamento de sistemas elétricos e no diagnóstico de problemas ocasionados por afundamentos de tensão.

Implementação e testes de campo de um ESP digital

Neste trabalho é apresentado o desenvolvimento e os resultados da implementação e testes em campo de um estabilizador de sistema de potência (ESP) projetado com técnica de controle digital para fins de amortecimento de modos de oscilação eletromecânica observáveis em sinais de potência elétrica medido em uma unidade hidro-geradora, de 350 MVA da Usina Hidrelétrica de Tucuruí. É apresentada e aplicada a metodologia de identificação de modelos paramétricos lineares do tipo auto regressivo com entradas exógenas (ARX), para estimação de modelos com capacidade de capturar a informação relevante (amortecimento e freqüência natural) dos modos eletromecânicos dominantes do sistema. De posse do modelo paramétrico ARX, é efetuada então a síntese da lei de controle digital amortecedor para o ESP, através da técnica de deslocamento radial dos pólos da função de transferência de malha fechada. Para a síntese da lei de controle digital, utilizou-se uma estrutura canônica do tipo RST. Para os testes de campo, a lei de controle amortecedor do ESP digital foi codificada em linguagem C e embarcada em um protótipo cujo hardware é baseado em microcontrolador modelo DSPIC 30F3014, o qual incorpora um grande número de periféricos para aquisição e comunicação de dados. Para avaliar o desempenho do ESP digital desenvolvido, testes experimentais foram realizados em uma unidade geradora de 350 MVA da casa de força número 1, da UHE de Tucuruí. O estabilizador desenvolvido atua através da modulação da referência de tensão do regulador automático de tensão da respectiva unidade geradora, de acordo com as oscilações observadas através da medida de potência elétrica no estator do gerador. Os resultados de testes de campo mostraram um excelente desempenho do ESP digital no amortecimento de um modo eletromecânico, de freqüência natural de aproximadamente 1,7 Hz, observado nos teste de campo realizado.

O estudo da influência do ciclo de partidas e paradas no estator de geradores síncronos

Os geradores síncronos se submetem a diferentes situações de paradas para desligamento forçado, desligamento por condições externas, desligamento programado para manutenção corretiva e preventiva. As diferentes variações de temperatura em que se submete o gerador síncrono, precisamente, ao seu estator, afeta o seu isolamento, que não acompanha essas constantes variações, devido ter um coeficiente de dilatação diferente. No presente trabalho busca-se identificar uma tendência de progressão do envelhecimento do isolamento de seus enrolamentos estatóricos em máquinas que sofrem repetitivas partidas e paradas e que ficam com sua reserva de energia desligada por muitas horas. O estudo foi aplicado nos cinco geradores síncronos da usina hidrelétrica de Balbina, pertencente ao grupo Eletrobrás, que apresenta tal problema. Esse estudo vem comprovar os ensaios de análise de descargas parciais dos geradores síncronos realizados periodicamente para diagnosticar problemas nos seus enrolamentos estatóricos que apresentam diferentes magnitudes e confirmar a hipótese de que máquinas que apresentam maiores magnitudes de descargas parciais em seus estatores são as que se submetem maiores quantidades de partidas e paradas e horas com energia de reserva desligada.

Abordagem probabilística para caracterização do sistema de marcação de sequenciamento multiplex na plataforma ABI SOLID

Os sequenciadores de nova geração como as plataformas Illumina e SOLiD geram uma grande quantidade de dados, comumente, acima de 10 Gigabytes de arquivos-texto. Particularmente, a plataforma SOLiD permite o sequenciamento de múltiplas amostras em uma única corrida (denominada de corrida multiplex) por meio de um sistema de marcação chamado Barcode. Esta funcionalidade requer um processo computacional para separação dos dados por amostra, pois, o sequenciador fornece a mistura de todas amostras em uma única saída. Este processo deve ser seguro a fim de evitar eventuais embaralhamentos que possam prejudicar as análises posteriores. Neste contexto, o presente trabalho propõe desenvolvimento de um modelo probabilístico capaz de caracterizar sistema de marcação utilizado em sequenciamentos multiplex. Os resultados obtidos corroboraram a suficiência do modelo obtido, o qual permite, dentre outras coisas, identificar faltas em algum passo do processo de sequenciamento; adaptar e desenvolver de novos protocolos para preparação de amostras, além de atribuir um Grau de Confiança aos dados gerados e guiar um processo de filtragem que respeite as características de cada sequenciamento, não descartando sequências úteis de forma arbitrária.

Influências de famílias wavelets e suas ordens no desempenho de um localizador de faltas em linhas aéreas de transmissão de energia elétrica

Essa dissertação tem por objetivo analisar a influência de famílias wavelets e suas ordens no desempenho de um algoritmo de localização de faltas a partir das ondas viajantes de dois terminais de uma linha de transmissão aérea. Tornou-se objetivo secundário a modelagem de um sistema elétrico de potência (SEP) para obtenção de um universo de faltas que validassem o localizador. Para isso, parte de um SEP da Eletrobrás-Eletronorte em 500/230 kV foi modelado no Alternative Transient Program (ATP) utilizando-se parâmetros reais. A Transformada Wavelet, via análise multiresolução (AMR), é empregada valendo-se de sua característica de localização temporal, permitindo caracterizações precisas de instantes de transitórios eletromagnéticos ocasionados por faltas, as quais geram ondas que ao se propagarem em direção aos terminais da linha contêm os tempos de propagação destas do local do defeito a tais terminais e podem ser convenientemente extraídos por tal transformada. Pela metodologia adotada no algoritmo, a diferença entre esses tempos determina com boa exatidão o local de ocorrência da falta sobre a linha. Entretanto, um dos agentes variantes do erro nessa estimação é a escolha da Wavelet usada na AMR dos sinais, sendo, portanto, a avaliação dessa escolha sobre o erro, objetivo principal do trabalho, justificada pela ainda inexistente fundamentação científica que garanta a escolha de uma wavelet ótima a uma certa aplicação. Dentre um leque de Wavelets discretas, obtiveram-se resultados adequados para 16 delas, havendo erros máximos inferiores aos 250 metros estipulados para a precisão. Duas Wavelets, a Db15 e a Sym17, sobressaíram-se ao errarem, respectivamente, 3,5 e 1,1 vezes menos que as demais. A metodologia empregada consta da: exportação dos dados das faltas do ATP para o MATLAB®; aplicação da transformação modal de Clarke; decomposição dos modos alfa e síntese dos níveis 1 de detalhes via AMR; cálculo de suas máximas magnitudes e determinação dos índices temporais; e por fim, a teoria das ondas viajantes equaciona e estima o local do defeito sobre a LT, sendo tudo isso programado no MATLAB e os erros de localização analisados estatisticamente no Microsoft Excell®. Ao final elaborou-se ainda uma GUI (Guide User Interface) para a Interface Homem-Máquina (IHM) do localizador, servindo também para análises gráficas de qualquer das contingências aplicadas ao SEP. Os resultados alcançados demonstram uma otimização de performance em razão da escolha da wavelet mais adequada ao algoritmo e norteiam para uma aplicação prática do localizador.

Estudos de estratégias de identificação paramétrica para detecção e diagnóstico de faltas em um processo industrial do tipo tanques comunicantes

Esta dissertação apresenta uma técnica para detecção e diagnósticos de faltas incipientes. Tais faltas provocam mudanças no comportamento do sistema sob investigação, o que se reflete em alterações nos valores dos parâmetros do seu modelo matemático representativo. Como plataforma de testes, foi elaborado um modelo de um sistema industrial em ambiente computacional Matlab/Simulink, o qual consiste em uma planta dinâmica composta de dois tanques comunicantes entre si. A modelagem dessa planta foi realizada através das equações físicas que descrevem a dinâmica do sistema. A falta, a que o sistema foi submetido, representa um estrangulamento gradual na tubulação de saída de um dos tanques. Esse estrangulamento provoca uma redução lenta, de até 20 %, na seção desse tubo. A técnica de detecção de falta foi realizada através da estimação em tempo real dos parâmetros de modelos Auto-regressivos com Entradas Exógenas (ARX) com estimadores Fuzzy e de Mínimos Quadrados Recursivos. Já, o diagnóstico do percentual de entupimento da tubulação foi obtido por um sistema fuzzy de rastreamento de parâmetro, realimentado pela integral do resíduo de detecção. Ao utilizar essa metodologia, foi possível detectar e diagnosticar a falta simulada em três pontos de operação diferentes do sistema. Em ambas as técnicas testadas, o método de MQR teve um bom desempenho, apenas para detectar a falta. Já, o método que utilizou estimação com supervisão fuzzy obteve melhor desempenho, em detectar e diagnosticar as faltas aplicadas ao sistema, constatando a proposta do trabalho.

Técnica de projeto de controlador fuzzy aplicada ao acionamento vetorial de motor de indução

Este trabalho investiga uma estratégia de controle fuzzy Takagi-Sugeno aplicada ao controle de velocidade do motor de indução. A estratégia implementa uma interpolação ponderada entre um conjunto de controladores locais previamente projetados. Ao ocorrer variações nas condições operacionais do motor de indução, os ganhos da lei de controle são ajustados automaticamente, de modo a manter satisfatório o desempenho do sistema de controle. Para o projeto do controlador fuzzy a representação em espaço de estados da planta foi considerada sob a forma de um sistema aumentado, incluindo-se uma nova variável de estado que, nesse caso, foi selecionada como sendo a integral do erro de velocidade. Tal formulação permitiu o projeto de controladores locais com a estrutura PI, através de realimentação completa de estados, com posicionamento de pólos. Como variáveis de operação para o chaveamento fuzzy dos controladores locais, foram selecionados as variáveis velocidade angular do rotor e a componente da corrente de estator responsável pelo torque elétrico do motor. Em seguida, a estabilidade do controlador fuzzy Takagi- Sugeno projetado foi comprovada através do critério de Lyapunov, para isso o problema de estabilidade foi escrito na forma de LMIs. O desempenho do controlador fuzzy Takagi-Sugeno foi avaliado através de estudos de simulação, e seus resultados comparados ao desempenho de um controlador PI convencional, para a regulação da velocidade do rotor. Os resultados obtidos nas simulações mostram que o emprego da estratégia proposta torna o sistema mais robusto a variações paramétricas no sistema de acionamento.

Tectonics and paleogeography of the Marajó Basin, northern Brazil

ABSTRACT: The Marajó Basin area presents geologic and geomorphologic features chiefly due to the Mesozoic extension and post-Miocene neotectonics. The extension event with an Early and a Late Cretaceous phases originated four sub-basins that constitutes the Marajó Basin, with a thick continental clastic sequence showing marine influence. NW and NNW normal faults and NE and ENE strike-slip faults controlled the basin geometry. The extension, related to the Equatorial Atlantic opening, propagated into the continent along crustal weakness zones of the Precambrian Tumucumaque, Amapá and Araguaia orogenic belts. The neotectonic event is a strike-slip regime which developed transtensional basins filled in by Upper Tertiary shallow marine (Pirabas Formation) and transitional sequences (Barreiras Group), followed by Quaternary fluvial deposits and transitional sequences derived from the Amazon and Tocantins rivers and the Marajoara estuary. The current landscape has a typical estuarine morphology. The coast morphology presents sea-cliffs on transitional Upper Tertiary sequences, while inwards dominate hills sustained by Mid-Pleistocene lateritic crust, with a flat erosive surface at 70 m. In the eastern Marajó Island several generations of paleochannels associated with fluvial-estuarine sequences are recognized, while a fluvial-marine plain is widespread on its western side.

Controladores fuzzy aplicados ao conversor de geradores de indução duplamente excitados em sistemas eólicos integrados a sistemas de potência

Este artigo apresenta novas estratégias de controle fuzzy aplicadas ao conversor estático interligado ao rotor dos geradores de indução duplamente excitados (DFIG), em esquemas eólicos de velocidade variável, integrados a redes elétricas reais. As estratégias de controle propostas são do tipo fuzzy "look-up-table" com supervisão. O desempenho de tais reguladores, baseados em inteligência computacional, é comparado com o obtido com controladores PI's a parâmetros fixos, na ocorrência de faltas no sistema elétrico de potência. Esses controladores fuzzy são destinados a garantir a continuidade da operação dos conversores, e melhorar o desempenho transitório do sistema, em relação aos controladores convencionais. Os resultados apresentados confirmam a eficácia das estratégias de controle adotadas. O modelo físico dos aerogeradores, consistindo de um grande número de turbinas eólicas, foi implementado através de uma ferramenta de simulação dinâmica, desenvolvida no ambiente computacional MATLABTM.

Extração de conhecimento em forma de regras difusas a partir de mapas auto-organizáveis de Kohonen: aplicação em diagnóstico de faltas incipientes em transformadores

Apesar das diversas vantagens oferecidas pelas redes neurais artificiais (RNAs), algumas limitações ainda impedem sua larga utilização, principalmente em aplicações que necessitem de tomada de decisões essenciais para garantir a segurança em ambientes como, por exemplo, em Sistemas de Energia. Uma das principais limitações das RNAs diz respeito à incapacidade que estas redes apresentam de explicar como chegam a determinadas decisões; explicação esta que seja humanamente compreensível. Desta forma, este trabalho propõe um método para extração de regras a partir do mapa auto-organizável de Kohonen, projetando um sistema de inferência difusa capaz de explicar as decisões/classificação obtidas através do mapa. A metodologia proposta é aplicada ao problema de diagnóstico de faltas incipientes em transformadores, em que se obtém um sistema classificatório eficiente e com capacidade de explicação em relação aos resultados obtidos, o que gera mais confiança aos especialistas da área na hora de tomar decisões.

Estudo experimental de técnicas de medição de descargas parciais nos enrolamentos estatóricos de hidrogeradores

Este trabalho tem como objetivo geral analisar experimentalmente técnicas para medição de descargas parciais nos enrolamentos estatóricos de hidrogeradores. Para atingir o objetivo proposto, primeiramente foi realizado um estudo bibliográfico sobre como é a estrutura dos enrolamentos do estator, com ênfase na isolação elétrica, assim como os conceitos de descargas parciais e técnicas empregadas para medição em hidrogeradores. Em seguida, foi montado em laboratório um arranjo experimental para o estudo das técnicas de medição de descargas parciais nos enrolamentos do estator utilizando seis acopladores diferentes, sendo quatro acopladores capacitivos convencionais de 80, 220, 500 e 1000 pF, e dois sistemas experimentais, uma antena SSC fabricada com microfita e um cabo coaxial RG 58, desencapado, utilizado como antena. A avaliação da sensibilidade destes seis sistemas foi divida em duas etapas: primeiro foi analisada a resposta dos acopladores à aplicação de um sinal conhecido de calibração e, em seguida, foi analisada a resposta dos mesmos durante ensaios com alta tensão. Para aplicação de um sinal conhecido de calibração os acopladores capacitivos apresentaram maior sensibilidade em relação à antena e ao cabo coaxial. Já a resposta durante os ensaios com alta tensão mostrou que a sensibilidade da antena é maior que a do método convencional de medição de descargas parciais em hidrogeradores, baseado em acopladores capacitivos de 80 e 220 pF. Os resultados obtidos em laboratório mostraram que a medição de descargas parciais através da antena pode ser realizada de forma satisfatória na avaliação da condição operacional dos enrolamentos do estator de hidrogeradores, sem perda significativa de informação e com aumento significativo de sensibilidade na quantificação do fenômeno, cerca de duas vezes maior que a sensibilidade do acoplador de 220 pF e quatro vezes a do acoplador de 80 pF.

Estudo da geometria e cinemática das rochas sedimentares arqueanas da mina do Igarapé do Azul – Carajás-PA

Mina de Manganês do Igarapé Azul posiciona-se geologicamente no interior do feixe da Falha Carajás, na porção central do Sistema Transcorrente de Carajás. O depósito do Manganês do Azul relaciona-se a rochas sedimentares pelíticas do Membro Azul, na base da Formação Águas Claras (Arqueano), em contato discordante, acima do Grupo Grão Pará (Nogueira et al, 1995). Três frentes de lavra a céu aberto estão atualmente em andamento na área: (1) Mina Principal (Mina 1), (2) Mina 2 e (3) Mina 3. Nestes locais encontram-se excelentes afloramentos de siltitos intercalados com argilitos e arenitos finos, intercalados com níveis manganesíferos. Essas rochas estão organizadas em conjuntos de dobras e falhas normais e inversas sob deformação heterogênea, particionada em diferentes escalas. As seções geológicas realizadas nas frentes de lavra mostram a predominância de siltitos intercalados com argilitos em contato com rochas pelíticas manganesíferas e minério (bióxido de Mn). Nessas rochas são comuns estruturas primárias tipo hummocky, estratificações cruzadas, e laminações plano-paralelas. O acamamento centimétrico a métrico (em média de 30 a 50 cm ) representa a principal estrutura primária, usada como marcador de deformação, observada nas rochas. A Mina do Igarapé Azul encontra-se dividida em dois blocos, separados por falha normal com rejeito de até dezenas de metros, com o bloco norte alto em relação ao bloco sul. O bloco sul encontra-se pouco deformado, apresentado uma regularidade no acamamento que mergulha com ângulos suaves para sul, colocando a camada de minério sucessivamente em níveis mais profundos na direção S. No bloco norte o acamamento apresenta um comportamento heterogêneo. A deformação é mais expressiva nessa região, estando o nível de minério deformado por dobras e falhas inversas. Além da cinemática vertical, as falhas apresentam deslocamento conjugado destral dando a essas feições um caráter oblíquo. Essa região pode ser definida como um corredor de deformação. O corredor observado no bloco norte, de acordo com os domínios principais separados pelas falhas anteriormente descritas, possui orientação NW-SE, com aproximadamente um quilômetro de extensão, sendo caracterizado por dobras assimétricas curvilineares com eixos de mergulhos suaves (10° a 25°) para NW e SE. Essas dobras são seccionadas por falhas normais sinuosas NW-SE e/ou E-W, com baixo ângulo de mergulho (em torno de 10° a 30°), subordinadas a transcorrências destrais, gerando em escala de detalhes feições como drag folds. Observam-se ainda falhas inversas retas e/ou sinuosas NW-SE e zonas de falhas sub-verticalizadas WNW -ESE. As dobras individuais nesta área são estruturas do tipo reversas, flexurais e com geometria en echelon com orientação semelhante às dobras curvilineares: eixos com baixo ângulo de mergulho (10° a 25°) e caimento para SE. O conjunto de feições anteriormente descrito desenha, em escala quilométrica, um antiforme aberto, provavelmente resultante da acomodação do acamamento em resposta a deformação dessas falhas. O paralelismo entre feições observadas na área da Mina do Igarapé Azul e os lineamentos maiores que desenham a Falha Carajás em planta sugere uma relação com dois importantes episódios deformacionais ocorridos durante a história tectônica da Falha Carajás. As falhas normais associadas a componente direcional destral, de maior expressividade da área da mina, estariam relacionadas ao episódio de transtensão destral responsável pela instalação da Falha Carajás anterior a 2.6 Ga (Pinheiro, 1997). As dobras, as falhas de cavalgamento e as zonas de falhas sub-verticalizadas estariam relacionados a deformações sob regime de transpressão sinistral, um segundo evento atuante na região, responsável pela reativação e inversão da maioria das estruturas próximas à zona da Falha Carajás (Pinheiro, 1997; Pinheiro & Holdsworth, 2000; Lima, 2002).

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.