Produção Científica

**Artigo em Revista**

DetecÃ§Ã£o de camadas delgadas usando sÃsmica de reflexÃ£o.A melhoria da resoluÃ§Ã£o vertical no dado sÃsmico tem se constituÃdo em um dos maiores desafios para o mÃ©todo de reflexÃ£o sÃsmica. O grande problema estÃ¡ relacionado Ã maior perda das amplitudes das altas freqÃ¼Ãªncias, quando comparada Ã s baixas, durante a propagaÃ§Ã£o da onda. Isto faz com que na composiÃ§Ã£o final do pulso sÃsmico a contribuiÃ§Ã£o das baixas freqÃ¼Ãªncias seja maior do que a das altas. Em princÃpio, isto nÃ£o deveria se constituir em obstÃ¡culo para a recuperaÃ§Ã£o plena do espectro, na medida em que o processamento sÃsmico dispÃµe de vÃ¡rias ferramentas para a equalizaÃ§Ã£o de todas as freqÃ¼Ãªncias presentes no volume sÃsmico. O que impede que essas aÃ§Ãµes logrem Ãªxito Ã© a presenÃ§a de um ruÃdo com comportamento randÃ´mico, o que faz com que exista uma freqÃ¼Ãªncia limite, freqÃ¼Ãªncia crÃtica, a partir da qual nÃ£o Ã© possÃvel recuperar o sinal de forma a efetivamente contribuir para o encurtamento do pulso no tempo. Tentativas de contornar este problema, a partir da atenuaÃ§Ã£o do ruÃdo aleatÃ³rio, tÃªm se mostrado insuficientes. A alternativa a estas tentativas Ã© o aumento da multiplicidade de empilhamento. Como a perda de amplitude com a freqÃ¼Ãªncia Ã© uma relaÃ§Ã£o exponencial, Ã© necessÃ¡ria multiplicidade extremamente elevada para que o ganho seja efetivo. A capacidade de utilizaÃ§Ã£o de altÃssimas multiplicidades Ã© exatamente a grande virtude da tÃ©cnica CRS (Common Reflection Surface). Neste trabalho propomos a combinaÃ§Ã£o do mÃ©todo CRS com o chamado balanceamento espectral (spectral whitening), para a recuperaÃ§Ã£o de altas freqÃ¼Ãªncias, com o objetivo de aumentar ganhos na resoluÃ§Ã£o vertical e, como conseqÃ¼Ãªncia, propiciar melhor discriminaÃ§Ã£o de reservatÃ³rios delgados. Os primeiros testes em dados sintÃ©ticos e reais apresentados neste trabalho sÃ£o bastante encorajadores e permitem concluir que a estratÃ©gia proposta tem bom potencial de aplicaÃ§Ã£o prÃ¡tica. |

**Artigo em Revista**

Intraplate earthquake swarm in Belo Jardim, NE Brazil: reactivation of a major Neoproterozoic shear zone (Pernambuco Lineament).Intraplate earthquakes in stable continental areas have been explained basically by reactivation of pre-existing zones of weakness, stress concentration, or both. Zones of weakness are usually identified as sites of the last major orogeny, provinces of recent alkaline intrusions, or stretched crust in ancient rifts. However, it is difficult to identify specific zones of weakness and intraplate fault zones are not always easily correlated with known geological features. Although Northeastern Brazil is one of the most seismically active areas in the country (magnitudes 5 roughly every 5 yr), with hypocentral depths shallower than âˆ¼10 km and seismic zones as long as 30â€“40 km, no clear relationship with the known surface geology can be usually established with confidence, and a clear identification of zones of weakness has not yet been possible. Here we present the first clear case of seismic activity occurring as reactivation of an old structure in Brazil: the Pernambuco Lineament, a major Neoproterozoic shear zone. The 2004 earthquake swarm of Belo Jardim (magnitudes up to 3.1) and the recurrent activities in the nearby towns of SÃ£o Caetano and Caruaru (magnitudes up to 4.0 and 3.8), show that the Pernambuco Lineament is a weak zone. A local seismic network showed that the Belo Jardim swarm of 2004 November occurred by normal faulting on a North dipping, Eâ€“W oriented fault plane in close agreement with the Eâ€“W trending structures within the Pernambuco Lineament. The Belo Jardim activity was concentrated in a 1.5 km (Eâ€“W) by 2 km (downdip) fault area, and average depth of 4.5 km. The nearby Caruaru activity occurs as both strike-slip and normal faulting, also consistent with local structures of the Pernambuco Lineament. The focal mechanisms of Belo Jardim, Caruaru and S. Caetano, indicate Eâ€“W compressional and Nâ€“S extensional principal stresses. The NS extension of this stress field is larger than that predicted by numerical models such as those of Coblentz & Richardson and we propose that additional factors such as flexural stresses from the nearby Sergipe-Alagoas marginal basin could also affect the current stress field in the Pernambuco Lineament. |

