Interpretação aeromagnética automática com uso da equação homogênea de Euler e sua aplicação na Bacia do Solimões

ABSTRACT: We present here a methodology for the rapid interpretation of aeromagnetic data in three dimensions. An estimation of the x, y and z coordinates of prismatic elements is obtained through the application of "Euler's Homogeneous equation" to the data. In this application, it is necessary to have only the total magnetic field and its derivatives. These components can be measured or calculated from the total field data. In the use of Euler's Homogeneous equation, the structural index, the coordinates of the corners of the prism and the depth to the top of the prism are unknown vectors. Inversion of the data by classical least-squares methods renders the problem ill-conditioned. However, the inverse problem can be stabilized by the introduction of both a priori information within the parameter vector together with a weighting matrix. The algorithm was tested with synthetic and real data in a low magnetic latitude region and the results were satisfactory. The applicability of the theorem and its ambiguity caused by the lack of information about the direction of total magnetization, inherent in all automatic methods, is also discussed. As an application, an area within the Solimões basin was chosen to test the method. Since 1977, the Solimões basin has become a center of exploration activity, motivated by the first discovery of gas bearing sandstones within the Monte Alegre formation. Since then, seismic investigations and drilling have been carried on in the region. A knowledge of basement structures is of great importance in the location of oil traps and understanding the tectonic history of this region. Through the application of this method a preliminary estimate of the areal distribution and depth of interbasement and sedimentary magnetic sources was obtained.

Inversão de dados sísmicos de reflexão a partir da curva do tempo de trânsito

A presente Dissertação de Mestrado tem como objetivo o estudo do problema de inversão sísmica baseada em refletores planos para arranjo fonte-comum (FC) e ponto-médiocomum (PMC). O modelo direto é descrito por camadas homogêneas, isotrópicas com interfaces plano-horizontais. O problema é relacionado ao empilhamento NMO baseado na otimização da função semblance, para seções PMC corrigidas de sobretempo normal (NMO). O estudo foi baseado em dois princípios. O primeiro princípio adotado foi de combinar dois grupos de métodos de inversão: um Método Global e um Método Local. O segundo princípio adotado foi o de cascata, segundo a teoria Wichert-Herglotz-Bateman, que estabelece que para conhecer uma camada inferior tem-se que conhecer primeiro a camada superior (dissecação). A aplicação do estudo é voltada à simulação sísmica de Bacia Sedimentar do Solimões e de Bacia Marinha para se obter uma distribuição local 1D de velocidades e espessuras para a subsuperfície em horizontes alvo. Sendo assim, limitamos a inversão entre 4 e 11 refletores, uma vez que na prática a indústria limita uma interpretação realizada apenas em número equivalente de 3 a 4 refletores principais. Ressalta-se que este modelo é aplicável como condição inicial ao imageamento de seções sísmicas em regiões geologicamente complexas com variação horizontal suave de velocidades. Os dados sintéticos foram gerados a partir dos modelos relacionados a informações geológicas, o que corresponde a uma forte informação a priori no modelo de inversão. Para a construção dos modelos relacionados aos projetos da Rede Risco Exploratório (FINEP) e de formação de recursos humanos da ANP em andamento, analisamos os seguintes assuntos relevantes: (1) Geologia de bacias sedimentares terrestre dos Solimões e ma rinha (estratigráfica, estrutural, tectônica e petrolífera); (2) Física da resolução vertical e horizontal; e (3) Discretização temporal-espacial no cubo de multi-cobertura. O processo de inversão é dependente do efeito da discretização tempo-espacial do campo de ondas, dos parâmetros físicos do levantamento sísmico, e da posterior reamostragem no cubo de cobertura múltipla. O modelo direto empregado corresponde ao caso do operador do empilhamento NMO (1D), considerando uma topografia de observação plana. O critério básico tomado como referência para a inversão e o ajuste de curvas é a norma 2 (quadrática). A inversão usando o presente modelo simples é computacionalmente atrativa por ser rápida, e conveniente por permitir que vários outros recursos possam ser incluídos com interpretação física lógica; por exemplo, a Zona de Fresnel Projetada (ZFP), cálculo direto da divergência esférica, inversão Dix, inversão linear por reparametrização, informações a priori, regularização. A ZFP mostra ser um conceito út il para estabelecer a abertura da janela espacial da inversão na seção tempo-distância, e representa a influência dos dados na resolução horizontal. A estimativa da ZFP indica uma abertura mínima com base num modelo adotado, e atualizável. A divergência esférica é uma função suave, e tem base física para ser usada na definição da matriz ponderação dos dados em métodos de inversão tomográfica. A necessidade de robustez na inversão pode ser analisada em seções sísmicas (FC, PMC) submetida a filtragens (freqüências de cantos: 5;15;75;85; banda-passante trapezoidal), onde se pode identificar, comparar e interpretar as informações contidas. A partir das seções, concluímos que os dados são contaminados com pontos isolados, o que propõe métodos na classe dos considerados robustos, tendo-se como referência a norma 2 (quadrados- mínimos) de ajuste de curvas. Os algoritmos foram desenvolvidos na linguagem de programação FORTRAN 90/95, usando o programa MATLAB para apresentação de resultados, e o sistema CWP/SU para modelagem sísmica sintética, marcação de eventos e apresentação de resultados.

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.

