Do seismic waves and fluid flow sense the same permeability?

We study the discrepancy between the effective flow permeability and the effective seismic permeability, that is, the effective permeability controlling seismic attenuation due to wave-induced fluid flow, in 2D rock samples having mesoscopic heterogeneities and in the presence of strong permeability fluctuations. In order to do so, we employ a numerical oscillatory compressibility test to determine attenuation and velocity dispersion due to wave-induced fluid flow in these kinds of media and compare the responses with those obtained by replacing the heterogeneous permeability field by constant values, including the average permeability as well as the effective flow permeability of the sample. The latter is estimated in a separate upscaling procedure by solving the steady-state flow equation in the rock sample under study. Numerical experiments let us verify that attenuation levels are less significant and the attenuation peak gets broader in the presence of such strong permeability fluctuations. Moreover, we observe that for very low frequencies the effective seismic permeability is similar to the effective flow permeability, while for very high frequencies it approaches the arithmetic average of the permeability field.

Do seismic waves sense fracture connectivity?

A defining characteristic of fractured rocks is their very high level of seismic attenuation, which so far has been assumed to be mainly due to wave-induced fluid flow (WIFF) between the fractures and the pore space of the embedding matrix. Using oscillatory compressibility simulations based on the quasi-static poroelastic equations, we show that another important, and as of yet undocumented, manifestation of WIFF is at play in the presence of fracture connectivity. This additional energy loss is predominantly due to fluid flow within the connected fractures and is sensitive to their lengths, permeabilities, and intersection angles. Correspondingly, it contains key information on the governing hydraulic properties of fractured rock masses and hence should be accounted for whenever realistic seismic models of such media are needed.

Seismic waves from an underground explosion in a salt bed : final report /

At head of title: "Project Gnome" ; "PNE-150F."

3D Simultaneous Inversion for Seismic Structure and Local Hypocenters in Germany under Consideration of Seismic Anisotropy in the Upper Mantel

The main task of this work has been to investigate the effects of anisotropy onto the propagation of seismic waves along the Upper Mantle below Germany and adjacent areas. Refraction- and reflexion seismic experiments proved the existence of Upper Mantle anisotropy and its influence onto the propagation of Pn-waves. By the 3D tomographic investigations that have been done here for the crust and the upper mantle, considering the influence of anisotropy, a gap for the investigations in Europe has been closed. These investigations have been done with the SSH-Inversionprogram of Prof. Dr. M. Koch, which is able to compute simultaneously the seismic structure and hypocenters. For the investigation, a dataset has been available with recordings between the years 1975 to 2003 with a total of 60249 P- and 54212 S-phase records of 10028 seismic events. At the beginning, a precise analysis of the residuals (RES, the difference between calculated and observed arrivaltime) has been done which confirmed the existence of anisotropy for Pn-phases. The recognized sinusoidal distribution has been compensated by an extension of the SSH-program by an ellipse with a slow and rectangular fast axis with azimuth to correct the Pn-velocities. The azimuth of the fast axis has been fixed by the application of the simultaneous inversion at 25° - 27° with a variation of the velocities at +- 2.5 about an average value at 8 km/s. This new value differs from the old one at 35°, recognized in the initial residual analysis. This depends on the new computed hypocenters together with the structure. The application of the elliptical correction has resulted in a better fit of the vertical layered 1D-Model, compared to the results of preceding seismological experiments and 1D and 2D investigations. The optimal result of the 1D-inversion has been used as initial starting model for the 3D-inversions to compute the three dimensional picture of the seismic structure of the Crust and Upper Mantle. The simultaneous inversion has showed an optimization of the relocalization of the hypocenters and the reconstruction of the seismic structure in comparison to the geology and tectonic, as described by other investigations. The investigations for the seismic structure and the relocalization have been confirmed by several different tests. First, synthetic traveltime data are computed with an anisotropic variation and inverted with and without anisotropic correction. Further, tests with randomly disturbed hypocenters and traveltime data have been proceeded to verify the influence of the initial values onto the relocalization accuracy and onto the seismic structure and to test for a further improvement by the application of the anisotropic correction. Finally, the results of the work have been applied onto the Waldkirch earthquake in 2004 to compare the isotropic and the anisotropic relocalization with the initial optimal one to verify whether there is some improvement.

