The use of shallow seismic reflection technique in near surface exploration of urban sites: An evaluation in the city of SãO Paulo, Brazil

In order to evaluate the use of shallow seismic technique to delineate geological and geotechnical features up to 40 meters depth in noisy urban areas covered with asphalt pavement, five survey lines were conducted in the metropolitan area of São Paulo City. The data were acquired using a 24-bit, 24-channel seismograph, 30 and 100 Hz geophones and a sledgehammer-plate system as seismic source. Seismic reflection data were recorded using a CMP (common mid point) acquisition method. The processing routine consisted of: prestack band-pass filtering (90-250 Hz); automatic gain control (AGC); muting (digital zeroin) of dead/noisy traces, ground roll, air-wave and refracted-wave; CMP sorting; velocity analyses; normal move-out corrections; residual static corrections; f-k filtering; CMP stacking. The near surface is geologically characterized by unconsolidated fill materials and Quaternary sediments with organic material overlying Tertiary sediments with the water table 2 to 5 m below the surface. The basement is composed of granite and gneiss. Reflections were observed from 40 milliseconds to 65 ms two-way traveltime and were related to the silt clay layer and fine sand layer contact of the Tertiary sediments and to the weathered basement. The CMP seismic-reflection technique has been shown to be useful for mapping the sedimentary layers and the bedrock of the São Paulo sedimentary basin for the purposes of shallow investigations related to engineering problems. In spite of the strong cultural noise observed in these urban areas and problems with planting geophones we verified that, with the proper equipment, field parameters and particularly great care in data collection and processing, we can overcome the adverse field conditions and to image reflections from layers as shallow as 20 meters.

Seismic activity triggered by water wells in the Parana Basin, Brazil

Triggered seismicity is commonly associated with deep water reservoirs or injection wells where water is injected at high pressure into the reservoir rock. However, earth tremors related solely to the opening of groundwater wells are extremely rare. Here we present a clear case of seismicity induced by pore-pressure changes following the drilling of water wells that exploit a confined aquifer in the intracratonic Parana Basin of southeastern Brazil. Since 2004, shallow seismic activity, with magnitudes up to 2.9 and intensities V MM, has been observed near deep wells (120-200 m) that were drilled in early 2003 near the town of Bebedouro. The wells were drilled for irrigation purposes, cross a sandstone layer about 60-80 m thick and extract water from a confined aquifer in fractured zones between basalt flow layers. Seismic activity, mainly event swarms, has occurred yearly since 2004, mostly during the rainy season when the wells are not pumped. During the dry season when the wells are pumped almost continuously, the activity is very low. A seismographic network, installed in March 2005, has located more than 2000 microearthquakes. The events are less than 1 km deep (mostly within the 0.5 km thick basalt layer) and cover an area roughly 1.5 km x 5 km across. The seismicity generally starts in a small area and expands to larger distances with an equivalent hydraulic diffusivity ranging from 0.06 to 0.6 m(2)/s. Geophysical and geothermal logging of several wells in the area showed that water from the shallow sandstone aquifer enters the well at the top and usually forms waterfalls. The waterfalls flow down the sides of the wells and feed the confined, fractured aquifer in the basalt layer at the bottom. Two seismic areas are observed: the main area surrounds several wells that are pumped continuously during the dry season, and a second area near another well (about 10 km from the first area) that is not used for irrigation and not pumped regularly. The main area displays cyclic annual activity, but the second area does not. We explain the earthquake swarms as being triggered by pore pressure diffusion in the fractured basalt layer due to additional pressure from the newly connected surface aquifer. This reaches critically prestressed areas up to a few kilometers away from the wells. During periods of continuous pumping, the reduction of pore pressure in the confined aquifer stops the seismic activity. Our study suggests that this kind of activity may be more common than previously thought and implies that many other cases of small tremors associated with the drilling of water wells may have gone unnoticed.

DIFFERENTIATING PATH AND SOURCE PROCESSES IN COMPLEX GEOLOGIC STRUCTURE THROUGH PASSIVE SEISMIC STUDIES: BERING GLACIER, ALASKA AND VILLARRICA VOLCANO, CHILE

