Estudo da estrutura subsuperficial da província Borborema com correlação de ruído sísmico

Ambient seismic noise has traditionally been considered as an unwanted perturbation in seismic data acquisition that "contaminates" the clean recording of earthquakes. Over the last decade, however, it has been demonstrated that consistent information about the subsurface structure can be extracted from cross-correlation of ambient seismic noise. In this context, the rules are reversed: the ambient seismic noise becomes the desired seismic signal, while earthquakes become the unwanted perturbation that needs to be removed. At periods lower than 30 s, the spectrum of ambient seismic noise is dominated by microseism, which originates from distant atmospheric perturbations over the oceans. The microsseism is the most continuous seismic signal and can be classified as primary – when observed in the range 10-20 s – and secondary – when observed in the range 5-10 s. The Green‘s function of the propagating medium between two receivers (seismic stations) can be reconstructed by cross-correlating seismic noise simultaneously recorded at the receivers. The reconstruction of the Green‘s function is generally proportional to the surface-wave portion of the seismic wavefield, as microsseismic energy travels mostly as surface-waves. In this work, 194 Green‘s functions obtained from stacking of one month of daily cross-correlations of ambient seismic noise recorded in the vertical component of several pairs of broadband seismic stations in Northeast Brazil are presented. The daily cross-correlations were stacked using a timefrequency, phase-weighted scheme that enhances weak coherent signals by reducing incoherent noise. The cross-correlations show that, as expected, the emerged signal is dominated by Rayleigh waves, with dispersion velocities being reliably measured for periods ranging between 5 and 20 s. Both permanent stations from a monitoring seismic network and temporary stations from past passive experiments in the region are considered, resulting in a combined network of 33 stations separated by distances between 60 and 1311 km, approximately. The Rayleigh-wave, dispersion velocity measurements are then used to develop tomographic images of group velocity variation for the Borborema Province of Northeast Brazil. The tomographic maps allow to satisfactorily map buried structural features in the region. At short periods (~5 s) the images reflect shallow crustal structure, clearly delineating intra-continental and marginal sedimentary basins, as well as portions of important shear zones traversing the Borborema Province. At longer periods (10 – 20 s) the images are sensitive to deeper structure in the upper crust, and most of the shallower anomalies fade away. Interestingly, some of them do persist. The deep anomalies do not correlate with either the location of Cenozoic volcanism and uplift - which marked the evolution of the Borborema Province in the Cenozoic - or available maps of surface heat-flow, and the origin of the deep anomalies remains enigmatic.

Análise sísmica usando transformada de Curvelet

Oil prospecting is one of most complex and important features of oil industry Direct prospecting methods like drilling well logs are very expensive, in consequence indirect methods are preferred. Among the indirect prospecting techniques the seismic imaging is a relevant method. Seismic method is based on artificial seismic waves that are generated, go through the geologic medium suffering diffraction and reflexion and return to the surface where they are recorded and analyzed to construct seismograms. However, the seismogram contains not only actual geologic information, but also noise, and one of the main components of the noise is the ground roll. Noise attenuation is essential for a good geologic interpretation of the seismogram. It is common to study seismograms by using time-frequency transformations that map the seismic signal into a frequency space where it is easier to remove or attenuate noise. After that, data is reconstructed in the original space in such a way that geologic structures are shown in more detail. In addition, the curvelet transform is a new and effective spectral transformation that have been used in the analysis of complex data. In this work, we employ the curvelet transform to represent geologic data using basis functions that are directional in space. This particular basis can represent more effectively two dimensional objects with contours and lines. The curvelet analysis maps real space into frequencies scales and angular sectors in such way that we can distinguish in detail the sub-spaces where is the noise and remove the coefficients corresponding to the undesired data. In this work we develop and apply the denoising analysis to remove the ground roll of seismograms. We apply this technique to a artificial seismogram and to a real one. In both cases we obtain a good noise attenuation

Anisotropia sísmica crustal na região de Cascavel - CE

