Produção Científica
Artigo em Revista
| Métodos estocásticos para modelagem de escoamento estacionário e transiente em meios porosos. Uma maneira de incorporar a incerteza das medidas de campo e a variabilidade espacial nas propriedades hidráulicas de aquÃferos e estabelecer distribuições de probabilidade para os parâmetros fÃsicos do meio. As propriedades estatÃsticas do potencial hidráulico são então determinadas pela solução numérica das equações diferenciais estocásticas que regem o regime de escoamento no aquÃfero. Neste trabalho descrevemos a utilização de dois métodos de elementos finitos estocásticos (o método de Monte Carlo e o método da Colocação) para estimar a média e a variância do potencial hidráulico para fluxo saturado em meio estatisticamente heterogêneo, supondo que o coeficiente de transmissividade hidráulica é descrito por um processo lognormal. Um dos fatores decisivos na precisão numérica dos métodos é o comprimento de correlação associado à transmissividade. Discutimos também algumas configurações de baterias de extração que foram recentemente propostas para uma exploração adequada do aquÃfero Recôncavo na bacia do rio Capivara (Bahia, Brasil), comparando dois modelos clássicos e introduzindo a aleatoriedade da transmissividade em um dos arranjos de poços que foram propostos. |
Artigo em Revista
| Influence of Sea Water Resistivity on MCSEM Data. O Marine Controlled Source ElectroMagnetic (MCSEM) é um método geofÃsico para a detecção de camadas resistivas contidas abaixo do assoalhooceânico. Neste trabalho nós modelamos os dados do MCSEM em um ambiente unidimensional, incluindo variações na resistividade da água do mar, na forma de camadas de resistividades uniformes. Estas variações na resistividade da água podem surgir devido à influência de correntes marinhas, gradientes de temperatura ou mudança na salinidade da água. Nós estudamos o efeito destas variações nos dados do método MCSEM. Nossos resultados mostram que a interpretação pode ser fortemente influenciada, principalmente quando são analisados os dados normalizados. Vimos que as mudanças na resistividade da água têm um efeito sobre os dados similar à quele de mudanças na profundidade da lâmina d'água, sendo que ambos influenciam na atenuação da chamada "air-wave". |
Artigo em Revista
| Estudo de formas implementacionais da equação da onda imagem para remigração na profundidade. Neste trabalho estudamos teórica e numericamente novas formas implementacionais da equação da onda imagem para remigração em profundidade. Trata-se de uma equação diferencial parcial de segunda ordem semelhante à equação da onda acústica. Determinamos, além da consistência e da estabilidade, as condições de dispersão e dissipação de seis esquemas de diferenças finitas para esta equação em diversas formas obtidas por mudanças de variáveis. Estas condições não podem ser satisfeitas simultaneamente, ie., um resultado preciso não é de fácil obtenção. Testes numéricos confirmam os resultados teóricos da estabilidade para três esquemas, mas falham nos outros. Os esquemas avançados na variável da velocidade e avançados ou atrasados em profundidade foram os que tiveram os resultados mais proveitosos, pois estes esquemas satisfazem as previsões teóricas quanto à sua estabilidade e podem ser utilizados para realizar a propagação da onda imagem, sendo o primeiro para velocidades crescentes e o segundo para velocidades decrescentes. |
Artigo em Revista
| Wave equation depth migration using complex Padé approximation. Propomos um novo método de migração em profundidade baseado na solução da equação da onda com densidade constante no domÃnio da freqüência. Uma aproximação de Padé complexa é usada para aproximar o operador de evolução aplicado na extrapolação do campo de ondas. Esse método reduz as imprecisões e instabilidades devido à s ondas evanescentes e produz imagens com menos ruÃdos numéricos que aquelas obtidas usando-se a aproximação de Padé real para o operador exponencial, principalmente em meios com fortes variações de velocidades. Testes em dados de afastamento nulo do modelo de sal SEG/EAGE e nos dados de tiro comum 2-D Marmousi foram realizados. Os resultados obtidos mostram que o método de migração proposto consegue lidar com fortes variações laterais e também tem uma boa resposta para refletores com mergulhos Ãngremes. Os resultados foram comparados à queles resultados obtidos com os métodos split-step Fourier (SSF), phase shift plus interpolarion (PSPI) e Fourier diferenças-finitas (FFD). |
Artigo em Revista
| Filtragem SVD aplicada à melhoria do rastreamento de horizontes sÃsmicos. No presente artigo apresentamos uma aplicação da filtragem SVD (Singular Value Decomposition ) para o mapeamento automático de horizontes sÃsmico. A filtragem SVD pode ser vista como um método de filtragem multicanal onde cada traço filtrado guarda certo grau de coerência com os traços imediatamente vizinhos. Esta filtragem preserva as relações de amplitude, fase e correlação espacial dos eventos sismicos, ao tempo em que permite eliminar o ruÃdo incoerente, normalmente associado aos últimos autovalores. A decomposição SVD é realizada sobre o subconjunto de traços vizinhos a cada traço da linha sÃsmica 2D ou de um volume 3D. O traço filtrado é obtido utilizando apenas alguns dos autovetores e autovalores associados. Ilustramos a aplicação do método sobre dados sÃsmicos terrestres. A melhoria da coerência dos eventos sÃsmicos permitiu maior robustez ao autotracking no mapeamento e interpretação automática dos horizontes sÃsmicos. A filtragem SVD e Ì computacionalmente eficiente e tem o mérito de melhorar significativamente a coerência, a consistência e a continuidade dos eventos de reflexão facilitando muito o “trabalhoâ€, do tracker na busca de padrões no processo de autotracking. |
Artigo em Revista
| Traveltime approximations for q-P waves in vertical transversely isotropy media. As exploration targets have become deeper, cable lengths have increased accordingly, making the conventional two term hyperbolic traveltime approximation produce increasingly erroneous traveltimes. To overcome this problem, many traveltime formulas were proposed in the literature that provide approximations of different quality. In this paper, we concentrate on simple traveltime approximations that depend on a single anisotropy parameter. We give an overview of a collection of such traveltime approximations found in the literature and compare their quality. Moreover, we propose some new single-parameter traveltime approximations based on the approximations found in the literature. The main advantage of our approximations is that some of them are rather simple analytic expressions that make them easy to use, while achieving the same quality as the better of the established formulas. |
Artigo em Revista
| Anisotropic complex Pade hybrid finite-difference depth migration. Standard real-valued finite-difference (FD) and Fourier finite-difference (FFD) migrations cannot handle evanescent waves correctly, which can lead to numerical instabilities in the presence of strong velocity variations. A possible solution to these problems is the complex Padé approximation, which avoids problems with evanescent waves by rotating the branch cut of the complex square root. We have applied this approximation to the acoustic wave equation for vertical transversely isotropic media to derive more stable FD and hybrid FD/FFD migrations for such media. Our analysis of the dispersion relation of the new method indicates that it should provide more stable migration results with fewer artifacts and higher accuracy at steep dips. Our studies lead to the conclusion that the rotation angle of the branch cut that should yield the most stable image is 60° for FD migration, as confirmed by numerical impulse responses and work with synthetic data. |
Artigo em Revista
| Including Lateral Velocity Variations into True-amplitude Common-Shot Wave-equation Migration. In heterogeneous media, standard one-way wave equations describe only the kinematic part of one-way wave propagation correctly. For a correct description of amplitudes, the one-way wave equations must be modified. In media with vertical velocity variations only, the resulting true-amplitude one-way wave equations can be solved analytically. In media with lateral velocity variations, these equations are much harder to solve and require sophisticated numerical techniques. We present an approach to circumvent these problems by implementing approximate solutions based on the one-dimensional analytic amplitude modifications. We use these approximations to show how to modify conventional migration methods such as split-step and Fourier finite-difference migrations in such a way that they more accurately handle migration amplitudes. Simple synthetic data examples in media with a constant vertical gradient demonstrate that the correction achieves the recovery of true migration amplitudes. Applications to the SEG/EAGE salt model and the Marmousi data show that the technique improves amplitude recovery in the migrated images in more realistic situations. |
Artigo em Revista
| Obliquity-correction imaging condition for reverse time migration. The quality of seismic images obtained by reverse time migration (RTM) strongly depends on the imaging condition. We propose a new imaging condition that is motivated by stationary phase analysis of the classical crosscorrelation imaging condition. Its implementation requires the Poynting vector of the source and receiver wavefields at the imaging point. An obliquity correction is added to compensate for the reflector dip effect on amplitudes of RTM. Numerical experiments show that using an imaging condition with obliquity compensation improves reverse time migration by reducing backscattering artifacts and improving illumination compensation. |
Artigo em Revista
| The role of fluids in triggering earthquakes: observations from reservoir induced seismicity in Brazil. We relocate precisely micro-earthquakes induced by the Açu reservoir in Brazil and observe seismicity migration consistent with pore-pressure diffusion on a single fault zone. Fluids are believed to play a major role in triggering tectonic earthquakes; reservoir induced seismicity provides a natural laboratory in which to investigate the spatio-temporal evolution and triggering of earthquakes caused by fluid and pore-pressure diffusion. Between 1994 and 1997, 267 earthquakes (ML≤ 2.1) were recorded and located beneath the Açu reservoir. The seismicity increased several months following annual water level peaks, implying that pore-pressure diffusion is the principal triggering mechanism. The small station spacing and very low-attenuation, Precambrian basement, rock enabled starting earthquake locations with uncertainties of only a few hundred metres. We relocate 155 earthquakes from the largest cluster at Açu using waveform cross-correlation to obtain groups of similar events. We use these groups to improve the pick accuracy (to subsample accuracy in 200 sample per second data), and then invert for more accurate hypocentral locations. Our uncertainties are on the order of 10 m, and our locations are more tightly clustered. We observe temporal migration of the earthquakes, both along strike, and to increasing depth. We observe a seismicity migration rate between 15 and 58 m d–1. The rate is highest during the time of peak seismicity rate, and there is some suggestion that the rate decreases with increasing depth. Peak depth in seismicity is reached 175 d after the water peak, that is 192 d after the water low, and the maximum depth then decreases at a similar rate to the rate of increase. Our observations are consistent with triggering by pore-pressure diffusion within a heterogeneous fault zone with an average hydraulic diffusivity of ∼0.06 m2 s–1 and fracture permeability of ∼6 × 10−16 m2. |
