We propose an efficient iterative data interpolation method using continuity along reflectors in seismic images via curvelet and discrete cosine transforms. The curvelet transform is a new multiscale transform that provides sparse representations for images that comprise smooth objects separated by piece-wise smooth discontinuities (e.g. seismic images). The advantage of using curvelets is that these frames are sparse for high-frequency caustic-free solutions of the wave-equation. Since we are dealing with less than ideal data (e.g. bandwidth-limited), we compliment the curvelet frames with the discrete cosine transform. The latter is motivated by the successful data continuation with the discrete Fourier transform. By choosing generic basis functions we circumvent the necessity to make parametric assumptions (e.g. through linear/parabolic Radon or demigration) regarding the shape of events in seismic data. Synthetic and real data examples demonstrate that our algorithm provides interpolated traces that accurately reproduce the wavelet shape as well as the AVO behavior along events in shot gathers.

A spherical wave AVO response is investigated by measuring ultrasonic reflection amplitudes from a water/Plexiglas interface. The experimental results show substantial deviation from the plane-wave reflection coefficients at large angles. However there is an excellent agreement between experimental data and full-wave numerical simulations performed with the reflectivity algorithm. By comparing the spherical-wave AVO response, modeled with different frequencies, to the plane-wave response, I show that the differences between the two are of such magnitude that three-term AVO inversion based on AVA curvature can be erroneous. I then propose an alternative approach to use critical angle information extracted from AVA curves, and show that this leads to a significant improvement of the estimation of elastic parameters. Azimuthal variation of the AVO response of a vertically fractured model also shows good agreement with anisotropic reflectivity simulations, especially in terms of extracted critical angles which indicated that (1) reflection measurements are consistent with the transmission measurements; (2) the anisotropic numerical simulation algorithm is capable of simulating subtle azimuthal variations with excellent accuracy; (3) the methodology of picking critical angles on seismograms using the inflection point is robust, even in the presence of random and/or systematic noise.

I conducted an amplitude vs. offset (AVO) analysis on newly acquired 3D seismic reflection data to detect elevated pore fluid content and pore fluid pressure along the Costa Rica convergent margin to address dewatering processes of subduction zone sediments. These data provide the highest quality 3D seismic data acquired to date along a convergent margin for detailed analysis of geophysical properties along the plate boundary fault interface. In 2011, a 55 km by 11 km 3D seismic reflection survey was completed using the R/V Marcus G. Langseth offshore western Costa Rica at the convergent margin of the Cocos and Caribbean plates. We applied pre-stack Kirchhoff time migration to a subset of these data across the frontal prism where amplitude versus offset (AVO) attributes were extracted along the decollement. When pore fluid pressure, l , exceeds 0.7, the pressure at which Poisson’s ratio begins to approach that of water, the AVO response of a fluid-filled, clay-rich decollement requires a high Poisson’s ratio and an excessively low seismic P-wave and S-wave velocity. Acute wedge taper, undercompacted subducted hemipelagic and pelagic sediments, and a smooth decollement in the northwest half of the survey correspond with decollement AVO response of relatively high values of Poisson’s ratio. These findings suggest increased pore fluid content and vertical containment of near-lithostatic pore fluid pressures within the decollement. In contrast, increased wedge taper angles, thin hemipelagic and pelagic sediments, and a rugose decollement beneath the southeastern frontal prism produce an AVO response interpreted as due to lower pore fluid contents and pressures. We propose that large-offset subducting basement normal faults in this area, as close as 20 m from the decollement, induce vertical fractures within the decollement that allow for fluid expulsion into the frontal prism and lower fluid pressure. Lateral variability of overpressure within the decollement shear zone of subduction margins is important in understanding the evolution of frontal prism strain accumulation and seismogenic rupture. / text