**Artigo em Revista**

Applications of Information Theory, Genetic Algorithms, and Neural Models to Predict Oil Flow.This work introduces a new information-theoretic methodology for choosing variables and their time lags in a prediction setting, particularly when neural networks are used in non-linear modeling. The first contribution of this work is the Cross Entropy Function (XEF) proposed to select input variables and their lags in order to compose the input vector of black-box prediction models. The proposed XEF method is more appropriate than the usually applied Cross Correlation Function (XCF) when the relationship among the input and output signals comes from a non-linear dynamic system. The second contribution is a method that minimizes the Joint Conditional Entropy (JCE) between the input and output variables by means of a Genetic Algorithm (GA). The aim is to take into account the dependence among the input variables when selecting the most appropriate set of inputs for a prediction problem. In short, theses methods can be used to assist the selection of input training data that have the necessary information to predict the target data. The proposed methods are applied to a petroleum engineering problem; predicting oil production. Experimental results obtained with a real-world dataset are presented demonstrating the feasibility and effectiveness of the method. |

**Artigo em Revista**

Evidence of hydraulic connectivity across deformation bands from field pumping tests: two examples from Tucano Basin, NE Brazil.It is assumed that deformation bands may compartmentalize aquifers or hydrocarbon reservoirs because these low-permeability structures may behave as barriers to fluid flow. To address the question whether there is, at a reservoir scale, hydraulic connectivity across a damage zone dominated by cataclastic deformation bands, we present results of two pumping tests carried out in a fluvial-deltaic phreatic sandstone aquifer from the Ilhas Group at Tucano Basin (NE Brazil) where intense concentration of deformation bands occurs. Both test sites are associated with macroscopic damage zones that are approximately 1 km in length and 15 m thick. In situ permeability measurements show values of 2000 mD for host rock and 0.1 mD for deformation bands. GPR profiles reveal good continuity of the primary sedimentary structures with almost no deformation band vertical offsets (less than 10 cm). The well locations for the pumping tests were chosen so that the damage zone is located between pumping and monitoring wells. Pumping tests in both cases revealed hydraulic connectivity across the damage zone since the observed stationary drawdown at monitoring wells was a considerable fraction of the drawdown observed in the pumped well. In one experiment using eight monitoring wells the drawdown cone is evolving through the damage zone instead of contouring it. Local deviations in the natural groundwater flow allow to say that the damage zone dimensions are quite large both in horizontal and vertical directions compared to the distances among the wells. The interpretation of all experimental results is that deformation bands do not fully compartmentalize the aquifer. Generalization of this result to hydrocarbon reservoirs has to take into account capillary effects which are not present in the studied case. |

**Artigo em Revista**

Ground-roll attenuation using a 2D time-derivative filter.We present a new filtering method for the attenuation of ground-roll. The method is based on the application of a bi-dimensional filter for obtaining the time-derivative of the seismograms. Before convolving the filter with the input data matrix, the normal moveout correction is applied to the seismograms with the purpose of flattening the reflections. The method can locally attenuate the amplitude of data of low frequency (in the ground-roll and stretch normal moveout region) and enhance flat events (reflections). The filtered seismograms can reveal horizontal or sub-horizontal reflections while vertical or sub-vertical events, associated with ground-roll, are attenuated. A regular set of samples around each neighbourhood data sample of the seismogram is used to estimate the time-derivative. A numerical approximation of the derivative is computed by taking the difference between the interpolated values calculated in both the positive and the negative neighbourhood of the desired position. The coefficients of the 2D time-derivative filter are obtained by taking the difference between two filters that interpolate at positive and negative times. Numerical results that use real seismic data show that the proposed method is effective and can reveal reflections masked by the ground-roll. Another benefit of the method is that the stretch mute, normally applied after the normal moveout correction, is unnecessary. The new filtering approach provides results of outstanding quality when compared to results obtained from the conventional FK filtering method. |

