Redução de ordem de modelos matemáticos lineares usando algoritmos genéticos

Nesta dissertação apresenta-se o problema de redução de ordem de modelos dinâmicos lineares, sob o ponto de vista de otimização via Algoritmos Genéticos. Uma função custo, obtida a partir da norma dos coeficientes do numerador da função de transferência do erro entre o modelo original e o reduzido, e minimizada por meio de um algoritmo genético, com consequente calculo dos parâmetros do modelo reduzido. O procedimento e aplicado em alguns exemplos que demonstram a validade da abordagem.

Conjuntos limite e bifurfações de campos de vetores suaves por partes no plano

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Metodologia de discretização espacial para o método radial de interpolação por pontos (RPIM) aplicada para solução numérica das equações de Maxwell

Neste trabalho, foi desenvolvido e implementado um método de discretização espacial baseado na lei de Coulomb para geração de pontos que possam ser usados em métodos meshless para solução das equações de Maxwell. Tal método aplica a lei de Coulomb para gerar o equilíbrio espacial necessário para gerar alta qualidade de discretização espacial para um domínio de análise. Este método é denominado aqui de CLDM (Coulomb Law Discretization Method ) e é aplicado a problemas bidimensionais. Utiliza-se o método RPIM (Radial Point Interpolation Method) com truncagem por UPML (Uniaxial Perfectlly Matched Layers) para solução das equações de Maxwell no domínio do tempo (modo TMz).

Treino sem consequências diferenciais: importância conceitual, metodológica e algumas implicações educacionais

Neste trabalho, retomamos os principais estudos sobre o treino de discriminações condicionais sem conseqüências diferenciais implementados em nosso laboratório com humanos desprovidos de história prévia de participação como sujeitos. Buscamos descrever os procedimentos utilizados e discutimos os resultados em torno de aspectos conceituais e metodológicos a respeito dos quais os estudos podem ser inter-relacionados - a questão do reforço (a temporalidade com que este segue o desempenho, enquanto uma característica envolvida em sua definição) e fontes alternativas de controle (o pareamento Modelo-Sc e Modelo-Si, bem como a reexposição às tentativas de treino). Finalmente, relacionamos o procedimento geral dos estudos com a postura do professor em sala de aula e sugerimos algumas implicações educacionais dele decorrentes, como podendo estar envolvidas na programação e implementação do ensino.

Discriminações condicionais sem conseqüências diferenciais em crianças: efeitos da história de treino precoce de simetria

O presente estudo avaliou a formação de classes de estímulos através de um treino de duas discriminações condicionais (AB e AC) com pareamento consistente modelo-comparação correta, sem conseqüências diferenciais imediatas e fading, aplicando-se testes de simetria e de equivalência após cada bloco de treino. Participaram quatro crianças do pré-escolar que foram expostas ao procedimento de ensino, em duas etapas: na Etapa 1, com estímulos usuais e na Etapa 2 com estímulos não usuais Para cada modelo, três estímulos de comparação foram apresentados simultaneamente. Cada modelo foi emparelhado consistentemente, com os estímulos de comparação, sendo que o estímulo de comparação correto e o modelo apareceram em fading ao longo do treino. Relações simétricas foram demonstradas com dois participantes na Etapa 1, mas não ocorreram relações emergentes na Etapa 2. Dois participantes transferiram o desempenho obtido de uma etapa para outra com o treino discriminativo, como uma espécie de learning set arbitrário. Os resultados indicam que a sequência de treino com estímulos usuais e não usuais e a natureza dos estímulos na Etapa 1 foram variáveis relevantes.

Linearização dos coeficientes de reflexão de ondas qP em meios anisotrópicos

Os coeficientes de reflexão de onda qP em uma interface plana que separa dois meios anisotrópicos dependem dos parâmetros elásticos dos meios envolvidos de maneira não linear. Aproximações lineares para a refletividade de ondas qP pressupondo fraco contraste de impedância e fraca anisotropia levam a uma forma simples para se fazer à análise de AVO/AVD. Neste trabalho a solução das equações de Zoeppritz é rescrita explicitamente em função de matrizes de impedância e polarização. Além disso, utilizando-se esta abordagem, é apresentada uma metodologia geral mais simples para se obter as formas linearizadas. Estas formas linearizadas para os coeficientes de reflexão da onda qP e das ondas convertidas para uma onda qP incidente apresentam resultados muito próximos da formulação exata para ângulos de incidência menores que 30° considerando-se um contraste de impedância moderado e anisotropia em limites geologicamente aceitáveis.