Cross-hole and seismic CPT tests in a tropical soil site

A large portion of Brazil is covered with tropical soils but literature about dynamic parameters of these soils is very limited. SCPT and cross-hole tests were carried out at an experimental research site inland in the state of São Paulo, Brazil. Shear wave velocities (VS) determined based on both tests are presented and compared. A good agreement was observed between both test results and the differences can be associated with soil variability, which was very sensitive to CPT tests. It was also observed that Go/q c ratio determined based on SCPT appears to be an interesting technique to help identify tropical soils. Copyright ASCE 2006.

Multiresolution spherical wavelet analysis in global seismic tomography

Every seismic event produces seismic waves which travel throughout the Earth. Seismology is the science of interpreting measurements to derive information about the structure of the Earth. Seismic tomography is the most powerful tool for determination of 3D structure of deep Earth's interiors. Tomographic models obtained at the global and regional scales are an underlying tool for determination of geodynamical state of the Earth, showing evident correlation with other geophysical and geological characteristics. The global tomographic images of the Earth can be written as a linear combinations of basis functions from a specifically chosen set, defining the model parameterization. A number of different parameterizations are commonly seen in literature: seismic velocities in the Earth have been expressed, for example, as combinations of spherical harmonics or by means of the simpler characteristic functions of discrete cells. With this work we are interested to focus our attention on this aspect, evaluating a new type of parameterization, performed by means of wavelet functions. It is known from the classical Fourier theory that a signal can be expressed as the sum of a, possibly infinite, series of sines and cosines. This sum is often referred as a Fourier expansion. The big disadvantage of a Fourier expansion is that it has only frequency resolution and no time resolution. The Wavelet Analysis (or Wavelet Transform) is probably the most recent solution to overcome the shortcomings of Fourier analysis. The fundamental idea behind this innovative analysis is to study signal according to scale. Wavelets, in fact, are mathematical functions that cut up data into different frequency components, and then study each component with resolution matched to its scale, so they are especially useful in the analysis of non stationary process that contains multi-scale features, discontinuities and sharp strike. Wavelets are essentially used in two ways when they are applied in geophysical process or signals studies: 1) as a basis for representation or characterization of process; 2) as an integration kernel for analysis to extract information about the process. These two types of applications of wavelets in geophysical field, are object of study of this work. At the beginning we use the wavelets as basis to represent and resolve the Tomographic Inverse Problem. After a briefly introduction to seismic tomography theory, we assess the power of wavelet analysis in the representation of two different type of synthetic models; then we apply it to real data, obtaining surface wave phase velocity maps and evaluating its abilities by means of comparison with an other type of parametrization (i.e., block parametrization). For the second type of wavelet application we analyze the ability of Continuous Wavelet Transform in the spectral analysis, starting again with some synthetic tests to evaluate its sensibility and capability and then apply the same analysis to real data to obtain Local Correlation Maps between different model at same depth or between different profiles of the same model.

Behaviour and applications of elastic waves in structures and metamaterials