Estudos isotópicos (Pb, O, H, S) em zonas alteradas e mineralizadas do depósito cupro-aurífero Visconde, Província Mineral de Carajás

O depósito cupro-aurífero Visconde está localizado na Província Mineral de Carajás, a cerca de 15 km a leste do depósito congênere de classe mundial Sossego. Encontra-se em uma zona de cisalhamento de direção WNW-ESE, que marca o contato das rochas metavulcanossedimentares da Bacia Carajás com o embasamento. Nessa zona ocorrem outros depósitos hidrotermais cupro-auríferos com características similares (Alvo 118, Cristalino, Jatobá, Bacaba, Bacuri, Castanha), que têm sido enquadrados na classe IOCG (Iron Oxide Copper-Gold), embora muitas dúvidas ainda existam quanto a sua gênese, principalmente no que diz respeito à idade da mineralização e fontes dos fluidos, ligantes e metais. O depósito Visconde está hospedado em rochas arqueanas variavelmente cisalhadas e alteradas hidrotermalmente, as principais sendo metavulcânicas félsicas (2968 ± 15 Ma), o Granito Serra Dourada (2860 ± 22 Ma) e gabros/dioritos. Elas registram diversos tipos de alteração hidrotermal com forte controle estrutural, destacando-se as alterações sódica (albita + escapolita) e sódico-cálcica (albita + actinolita ± turmalina ± quartzo ± magnetita ± escapolita), mais precoces, que promoveram a substituição ubíqua de minerais primários das rochas e a disseminação de calcopirita, pirita, molibdenita e pentlandita. Dados isotópicos de oxigênio e hidrogênio de minerais representativos desses tipos de alteração mostram que os fluidos hidrotermais foram quentes (410 – 355°C) e ricos em ¹⁸O (δ¹⁸O_H2O= +4,2 a 9,4‰). Sobreveio a alteração potássica, caracterizada pela intensa biotitização das rochas, a qual ocorreu concomitantemente ao desenvolvimento de foliação milonítica, notavelmente desenhada pela orientação de palhetas de biotita, que precipitaram de fluidos com assinatura isotópica de oxigênio similar à dos estágios anteriores (δ¹⁸O_H2O entre +4,8 e +7,2‰, a 355°C). Microclina e alanita são outras fases características desse estágio, além da calcopirita precipitada nos planos da foliação. A temperaturas mais baixas (230 ± 11°C), fluidos empobrecidos em ¹⁸O (δ¹⁸O_H2O = -1,3 a +3,7‰) geraram associações de minerais cálcico-magnesianos (albita + epidoto + clorita ± calcita ± actinolita) que são contemporâneas à mineralização. Valores de δ¹⁸D_H2O e δO_H2O indicam que os fluidos hidrotermais foram inicialmente formados por águas metamórficas e formacionais, a que se misturou alguma água de fonte magmática. Nos estágios tardios, houve considerável influxo de águas superficiais. Diluição e queda da temperatura provocaram a precipitação de abundantes sulfetos (calcopirita ± bornita ± calcocita ± digenita), os quais se concentraram principalmente em brechas tectônicas - os principais corpos de minério - que chegam a conter até cerca de 60% de sulfetos. Veios constituídos por minerais sódico-cálcicos também apresentam comumente sulfetos. A associação de minerais de minério e ganga indica uma assinatura de Cu-Au- Fe-Ni-ETRL-B-P para a mineralização. Os valores de δ³⁴S (-1,2 a +3,4‰) de sulfetos sugerem enxofre de origem magmática (proveniente da exsolução de magmas ou da dissolução de sulfetos das rochas ígneas pré-existentes) e precipitação em condições levemente oxidantes. Datação do minério por lixiviação e dissolução total de Pb em calcopirita forneceu idades de 2736 ± 100 Ma e 2729 ± 150 Ma, que indicam ser a mineralização neoarqueana e, a despeito dos altos erros, permite descartar um evento mineralizador paleoproterozoico. A idade de 2746 ± 7 Ma (MSDW=4,9; evaporação de Pb em zircão), obtida em um corpo granítico não mineralizado (correlacionado à Suíte Planalto) que ocorre na área do depósito, foi interpretada como a idade mínima da mineralização. Assim, a formação do depósito Visconde teria relação com o evento transpressivo ocorrido entre 2,76 e 2,74 Ga, reponsável pela inversão da Bacia Carajás e pela geração de magmatismo granítico nos domínios Carajás e de Transição. Esse evento teria desencadeado reações de devolatilização em rochas do Supergrupo Itacaiúnas, ou mesmo, provocado a expulsão de fluidos conatos salinos aprisionados em seus intertícios. Esses fluidos teriam migrado pelas zonas de cisalhamento e reagido com as rochas (da bacia e do embasamento) pelas quais se movimentaram durante a fase dúctil. As concentrações subeconômicas do depósito Visconde devem ser resultado da ausência de grandes estruturas que teriam favorecido maior influxo de fluidos superficiais, tal como ocorreu na formação dos depósitos Sossego e Alvo 118.