Desempenho e características de carcaças de cordeiros suplementados com diferentes níveis de resíduo de biodiesel

Estudou-se a suplementação de dietas de cordeiros confinados com resíduo de biodiesel, utilizando-se 25 ovinos, SRD, machos, castrados, com média de peso de 20kg. O delineamento experimental foi de blocos ao acaso, com cinco repetições por tratamento e esquema de parcelas subdivididas para os dados de biometria. O experimento teve a duração de 70 dias. As dietas eram compostas de 34% de volumoso e 61% de concentrado, à base de milho, soja, minerais e 5% de óleo de dendê ou resíduo de biodiesel, em percentagens crescentes – zero, 25, 50, 75 e 100%. Os animais foram abatidos ao final do experimento para avaliação da carcaça. Foram observados efeitos significativos (P<0,05) de percentagem de inclusão do resíduo, para as variáveis de desempenho e biometria demonstrados nos efeitos lineares crescentes das equações de regressão, assim como para peso da carcaça quente (PCQ) –, com médias de 14; 15,92; 16,14; 16,42 e 18,02% – e peso da carcaça fria (PCF) –, com médias de 13,12; 14,78; 15,06; 15,70 e 17,25% –, para dietas com, respectivamente, 0, 25, 50, 100% de resíduo de biodiesel de dendê. A utilização de resíduo de biodiesel de dendê na alimentação de cordeiros em crescimento é alternativa para aumentar a densidade energética de suas dietas.

Sólitons latentes, cordas negras, membranas negras e equações de estado em modelos de Kaluza-Klein

Neste trabalho investigamos soluções solitônicas em modelos de Kaluza-Klein com um número arbitrário de espaços internos toroidais, que descrevem o campo gravitacional de um objeto massivo compacto. Cada toro d_i-dimensional possui um fator de escala independente C_i, i = 1, ..., N, que é caracterizado pelo parâmetro ᵞi. Destacamos a solução fisicamente interessante correspondente à massa puntual. Para a solução geral obtemos equações de estado nos espaços externo e interno. Estas equações demonstram que a massa pontual solitônica possui equações de estado tipo poeira em todos os espaços. Obtemos também os parâmetros pósnewtonianos que nos possibilitam encontrar as fórmulas da precessão do periélio, do desvio da luz e do atraso no tempo de ecos de radar. Além disso, os experimentos gravitacionais levam a uma forte limitação nos parâmetros do modelo: T = ƩNi=1 diYi = −(2, 1±2, 3)×10⁻⁵. A solução para massa pontual com Y1 = . . . = YN = (1+ƩNi=1 di)⁻¹ contradiz esta restrição. A imposição T = 0 satisfaz essa limitação experimental e define uma nova classe de soluções que são indistinguíveis para a relatividade geral. Chamamos estas soluções de sólitons latentes. Cordas negras e membranas negras com Yi = 0 pertencem a esta classe. Além disso, a condição de estabilidade dos espaços internos destaca cordas/membranas negras de sólitons latentes, conduzindo exclusivamente para as equações de estado de corda/membrana negra p_i = −ε/2, i = 1, . . . ,N, nos espaços internos e ao número de dimensões externas d₀ = 3. As investigações do fluido perfeito multidimensional estático e esfericamente simétrico com equação de estado tipo poeira no espaço externo confirmam os resultados acima.

Estimativa de parâmetros elásticos em meios anisotrópicos

