Effects of gender on stroke rates, critical speed and velocity of a 30-min swim in young swimmers

Our objective was to analyze the effect of gender on the relationship between stroke rates corresponding to critical speed (SRCS) and maximal speed of 30 min (SRS30) in young swimmers. Twenty two males (GM1) (Age = 15.4 ± 2.1 yr., Body mass = 63.7 ± 12.9 kg, Stature = 1.73 ± 0.09 m) and fourteen female (GF) swimmers (Age = 15.1 ± 1.6 yr., Body mass = 58.3 ± 8.8 kg, Stature = 1.65 ± 0.06 m) were studied. A subset of males (GM2) was matched to the GF by their velocity for a 30 min swim (S30). The critical speed (CS) was determined through the slope of the linear regression line between the distances (200 and 400 m) and participant's respective times. CS was significantly higher than S30 in males (GM1 - 1.25 and 1.16 and GM2 - 1.21 and 1.12 m·s-1) and females (GF - 1.15 and 1.11 m·s-1). There was no significant difference between SRCS and SRS30 in males (GM1 - 34.16 and 32.32 and GM2 - 34.67 and 32.46 cycle·s-1, respectively) and females (GF - 34.18 and 33.67 cycle·s-1-1, respectively). There was a significant correlation between CS and S30 (GM1 - r = 0.89, GF - r = 0.94 and GM2 - r = 0.90) and between SRCS and SRS30 (GM1 - r = 0.89, GF - r = 0.80 and GM2 - r = 0.88). Thus, the relationship between SRCS and SRS30 is not influenced by gender, in swimmers with similar and different aerobic capacity levels. ©Journal of Sports Science and Medicine (2007).

Determination of critical velocity and anaerobic capacity of running rats

The aim of this study was to validate a non-invasive protocol to determine aerobic and anaerobic capacity of treadmill running rats. Thirteen male Wistar rats (90 days old) were submitted to 4 exercise tests, consisting of running at 25, 30, 35 and 40 m min-1, continuously until exhaustion. For the critical velocity (CV) and anaerobic running capacity (ARC) estimations, the hyperbolic curve (velocity versus time to exhaustion (tlim)) was linearized to V= CV+ARC/tlim, where the CV and ARC were linear and slope coefficients, respectively. In order to verify if the CV was the maximal aerobic intensity, the rats were submitted to the maximal lactate steady state test (MLSS) composed of three 25-minute tests of continuous running trials at 15, 20 and 25 m min-1, with blood collection every 5 minutes. The CV was obtained at 22.8±0.7 m min-1 and the ARC, at 26.80±2.77 m. The MLSS was observed at 20m min-1, with blood lactate 3.84 ± 0.31 mmol L-1. There was a progressive increase in lactate concentration at 25 m min-1. The CV and MLSS were different, but presented a high and significant correlation (r=0.81). These results indicate that the non-invasive protocol can be used for physical evaluation of aerobic running rats, but the ARC should still be further investigated.

Effects of 12-week overground walking training at ventilatory threshold velocity in type 2 diabetic women

This study analyzed the effects of overground walking training at ventilatory threshold (VT) velocity on glycaemic control, body composition, physical fitness and lipid profile in DM2 women. Nineteen sedentary patients were randomly assigned to a control group (CG; n=10, 55.9±2.2 years) or a trained group (TG; n=9, 53.4±2.3 years). Both groups were subjected to anthropometric measures, a 12-h fasting blood sampling and a graded treadmill exercise test at baseline and after a 12-week period, during which TG followed a training program involving overground walking at VT velocity for 20-60min/session three times/week. Significant group×time interactions (P<0.05) in glycated hemoglobin (HbA1c), body mass, body mass index (BMI), peak oxygen uptake (VO 2peak) and exercise duration were observed as effects of training exercise, whereas intervention did not induced significant changes (P>0.05) in fasting blood glucose, submaximal fitness parameters and lipid profile. Our results suggest that overground walking training at VT velocity improves long term glycaemic control, body composition and exercise capacity, attesting for the relevance of this parameter as an effective strategy for the exercise intensity prescription in DM2 population. © 2011 Elsevier B.V.

Determination of minimum fluidization velocity for gas-solid beds by experimental data and numerical simulations