The Marcellus Shale is an important resource play prevalent in several states in the eastern United States. The productive zone of the Marcellus Shale has variations in rock properties such as clay content, kerogen content and pore aspect ratio, and these variations may strongly effect elastic anisotropy. The objective of this study is to characterize surface seismic sensitivity for variations in anisotropic parameters relating to kerogen content and aspect ratio of kerogen saturated pores. The recognized sensitivity may aid to characterize these reservoir from surface seismic observations for exploration and production of hydrocarbon. In this study, I performed VTI anisotropic modeling based on geophysical wireline log data from Harrison County, WV. The wireline log data includes spectral gamma, density, resistivity, neutron porosity, monopole and dipole sonic logs. Borehole log data were analyzed to characterize the Marcellus Shale interval, and quantify petrophysical properties such as clay content, kerogen content and porosity. A rock physics model was employed to build link between petrophysical properties and elastic constants. The rock physics model utilized differential effective medium (DEM) theory, bounds and mixing laws and fluid substitution equations in a model scheme to compute elastic constants for known variations in matrix composition, kerogen content and pore shape distribution. The seismic simulations were conducted applying a vertical impulse source and three component receivers. The anisotropic effect to angular amplitude variations for PP, PS and SS reflections were found to be dominantly controlled by the Thomsen Ɛ parameter, characterizing seismic velocity variations with propagation direction. These anisotropic effect to PP data can be seen at large offset (>15o incidence angle). The most sensitive portion of PS reflections was observed at mid offset (15o-30o). I also analyzed seismic sensitivity for variations in kerogen content and aspect ratio of structural kerogen. Elastic constants were computed for 5%, 10%, 20% and 30% kerogen content from rock physics model and provided to the seismic model. For both kerogen content and aspect ratio model, PP amplitudes varies significantly at zero to near offset while PS amplitude varied at mid offsets (12 to 30 degree angle of incidences). / text

We present a robust iterative sparseness-constrained interpolation algorithm using 2/3D curvelet frames and Fourier-like transforms that exploits continuity along reflectors in seismic data. By choosing generic transforms, we circumvent the necessity to make parametric assumptions (e.g. through linear/parabolic Radon or demigration) regarding the shape of events in seismic data. Simulation and real data examples for data with moderately sized gaps demonstrate that our algorithm provides interpolated traces that accurately reproduce the wavelet shape as well as the AVO behavior. Our method also shows good results for de-aliasing judged by the behavior of the (f-k)-spectrum before and after regularization.

Neste trabalho avaliou-se a potencialidade do uso da sísmica rasa de reflexão multicomponente para investigação geológica-geotécnica. Foram abordados vários aspectos relacionados à aquisição dos dados sísmicos de reflexão multicomponente, com o objetivo de entender as vantagens e limitações do método para aplicação em investigações de subsuperfície rasa. Os ensaios de campo foram realizados em duas áreas, ambas em terrenos da Bacia Sedimentar de São Paulo, em área urbana da cidade de São Paulo. Para a interpretação dos dados sísmicos multicomponente foram investigados procedimentos para a inversão não-linear dos tempos de trânsito e das amplitudes. O testes realizados orientaram a escolha da aproximação não-hiperbólica mais apropriada para calculo dos tempos de trânsito visando à análise de velocidades do pacote acima do refletor. O estudo numérico desenvolvido para a inversão das amplitudes mostrou a viabilidade da estimativa das velocidades e densidades, acima e abaixo do refletor, empregando-se as equações de Zoeppritz para as ondas refletidas PP, PSv, SvP e SvSv, antes e depois do ângulo crítico. Dada a complexidade da inversão nao-linear das amplitudes, se fez necessário elaborar uma estratégia estocástica de otimização e desenvolver uma nova abordagem para análise da função objetivo multi-dimensional, garantindo confiabilidade ao resultado da inversão não-linear. Os resultados deste trabalho mostraram o potencial da sísmica de reflexão rasa multicomponente para caracterização geológica-geotécnica, possibilitando um melhor entendimento das camadas superficiais. / This thesis aims to evaluate the useful of the multicomponent seismic methods for shallow investigations, mainly its potential for the geotechnical and geological characterization of the nearsurface. Several aspects regarding the acquisition and processing data of multicomponent seismic data are discussed. They were based on data set acquired in the urban area of Sao Paulo city, Brazil. Two different areas were investigated. Both located in sedimentary terrains belonging to the Sao Paulo Sedimentary Basin. We present a non-linear travel time and seismic amplitude inversion scheme to quantitative interpretation of multicomponent seismic data. Several tests were performed to guide the choice of non-hyperbolic equation more suitable for travel time inversion aiming the velocity analysis above the reflector. A numerical experiment developed to solve the nonlinear inversion of seismic amplitudes showed the feasibility to estimate seismic interval velocities and layer densities above and below the reflector using the exact Zoeppritz equations for PP, PSv, SvP e SvSv reflected waves, before and after critical angle. Due to the apparent complexity of the nonlinear seismic amplitude inversion, it was necessary elaborate the strategy for stochastic optimization and develop a new approach to analyze the multi-dimensional objective function, with different implications for the accuracy and efficiency of the non-linear inversion. The study show the benefits of using the multicomponent seismic method for shallow geological-geotechnical characterization, improving the nearsurface understanding, once allows an integrated analyzes of a more complete record of the wave field.