Inversion method for obtaining electrical parameters for soil in the rural region in Brazil

This paper presents an inversion methodology through weighted least squares to obtain the electrical parameters for the soil of a typical mid-western region in Brazil using the model based in the formalism of the parabolic equations to calculate the electric field intensity received. To validate this methodology, the results of the radio signal measurement campaign conducted in six radial routes leaving the city of Brasilia, Federal District, where the transmitter was located, were used. The measurements were compared to computer simulations and, thus, the optimal values for the electric conductivity and relative permittivity for the soil of the region could be estimated. Finally, a quantitative analysis of these parameters was performed with the values found in the literature, which demonstrated the effectiveness of the proposed methodology.

Aeroelastic stability analysis using linear matrix inequalities

The present work describes an alternative methodology for identification of aeroelastic stability in a range of varying parameters. Analysis is performed in time domain based on Lyapunov stability and solved by convex optimization algorithms. The theory is outlined and simulations are carried out on a benchmark system to illustrate the method. The classical methodology with the analysis of the system's eigenvalues is presented for comparing the results and validating the approach. The aeroelastic model is represented in state space format and the unsteady aerodynamic forces are written in time domain using rational function approximation. The problem is formulated as a polytopic differential inclusion system and the conceptual idea can be used in two different applications. In the first application the method verifies the aeroelastic stability in a range of air density (or its equivalent altitude range). In the second one, the stability is verified for a rage of velocities. These analyses are in contrast to the classical discrete analysis performed at fixed air density/velocity values. It is shown that this method is efficient to identify stability regions in the flight envelope and it offers promise for robust flutter identification.

Measurement of the t-channel single-top-quark production cross section and of the |V (tb)| CKM matrix element in pp collisions at root s=8 TeV

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Hydroxyapatite/beta-tricalcium phosphate and enamel matrix derivative for treatment of proximal class II furcation defects: a randomized clinical trial

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

Mesenchymal stem cells modulate release of matrix proteins from tendon surfaces in vitro: a potential beneficial therapeutic effect

Aim: Injury of tendons contained within a synovial environment, such as joint, bursa or tendon sheath, frequently fails to heal and releases matrix proteins into the synovial fluid, driving inflammation. This study investigated the effectiveness of cells to seal tendon surfaces and provoke matrix synthesis as a possible effective injectable therapy. Materials & methods: Equine flexor tendon explants were cultured overnight in suspensions of bone marrow and synovium-derived mesenchymal stems cells and, as controls, two sources of fibroblasts, derived from tendon and skin, which adhered to the explants. Release of the most abundant tendon extracellular matrix proteins into the media was assayed, along with specific matrix proteins synthesis by real-time PCR. Results: Release of extracellular matrix proteins was influenced by the coating cell type. Fibroblasts from skin and tendon appeared less capable of preventing the release of matrix proteins than mesenchymal stems cells. Conclusion: The source of cell is an important consideration for cell therapy.

Polymer Matrix-Based Piezoelectric Composite for Structural Health Monitoring

Structural health monitoring (SHM) refers to the procedure of assessing the structure conditions continuously so it is an alternative to conventional nondestructive evaluation (NDE) techniques [1]. With the growing developments in sensor technology acoustic emission (AE) technology has been attracting attention in SHM applications. AE are characterized by waves produced by the sudden internal stress redistribution caused by the changes in the internal structure, such as fatigue, crack growth, corrosion, etc. Piezoelectric materials such as Lead Zirconate Titanate (PZT) ceramic have been widely used as sensor due to its high electromechanical coupling factor and piezoelectric d coefficients. Because of the poor mechanical characteristic and the lack in the formability of the ceramic, polymer matrix-based piezoelectric composites have been studied in the last decade in order to obtain better properties in comparison with a single phase material. In this study a composite film made of polyurethane (PU) and PZT ceramic particles partially recovered with polyaniline (PAni) was characterized and used as sensor for AE detection. Preliminary results indicate that the presence of a semiconductor polymer (PAni) recovering the ceramic particles, make the poling process easier and less time consuming. Also, it is possible to observe that there is a great potential to use such type of composite as sensor for structure health monitoring.

Urethral striated muscle and extracellular matrix morphological characteristics among mildly diabetic pregnant rats: translational approach

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

Effects of antiproteolytic agents on corneal epithelial viability and matrix metalloproteinase-2 and metalloproteinase-9 activity in alkali-burned corneas of rats

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Inactivation of Matrix-bound Matrix Metalloproteinases by Cross-linking Agents in Acid-etched Dentin

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stabilization of dentin matrix after cross-linking treatments, in vitro

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)