The present thesis focuses on elastic waves behaviour in ordinary structures as well as in acousto-elastic metamaterials via numerical and experimental applications. After a brief introduction on the behaviour of elastic guided waves in the framework of non-destructive evaluation (NDE) and structural health monitoring (SHM) and on the study of elastic waves propagation in acousto-elastic metamaterials, dispersion curves for thin-walled beams and arbitrary cross-section waveguides are extracted via Semi-Analytical Finite Element (SAFE) methods. Thus, a novel strategy tackling signal dispersion to locate defects in irregular waveguides is proposed and numerically validated. Finally, a time-reversal and laser-vibrometry based procedure for impact location is numerically and experimentally tested. In the second part, an introduction and a brief review of the basic definitions necessary to describe acousto-elastic metamaterials is provided. A numerical approach to extract dispersion properties in such structures is highlighted. Afterwards, solid-solid and solid-fluid phononic systems are discussed via numerical applications. In particular, band structures and transmission power spectra are predicted for 1P-2D, 2P-2D and 2P-3D phononic systems. In addition, attenuation bands in the ultrasonic as well as in the sonic frequency regimes are experimentally investigated. In the experimental validation, PZTs in a pitch-catch configuration and laser vibrometric measurements are performed on a PVC phononic plate in the ultrasonic frequency range and sound insulation index is computed for a 2P-3D phononic barrier in the sonic frequency range. In both cases the numerical-experimental results comparison confirms the existence of the numerical predicted band-gaps. Finally, the feasibility of an innovative passive isolation strategy based on giant elastic metamaterials is numerically proved to be practical for civil structures. In particular, attenuation of seismic waves is demonstrated via finite elements analyses. Further, a parametric study shows that depending on the soil properties, such an earthquake-proof barrier could lead to significant reduction of the superstructure displacement.

Seismic data from Colle Gnifetti glacie saddle around borehole KCI, Swiss Alps

Two seismic surveys were carried out on the high-altitude glacier saddle, Colle Gnifetti, Monte Rosa, Italy/Switzerland. Explosive and vibroseismic sources were tested to explore the best way to generate seismic waves to deduce shallow and intermediate properties (<100 m) of firn and ice. The explosive source (SISSY) excites strong surface and diving waves, degrading data quality for processing; no englacial reflections besides the noisy bed reflector are visible. However, the strong diving waves are analyzed to derive the density distribution of the firn pack, yielding results similar to a nearby ice core. The vibrator source (ElViS), used in both P- and SH-wave modes, produces detectable laterally coherent reflections within the firn and ice column. We compare these with ice-core and radar data. The SH-wave data are particularly useful in providing detailed, high-resolution information on firn and ice stratigraphy. Our analyses demonstrate the potential of seismic methods to determine physical properties of firn and ice, particularly density and potentially also crystal-orientation fabric.

Seismic modelling and ava analysis of hidrocarbon traps

In this study, forward seismic modelling of four geological models with Hydrocarbon (HC) traps were performed by ray tracing method to produce synthetic seismogram of each model. The idea is to identify the Hydrocarbon Indicators (HCI‟s) such as bright spot, flat spot, dim spot and Bottom Simulating Reflector (BSR) in the synthethic seismogram. The modelling was performed in DISCO/FOCUS 5.0 seismic data processing programme. Strong positive and negative reflection amplitudes and some artifact reflection horizons were observed on produced seismograms due to rapid changes in subsurface velocity and geometry respectively Additionally, Amplitude-versus-angle (AVA) curves of each HCIs was calculated by the Crewes Zoeppritz Explorer programme. AVA curves show that how the reflection coefficients change with the density and the P and S wave velocities of each layer such as oil, gas, gas hydrate or water saturated sediments. Due to AVA curves, an increase in reflection amplitude with incident angle of seismic waves corresponds to an indicator of a hydrocarbon reservoir

Seismic velocity study /

"Issued: May 15, 1963"--Cover ; "December 1962"--Title page.

Estudo do Ruído Sísmico no Arquipélago de Cabo Verde

Traditionally the study of Earth's internal structure has been based on analysis of seismic waves recorded during the occurrence of earthquakes. Recently, some researchers have shown to be possible to use ambient seismic noise (generally generated by ocean storms) to obtain dispersion curves with satisfactory results. Moreover, the use of ambient seismic noise measurement has advantages because they are independent of the earthquake-station geometry. This study aims to determine the dispersion curves of Rayleigh waves from the analysis of ambient seismic noise records, obtained in Cape Verde archipelago, whose origin and structure has attracted interest from several researchers. Analyzing the data collected between November 2007 and September 2008, Green's functions were determined. The analysis of these functions allowed to construct energy diagrams and to measure dispersion curves of group velocities between pairs of stations.

