· · · Institutional Repository

For shallow subsoil, the estimates of in-situ porosity and water saturation are important, but until now it has been difficult to estimate them reliably. We relate seismic and GPR reflection coefficients to porosity and water saturation using a shared earth model. Using this model, we propose a method to integrate quantitatively seismic and GPR angle-dependent reflection coefficients. The new approach has been tested through numerical simulations. The results clearly show that from either seismic or GPR data alone it is impossible to obtain unique estimates for porosity and water saturation; however, a correct integration of those two data types leads to unique and stable estimates at a subsoil layer boundary. Potential applications of this approach exist in hydrogeology and environmental, agricultural and geotechnical engineering.

Small-throw seismogenic fault segments hidden in the Holocene sediments are crucial but difficult targets in seismic exploration. We report here the detection of the deformation pattern and a concealed fault segment in the unconsolidated sediments at Vila Franca Xira, Portugal, through identification in shear wave reflection data of multiple signatures of ductile deformation associated with faulting. We find step-like changes in the stacking velocity along a shallow subsoil layer boundary, indicating synsedimentary faulting. We also recognize a consistent distortion in the moveout of the reflection events in the raw shear wave data. Synthetic modeling of seismic data helps in interpreting these observations and identifying backscattered energy from a steeply dipping shallow fault zone. Prior to this finding, there was no evidence for Holocene activity of this fault, although the fault is considered to be the most probable source for the disastrous 1531 earthquake.

A reliable estimate of the in‐situ permeability of a porous layer in the subsurface is extremely difficult to obtain. We have observed that at the field seismic frequency band the poroelastic behavior for different seismic wavetypes can differ in such a way that their combination gives unique estimates of in‐situ permeability and porosity simultaneously. This is utilized in the integration of angle‐ and frequency‐dependent poroelastic reflection coefficients in a cost function. Realistic numerical simulations show that the estimated values of permeability and porosity are robust against uncertainties in the employed poroelastic mechanism and in the data. Potential applications of this approach exist in hydrocarbon exploration, hydrogeology, and geotechnical engineering.

A careful look into the pertinent models of poroelasticity reveals that in water-saturated sediments or soils, the seismic (P and S wave) velocity dispersion and attenuation in the low field-seismic frequency band (20–200 Hz) have a contrasting behavior in the porosity-permeability domain. Taking advantage of this nearly orthogonal behavior, a new approach has been proposed, which leads to unique estimates of both porosity and permeability simultaneously. Through realistic numerical tests, the effect of maximum frequency content in data and the integration of P and S waves on the accuracy and robustness of the estimates are demonstrated.

A careful look into the pertinent models of poroelasticity reveals that in water-saturated sediments or soils, the seismic (P and S wave) velocity dispersion and attenuation in the low field-seismic frequency band (20-200 Hz) have a contrasting behaviour in the porosity-permeability domain.Taking advantage of this nearly orthogonal behaviour, a new approach has been proposed, which leads to unique estimates of both porosity and permeability simultaneously. Through realistic numerical tests, the effect of maximum frequency content in data and the integration of P and S waves on the accuracy and robustness of the estimates are demonstrated.

Application of seismic interferometry to records from receivers at the Earth’s surface from sources in wells retrieves the reflection response measured at the receivers as if from virtual sources located also at the surface. When the wavefields experience intrinsic losses during propagation, non-physical arrivals (ghosts) would appear in the retrieved result. These ghosts appear due to waves that reflect inside a subsurface layer. Thus, a ghost contains information about the seismic properties of the specific layer. We show how such ghosts can be used to monitor layer-specific changes in the velocity and intrinsic losses in the subsurface. We show how to identify the ghosts using numerical-modelling results from a vertical well, and how to estimate the layer-specific velocity and quality-factor changes using numerical-modelling results from a horizontal well as well as ultrasonic S-wave laboratory data.

Seismic waves converted from compressional to shear mode in the shallow subsurface can be useful not only for obtaining shear-wave velocity information but also for improved processing of deeper reflection data. These waves generated at deep seas have been used successfully in hydrocarbon exploration; however, acquisition of good-quality converted-wave data in shallow marine environments remains challenging. We have looked into this problem through field experiments and synthetic modeling. A high-resolution seismic survey was conducted in a shallow-water canal using different types of seismic sources; data were recorded with a four-component water-bottom cable. Observed events in the field data were validated through modeling studies. Compressional waves converted to shear waves at the water bot-tom and at shallow reflectors were identified. The shear waves showed distinct linear polarization in the horizontal plane and low velocities in the marine sediments. Modeling results indicated the presence of a nongeometric shear-wave arrival excited only when the dominant wavelength exceeded the height of the source with respect to the water/sediment interface, as observed in air-gun data. This type of shear wave has a traveltime that corresponds to the raypath originating not at the source but at the interface directly below the source. Thus, these shear waves, excited by the source/water-bottom coupled system, kinematically behave as if they were generated by an S-wave source placed at the water bottom. In a shallow-water environment, the condition appears to be favorable for exciting such shear waves with nongeometric arrivals. These waves can provide useful information of shear-wave velocity in the sediments.

A digital 3C array seismic cone penetrometer has been developed for multidisciplinary geophysical and geotechnical applications. Seven digital triaxial microelectromechanical system accelerometers are installed at 0.25-m intervals to make a 1.5-m-long downhole seismic array. The accelerometers have a flat response up to 2 kHz. The seismic array is attached to a class 1 digital seismic cone, which measures cone tip resistance, sleeve friction, pore-pressure, and inclination. The downhole 3C array can be used together with impulsive seismic sources and/or high-frequency vibrators that are suitable for high-resolution shallow applications. Results from two field experiments showed that a good data quality, including a constant source function within an array, and a dense depth-sampling allowed robust estimation of seismic velocity profiles in the shallow subsoil. Using horizontal and vertical seismic sources, downhole 9C seismic array data can be easily acquired. The quality of the shear-wave data is much superior when the surface seismic source is a controlled, high-frequency vibrator in stead of a traditional sledge hammer. A remarkable correlation in depth, in a fine scale, between low-strain seismic shear wave velocity and high-strain cone tip resistance could be observed. The array measurements of the full-elastic wavefield and the broad spectral bandwidth are useful in investigating frequency-dependent seismic wave propagation in the porous near-surface soil layers, which is informative of the in situ fluid-flow properties. Stable estimates of dispersive seismic velocity and attenuation can be obtained.

The definition of appropriate places for the development of paleoseismic studies is extremely important in earthquake engineering site investigations. The seismic reflection method is routinely used to locate shallow fault segments where these do not outcrop, like in low slip-rate areas where faults identified in Quaternary sediments have vertical throws less than 2 m. The Lower Tagus Valley (Portugal), covered by 50 m of alluvium sediments, is one of these areas. To find fault segments in this geological environment is a daunting task. Any displacement in the seismic data can be a velocity anomalies and/or statics effect. To illuminate the presence of a fault in Holocene sediments, we acquired an S-wave seismic reflection profile at V. F. de Xira and modelled the response of a fault segment in order to recognise it in our data. The signature of a fault segment can be a change in the shape/attitude of the reflection hyperbolae but reflections from the fault plane are rather weak. These are controlled by fault width and velocity contrast between damaged zone-adjacent sediments. The resemblance between the modelled and field data of V. F. Xira supports the conclusion that the fault affects the Holocene alluvium and is still active.