Amine Dhemaied

Seismic Surface-wave Prospecting Methods for Sinkhole Hazard Assessment along the Dead Sea Shoreline

The Dead Sea (DS) coastal areas have been dramatically hit by sinkhole occurrences since around 1990. It has been shown that the sinkholes along both Israeli and Jordanian shorelines are linked to evaporate karst cavities that are formed by slow salt dissolution. Both the timing and location of sinkholes suggest that: 1) the salt weakens as the result of unsaturated water circulation, thus enhancing the karstification process; and 2) sinkholes appear to be related to the decompaction of the sediments above karstified zones. The location, depth, thickness and weakening of salt layers along the DS shorelines, as well as the thickness and mechanical properties of the upper sedimentary deposits, are thus considered as controlling factors of this on-going process. The knowledge of shear-wave velocities (Vs) should add valuable insights on mechanical properties of both the salt and its overburden. We have suggested Vs estimation using surface-wave prospecting methods, based on surface-wave dispersion measurements and inversion. Two approaches have been used. Along the Israeli shoreline, Vs mapping has been performed to discriminate weak and hard zones within salt layers, after calibration of inverted Vs near boreholes. It has been shown that there is a Vs increase in the DS direction. Initially examined weak zones, located near the salt edge, associated with karstified salt, are characterized by Vs values of 760–1,050 m/s, and extend 60–100 m from the salt edge in the DS direction. Hard salt zones with velocity Vs values greater than 1,500 m/s are located at distances of more than 100–220 m from the salt edge. Finally, transition zones (1,050 , Vs , 1,500 m/s) have a 40–160 m spread. On a Jordanian site, roll-along acquisition and dispersion stacking has been performed to achieve multi-modal dispersion measurements along linear profiles. Inverted pseudo-2-D Vs sections present low Vs anomalies in the vicinity of existing sinkholes and made it possible to detect decompacted sediments associated with potential sinkhole occurrences. Moreover, Vs profiles showed a high velocity unit at 40–50 m depth that can be interpreted as a salt layer.

2D characterization of near-surface V P/V S: surface-wave dispersion inversion versus refraction tomography

The joint study of pressure (P-) and shear (S-) wave velocities (Vp and Vs ), as well as their ratio (Vp /Vs), has been used for many years at large scales but remains marginal in near-surface applications. For these applications, and are generally retrieved with seismic refraction tomography combining P and SH (shear-horizontal) waves, thus requiring two separate acquisitions. Surface-wave prospecting methods are proposed here as an alternative to SH-wave tomography in order to retrieve pseudo-2D Vs sections from typical P-wave shot gathers and assess the applicability of combined P-wave refraction tomography and surface-wave dispersion analysis to estimate Vp/Vs ratio. We carried out a simultaneous P- and surface-wave survey on a well-characterized granite-micaschists contact at Ploemeur hydrological observatory (France), supplemented with an SH-wave acquisition along the same line in order to compare Vs results obtained from SH-wave refraction tomography and surface-wave profiling. Travel-time tomography was performed with P- and SH- wave first arrivals observed along the line to retrieve Vtomo p and Vtomo s models. Windowing and stacking techniques were then used to extract evenly spaced dispersion data from P-wave shot gathers along the line. Successive 1D Monte Carlo inversions of these dispersion data were performed using fixed Vp values extracted from Vtomo p the model and no lateral constraints between two adjacent 1D inversions. The resulting 1D Vsw s models were then assembled to create a pseudo-2D Vsw s section, which appears to be correctly matching the general features observed on the section. If the pseudo-section is characterized by strong velocity incertainties in the deepest layers, it provides a more detailed description of the lateral variations in the shallow layers. Theoretical dispersion curves were also computed along the line with both and models. While the dispersion curves computed from models provide results consistent with the coherent maxima observed on dispersion images, dispersion curves computed from models are generally not fitting the observed propagation modes at low frequency. Surface-wave analysis could therefore improve models both in terms of reliability and ability to describe lateral variations. Finally, we were able to compute / sections from both and models. The two sections present similar features, but the section obtained from shows a higher lateral resolution and is consistent with the features observed on electrical resistivity tomography, thus validating our approach for retrieving Vp/Vs ratio from combined P-wave tomography and surface-wave profiling.

Laser-Doppler acoustic probing of granular media with in-depth property gradient and varying pore pressures

Non-contacting ultrasonic techniques recently proved to be efficient in the physical modeling of seismic-wave propagation at various application scales, as for instance in the context of geological analogue and seismic modeling. An innovative experimental set-up is proposed here to perform laser-Doppler acoustic probing of unconsolidated granular media with varying pore pressures. The preliminary experiments presented here provide reproducible results and exploitable data, thus validating both the proposed medium preparation and pressure gradient generation procedure.