This dissertation presents a study on crustal seismic anisotropy in Cascavel - CE. The earthquake data employed here are from the Seismological Laboratory at Universidade Federal do Rio Grande do Norte (UFRN) and were colected from 29 September 1997 to 05 march 1998 using six three-component digital seismographic stations. In general, the cause of the observed seismic anisotropy in many regions of the world is interpreted in terms of fluid-filled stress aligned microcracks in the rockmass (EDA). In other words, the polarisation direction of the faster shear-wave splitting is parallel to SHmax. However, other researches on seismic anisotropy carried out in NE Brazil have shown a remarkable consistency of the faster shear-wave polarisation direction with the direction of the Precambrian fabric. The present work is another case study that is used to investigate this issue. In order to map the Precambrian fabric we used aeromagnetic data, since the study area is mostly covered with sediments (up to 50m thick) and in-situ field mapping would be very difficult to be carried out. According to the results from the present research, the observations of the faster shear-wave polarisation directions in two seismographic stations in Cascavel region are best explained in the framework of EDA. For the remaining two stations, the observed anisotropy may have two interpretions: (i) - 90_ flips of the direction of polarisation of the faster shear-wave, since that the event-to-station ray path would be through the fracture zone and hence would travel under a higher pore pressure and (ii) - the observed seismic anisotropy would agree with the direction on the ductile Precambrian fabric

Estudo da atividade sísmica em São Caetano-PE em 2007

In this dissertation we studied the seismic activity in the São Caetano county, Pernambuco State, Northeastern Brazil, located near the Pernambuco Lineament. The Pernambuco Lineament is a one of Neoproterozoic continental-scale shear zones that deforms the Borborema province. The seismicity estudied occurred in a NE trending branch of Pernambuco Lineament. The seismic activity in São Caetano started in 2006 and in May 20th, 2006 a 4,0 mb earthquake hit there. This was the largest earthquake ever reported in Pernambuco State. This dissertation is the result of a campaign done in the period from Februay 1th 2007 to July 31 th 2007. In this campaign up to nine three-component digital seismographic stations were deployed and the collected data was used to determine hypocenters and focal mechanism. A total of 214 earthquakes, recorded at least by three stations, were analyzed. To determine hypocenters and time origin the HYPO71 program was used assuming a half-space model with parameters : VP (P-wave velocity) equal to 5.90 km/s and the ratio VP/VS 1.70, where VS is the S-wave velocity. The earthquakes hypocentral distribution was approximately 4 km long and agrees with the NE-SW direction of the Pernambuco Lineamento branch. Hypocentres depth range from 2 to 8 km. The composed focal mechanism was made from a group of 14 selected earthquakes. We try firstly to find the fault plane solution matching the polarity distribution at stations, using the FPFIT program. The result was 43 deg ± 15 deg for strike, 59 deg ± 9 deg for dip and -142 deg ± 15 deg for rake. We also fitted a plane using the hypocentral distribution to obtain the dip and azimuth of the hypocentral distribution. The results obtained by this fit were 58 deg for the azimuth, 55 deg for the dip and -155 deg for rake. This result shows a mechanism of a strike-slip dextral fault with a normal component. This dissertation shows, once more, that there is a good correlation between the seismic activity and geological features in the region near the Pernambuco Lineament and its NE branches

Aquisição, processamento e análise de dados de sísmica de alta resolução na plataforma continental norte do Rio Grande do Norte: Vale inciso do rio açu