As medidas de amplitude, polarização e vagarosidade contem informações sobre o meio onde a propagação de onda ocorre. Esta tese investiga esses dados com objetivo de estimar as propriedades elásticas deste meio. Coeficientes de reflexão podem ser estimados das amplitudes dos dados e dependem de forma não linear dos contrastes dos parâmetros elásticos e do contraste de densidade entre os meios separados por uma interface. Quando o contraste de impedância é fraco, as aproximações lineares para a refletividade qP são mais convenientes para inversão da densidade e dos parâmetros elásticos usando as análises de amplitude versus ângulo de incidência (AVO) e amplitude versus a direção do plano de incidência (AVD). Escrevendo as equações de Zoepprittz de forma separada nos permite escrever uma solução destas equações em termos das matrizes de impedância e polarização. Usando esta solução são determinadas aproximações lineares para a refletividade da onda qP considerando fraco contraste de impedância, fraca anisotropia mas com classe de simetria de arbitrária. As linearizações são avaliadas para diferentes geometrias de aquisição e várias escolhas do meio de referência. Estas aproximações apresentam bom desempenho comparado com o valor exato do coeficiente de reflexão da onda qP e de suas ondas convertidas para incidências de até 30° e meios que obedecem à hipótese de fraca anisotropia. Um conjunto de fraturas orientado é representado efetivamente por um meio transversalmente isotrópico (TI), as aproximações lineares da refletividade da onda qP podem ser usadas para estimar a orientação de fratura. Partindo deste pressuposto este problema consiste em estimar a orientação do eixo de simetria a partir de dados de refletividade de onda qP. Este trabalho mostra que são necessários múltiplos azimutes e múltiplas incidências para se obter uma estimativa estável. Também é mostrado que apenas os coeficientes das ondas qS e qT são sensíveis ao mergulho da fratura. Foi investigada a estimativa da anisotropia local através de dados de VSP multiazimutal dos vetores de polarização e vagarosidade. Foram usadas medidas da componente vertical do vetor de vagarosidade e o vetor de polarização de ondas qP diretas e refletidas. O esquema de inversão é validado através de exemplos sintéticos considerando diferentes escolhas do vetor normal à frente de onda no meio de referência, meios de referências e geometria de aquisição. Esta análise mostra que somente um subgrupo dos parâmetros elástico pode ser estimado. Uma importante aplicação desta metodologia é o seu potencial para a determinação de classes de anisotropia. A aplicação desta metodologia aos dados do mar de Java mostra que os modelos isotrópicos e TIV são inadequados para o ajuste desses dados.

Processamento interpretativo de dados magnetométricos e inversão de dados gravimétricos aplicados à prospecção de hidrocarbonetos

Apresentamos dois métodos de interpretação de dados de campos potenciais, aplicados à prospecção de hidrocarbonetos. O primeiro emprega dados aeromagnéticos para estimar o limite, no plano horizontal, entre a crosta continental e a crosta oceânica. Este método baseia-se na existência de feições geológicas magnéticas exclusivas da crosta continental, de modo que as estimativas das extremidades destas feições são usadas como estimativas dos limites da crosta continental. Para tanto, o sinal da anomalia aeromagnética na região da plataforma, do talude e da elevação continental é amplificado através do operador de continuação analítica para baixo usando duas implementações: o princípio da camada equivalente e a condição de fronteira de Dirichlet. A maior carga computacional no cálculo do campo continuado para baixo reside na resolução de um sistema de equações lineares de grande porte. Este esforço computacional é minimizado através do processamento por janelas e do emprego do método do gradiente conjugado na resolução do sistema de equações. Como a operação de continuação para baixo é instável, estabilizamos a solução através do funcional estabilizador de primeira ordem de Tikhonov. Testes em dados aeromagnéticos sintéticos contaminados com ruído pseudo-aleatório Gaussiano mostraram a eficiência de ambas as implementações para realçar os finais das feições magnéticas exclusivas da crosta continental, permitindo o delineamento do limite desta com a crosta oceânica. Aplicamos a metodologia em suas duas implementações a dados aeromagnéticos reais de duas regiões da costa brasileira: Foz do Amazonas e Bacia do Jequitinhonha. O segundo método delineia, simultaneamente, a topografia do embasamento de uma bacia sedimentar e a geometria de estruturas salinas contidas no pacote sedimentar. Os modelos interpretativos consistem de um conjunto de prismas bidimensionais verticais justapostos, para o pacote sedimentar e de prismas bidimensionais com seções verticais poligonais para as estruturas salinas. Estabilizamos a solução, incorporando características geométricas do relevo do embasamento e das estruturas salinas compatíveis com o ambiente geológico através dos estabilizadores da suavidade global, suavidade ponderada e da concentração de massa ao longo de direções preferenciais, além de vínculos de desigualdade nos parâmetros. Aplicamos o método a dados gravimétricos sintéticos produzidos por fontes 2D simulando bacias sedimentares intracratônicas e marginais apresentando densidade do pacote sedimentar variando com a profundidade segundo uma lei hiperbólica e abrigando domos e almofadas salinas. Os resultados mostraram que o método apresenta potencial para delinear, simultaneamente, as geometrias tanto de almofadas e domos salinos, como de relevos descontínuos do embasamento. Aplicamos o método, também, a dados reais ao longo de dois perfis gravimétricos sobre as Bacias de Campos e do Jequitinhonha e obtivemos interpretações compatíveis com a geologia da área.

