Donatienne Leparoux

Surface-wave inversion limitations from laser-Doppler physical modeling

Surface-wave dispersion inversion is growing in popularity for geotechnical applications, due to its noninvasive character, relative straightforward field procedures and interpretation, especially when the subsurface structure is locally assumed to be one-dimensional (1D). Here, laser-Doppler physical modeling of surface-wave propagation is used to address issues of surface-wave depth penetration, the presence of dipping layers, and the associated limitations and systematic errors propagated in conventional 1D surface-wave inversion. Flat-layered models show that, with an active source and linear spread, the maximum resolvable wavelength of the Rayleigh-wave fundamental mode is on the order of 40% of the spread length. Linearised inversions confirm the rule of thumb that the depth penetration is 20–25% of the spread length, and that correct a priori layer interface depths from refraction analysis allow more accurate results. However, even under optimal conditions, failing to account for a dominant higher mode at low frequency when a stiff shallow layer is present, causes an overestimate of deeper layer shear-wave velocity. Moreover, a layer dip of only a few degrees can significantly bias the surface-wave inversion. If the incorrect a priori information from a single-shot refraction analysis is incorporated in the inverse problem, estimated interface depth depends on the shot position and deeper layer shear-wave velocity is underestimated. Even if correct a priori constraints are used, an underestimate of half-space shear-wave velocity of up to 25% remains.

Small scale adaptation of the seismic full waveform inversion method ‐ Application to civil engineering applications

Full Waveform Inversion (FWI) is a very general multi‐parameters quantitative imaging method originally developed to obtain high resolution images of velocities and attenuation in the natural underground media. FWI promises interesting performances for civil engineering applications like cavity detection. Performances of the FWI method in seismic exploration are difficult to quantify because in real experimentations, the properties of the media are not known. Furthermore, characteristics of the source and coupling of receivers are not controlled. In order to quantify the performances of FWI and its adaptability to near subsurface applications, small scale physicals models are realized and a non‐contact laser ultrasonic laboratory is dedicated to simulate seismic reflection measurement configurations at ultrasonic scale. Seismograms well reproduce real scale data in terms of waveforms and repartition of energy. A FWI algorithm developed by Geosciences Azur and applied in geological investigation is used. T...

Seismic dispersion analysis feasibility for the subgrade investigation: measurement, experimental and numerical modeling

In the construction of roads or railways, the capping layer is the last layer of the earthworks phase. This layer can be made of non bound aggregates or by a treatment of a soil with lime and/or hydraulic binder, such as cement or hydraulic road binders. In the latter case, the in situ testing of the capping layer performances should encompass treated soil modulus measurements, but sampling such materials is not often satisfactory, because of the risk of material degradation by the sampling itself. Consequently, a non destructive method, aiming at measuring the modulus of the treated materials, could be very useful. For this reason, we propose to study the feasibility of the seismic surface or guided waves dispersion analysis in order to recover the depth and the S wave velocity of the subgrade. Previous works provided results and analysis of the dispersion curves concerning the pavement auscultation (Ryden et al, 2004). However, in these cases, the subgrade was an underlying layer in the global zone of interest that includes the upper pavement layers where the measurement surface is the thin asphalt layer. In the present study, we focus on the subgrade layer in the case of under construction roads, before the shallower pavement layers are built because it should help to qualify the project acceptance concerning this earthworks phase. In this context, the issue deals with a two layers medium case where the investigated subgrade, whom the top is the measurement surface, lays above a low velocity zone, i.e. the natural soil. As described by Ryden et al. (2004), the resulting dispersion curves in the case of a high velocity upper layer should be typical of Lamb waves dispersion curves and could bring out higher modes that could be difficult to pick. In this case, they advocated the entire dispersion diagram inversion to avoid any subjective picking in the data (Ryden et al, 2006). However, the treated soil can contain heterogeneities unfavourable to the assumption of homogeneous layers presupposed to the dispersion diagram calculation. Thus a first feasibility stage, i.e. a field experimental data acquisition, was conducted to define the ability of seismic data to provide a coherent dispersion diagram in the spectral content required. The dispersion curve have been extracted and inverted with an iterative weighted least squares local minimization method (Hermann, 2002). In order to consider the possibility of inverting the entire dispersion diagram, a second feasibility stage consisted in analysing all the events that possibly occur in the dispersion diagram. For that, the measurement experience was reproduced at reduced scale in laboratory as a perfectly controlled experimental modelling approach. These data and more precisely the dispersion diagram is compared in one hand to the theoretical curves associated to the leakage attenuations and in an other hand to the theoretical dispersion diagram numerically calculated with an original method taking into account the source effects.