Measuring shallow seismic sources provides a way to reveal processes that cannot be directly observed, but the correct interpretation and value of these signals depend on the ability to distinguish source from propagation effects. Furthermore, seismic signals produced by a resonating source can look almost identical to those produced by impulsive sources, but modified along the path. Distinguishing these two phenomena can be accomplished by examining the wavefield with small aperture arrays or by recording seismicity near to the source when possible. We examine source and path effects in two different environments: Bering Glacier, Alaska and Villarrica Volcano, Chile. Using three 3-element seismic arrays near the terminus of the Bering Glacier, we have identified and located both terminus calving and iceberg breakup events. We show that automated array analysis provided a robust way to locate icequake events using P waves. This analysis also showed that arrivals within the long-period codas were incoherent within the small aperture arrays, demonstrating that these codas previously attributed to crack resonance were in fact a result of a complicated path rather than a source effect. At Villarrica Volcano, seismometers deployed from near the vent to ~10 km revealed that a several cycle long-period source signal recorded at the vent appeared elongated in the far-field. We used data collected from the stations nearest to the vent to invert for the repetitive seismic source, and found it corresponded to a shallow force within the lava lake oriented N75°E and dipping 7° from horizontal. We also used this repetitive signal to search the data for additional seismic and infrasonic properties which included calculating seismic-acoustic delay times, volcano acoustic-seismic ratios and energies, event frequency, and real-time seismic amplitude measurements. These calculations revealed lava lake level and activity fluctuations consistent with lava lake level changes inferred from the persistent infrasonic tremor.

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

Aplicação de softwares livres no Processamento de dados sísmicos de alta resolução

The increasing use of shallow seismic methods of high resolution, for investigations of geological problems, environmental or industrial, has impelled the development of techniques, flows and computational algorithms. The practice of applying techniques for processing this data, until recently it wasn t used and the interpretation of the data was made as they were acquired. In order to facilitate and contribute to the improvement of the practices adopted, was developed a free graphical application and open source, called OpenSeismic which is based on free software Seismic Un*x, widely used in the treatment of conventional seismic data used in the exploration of hydrocarbon reservoirs. The data used to validate the initiative were marine seismic data of high resolution, acquired by the laboratory of Geology and Marine Geophysics and Environmental Monitoring - GGEMMA, of the Federal University of Rio Grande do Norte UFRN, for the SISPLAT Project, located at the region of paleo-valley of the Rio Acu. These data were submitted to the processing flow developed by Gomes (2009), using the free software developed in this work, the OpenSeismic, as well other free software, the Seismic Un*x and the commercial software ProMAX, where despite its peculiarities has presented similar results

Obtencao e interpretacao de dados de Reflexao Sismica Superficial (Rasa) - exemplos e aplicacoes

Successful application of shallow seismic reflection method is related directly to the ability of the ground to transmit high-frequency seismic energy from a seismic source. The dominant frequencies of reflection data are in the range of 50 to 100 Hz and depend on the surface materials and water table. A spread of geophones send the siemsic signal to be recorded on a 12, 24 or 48 channel portable seismograph, using a single high-frequency geophone per channel. Although it is possible to find seismographs with display and processing unities, it is also possible to transfer digitalized data to a personal computer to be processed and interpreted by using specific programs. Some results of two recent field studies to specify the underground structure in the tidal flats area of Baie St. Paul (Quebec-Canada) and in glacial terrains, in the Waterloo (Ontario-Canada) region. -from English summary

TDEM survey in an area of seismicity induced by water wells in Parana sedimentary basin, Northern Sao Paulo State, Brazil

This article presents TDEM results from an area with recent induced shallow seismicity. The purpose was to do a geoelectrical mapping of sedimentary and fractured basaltic aquifers for better understanding of the hydrogeologic setting. The study area is in the Parana basin where flood basalts are overlain by sedimentary units near the city of Bebedouro, northern Sao Paulo State, Brazil. 86 TDEM soundings were acquired in an area of 90 km(2) in the Andes and Botafogo study areas. The soundings were chosen next to wells for calibration, and also along profiles crossing the seismically active areas. 1-D interpretation results showed the general geoelectrical stratigraphy of this part of the Parana basin. The upper geoelectrical layer is the shallow sedimentary aquifer (Adamantina formation) with less than 80 m thickness. The second geoelectrical layer contains the upper basalts of the Serra Geral formation at about 60-80 m depths. A saturated fractured basalt zone between 100 and 300 m depths was identifiable on various TDEM soundings. This depth range corresponds to the range of hypocentral depths for more than 3000 micro-earthquakes in this area. The lower basalt layer was estimated to lie between 400 and 650 m depth. The deepest geoelectrical layer detected by various TDEM soundings corresponds to the Botucatu sandstone (Guarani aquifer). Results suggest that the high-discharge wells are located in the fractured zone in the middle basalt of the Serra Geral formation. There is a good correlation between seismically active areas, high discharge wells (>190 m(3)/h), and fracture zones in the middle basalt. The results reinforce the hypothesis that the shallow seismic activity in the Bebedouro region is being triggered by high rates of groundwater withdrawal. (C) 2012 Elsevier B.V. All rights reserved.