On the modern Continental Shelf to the north of Rio Grande do Norte state (NE Brazil) is located a paleo-valley, submerged during the last glacial sea-level lowstand, that marks continuation of the most important river of this area (Açu River). Despite the high level of exploration activity of oil industry, there is few information about shallow stratigraphy. Aiming to fill this gap, situated on the Neogene, was worked a marine seismic investigation, the development of a processing flow for high resolution data seismic, and the recognition of the main feature morphology of the study area: the incised valley of the River Açu. The acquisition of shallow seismic data was undertaken in conjunction with the laboratory of Marine Geology/Geophysics and Environmental Monitoring - GGEMMA of Federal University of Rio Grande do Norte UFRN, in SISPLAT project, where the geomorphological structure of the Rio paleovale Açu was the target of the investigation survey. The acquisition of geophysical data has been over the longitudinal and transverse sections, which were subsequently submitted to the processing, hitherto little-used and / or few addressed in the literature, which provided a much higher quality result with the raw data. Once proposed for the flow data was developed and applied to the data of X-Star (acoustic sensor), using available resources of the program ReflexW 4.5 A surface fluvial architecture has been constructed from the bathymetric data and remote sensing image fused and draped over Digital Elevation Models to create three-dimensional (3D) perspective views that are used to analyze the 3D geometry geological features and provide the mapping morphologically defined. The results are expressed in the analysis of seismic sections that extend over the region of the continental shelf and upper slope from mouth of the Açu River to the shelf edge, providing the identification / quantification of geometrical features such as depth, thickness, horizons and units seismic stratigraphyc area, with emphasis has been placed on the palaeoenvironmental interpretation of discordance limit and fill sediment of the incised valley, control by structural elements, and marked by the influence of changes in the sea level. The interpretation of the evolution of this river is worth can bring information to enable more precise descriptions and interpretations, which describes the palaeoenvironmental controls influencing incised valley evolution and preservation to provide a better comprehensive understanding of this reservoir analog system

Interpretação sísmica e análise de atributos - área do levantamento sísmico 3D de siririzinho (Bacia Sergipe-Alagoas, NE do Brasil)

The main objective of the present thesis was the seismic interpretation and seismic attribute analysis of the 3D seismic data from the Siririzinho high, located in the Sergipe Sub-basin (southern portion of Sergipe-Alagoas Basin). This study has enabled a better understanding of the stratigraphy and structure that the Siririzinho high experienced during its development. In a first analysis, we used two types of filters: the dip-steered median filter, was used to remove random noise and increase the lateral continuity of reflections, and fault-enhancement filter was applied to enhance the reflection discontinuities. After this filtering step similarity and curvature attributes were applied in order to identify and enhance the distribution of faults and fractures. The use of attributes and filtering greatly contributed to the identification and enhancement of continuity of faults. Besides the application of typical attributes (similarity and curvature) neural network and fingerprint techniques were also used, which generate meta-attributes, also aiming to highlight the faults; however, the results were not satisfactory. In a subsequent step, well log and seismic data analysis were performed, which allowed the understanding of the distribution and arrangement of sequences that occur in the Siririzinho high, as well as an understanding of how these units are affected by main structures in the region. The Siririzinho high comprises an elongated structure elongated in the NS direction, capped by four seismo-sequences (informally named, from bottom to top, the sequences I to IV, plus the top of the basement). It was possible to recognize the main NS-oriented faults, which especially affect the sequences I and II, and faults oriented NE-SW, that reach the younger sequences, III and IV. Finally, with the interpretation of seismic horizons corresponding to each of these sequences, it was possible to define a better understanding of geometry, deposition and structural relations in the area.

Sismicidade e esforços tectônicos na zona sísmica Acaraú, Nordeste do Brasil

Inside the Borborema Province the Northwestern Ceará (NC) is one of the most seismic active regions. There are reports of an earthquake occurred in 1810 in the Granja town. On January, 2008 the seismic activity in NC has increased and it was deployed a seismographic network with 11 digital stations. In 2009, another earthquake sequence began and it was deployed another seismographic network in the Santana do Acaraú town with 6 stations. This thesis presents the results obtained by analyzing the data recorded in these two networks. The epicentral areas are located near the northeastern part of the Transbrasiliano Lineament, a shear zone with NE-SW-trending that cuts the study area. The hypocenters are located between 1km and 8km. The strike-slip focal mechanisms were found, which is predominant in the Borborema Province. An integration of seismological, geological and geophysical data was performed and it show that the seismogenic faults found are oriented in the same direction to the local brittle structures observed in field and magnetic lineaments. The SHmax (maximum compressional stress) direction in NC was estimated using an inversion of seven focal mechanisms. The horizontal maximum compression stress (σ1 = 300°) with orientation NW-SE and extension (σ3 = 210°) with NE-SW and σ2 vertical. These results are consistent with results of previous studies. The seismic activity recorded in NC is not related to a possible reactivation of the Transbrasiliano Lineament, by now.