The purpose of this work is to predict the minimum fluidization velocity Umf in a gas-solid fluidized bed. The study was carried out with an experimental apparatus for sand particles with diameters between 310μm and 590μm, and density of 2,590kg/m3. The experimental results were compared with numerical simulations developed in MFIX (Multiphase Flow with Interphase eXchange) open source code [1], for three different sizes of particles: 310mum, 450μm and 590μm. A homogeneous mixture with the three kinds of particles was also studied. The influence of the particle diameter was presented and discussed. The Ergun equation was also used to describe the minimum fluidization velocity. The experimental data presented a good agreement with Ergun equation and numerical simulations. Copyright © 2011 by ASME.

Anaerobic running capacity determined from the critical velocity model is not significantly associated with maximal accumulated oxygen deficit in army runners

Purpose: The aim of this study was to verify whether there is an association between anaerobic running capacity (ARC) values, estimated from two-parameter models, and maximal accumulated oxygen deficit (MAOD) in army runners. Methods: Eleven, trained, middle distance runners who are members of the armed forces were recruited for the study (20 ± 1 years). They performed a critical velocity test (CV) for ARC estimation using three mathematical models and an MAOD test, both tests were applied on a motorized treadmill. Results: The MAOD was 61.6 ± 5.2 mL/kg (4.1 ± 0.3 L). The ARC values were 240.4 ± 18.6 m from the linear velocity-inverse time model, 254.0 ± 13.0 m from the linear distance-time model, and 275.2 ± 9.1 m from the hyperbolic time-velocity relationship (nonlinear 2-parameter model), whereas critical velocity values were 3.91 ± 0.07 m/s, 3.86 ± 0.08 m/s and 3.80 ± 0.09 m/s, respectively. There were differences (P < 0.05) for both the ARC and the CV values when compared between velocity-inverse time linear and nonlinear 2-parameter mathematical models. The different values of ARC did not significantly correlate with MAOD. Conclusion: In conclusion, estimated ARC did not correlate with MAOD, and should not be considered as an anaerobic measure of capacity for treadmill running. © 2013 Elsevier Masson SAS. All rights reserved.

Are early and late rate of force development differently influenced by fast-velocity resistance training?

This study examined the effect of fast-velocity concentric isokinetic resistance training (FV) on the rate of force development (RFD) at early (<100 ms) and late phases (>100 ms) of rising muscle force. Nine men participated in a 6-week resistance training intervention for the lower body, and nine matched subjects participated as controls (CON). During concentric isokinetic (180°s-1) knee extension training, subjects were instructed to do each contraction 'as fast and forcefully as possible'. Maximal muscle strength (MVC) and RFD (0-10, 0-20, ..., 0-250 ms from the onset of contraction) were measured during maximal voluntary isometric contraction of the knee extensors (KE). There were no significant changes in MVC of KE in both groups after intervention (FV = 314·2 ± 101·1 versus 338·7 ± 88·0 N{bullet operator}m, P>0·05; CON = 293·3 ± 94·8 versus 280·0 ± 72·2 N{bullet operator}m, P>0·05). The RFD increased 39-71% at time intervals up to 90 ms from the onset of the contraction (P<0·05), whereas no change occurred at later time intervals. Similarly, relative RFD (i.e.%MVC{bullet operator}s-1) (RFDr) increased 33-56% at time intervals up to 70 ms from the onset of the contraction (P<0·05). It can be concluded that a short period of resistance training performed with concentric fast-velocity isokinetic muscle contractions is able to enhance RFD and RFDr obtained at the early phase of rising muscle force. © 2013 The Authors Clinical Physiology and Functional Imaging © 2013 Scandinavian Society of Clinical Physiology and Nuclear Medicine.

Strain rate of the South American lithospheric plate by SIRGAS-CON geodetic observations