Neste trabalho avaliou-se a potencialidade do uso da sísmica rasa de reflexão multicomponente para investigação geológica-geotécnica. Foram abordados vários aspectos relacionados à aquisição dos dados sísmicos de reflexão multicomponente, com o objetivo de entender as vantagens e limitações do método para aplicação em investigações de subsuperfície rasa. Os ensaios de campo foram realizados em duas áreas, ambas em terrenos da Bacia Sedimentar de São Paulo, em área urbana da cidade de São Paulo. Para a interpretação dos dados sísmicos multicomponente foram investigados procedimentos para a inversão não-linear dos tempos de trânsito e das amplitudes. O testes realizados orientaram a escolha da aproximação não-hiperbólica mais apropriada para calculo dos tempos de trânsito visando à análise de velocidades do pacote acima do refletor. O estudo numérico desenvolvido para a inversão das amplitudes mostrou a viabilidade da estimativa das velocidades e densidades, acima e abaixo do refletor, empregando-se as equações de Zoeppritz para as ondas refletidas PP, PSv, SvP e SvSv, antes e depois do ângulo crítico. Dada a complexidade da inversão nao-linear das amplitudes, se fez necessário elaborar uma estratégia estocástica de otimização e desenvolver uma nova abordagem para análise da função objetivo multi-dimensional, garantindo confiabilidade ao resultado da inversão não-linear. Os resultados deste trabalho mostraram o potencial da sísmica de reflexão rasa multicomponente para caracterização geológica-geotécnica, possibilitando um melhor entendimento das camadas superficiais. / This thesis aims to evaluate the useful of the multicomponent seismic methods for shallow investigations, mainly its potential for the geotechnical and geological characterization of the nearsurface. Several aspects regarding the acquisition and processing data of multicomponent seismic data are discussed. They were based on data set acquired in the urban area of Sao Paulo city, Brazil. Two different areas were investigated. Both located in sedimentary terrains belonging to the Sao Paulo Sedimentary Basin. We present a non-linear travel time and seismic amplitude inversion scheme to quantitative interpretation of multicomponent seismic data. Several tests were performed to guide the choice of non-hyperbolic equation more suitable for travel time inversion aiming the velocity analysis above the reflector. A numerical experiment developed to solve the nonlinear inversion of seismic amplitudes showed the feasibility to estimate seismic interval velocities and layer densities above and below the reflector using the exact Zoeppritz equations for PP, PSv, SvP e SvSv reflected waves, before and after critical angle. Due to the apparent complexity of the nonlinear seismic amplitude inversion, it was necessary elaborate the strategy for stochastic optimization and develop a new approach to analyze the multi-dimensional objective function, with different implications for the accuracy and efficiency of the non-linear inversion. The study show the benefits of using the multicomponent seismic method for shallow geological-geotechnical characterization, improving the nearsurface understanding, once allows an integrated analyzes of a more complete record of the wave field.

En exploration pétrolière, la technique d'imagerie la plus couramment utilisée reste l'imagerie par sismique réflexion. Elle permet d'obtenir une première représentation des structures géologiques mais aussi une estimation du coefficient de réflexion des interfaces du sous-sol. Les variations de ce coefficient dépendent des paramètres pétrophysiques des milieux situés de part et d'autre de la discontinuité donnant naissance à la réflexion (densité et vitesses des ondes de propagation et de cisaillement). Une analyse de ces variations permet ainsi d'évaluer les paramètres élastiques qui sont des données indispensables et complémentaires à l'interprétation des réservoirs. Les traitements classiques consistent à étudier ces variations directement sur des données enregistrées en fonction du déport et du temps de propagation. L'analyse s'effectue par inversion d'une approximation linéaire des équations de Zoeppritz régissant les variations théoriques de l'amplitude. Nous proposons ici d'étudier les variations de l'amplitude après la phase d'imagerie, directement dans le domaine profondeur. Pour réaliser cette analyse, nous proposons un traitement spécifique après l'imagerie et avant l'analyse de l'amplitude pour corriger les impacts liés à la qualité du modèle de vitesse, à l'interférence des ondelettes entre deux réflexions proches et à l'étirement de l'ondelette. Ce traitement a aussi permis l'obtention d'un squelette du sous-sol conduisant à une première interprétation en profondeur des réflecteurs majeurs.