Medium change monitoring using ambient seismic noise and coda wave interferometry: examples from intraplate NE Brazil and the Mid-Atlantic Ridge

This thesis presents and discusses the results of ambient seismic noise correlation for two different environments: intraplate and Mid-Atlantic Ridge. The coda wave interferometry method has also been tested for the intraplate data. Ambient noise correlation is a method that allows to retrieve the structural response between two receivers from ambient noise records, as if one of the station was a virtual source. It has been largely used in seismology to image the subsurface and to monitor structural changes associated mostly with volcanic eruptions and large earthquakes. In the intraplate study, we were able to detect localized structural changes related to a small earthquake swarm, which main event is mR 3.7, North-East of Brazil. We also showed that the 1-bit normalization and spectral whitening result on the loss of waveform details and that the phase auto-correlation, which is amplitude unbiased, seems to be more sensitive and robust for our analysis of a small earthquake swarm. The analysis of 6 months of data using cross-correlations detect clear medium changes soon after the main event while the auto-correlations detect changes essentially after 1 month. It could be explained by fluid pressure redistribution which can be initiated by hydromechanical changes and opened path ways to shallower depth levels due to later occurring earthquakes. In the Mid-Atlantic Ridge study, we investigate structural changes associated with a mb 4.9 earthquake in the region of the Saint Paul transform fault. The data have been recorded by a single broadband seismic station located at less than 200 km from the Mid-Atlantic ridge. The results of the phase auto-correlation for a 5-month period, show a strong co-seismic medium change followed by a relatively fast post-seismic recovery. This medium change is likely related to the damages caused by the earthquake’s ground shaking. The healing process (filling of the new cracks) that lasted 60 days can be decomposed in two phases, a fast recovery (70% in ~30 days) in the early post-seismic stage and a relatively slow recovery later (30% in ~30 days). In the coda wave interferometry study, we monitor temporal changes of the subsurface caused by the small intraplate earthquake swarm mentioned previously. The method was first validated with synthetics data. We were able to detect a change of 2.5% in the source position and a 15% decrease of the scatterers’ amount. Then, from the real data, we observed a rapid decorrelation of the seismic coda after the mR 3.7 seismic event. This indicates a rapid change of the subsurface in the fault’s region induced by the earthquake.

Medium change monitoring using ambient seismic noise and coda wave interferometry: examples from intraplate NE Brazil and the Mid-Atlantic Ridge

This thesis presents and discusses the results of ambient seismic noise correlation for two different environments: intraplate and Mid-Atlantic Ridge. The coda wave interferometry method has also been tested for the intraplate data. Ambient noise correlation is a method that allows to retrieve the structural response between two receivers from ambient noise records, as if one of the station was a virtual source. It has been largely used in seismology to image the subsurface and to monitor structural changes associated mostly with volcanic eruptions and large earthquakes. In the intraplate study, we were able to detect localized structural changes related to a small earthquake swarm, which main event is mR 3.7, North-East of Brazil. We also showed that the 1-bit normalization and spectral whitening result on the loss of waveform details and that the phase auto-correlation, which is amplitude unbiased, seems to be more sensitive and robust for our analysis of a small earthquake swarm. The analysis of 6 months of data using cross-correlations detect clear medium changes soon after the main event while the auto-correlations detect changes essentially after 1 month. It could be explained by fluid pressure redistribution which can be initiated by hydromechanical changes and opened path ways to shallower depth levels due to later occurring earthquakes. In the Mid-Atlantic Ridge study, we investigate structural changes associated with a mb 4.9 earthquake in the region of the Saint Paul transform fault. The data have been recorded by a single broadband seismic station located at less than 200 km from the Mid-Atlantic ridge. The results of the phase auto-correlation for a 5-month period, show a strong co-seismic medium change followed by a relatively fast post-seismic recovery. This medium change is likely related to the damages caused by the earthquake’s ground shaking. The healing process (filling of the new cracks) that lasted 60 days can be decomposed in two phases, a fast recovery (70% in ~30 days) in the early post-seismic stage and a relatively slow recovery later (30% in ~30 days). In the coda wave interferometry study, we monitor temporal changes of the subsurface caused by the small intraplate earthquake swarm mentioned previously. The method was first validated with synthetics data. We were able to detect a change of 2.5% in the source position and a 15% decrease of the scatterers’ amount. Then, from the real data, we observed a rapid decorrelation of the seismic coda after the mR 3.7 seismic event. This indicates a rapid change of the subsurface in the fault’s region induced by the earthquake.

