Shear-wave splitting and azimuthal variations of seismic attributes are two major anisotropic effects induced by vertically-aligned fractures. They both have influences on seismic data processing and interpretation, and provide information on fracture properties. Azimuthal variations in P-wave data have been intensively studied to improve imaging and obtain fracture parameters. However, azimuthal variations in PS-converted wave seismic data, particularly the velocity variation in PS-converted wave data, have not been well studied. Shear-wave splitting has been frequently used to estimate fracture directions and densities. However, its influence on the azimuthal variations of PS-converted wave data has also lacked a proper analysis. In this thesis, I analyse the anisotropic behaviour of PS-converted wave seismic data in the presence of azimuthal anisotropy, which includes the azimuthal variation of the PSconverted wave and PS-converted wave splitting. First, I demonstrate the robustness of PS-converted wave splitting for fracture characterisation. PS-converted wave seismic data is also influenced by the splitting effect due to its upgoing shear-wave leg. This important feature enables the application of shear-wave splitting analysis to PS-converted wave seismic data. I use synthetic data to show the necessity for separation of the split PS-converted waves. Then I apply the PS-converted wave splitting analysis to Sanhu 3D3C land seismic data. By separation of the fast and slow PS-converted waves and compensation for the time delays, the imaging quality has been improved. Dominant fracture properties obtained from the splitting analysis show a good correlation with the stress-field data. However, this work is accomplished by assuming only one set of vertical fractures in processing a given time window. In future work a specific layer-stripping algorithm could be constructed and applied. . Second, I study azimuthal variations of velocities in PS-converted wave seismic data. It involves two major parts: analysing azimuthal variations of NMO velocities to improve imaging, and examining the sensitivity of azimuthal variations to different fluid saturations. For a layer with HTI anisotropy induced by a set of vertical fractures, seismologists usually analyse the azimuthal behaviour exhibited on the radial and transverse components, on which PS-converted wave data are recorded. However, PS-converted waves also undergo shear-wave splitting, which complicates the azimuthal variations of PS-converted wave data. I demonstrate that it is essential to separate the fast P-SV1 wave from the slow P-SV2 wave, before applying any azimuthal analysis. I derive an equation describing the azimuthal variation in PSconverted wave NMO velocities, which shows the variation can be approximated into an ellipse. Based on this theory, I build a workflow to analyse the azimuthal variations of velocities in PS-converted wave data and apply this workflow to synthetic data. The imaging quality can be improved by using this workflow. Different fluid saturations in fractures have different influences on the azimuthal variations of both P-wave and PS-converted wave data. I perform a numerical study to understand how dry or water-saturated fractures control the azimuthal variations. Through theoretical and synthetic studies, I find that the azimuthal variation of velocities in PS-converted wave data is sensitive to different fluid saturations. By analysing the azimuthal variation, the fracture properties can also be estimated, but results are not as robust as those from PS-converted wave splitting analysis. I find that azimuthal variations of fast P-SV1 and slow P-SV2 waves show in-phase characteristics in dry fractures, but exhibit out-of-phase characteristics in water-saturated fractures. This important feature could open a new application for using PS-converted wave seismic data to distinguish oil-filled fractures from gas-filled fractures. In cases where multiple HTI layers are involved, I have developed a specific layer-stripping method to analyse both azimuthal variations and splitting effects of PS-converted waves. By applying this method to synthetic data, the fracture properties of each HTI layer can be estimated. The analysis of azimuthal variations in PS-converted wave velocities is applied to Daqing 3D3C land data. By using azimuthal velocity models in the PS-converted wave seismic data processing, the imaging quality is improved, especially in the anticline area where intensive fractures are likely to be developed. Furthermore, all fracture information obtained from analysis of azimuthal variations and splitting effects is compared with the stress-field data. The results from splitting analysis show a better correlation with the stress-field study. Finally, it is important to conclude that the analysis of PS-converted wave splitting is a robust method to estimate fracture directions and densities. However, it is not sensitive to different fluid saturations, which limits its application to fractured reservoir characterisation. Azimuthal variations of PS-converted wave seismic data can be analysed to improve imaging quality. Moreover their sensitivity to fluid saturations may provide a new way to discriminate between oil-filled and gas-filled fractures. However, the analysis of azimuthal variations is not as robust as the analysis of splitting effects, and it may require appropriate calibration with other fracture characterisation methods.

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.