Análise das aproximações RPP e RSP para meios isotrópicos

Neste trabalho são apresentadas aproximações lineares e quadráticas das equações de Zoeppritz para a obtenção dos coeficientes de reflexão e transmissão de eventos P-P e P-S em função dos ângulos de incidência e da média angular, bem como a análise de inversão linear de AVO, considerando os eventos de reflexão P-P e P-S dissociados e combinados. O uso das chamadas aproximações pseudo-quadráticas foi aplicado para obtenção de aproximações quadráticas apenas para eventos-PP, em torno dos contrastes médios das velocidades de ondas compressionais e cisalhantes e da razão Vs/Vp. Os resultados das aproximações desenvolvidas neste trabalho mostram que as aproximações quadráticas são mais precisas que as lineares, nas duas versões angulares. As comparações entre as aproximações em termos do ângulo de incidência e da média angular mostram que as aproximações quadráticas são equivalentes dentro do limite angular de [0º, 30º]. Por outro lado, as aproximações lineares em função do ângulo de incidência mostram-se mais precisas que as aproximações lineares em função da média angular. Na inversão linear, fez-se análises de sensibilidade e de ambigüidade e observou-se que, nos caso de eventos de reflexão P-P e P-S dissociados, apenas um parâmetro pode ser estimado e que a combinação destes eventos consegue estabilizar a inversão, permitindo a estimativa de dois dentre os parâmetros físicos dos meios (contrastes de impedância, de velocidade de onda P e de módulo de cisalhamento).

Estudo do comportamento das feições lineares observado em dados aeromagnéticos da Bacia do Paraná

O processamento dos dados aeromagnéticos em imagens revelou-se como um bom método de interpretação, que aqui denominamos de método das imagens sombreadas. A aplicação deste método foi feita em dados aeromagnéticos de uma área da Bacia do Paraná, com o objetivo de realçar os lineamentos que não são evidentes nos mapas de contorno. Para subsidiar os nossos estudos, utilizamos imagens de LANDSAT porque muitos lineamentos presentes na bacia não possuem resposta ao levantamento aeromagnético. Com a integração destas duas imagens definimos duas direções dominantes para os lineamentos nos seguintes intervalos: N40° — 60°W e N40° — 65°E. Com a definição deste padrão, obtivemos informações do arcabouço estrutural e, consequentemente, uma idéia da evolução tectônica da bacia. A outra parte desenvolvida neste trabalho foi o cálculo do vetor de magnetização total nos principais lineamentos da bacia, e eles apresentaram um total de 5 inclinações de magnetização diferentes. Essas inclinações foram explicadas pela superposição do efeito da magnetização induzida, efeito de reversão do campo geomagnético, e o efeito da magnetização anisotrópica causada pela desmagnetização nos diques.

Estudo das feições lineares observadas em dados aeromagnéticos da Bacia do Solimões a partir de imagens digitais de sombreamento do relevo anômalo total