A multiyear solution of the SIRGAS-CON network was used to estimate the strain rates of the earth surface from the changing directions of the velocity vectors of 140 geodetic points located in the South American plate. The strain rate was determined by the finite element method using Delaunay triangulation points that formed sub-networks; each sub-network was considered a solid and homogeneous body. The results showed that strain rates vary along the South American plate and are more significant on the western portion of the plate, as expected, since this region is close to the subduction zone of the Nazca plate beneath the South American plate. After using Euler vectors to infer Nazca plate movement and to orient the velocity vectors of the South American plate, it was possible to estimate the convergence and accommodation rates of the Nazca and South American plates, respectively. Strain rate estimates permitted determination of predominant contraction and/or extension regions and to establish that contraction regions coincide with locations with most of the high magnitude seismic events. Some areas with extension and contraction strains were found to the east within the stable South American plate, which may result from different stresses associated with different geological characteristics. These results suggest that major movements detected on the surface near the Nazca plate occur in regions with more heterogeneous geological structures and multiple rupture events. Most seismic events in the South American plate are concentrated in areas with predominant contraction strain rates oriented northeast-southwest; significant amounts of elastic strain can be accumulated on geological structures away from the plate boundary faults; and, behavior of contractions and extensions is similar to what has been found in seismological studies. © 2013 Elsevier Ltd.

Desenvolvimento de sistema para sísmica up-hole em conjunto com o ensaio SPT

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

Migração em profundidade pré-empilhamento utilizando os atributos cinemáticos do empilhamento por superfície de reflexão comum

O empilhamento por superfície de reflexão comum (ou empilhamento SRC), conhecido como empilhamento CRS, do inglês Commom reflection surface, constitui-se em um novo método para o processamento sísmico na simulação de seções afastamento nulo (AN) e afastamento comum (AC). Este método é baseado em uma aproximação paraxial hiperbólica de segunda ordem dos tempos de trânsito de reflexão na vizinhança de um raio central. Para a simulação de seção AN, o raio central é um raio normal, enquanto que para a simulação de uma seção AC o raio central é um raio de afastamento finito. Em adição à seção AN, o método de empilhamento SRC também fornece estimativas dos atributos cinemáticos do campo de onda, sendo aplicados, por exemplo, na determinação (por um processo de inversão) da velocidade intervalar, no cálculo do espalhamento geométrico, na estimativa da zona de Fresnel, e também na simulação de eventos de tempos de difrações, este último tendo uma grande importância para a migração pré-empilhamento. Neste trabalho é proposta uma nova estratégia para fazer uma migração em profundidade pré-empilhamento, que usa os atributos cinemáticos do campo de onda derivados do empilhamento SRC, conhecido por método CRS-PSDM, do inglês CRS based pre-stack depth migration. O método CRS-PSDM usa os resultados obtidos do método SRC, isto é, as seções dos atributos cinemáticos do campo de onda, para construir uma superfície de tempos de trânsito de empilhamento, ao longo da qual as amplitudes do dado sísmico de múltipla cobertura são somadas, sendo o resultado da soma atribuído a um dado ponto em profundidade, na zona alvo de migração que é definida por uma malha regular. Similarmente ao método convencional de migração tipo Kirchhoff (K-PSDM), o método CRS-PSDM precisa de um modelo de velocidade de migração. Contrário ao método K-PSDM, o método CRS-PSDM necessita apenas computar os tempos de trânsito afastamento nulo, ao seja, ao longo de um único raio ligando o ponto considerado em profundidade a uma dada posição de fonte e receptor coincidentes na superfície. O resultado final deste procedimento é uma imagem sísmica em profundidade dos refletores a partir do dado de múltipla cobertura.

Modelagem sísmica via métodos das diferenças finitas: caso da bacia do Amazonas

Este trabalho tem por objetivo apresentar os resultados da modelagem sísmica em meios com fortes descontinuidades de propriedades físicas, com ênfase na existência de difrações e múltiplas reflexões, tendo a Bacia do Amazonas como referência à modelagem. As condições de estabilidade e de fronteiras utilizadas no cálculo do campo de ondas sísmicas foram analisadas numericamente pelo método das diferenças finitas, visando melhor compreensão e controle da interpretação de dados sísmicos. A geologia da Bacia do Amazonas é constituída por rochas sedimentares depositadas desde o Ordoviciano até o Recente que atingem espessuras da ordem de 5 km. Os corpos de diabásio, presentes entre os sedimentos paleozóicos, estão dispostos na forma de soleiras, alcançam espessuras de centenas de metros e perfazem um volume total de aproximadamente 90000 Km³. A ocorrência de tais estruturas é responsável pela existência de reflexões múltiplas durante a propagação da onda sísmica o que impossibilita melhor interpretação dos horizontes refletores que se encontram abaixo destas soleiras. Para representar situações geológicas desse tipo foram usados um modelo (sintético) acústico de velocidades e um código computacional elaborado via método das diferenças finitas com aproximação de quarta ordem no espaço e no tempo da equação da onda. A aplicação dos métodos de diferenças finitas para o estudo de propagação de ondas sísmicas melhorou a compreensão sobre a propagação em meios onde existem heterogeneidades significativas, tendo como resultado boa resolução na interpretação dos eventos de reflexão sísmica em áreas de interesse. Como resultado dos experimentos numéricos realizados em meio de geologia complexa, foi observada a influência significativa das reflexões múltiplas devido à camada de alta velocidade, isto provocou maior perda de energia e dificultou a interpretação dos alvos. Por esta razão recomenda-se a integração de dados de superfície com os de poço, com o objetivo de obter melhor imagem dos alvos abaixo das soleiras de diabásio.