Maps and acoustic profiles of free gas distribution, Western Baltic Sea and Kattgat-Skagerrak area

A shallow gas depth-contour map covering the Skagerrak-western Baltic Sea region has been constructed using a relatively dense grid of existing shallow seismic lines. The digital map is stored as an ESRI shape file in order to facilitate comparison with other data from the region. Free gas usually occurs in mud and sandy mud but is observed only when sediment thickness exceeds a certain threshold value, depending on the water depth of the area in question. Gassy sediments exist at all water depths from approx. 20 m in the coastal waters of the Kattegat to 360 m in the Skagerrak. In spite of the large difference in water depths, the depth of free gas below seabed varies only little within the region, indicating a relatively fast movement of methane in the gas phase towards the seabed compared to the rate of diffusion of dissolved methane. Seeps of old microbial methane occur in the northern Kattegat where a relatively thin cover of sandy sediments exists over shallow, glacially deformed Pleistocene marine sediments. Previous estimates of total methane escape from the area may be correct but the extrapolation of local methane seepage rate data to much larger areas on the continental shelf is probably not justified. Preliminary data on porewater chemistry were compared with the free gas depth contours in the Aarhus Bay area, which occasionally suffers from oxygen deficiency, in order to examine if acoustic gas mapping may be used for monitoring the condition of the bay.

Caracterização morfológica-sedimentar do vale inciso Apodi-Mossoró e plataforma continental adjacente - Bacia Potiguar Offshore

Systems of incised valleys have been studied in different continental shelves, including the Brazilian continental margin. The interest to characterize this feature is given by the information that it can provide variations on sea level, as well as the ability to host economically significant quantities of hydrocarbons in reservoirs located in deposits filling of the incised valleys. This thesis has the overall objective to characterize the morphology and sedimentary cover of the incised valley Apodi-Mossoró, located in the Northern Continental shelf of Rio Grande do Norte state, adjacent to Areia Branca city. The methodology included the integration of satellite imagery, bathymetric data, sedimentological data, shallow seismic, and the identification of foraminifera. The results indicate that the ApodiMossró incised valley is currently formed by two channels, shallow channel and deep channel, which have distinct morphological and sedimentological characteristics. The deep channel has connection with one of the heads of the Apodi Canyon, located in the slope area. The acquisition, processing and interpretation of shallow seismic data allowed the recognition of the depositional surface, erosional surface, discordance, and sismofaceis. The erosional surface mapped from shallow seismic sections is possibly a indicative of an ancient surface of valley incision, where it would probably be associated with the limit Pleistocene/Holocene. Different sismofaceis were identified and reflect the rise in sea level with standards sometimes agradacional, sometimes progradational. The thickness of sediments on this surface was estimated at a maximum of 22m thick in the central portion of the incised valley. Statistically, there are differences between the adjacent continental shelf and channels, and between these channels, for the content of calcium carbonate, organic matter, sand and mud perceptual, except for the gravel grain size. The analysis of living and dead foraminifera showed the presence of fifty species distributed in regards to morphology, depth and type of sediment. Four type of seismic echocharacteres were identified and mapped, as well as their bedforms, indicating different sedimentary processes along the incised valley. The integration of results suggests an activation of the Apodi-Mossoró incised valley in the Late Pleistocene.

Segmentation and kinematics of the North America-Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high-resolution multibeam echo-sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate-boundary structures are a series of strike-slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre-existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike-slip regime. Along the most recent trace of the SOFZ, we measured a strike-slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS-derived motion of 9.8 ± 2 mm a−1 has remained stable during the entire Quaternary.

Bandwidth requirements for shallow, high-resolution seismic reflection

The optimum bandwidth for shallow, high-resolution seismic reflection differs from that required for conventional petroleum reflection. An understanding of this issue is essential for correct choice of acquisition instrumentation. Numerical modelling of simple Bowen Basin coal structures illustrates that, for high-resolution imaging, it is important to accurately record all frequencies up to the limit imposed by earth scattering. On the contrary, the seismic image is much less dependent on frequencies at the lower end of the spectrum. These quantitative observations support the use of specialised high-frequency geophones for high-resolution seismic imaging. Synthetic seismic inversion trials demonstrate that, irrespective of the bandwidth of the seismic data, additional low-frequency impedance control is essential for accurate inversion. Inversion provides no compelling argument for the use of conventional petroleum geophones in the high-resolution arena.

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

Estimation of inelastic seismic material properties of a surgical sea-bed from multi-component marine seismic data