The CO2SINK (and CO2MAN) project is the first onshore CO2 storage project in Europe. The research site is located near the town of Ketzin, close to Potsdam in Germany. Injection started in June 2008, with a planned injection target of 100,000 tonnes of CO2. In February 2011, around 45, 000 tons of CO2 had been injected into the saline aquifer at an approximate depth of 630 m. This thesis focuses on time-lapse analysis of borehole seismic data, surface seismic data and reservoir characterization at the Ketzin site. Baseline Moving Source Profiling (MSP) data were acquired in the borehole Ketzin 202/2007 (OW2), along seven lines in 2007. The zero-offset Vertical Seismic Profile (VSP) data were acquired in the same borehole. The main objective of the VSP and MSP survey was to generate high-resolution seismic images around the borehole. After modeling and data processing, the sandy layers within the Stuttgart Formation can potentially be imaged in the VSP and MSP data whereas reflections from these layers are not as clearly observed in the 3D surface seismic data. 2D and pseudo-3D time-lapse seismic surveys were conducted at the Ketzin site. Interpretation of 2D baseline and repeat stacks shows that no CO2 leakage related time lapse signature is observable where the 2D lines allow monitoring of the reservoir. This is consistent with the time-lapse results of the 3D surveys showing an increase in reflection amplitude just centered around the injection well. The results from the pseudo-3D surveys are also consistent with the 3D seismic time-lapse studies and show that the sparse pseudo-3D geometry can be used to qualitatively map the CO2 in the reservoir with significantly less effect than the full 3D surveying. The 2nd pseudo-3D repeat survey indicates preferential migration of the CO2 to the west. There are no indications of migration into the caprock on either of the repeat surveys. Amplitude Versus Offset (AVO) analysis was performed on both 2D and 3D repeat surveys. A Class 3 AVO anomaly is clearly observed on the 3D repeat data and matches the synthetic modeling well. No AVO anomaly was observed on the 2D repeat data, which was anticipated, but the result shows signs of a pressure response at the reservoir level in the data. Reflection coefficients were calculated using surface seismic data (3D and pseudo-3D) at the site. Pre-injection calculations agree well with calculations from logging data. Post-injection calculations are in general agreement with the seismic modeling, but generally show higher amplitudes than those expected. The full 3D data show a better image of the reflection coefficients before and after injection than the pseudo-3D data and can potentially be used to make quantitative calculations of CO2 volumes. The pseudo-3D data only provide qualitative information.

I here present a combined AVO analysis of P-P and P-S reflection data whose objective is to improve the identification of lithology by estimating the specific values of Poisson's ratio, [sigma], for each rock formation in a given geological model, rather than a contrast between formations. Limited knowledge on the elastic parameters of a given rock formation and difficulty regarding the availability and processing of P-S data constitute hindrances of lithology identification. Considering that ocean bottom seismology (OBS) has aided in solving the problem of P-S data availability, limited information on elastic parameters is still a challenge, and the focus of this thesis. The present analysis is based on Zoeppritz' solution for the P-P and P-S reflection coefficients, RPP and RPS, with a slight modification. We used the normalized P-S reflection coefficient; i.e., R'PS = RPS / sin [theta] for [theta] > 0, instead of RPS, where [theta] is the incident angle. By normalizing RPS, we avoid dealing with the absence of converted S-waves at small incident angles and enhance the similar linear behavior of the P-P and normalized P-S reflection coefficients at small angles of incidence. We have used the linearity of RPP and R'PS at angles smaller than 35 degrees to simultaneously estimate the average VP/VS ratio, the contrasts of P- and S-wave velocities, and the contrast of density. Using this information, we solve for Poisson's ratio of each formation, which may enable lithology discrimination. The feasibility of this analysis was demonstrated using nonlinear synthetic data (i.e., finite-difference data). The results in estimating Poisson's ratio yielded less than 5 percent error. We generalize this new combined P-P and P-S AVO analysis for dipping interfaces. Similarly to the nondipping interface case, our derivations show that the amplitude variation with offset (AVO) of P-P and P-S for a dipping interface can be cast into intercepts and gradients. However, these intercepts and gradients depend on the angle of the dipping interface. Therefore, we further generalize our analysis by including a migration step that allows us to find the dipping angle. Because seismic data is not available in terms of RPP and R'PS, this process includes recovery of reflection coefficients after migrating the data and correcting for geometrical spreading, as done by Ikelle et al. (1986 and 1988). The combination of all of these steps, namely geometrical-spreading correction, migration, and AVO analysis, is another novelty of this thesis, which leads to finding the specific values of Poisson's ratio of each rock formation directly from the seismic data.