Migração (2,5-D) com amplitudes verdadeiras em meios com gradiente constante de velocidade

A migração com amplitudes verdadeiras de dados de reflexão sísmica, em profundidade ou em tempo, possibilita que seja obtida uma medida dos coeficientes de reflexão dos chamados eventos de reflexão primária. Estes eventos são constituídos, por exemplo, pelas reflexões de ondas longitudinais P-P em refletores de curvaturas arbitrárias e suaves. Um dos métodos mais conhecido é o chamado migração de Kirchhoff, através do qual a imagem sísmica é produzida pela integração do campo de ondas sísmicas, utilizando-se superfícies de difrações, denominadas de Superfícies de Huygens. A fim de se obter uma estimativa dos coeficientes de reflexão durante a migração, isto é a correção do efeito do espalhamento geométrico, utiliza-se uma função peso no operador integral de migração. A obtenção desta função peso é feita pela solução assintótica da integral em pontos estacionários. Tanto no cálculo dos tempos de trânsito como na determinação da função peso, necessita-se do traçamento de raios, o que torna a migração em situações de forte heterogeneidade da propriedade física um processo com alto custo computacional. Neste trabalho é apresentado um algoritmo de migração em profundidade com amplitudes verdadeiras, para o caso em que se tem uma fonte sísmica pontual, sendo o modelo de velocidades em subsuperfície representado por uma função que varia em duas dimensões, e constante na terceira dimensão. Esta situação, conhecida como modelo dois-e-meio dimensional (2,5-D), possui características típicas de muitas situações de interesse na exploração do petróleo, como é o caso da aquisição de dados sísmicos 2-D com receptores ao longo de uma linha sísmica e fonte sísmica 3-D. Em particular, é dada ênfase ao caso em que a velocidade de propagação da onda sísmica varia linearmente com a profundidade. Outro tópico de grande importância abordado nesse trabalho diz respeito ao método de inversão sísmica denominado empilhamento duplo de difrações. Através do quociente de dois empilhamentos com pesos apropriados, pode-se determinar propriedades físicas e parâmetros geométricos relacionados com a trajetória do raio refletido, os quais podem ser utilizados a posteriori no processamento dos dados sísmicos, visando por exemplo, a análise de amplitudes.

Determinação das velocidades intervalares usando a teoria paraxial do raio: aproximação de segunda ordem dos tempos de trânsito

Neste trabalho foi desenvolvido um método de solução ao problema inverso para modelos sísmicos compostos por camadas homogêneas e isotrópicas separadas por superfícies suaves, que determina as velocidades intervalares em profundidade e calcula a geometria das interfaces. O tempo de trânsito é expresso como uma função de parâmetros referidos a um sistema de coordenadas fixo no raio central, que é determinada numericamente na superfície superior do modelo. Essa função é posteriormente calculada na interface anterior que limita a camada não conhecida, através de um processo que determina a função característica em profundidade. A partir da função avaliada na interface anterior se calculam sua velocidade intervalar e a geometria da superfície posterior onde tem lugar a reflexão do raio. O procedimento se repete de uma forma recursiva nas camadas mais profundas obtendo assim a solução completa do modelo, não precisando em nenhum passo informação diferente à das camadas superiores. O método foi expresso num algoritmo e se desenvolveram programas de computador, os quais foram testados com dados sintéticos de modelos que representam feições estruturais comuns nas seções geológicas, fornecendo as velocidades em profundidade e permitindo a reconstrução das interfaces. Uma análise de sensibilidade sobre os programas mostrou que a determinação da função característica e a estimação das velocidades intervalares e geometria das interfaces são feitos por métodos considerados estáveis. O intervalo empírico de aplicabilidade das correções dinâmicas hiperbólicas foi tomado como uma estimativa da ordem de magnitude do intervalo válido para a aplicação do método.