Filtro kl/dvs para atenuação do ruído "ground roll"

In geophysics there are several steps in the study of the Earth, one of them is the processing of seismic records. These records are obtained through observations made on the earth surface and are useful for information about the structure and composition of the inaccessible parts in great depths. Most of the tools and techniques developed for such studies has been applied in academic projects. The big problem is that the seismic processing power unwanted, recorded by receivers that do not bring any kind of information related to the reflectors can mask the information and/or generate erroneous information from the subsurface. This energy is known as unwanted seismic noise. To reduce the noise and improve a signal indicating a reflection, without losing desirable signals is sometimes a problem of difficult solution. The project aims to get rid of the ground roll noise, which shows a pattern characterized by low frequency, low rate of decay, low velocity and high amplituds. The Karhunen-Loève Transform is a great tool for identification of patterns based on the eigenvalues and eigenvectors. Together with the Karhunen-Loève Transform we will be using the Singular Value Decomposition, since it is a great mathematical technique for manipulating data

Análise da trajetória de poços horizontais em reservatório de óleo pesado

The application of thermal methods, to increase the recovery of heavy oil in mature fields through drainage with multilateral and horizontal wells, has been thoroughly studied, theorically, experimentally, testing new tools and methods. The continuous injection of steam, through a steam injector well and a horizontal producer well in order to improve horizontal sweep of the fluid reservoir, it is an efficient method. Starting from an heterogeneous model, geologically characterized, modeling geostatistics, set history and identification of the best path of permeability, with seismic 3D, has been dubbed a studying model. It was studied horizontal wells in various directions in relation to the steam and the channel of higher permeability, in eight different depths. Into in the same area were studied, the sensitivity of the trajectories of horizontal wells, according to the depth of navigation. With the purpose of obtaining the highest output of oil to a particular flow, quality, temperature and time for the injection of steam. The wells studied showed a significant improvement in the cumulative oil recovery in one of the paths by promoting an alternative to application in mature fields or under development fields with heavy oil

Remoção de ruídos sísmicos utilizando transformada de wavelet 1D e 2D com software em desenvolvimento

In the Hydrocarbon exploration activities, the great enigma is the location of the deposits. Great efforts are undertaken in an attempt to better identify them, locate them and at the same time, enhance cost-effectiveness relationship of extraction of oil. Seismic methods are the most widely used because they are indirect, i.e., probing the subsurface layers without invading them. Seismogram is the representation of the Earth s interior and its structures through a conveniently disposed arrangement of the data obtained by seismic reflection. A major problem in this representation is the intensity and variety of present noise in the seismogram, as the surface bearing noise that contaminates the relevant signals, and may mask the desired information, brought by waves scattered in deeper regions of the geological layers. It was developed a tool to suppress these noises based on wavelet transform 1D and 2D. The Java language program makes the separation of seismic images considering the directions (horizontal, vertical, mixed or local) and bands of wavelengths that form these images, using the Daubechies Wavelets, Auto-resolution and Tensor Product of wavelet bases. Besides, it was developed the option in a single image, using the tensor product of two-dimensional wavelets or one-wavelet tensor product by identities. In the latter case, we have the wavelet decomposition in a two dimensional signal in a single direction. This decomposition has allowed to lengthen a certain direction the two-dimensional Wavelets, correcting the effects of scales by applying Auto-resolutions. In other words, it has been improved the treatment of a seismic image using 1D wavelet and 2D wavelet at different stages of Auto-resolution. It was also implemented improvements in the display of images associated with breakdowns in each Auto-resolution, facilitating the choices of images with the signals of interest for image reconstruction without noise. The program was tested with real data and the results were good

Subsídios para o projeto de estruturas sismo resistentes