Seismic-Wave Propagation Modeling in Viscoelastic Media Using the Auxiliary Differential Equation Method

In many seismic applications, a constant quality factor is used to describe the constitutive laws of viscoelastic materials, characterized by frequency-independent attenuation characteristics. In such cases, the frequency dependence of the medium’s properties is not taken into account. To overcome this drawback, we proposed an elegant finite difference time domain implementation, with an auxiliary differential equation technique to explicitly solve any stress-strain relation. This scheme is inherited from the formalism of electromagnetism and is based on the separation of the propagation equations from the constitutive law defined in the stress-strain equation. The conventional assumption of a constant quality factor assumption can then be easily avoided in the modeling of seismic-wave propagation in viscoelastic media. We developed such a method and simulated synthetic traces over a simple, 2D viscoelastic homogeneous medium using a Zener model. Wave propagation phase velocities were estimated by means of dispersion analysis and appeared to match theoretical values over a reasonable frequency range. We also measured the material’s attenuation behavior by studying the quality factor, thereby validating our approach.

Contribution of Seismic Methods to Hydrogeophysics

The characterisation and monitoring of aquifer systems mainly rely on piezometric and log data. Delineating spatial variations of lithology between piezometers is a delicate task, which inevitably generates errors possibly propagating into hydrogeological models. Seismic methods have been proposed to: (i) improve the low spatial resolution of borehole data, (ii) provide a characterisation of the subsurface geometry, and (iii) estimate the physical parameters of the medium influenced by the presence of water and the associated flow and transport processes. The joint study of pressure (P-) and shear (S-) wave seismic velocities (VP and VS, respectively), whose evolution is strongly decoupled in the presence of fluid, has been proposed through the estimation of the VP/VS ratio and Poissons ratio. A specific methodology has been developed for the combined exploitation of P- and surface waves present on single seismic records. The use of this methodology in several geological and hydrogeological contexts allowed for estimating VP/VS ratio lateral and temporal variations in good agreement with a priori geological information and existing geophysical and piezometric data. Laser-based ultrasonic techniques were also proposed to put these processing techniques in practice on perfectly controlled physical models and study elastic wave propagation in partially saturated porous media.

Seismic Surface-wave Analysis for Railway Platform Auscultation

The renewal of existent railways requires the characterisation of the mechanical properties of railway platforms (RP), thus raising the need to select appropriate maintenance actions. Conventional techniques (geotechnical soundings, coring) remain local, destructive, expensive and with low yields. Using non-destructive investigation techniques for local diagnosis and monitoring thus appears of great interest for enhancing RP control. Seismic surface-wave methods have been proposed to estimate in situ mechanical parameters of the superficial layers below railways. In this context, a joint geotechnical and seismic survey was carried out along the Northern Europe high-speed line (LGV) in order to precisely determine the origins of a phenomenon affecting the geometry of the track. Strong a priori knowledge of the RP structure allowed for inverting dispersion measurements for 1D VS models along the track. The results showed a contrast of VS in the loess lying below the RP, between areas where the phenomenon was observed and those it was not. This contrast was confirmed by Bender Elements measurements of VS performed on core drilling samples, and corresponded to the lateral variations observed along the track. These results encourage considering dispersion measurements as an appropriate tool of RP monitoring.

Surface-wave Analyses in Unconsolidated Granular Models with Increasing Degrees of Complexity

Using micrometric glass beads, we build small scale physical models with increasing degrees of complexity in order to address theoretical and methodological issues of seismic methods (velocity gradients, lateral variations, pore overpressure, etc.). We simulate seismic records at the surface of the laboratory models thanks to a mechanical source and a laser-Doppler vibrometer. From recorded seismograms, we are able to invert surface-wave dispersion for one or two-dimensional velocity structures. These experiments are for instance used as benchmarks for processing and inversion techniques, enable the validation of numerical methods, or make it possible to study issues related to pore fluids.

Surface-wave Dispersion Stacking on a Granite- micaschists Contact at Ploemeur Hydrological Observatory, France

In the context of a geophysical survey at the Ploemeur hydrological observatory (France), we performed surface-wave profiling for the characterisation of shallow subsurface Shear-wave velocities. Since we anticipated lateral variations but needed great investigation depth, we deployed multifold acquisition geometries and used roll-along dispersion stacking to enable efficient measurements of multi-modal dispersion data. Several offset moving windows have been tested. Represented as pseudo-sections, the phase velocities extracted using a 12-trace window clearly showed three areas coherent with field observation and interestingly consistent with electrical conductivities and P-wave first arrival times. This cross-quality control has been of great help in the choice of the moving window size and revealed itself to be a rewarding step prior to the inversion process.