Empilhamento pelo método superfície de reflexão comum 2-D com topografia e introdução ao caso 3-D

O método de empilhamento sísmico CRS simula seções sísmicas ZO a partir de dados de cobertura múltipla, independente do macro-modelo de velocidades. Para meios 2-D, a função tempo de trânsito de empilhamento depende de três parâmetros, a saber: do ângulo de emergência do raio de reflexão normal (em relação à normal da superfície) e das curvaturas das frentes de onda relacionadas às ondas hipotéticas, denominadas NIP e Normal. O empilhamento CRS consiste na soma das amplitudes dos traços sísmicos em dados de múltipla cobertura, ao longo da superfície definida pela função tempo de trânsito do empilhamento CRS, que melhor se ajusta aos dados. O resultado do empilhamento CRS é assinalado a pontos de uma malha pré-definida na seção ZO. Como resultado tem-se a simulação de uma seção sísmica ZO. Isto significa que para cada ponto da seção ZO deve-se estimar o trio de parâmetros ótimos que produz a máxima coerência entre os eventos de reflexão sísmica. Nesta Tese apresenta-se fórmulas para o método CRS 2-D e para a velocidade NMO, que consideram a topografia da superfície de medição. O algoritmo é baseado na estratégia de otimização dos parâmetros de fórmula CRS através de um processo em três etapas: 1) Busca dos parâmetros, o ângulo de emergência e a curvatura da onda NIP, aplicando uma otimização global, 2) busca de um parâmetro, a curvatura da onda N, aplicando uma otimização global, e 3) busca de três parâmetros aplicando uma otimização local para refinar os parâmetros estimados nas etapas anteriores. Na primeira e segunda etapas é usado o algoritmo Simulated Annealing (SA) e na terceira etapa é usado o algoritmo Variable Metric (VM). Para o caso de uma superfície de medição com variações topográficas suaves, foi considerada a curvatura desta superfície no algoritmo do método de empilhamento CRS 2-D, com aplicação a dados sintéticos. O resultado foi uma seção ZO simulada, de alta qualidade ao ser comparada com a seção ZO obtida por modelamento direto, com uma alta razão sinal-ruído, além da estimativa do trio de parâmetros da função tempo de trânsito. Foi realizada uma nálise de sensibilidade para a nova função de tempo de trânsito CRS em relação à curvatura da superfície de medição. Os resultados demonstraram que a função tempo de trânsito CRS é mais sensível nos pontos-médios afastados do ponto central e para grandes afastamentos. As expressões da velocidade NMO apresentadas foram aplicadas para estimar as velocidades e as profundidades dos refletores para um modelo 2-D com topografia suave. Para a inversão destas velocidades e profundidades dos refletores, foi considerado o algoritmo de inversão tipo Dix. A velocidade NMO para uma superfície de medição curva, permite estimar muito melhor estas velocidades e profundidades dos refletores, que as velocidades NMO referidas as superfícies planas. Também apresenta-se uma abordagem do empilhamento CRS no caso 3-D. neste caso a função tempo de trânsito depende de oito parâmetros. São abordadas cinco estratégias de busca destes parâmetros. A combinação de duas destas estratégias (estratégias das três aproximações dos tempos de trânsito e a estratégia das configurações e curvaturas arbitrárias) foi aplicada exitosamente no empilhamento CRS 3-D de dados sintéticos e reais.

Empilhamento sísmico por superfície de reflexão comum: um novo algoritmo usando otimização global e local