Through this research is detailed the Brazilian seismic code focused on concrete projects design related to seismic engineering. At the beginning of the research is debated the fundaments of the seismic effects, the influence factors to the development of seismic effects and also relates the main data registration happened in Brazil. The second step is study the Brazilian seismic code explaining all the concepts related to it and does a compilation to the most important international seismic code. At this research is developed the designing of a building submitted to horizontal equivalent seismic forces and the modal process based on the answer spectrum based on the brazilin seismic code. It was also developed the design of a commercial building submitted to seismic loads based on the Brazilian code answer spectrum and compared to the same building submitted to wind loads.The research also focus on projects conception and detailing of seismic engineering Project design. At the study of seismic engineering it was concluded that seismic effects require special focus on concrete structures design, proving that is the essential consider the seismic effects

Escalabilidade Paralela de um Algoritmo de Migração Reversa no Tempo (RTM) Pré-empilhamento

The seismic method is of extreme importance in geophysics. Mainly associated with oil exploration, this line of research focuses most of all investment in this area. The acquisition, processing and interpretation of seismic data are the parts that instantiate a seismic study. Seismic processing in particular is focused on the imaging that represents the geological structures in subsurface. Seismic processing has evolved significantly in recent decades due to the demands of the oil industry, and also due to the technological advances of hardware that achieved higher storage and digital information processing capabilities, which enabled the development of more sophisticated processing algorithms such as the ones that use of parallel architectures. One of the most important steps in seismic processing is imaging. Migration of seismic data is one of the techniques used for imaging, with the goal of obtaining a seismic section image that represents the geological structures the most accurately and faithfully as possible. The result of migration is a 2D or 3D image which it is possible to identify faults and salt domes among other structures of interest, such as potential hydrocarbon reservoirs. However, a migration fulfilled with quality and accuracy may be a long time consuming process, due to the mathematical algorithm heuristics and the extensive amount of data inputs and outputs involved in this process, which may take days, weeks and even months of uninterrupted execution on the supercomputers, representing large computational and financial costs, that could derail the implementation of these methods. Aiming at performance improvement, this work conducted the core parallelization of a Reverse Time Migration (RTM) algorithm, using the parallel programming model Open Multi-Processing (OpenMP), due to the large computational effort required by this migration technique. Furthermore, analyzes such as speedup, efficiency were performed, and ultimately, the identification of the algorithmic scalability degree with respect to the technological advancement expected by future processors

Morfologia em alta resolução espacial de estrelas de tipo Mira e R Corenae Borealis

Mira and R Coronae Borealis (R CrB) variable stars are evolved objects surrounded by circumstellar envelopes (CSE) composed of the ejected stellar material. We present a detailed high-spatial resolution morfological study of the CSE of three stars: IRC+10216, the closest and more studied Carbon-Rich Mira; o Ceti, the prototype of the Mira class; and RY Sagitarii (RY Sgr), the brightest R CrB variable of the south hemisphere. JHKL near-infrared adaptive optics images of IRC+10216 with high dynamic range and Vband images with high angular resolution and high depth, collected with the VLT/NACO and VLT/FORS1 instruments, were analyzed. NACO images of o Ceti were also analyzed. Interferometric observations of RY Sgr collected with the VLTI/MIDI instrument allowed us to explore its CSE innermost regions (»20 40 mas). The CSE of IRC+10216 exhibit, in near-infrared, clumps with more complex relative displacements than proposed in previous studies. In V-band, the majority of the non-concentric shells, located in the outer CSE layers, seem to be composed of thinner elongated shells. In a global view, the morphological connection between the shells and the bipolar core of the nebulae, located in the outer layers, together with the clumps, located in the innermost regions, has a difficult interpretation. In the CSE of o Ceti, preliminar results would be indicating the presence of possible clumps. In the innermost regions (.110 UA) of the CSE of RY Sgr, two clouds were detected in different epochs, embedded in a variable gaussian envelope. Based on a rigorous verification, the first cloud was located at »100 R¤ (or »30 AU) from the centre, toward the east-north-east direction (modulo 180o) and the second one was almost at a perpendicular direction, having aproximately 2£ the distance of the first cloud. This study introduces new constraints to the mass-loss history of these kind of variables and to the morphology of their innermost CSE regions