A Bacia do Solimões é caracterizada por um padrão magnético complexo, expresso ao longo de toda sua extensão, e produzido pela superposição de anomalias magnéticas que se apresentam sob a forma de feições lineares. Tais feições refletem as diferentes atividades tectônicas que atuaram na Região Amazônica através do Pré-Cambriano e Fanerozóico. Neste trabalho, foi aplicado nos dados aeromagnéticos da Bacia do Solimões um método de processamento de imagens digitais de sombreamento do relevo magnético anômalo total que, dadas suas características metodológicas, permitindo utilizá-lo como um filtro direcional, possibilitou a definição de aspectos relevantes no contexto das relações entre os lineamentos magnéticos E/W, NE/SW, e NW/SE. Nesse sentido, foram identificados padrões de lineamentos magnéticos que refletem a existência de zonas de cisalhamento transcorrentes dextrais, orientadas preferencialmente na direção E/W. A interação entre os diversos segmentos transcorrentes promoveu o desenvolvimento de um regime predominantemente transpressivo, representado por falhas reversas associadas aos lineamentos magnéticos E/W e NE/SW; e duplexes direcionais formando falhas em flor positivas e negativas associadas, respectivamente, aos lineamentos magnéticos orientados na direção N70-80E e N70-80W. A análise quantitativa permitiu explicar dois aspectos importantes em relação às feições lineares observadas nas imagens digitais. O primeiro mostra, através de modelamentos baseados na superposição de prismas bidimensionais, que estas feições lineares magnéticas podem ser explicadas pela superposição de fontes profundas intraembasamento altamente magnéticas, e fontes rasas de alta frequência, sendo estas associadas a falhas reversas ao longo dos níveis de diabásio, presentes em forma de intrusões nos sedimentos paleozóicos da Bacia do Solimões. O segundo aspecto, baseado na utilização dos conceitos da cross-covariância, constata a presença de 'offsets' dextrais, associados aos lineamentos magnéticos NE—SW, ao longo da direção E—W. Este fato mostra, quantitativamente, que o padrão magnético desta região pode ser explicado pela presença de zonas de cisalhamento transcorrestes dextrais, cujos processos tectônicos associados foram fortemente condicionados por zonas de fraquezas pré-existentes (Pré-Cambrianas, Paleozóicas) durante o Mesozóico e Cenozóico.

Uso das rotações de givens modificadas como um método direto para obtenção e atualização das soluções em sistemas com acumulação seqüencial de dados

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Inversão gravimétrica rápida do relevo do embasamento aplicando o funcional da variação total

Apresentamos um novo método de inversão linear bidimensional de dados gravimétricos produzidos por bacias sedimentares com relevo do embasamento descontínuo. O método desenvolvido utiliza um modelo interpretativo formado por um conjunto de fitas horizontais bidimensionais justapostas cujas espessuras são os parâmetros a serem estimados. O contraste de densidade entre o embasamento e os sedimentos é presumido constante e conhecido. As estimativas das espessuras foram estabilizadas com o funcional da Variação Total (VT) que permite soluções apresentando descontinuidades locais no relevo do embasamento. As estimativas do relevo são obtidas através da resolução de um sistema de equações lineares, resolvido na norma L₁. Como métodos lineares subestimam as estimativas de profundidade do embasamento de bacias maiores que cerca de 500 m, amplificamos as estimativas de profundidade através da modificação da matriz associada ao modelo interpretativo de fitas. As estimativas obtidas através deste procedimento são em geral ligeiramente superestimadas. Desse modo, elas são corrigidas através de uma correção definida pela expressão da placa Bouguer. Testes em dados sintéticos e reais produziram resultados comparáveis aos produzidos pelo método não linear, mas exigiram menor tempo computacional. A razão R entre os tempos exigidos pelo método não linear e o método proposto cresce com o número de observações e parâmetros. Por exemplo, para 60 observações e 60 parâmetros, R é igual a 4, enquanto para 2500 observações e 2500 parâmetros R cresce para 16,8. O método proposto e o método de inversão não linear foram aplicados também em dados reais do Steptoe Valley, Nevada, Estados Unidos, e da ponte do POEMA, no Campus do Guamá em Belém, produzindo soluções similares às obtidas com o método não linear exigindo menor tempo computacional.