To date, state-of-the-art seismic material parameter estimates from multi-component sea-bed seismic data are based on the assumption that the sea-bed consists of a fully elastic half-space. In reality, however, the shallow sea-bed generally consists of soft, unconsolidated sediments that are characterized by strong to very strong seismic attenuation. To explore the potential implications, we apply a state-of-the-art elastic decomposition algorithm to synthetic data for a range of canonical sea-bed models consisting of a viscoelastic half-space of varying attenuation. We find that in the presence of strong seismic attenuation, as quantified by Q-values of 10 or less, significant errors arise in the conventional elastic estimation of seismic properties. Tests on synthetic data indicate that these errors can be largely avoided by accounting for the inherent attenuation of the seafloor when estimating the seismic parameters. This can be achieved by replacing the real-valued expressions for the elastic moduli in the governing equations in the parameter estimation by their complex-valued viscoelastic equivalents. The practical application of our parameter procedure yields realistic estimates of the elastic seismic material properties of the shallow sea-bed, while the corresponding Q-estimates seem to be biased towards too low values, particularly for S-waves. Given that the estimation of inelastic material parameters is notoriously difficult, particularly in the immediate vicinity of the sea-bed, this is expected to be of interest and importance for civil and ocean engineering purposes.

Basement depths in the Parecis Basin (Amazon) with receiver functions from small local earthquakes in the Porto dos Gauchos seismic zone

Receiver functions from small local earthquakes were used to determine sediment thicknesses in Porto dos Gauchos seismic zone (PGSZ), Parecis basin, Amazonian craton, Brazil. The high velocity contrast between basement and sediments (P-wave velocities of 6.1 and 3.0 km/s, respectively) favors the generation of clear P-to-S converted phases (Ps) seen in the radial component, and also S-to-P conversions (Sp) seen in the vertical component. A reference 10 velocity model determined with shallow refraction experiment in PGSZ helped to convert Ps P time differences to basement depths at 15 stations deployed for aftershocks studies. The results of receiver function integrated with the shallow refraction reveal that the basement depths in the PGSZ increases from the basin border in the north up to about 600 m depth in the south. The basement topography, however, does not vary smoothly and a basement high with a steep topography was detected near the epicentral area. A 400 m elevation difference within 1.7 km distance suggests a possible border fault of a buried graben. This feature seems to be oriented roughly WSW-ENE and could indicate basement structures related to the seismicity of the Porto dos Gauchos Seismic Zone. (C) 2011 Elsevier Ltd. All rights reserved.

The Oligocene-Neogene deep-sea Columbia Channel system in the South Brazilian Basin: Seismic stratigraphy and environmental changes

The Columbia Channel (CCS) system is a depositional system located in the South Brazilian Basin, south of the Vitoria-Trindade volcanic chain. It lies in a WNW-ESE direction on the continental rise and abyssal plain, at a depth of between 4200 and 5200 m. It is formed by two depocenters elongated respectively south and north of the channel that show different sediment patterns. The area is swept by a deep western boundary current formed by AABW. The system has been previously interpreted has a mixed turbidite-contourite system. More detailed study of seismic data permits a more precise definition of the modern channel morphology, the system stratigraphy as well as the sedimentary processes and control. The modern CCS presents active erosion and/or transport along the channel. The ancient Oligo-Neogene system overlies a ""upper Cretaceous-Paleogene"" sedimentary substratum (Unit U1) bounded at the top by a major erosive ""late Eocene-early Oligocene"" discordance (D2). This ancient system is subdivided into 2 seismic units (U2 and U3). The thick basal U2 unit constitutes the larger part of the system. It consists of three subunits bounded by unconformities: D3 (""Oligocene-Miocene boundary""), D4 (""late Miocene"") and D5 (""late Pliocene""). The subunits have a fairly tabular geometry in the shallow NW depocenter associated with predominant turbidite deposits. They present a mounded shape in the deep NE depocenter, and are interpreted as forming a contourite drift. South of the channel, the deposits are interpreted as a contourite sheet drift. The surficial U3 unit forms a thin carpet of deposits. The beginning of the channel occurs at the end of U1 and during the formation of D2. Its location seems to have been determined by active faults. The channel has been active throughout the late Oligocene and Neogene and its depth increased continuously as a consequence of erosion of the channel floor and deposit aggradation along its margins. Such a mixed turbidite-contourite system (or fan drift) is characterized by frequent, rapid lateral facies variations and by unconformities that cross the whole system and are associated with increased AABW circulation. (C) 2009 Elsevier B.V. All rights reserved.