O método de empilhamento sísmico por Superfície de Reflexão Comum (ou empilhamento SRC) produz a simulação de seções com afastamento nulo (NA) a partir dos dados de cobertura múltipla. Para meios 2D, o operador de empilhamento SRC depende de três parâmetros que são: o ângulo de emergência do raio central com fonte-receptor nulo (β₀), o raio de curvatura da onda ponto de incidência normal (R_NIP) e o raio de curvatura da onda normal (R_N). O problema crucial para a implementação do método de empilhamento SRC consiste na determinação, a partir dos dados sísmicos, dos três parâmetros ótimos associados a cada ponto de amostragem da seção AN a ser simulada. No presente trabalho foi desenvolvido uma nova sequência de processamento para a simulação de seções AN por meio do método de empilhamento SRC. Neste novo algoritmo, a determinação dos três parâmetros ótimos que definem o operador de empilhamento SRC é realizada em três etapas: na primeira etapa são estimados dois parâmetros (β°₀ e R°_NIP) por meio de uma busca global bidimensional nos dados de cobertura múltipla. Na segunda etapa é usado o valor de β°₀ estimado para determinar-se o terceiro parâmetro (R°_N) através de uma busca global unidimensional na seção AN resultante da primeira etapa. Em ambas etapas as buscas globais são realizadas aplicando o método de otimização Simulated Annealing (SA). Na terceira etapa são determinados os três parâmetros finais (β₀, R_NIP e R_N) através uma busca local tridimensional aplicando o método de otimização Variable Metric (VM) nos dados de cobertura múltipla. Nesta última etapa é usado o trio de parâmetros (β°₀, R°_NIP, R°_N) estimado nas duas etapas anteriores como aproximação inicial. Com o propósito de simular corretamente os eventos com mergulhos conflitantes, este novo algoritmo prevê a determinação de dois trios de parâmetros associados a pontos de amostragem da seção AN onde há intersecção de eventos. Em outras palavras, nos pontos da seção AN onde dois eventos sísmicos se cruzam são determinados dois trios de parâmetros SRC, os quais serão usados conjuntamente na simulação dos eventos com mergulhos conflitantes. Para avaliar a precisão e eficiência do novo algoritmo, este foi aplicado em dados sintéticos de dois modelos: um com interfaces contínuas e outro com uma interface descontinua. As seções AN simuladas têm elevada razão sinal-ruído e mostram uma clara definição dos eventos refletidos e difratados. A comparação das seções AN simuladas com as suas similares obtidas por modelamento direto mostra uma correta simulação de reflexões e difrações. Além disso, a comparação dos valores dos três parâmetros otimizados com os seus correspondentes valores exatos calculados por modelamento direto revela também um alto grau de precisão. Usando a aproximação hiperbólica dos tempos de trânsito, porém sob a condição de R_NIP = R_N, foi desenvolvido um novo algoritmo para a simulação de seções AN contendo predominantemente campos de ondas difratados. De forma similar ao algoritmo de empilhamento SRC, este algoritmo denominado empilhamento por Superfícies de Difração Comum (SDC) também usa os métodos de otimização SA e VM para determinar a dupla de parâmetros ótimos (β₀, R_NIP) que definem o melhor operador de empilhamento SDC. Na primeira etapa utiliza-se o método de otimização SA para determinar os parâmetros iniciais β°₀ e R°_NIP usando o operador de empilhamento com grande abertura. Na segunda etapa, usando os valores estimados de β°₀ e R°_NIP, são melhorados as estimativas do parâmetro R_NIP por meio da aplicação do algoritmo VM na seção AN resultante da primeira etapa. Na terceira etapa são determinados os melhores valores de β°₀ e R°_NIP por meio da aplicação do algoritmo VM nos dados de cobertura múltipla. Vale salientar que a aparente repetição de processos tem como efeito a atenuação progressiva dos eventos refletidos. A aplicação do algoritmo de empilhamento SDC em dados sintéticos contendo campos de ondas refletidos e difratados, produz como resultado principal uma seção AN simulada contendo eventos difratados claramente definidos. Como uma aplicação direta deste resultado na interpretação de dados sísmicos, a migração pós-empilhamento em profundidade da seção AN simulada produz uma seção com a localização correta dos pontos difratores associados às descontinuidades do modelo.

Imageamento homeomórfico de refletores sísmicos