**Artigo em Revista**

Time evolution of the wave equation using rapid expansion method.Forward modeling of seismic data and reverse time migration are based on the time evolution of wavefields. For the case of spatially varying velocity, we have worked on two approaches to evaluate the time evolution of seismic wavefields. An exact solution for the constant-velocity acoustic wave equation can be used to simulate the pressure response at any time. For a spatially varying velocity, a one-step method can be developed where no intermediate time responses are required. Using this approach, we have solved for the pressure response at intermediate times and have developed a recursive solution. The solution has a very high degree of accuracy and can be reduced to various finite-difference time-derivative methods, depending on the approximations used. Although the two approaches are closely related, each has advantages, depending on the problem being solved. |

**Artigo em Revista**

Trigonal meshes in diffraction tomography with optimum regularization: an application for carbon sequestration monitoring. |

**Artigo em Revista**

Migration velocity analysis by double path-integral migration.The idea of path-integral imaging is to sum over the migrated images obtained for a set of migration velocity models. Those velocities where common-image gathers align horizontally are stationary, thus favoring these images in the overall stack. In this way, the overall image forms with no knowledge of the true velocity model. However, the velocity information associated with the final image can be determined in the process. By executing the path-integral imaging twice, weighting one of the stacks with the velocity value, the stationary velocities that produce the final image can then be extracted by a division of the two images. A numerical example demonstrates that quantitative information about the migration velocity model can be determined by double path-integral migration. |

**Artigo em Revista**

On the estimation of local slopes.Current time-processing algorithms often are based on one-parameter or multiparameter coherency analysis (semblance) schemes applied to the data. Such procedures, besides being computationally expensive, lead to significant uncertainties in the searched parameters. Conventional semblance methods can be avoided for a number of imaging tasks if local slopes can be extracted directly from prestack dataâ€”for example, by filtering schemes. Although the idea is not new, it has revived for various purposes, such as velocity analysis, Ï„-p imaging, migration to zero offset, and time migration. We propose a simple, straightforward correction to linear plane-wave destructors based on the observation that in addition to the local slope, its inverse can be extracted from the data in a fully analogous way. Combining the information of both extractions yields a simple yet effective correction to the local slopes. The naive application of linear plane-wave destructors with our correction produces high-quality results, even with a high noise level and interfering events. |

**Artigo em Revista**

Source parameters and rupture velocity of small? 2.1 reservoir induced earthquakes.We calculate stress drop and rupture speed for ML â‰¤ 2.1 shallow reservoir induced earthquakes and ï¬nd them to be similar to those of large, natural earthquakes. Previous studies have suggested that hydrofractures, mining and reservoir-induced earthquakes have lower average stress drop than natural tectonic earthquakes. This difference might result from the different tectonic setting or the shallower hypocentral depths of induced earthquakes. Alternatively, difï¬culties in correcting for attenuation and site effects in earlier studies may lead to underestimation of stress drop. In addition, most studies assume the rupture velocity of small reservoir induced earthquakes to be the same as for the large earthquakes. We analyse a set of 101 ML â‰¤ 2.1 earthquakes induced by changing water level in the AcÂ¸u Reservoir, NE Brazil. The earthquakes are shallow, (depth <5 km) and the region has negligible natural seismicity. We use three different approaches to calculate the source parameters of the six largest (1.9 â‰¤ ML â‰¤ 2.1) earthquakes. We model the individual spectra to ï¬nd corner frequency, frequency-independent Q, and long period amplitude. We use collocated small earthquakes as empirical Greenâ€™s functions to calculate the spectral ratios, and determine the relative source time functions. Estimates of the source duration and corner frequency imply stress drops in the range of 26â€“179 MPa. These are similar to, or higher than tectonic earthquakes suggesting that the shallow hypocentral depth and the presence of water do not affect stress drop. We observe clear directivity for one of the earthquakes, and use the azimuthal variation in pulse width to estimate a rupture velocity of â‰¥0.6Î². |