Localization and characterization of an active fault in an urbanized area in central Guatemala by means of geoelectrical imaging

The Polochic and Motagua faults define the active plate boundary between the North American and Caribbean plates in central Guatemala. A splay of the Polochic Fault traverses the rapidly growing city of San Miguel Uspantan that is periodically affected by destructive earthquakes. This fault splay was located using a 2D electrical resistivity tomography (ERT) survey that also characterized the fault damage zone and evaluated the thickness and nature of recent deposits upon which most of the city is built. ERT images show the fault as a similar to 50 m wide, near-vertical low-resistivity anomaly, bounded within a few meters by high resistivity anomalies. Forward modeling reproduces the key aspects of the observed electrical resistivity data with remarkable fidelity thus defining the overall location, geometry, and internal structure of the fault zone as well as the affected lithologies. Our results indicate that the city is constructed on a similar to 20 m thick surficial layer consisting of poorly consolidated, highly porous, water-logged pumice. This soft layer is likely to amplify seismic waves and to liquefy upon moderate to strong ground shaking. The electrical conductivity as well as the major element chemistry of the groundwater provides evidence to suggest that the local aquifer might, at least in part, be fed by water rising along the fault. Therefore, the potential threat posed by this fault splay may not be limited to its seismic activity per se, but could be compounded its potential propensity to enhance seismic site effects by injecting water into the soft surficial sediments. The results of this study provide the basis for a rigorous analysis of seismic hazard and sustainable development of San Miguel Uspantan and illustrate the potential of ERT surveying for paleoseismic studies.

Sediment deformation in Miocene and post-Miocene strata, Northeastern Brazil: Evidence for paleoseismicity in a passive margin

Stratigraphic intervals characterized by varied and complex styles of soft-sediment deformation structures are well preserved in Miocene and Late Pleistocene to Holocene deposits of a sedimentary basin located in Northeastern Brazil. The Miocene strata, represented by the Barreiras Formation, record only brittle structures, including numerous faults and fractures with straight and high angle-dipping planes that are often filled with sands derived from overlying beds. Folds consisting of broad anticlines and synclines are also present in this unit. The late Pleistocene to Holocene deposits, named Post-Barreiras Sediments, contain an indurated sandy package with a large variety of ductile and brittle deformation structures (i.e., massive sandstones with isolated sand fragments and breccias, undulatory strata, sand dykes and diapirs, sinks and bowls, pebbly pockets, plunged sediment mixtures, fitted sand masses, cone-shaped cracks, fault grading and sedimentary enclaves). These features, confined to sharp-based stratigraphic horizons that progressively grade downward into undisturbed deposits, are related to seismic shocks of high surface-wave magnitude (i.e., Ms>5 or 6). Amalgamated seismites suggest that previously formed seismites were affected by subsequent seismic-wave propagation. Seismic waves caused by activity along one, or most likely, several tectonic structures would have propagated throughout the depositional environment, producing laterally extensive seismites. The close proximity to earthquake epicenters would have promoted pervasive re-sedimentation due to pore overpressure, resulting high volumes of massive sandstones and breccia. The similarity between deposits with correlatable strata from many other areas along the Brazilian coast allows raise the hypothesis that the seismic episodes might have affected sedimentation patterns in a large (i.e., extension of several hundreds of kilometers) geographic area. Thus, the modern seismicity recorded along Northeastern Brazil was recurrent during the Quaternary and, perhaps, also in the Pliocene. The estimated high magnitude of the seismic events and the great regional extent of the affected area demonstrate that the Brazilian coast experienced tectonic stress through the last geological episodes of its evolution, which would have favored sediment accumulation and penecontemporaneous re-sedimentation. This geological context is unexpected in a passive margin, inducing to revisit the debate on how active is a passive margin. (C) 2010 Elsevier B.V. All rights reserved.