Neste trabalho é apresentada uma nova técnica para a realização do empilhamento sísmico, aplicada ao problema do imageamento de refletores fixos em um meio bidimensional, suavemente heterogêneo, isotrópico, a partir de dados de reflexão. Esta nova técnica chamada de imageamento homeomórfico tem como base a aproximação geométrica do raio e propriedades topológicas dos refletores. São utilizados, portanto, os conceitos de frente de onda, ângulo de incidência, raio de curvatura da frente de onda, cáustica e definição da trajetória do raio; de tal modo que a imagem obtida mantém relações de homeomorfismo com o objeto que se deseja imagear. O empilhamento sísmico é feito, nesta nova técnica de imageamento, aplicando-se uma correção local do tempo, ∆ t, ao tempo de trânsito, t, do raio que parte da fonte sísmica localizada em x_o, reflete-se em um ponto de reflexão, C_o, sendo registrado como uma reflexão primária em um geofone localizado em x_g, em relação ao tempo de referência t_o no sismograma, correspondente ao tempo de trânsito de um raio central. A fórmula utilizada nesta correção temporal tem como parâmetros o raio de curvatura R_o, o ângulo de emergência β_o da frente de onda, no instante em que a mesma atinge a superfície de observação, e a velocidade v_o considerada constante nas proximidades da linha sísmica. Considerando-se uma aproximação geométrica seguido um círculo para a frente de onda, pode-se estabelecer diferentes métodos de imageamento homeomórfico dependendo da configuração de processamento. Sendo assim tem-se: 1) Método Elemento de Fonte (Receptor) Comum (EF(R)C). Utiliza-se uma configuração onde se tem um conjunto de sismogramas relacionado com uma única fonte (receptor), e considera-se uma frente de onda real (de reflexão); 2) Método Elemento de Reflexão Comum (ERC). Utiliza-se uma configuração onde um conjunto de sismogramas é relacionado com um único ponto de reflexão, e considera-se uma frente de onda hipoteticamente originada neste ponto; 3) Método Elemento de Evoluta Comum (EEC). Utiliza-se uma configuração onde cada sismograma está relacionado com um par de fonte e geofone coincidentemente posicionados na linha sísmica, e considera-se uma frente de onda hipoteticamente originada no centro de curvatura do refletor. Em cada um desses métodos tem-se como resultados uma seção sísmica empilhada, u(x_o, t_o); e outras duas seções denominadas de radiusgrama, R_o (x_o, t_o), e angulograma, β_o(x_o, t_o), onde estão os valores de raios de curvatura e ângulos de emergência da frente de onda considerada no instante em que a mesma atinge a superfície de observação, respectivamente. No caso do método denominado elemento refletor comum (ERC), a seção sísmica resultante do empilhamento corresponde a seção afastamento nulo. Pode-se mostrar que o sinal sísmico não sofre efeitos de alongamento como consequência da correção temporal, nem tão pouco apresenta problemas de dispersão de pontos de reflexão como consequência da inclinação do refletor, ao contrário do que acontece com as técnicas de empilhamento que tem por base a correção NMO. Além disto, por não necessitar de um macro modelo de velocidades a técnica de imageamento homeomórfico, de um modo geral, pode também ser aplicada a modelos heterogêneos, sem perder o rigor em sua formulação. Aqui também são apresentados exemplos de aplicação dos métodos elemento de fonte comum (EFC) (KEYDAR, 1993), e elemento refletor comum (ERC) (STEENTOFT, 1993), ambos os casos com dados sintéticos. No primeiro caso, (EFC), onde o empilhamento é feito tendo como referência um raio central arbitrário, pode-se observar um alto nível de exatidão no imageamento obtido, além do que é dada uma interpretação para as seções de radiusgrama e angulograma, de modo a se caracterizar aspectos geométricos do model geofísico em questão. No segundo caso, (ERC), o método é aplicado a série de dados Marmousi, gerados pelo método das diferenças finitas, e o resultado é comparado com aquele obtido por métodos convecionais (NMO/DMO) aplicados aos mesmos dados. Como consequência, observa-se que através do método ERC pode-se melhor detectar a continuidade de refletores, enquanto que através dos métodos convencionais caracterizam-se melhor a ocorrência de difrações. Por sua vez, as seções de radiusgrama e angulograma, no método (ERC), apresentam um baixo poder de resolução nas regiões do modelo onde se tem um alto grau de complexidade das estruturas. Finalmente, apresenta-se uma formulação unificada que abrange os diferentes métodos de imageamento homeomórfico citados anteriormente, e também situações mais gerais onde a frente de onda não se aproxima a um círculo, mas a uma curva quadrática qualquer.