Joint inversion of seismic and electrical data in saturated porous media

Joint inversion strategies and physical constraints on model parameters may be used to mitigate equivalence problems caused by solution non-uniqueness. This strategy is quite a common practice in exploration geophysics, where dedicated rock physical studies are usually carried out, while it is not so frequent in near surface geophysics. We use porosity as a constraint among seismic wave velocities and electrical resistivity in a deterministic joint inversion algorithm for surface wave dispersion, P-wave traveltimes and apparent resistivity from vertical electrical sounding. These data are often available for near surface characterization. We show that the physical constraint among model parameters leads to internally consistent geophysical models in which solution non-uniqueness is mitigated. Moreover, an estimate of soil porosity is obtained as a relevant side product of the procedure. In particular, we consider a clean sand deposit and hence the appropriate formulations for the computation of porosity from seismic velocities and resistivity are implemented in the algorithm. We first demonstrate how the non-uniqueness of the solution is reduced in a synthetic case and then we applied the algorithm to a real-case study. The algorithm is here developed for one-dimensional condition and for granular soils to better investigate the physical constraint only, but it can be extended to the two-dimensional or three-dimensional case as well as to other materials with the adoption of proper rock physical relationships

The Polito Surface Wave flat-file Database (PSWD): statistical properties of test results and some inter-method comparisons

The compilation and maintenance of experimental databases are of crucial importance in all research fields, allowing for researchers to develop and test new methodologies. In this work, we present a flat-file database of experimental dispersion curves and shear wave velocity profiles, mainly from active surface wave testing, but including also data from passive surface wave testing and invasive methods. The Polito Surface Wave flat-file Database (PSWD) is a gathering of experimental measurements collected within the past 25 years at different Italian sites. Discussion on the database content is reported in this paper to evaluate some statistical properties of surface wave test results. Comparisons with other methods for shear wave velocity measurements are also considered. The main novelty of this work is the homogeneity of the PSWD in terms of processing and interpretation methods. A common processing strategy and a new inversion approach were applied to all the data in the PSWD to guarantee consistency. The PSWD can be useful for further correlation studies and is made available as a reference benchmark for the validation and verification of novel interpretation procedures by other researchers

Illumination diagnosis for retrieval of reflections from ambient-noise seismic data in the Siilinjärvi mining site, Finland

Reflection seismic methods are becoming popular in mineral exploration, because they allow high-resolution delineation of the exploration targets, even at great depths. Seismic interferometry can be used to retrieve reflections from passive seismic data, removing the need for active seismic sources and, therefore, reducing the cost and environmental impact of exploration. The retrieval of reflections can be challenging, since passive seismic records are typically dominated by surface waves. Therefore, illumination diagnosis, a method which allows the isolation of the portions of the passive data where body-wave signals are stronger, can be a valuable step that improves the quality of the reflections retrieved from seismic interferometry and reduces the overall computational cost of the processing stage. Here, we validate the performance of the method to effectively isolate the portions of the passive data dominated by body waves and apply it on an ambient-noise seismic dataset acquired in the Siilinjärvi mining site in Finland

Challenges in shallow target reconstruction by 3D elastic full-waveform inversion - Which initial model?

Elastic full-waveform inversion (FWI) is a powerful tool for high-resolution subsurface multiparameter characterization. However, 3D FWI applied to land data for near-surface applications is particularly challenging because the seismograms are dominated by highly energetic, dispersive, and complex-scattered surface waves (SWs). In these conditions, a successful deterministic FWI scheme requires an accurate initial model. Our study, primarily focused on field data analysis for 3D applications, aims at enhancing the resolution in the imaging of complex shallow targets, by integrating devoted SW analysis techniques with a 3D spectral-element-based elastic FWI. From dispersion curves, extracted from seismic data recorded over a sharp-interface shallow target, we build different initial S-wave (VS) and P-wave (VP) velocity models (laterally homogeneous and laterally variable), using a specific data transform. Starting from these models, we carry out 3D FWI tests on synthetic and field data, using a relatively straightforward inversion scheme. The field data processing before FWI consists of band-pass filtering and muting of noisy traces. During FWI, a weighting function is applied to the far-offset traces. We test 2D and 3D acquisition layouts, with different positions of the sources and variable offsets. The 3D FWI workflow enriches the overall content of the initial models, allowing a reliable reconstruction of the shallow target, especially when using laterally variable initial models. Moreover, a 3D acquisition layout guarantees a better reconstruction of the target's shape and lateral extension. In addition, the integration of model-oriented (preliminary monoparametric FWI) and data-oriented (time windowing) strategies into the main optimization scheme has produced further improvement of the FWI results

Application of Surface wave methods for seismic site characterization

Surface-wave dispersion analysis is widely used in geophysics to infer a shear wave velocity model of the subsoil for a wide variety of applications. A shear-wave velocity model is obtained from the solution of an inverse problem based on the surface wave dispersive propagation in vertically heterogeneous media. The analysis can be based either on active source measurements or on seismic noise recordings. This paper discusses the most typical choices for collection and interpretation of experimental data, providing a state of the art on the different steps involved in surface wave surveys. In particular, the different strategies for processing experimental data and to solve the inverse problem are presented, along with their advantages and disadvantages. Also, some issues related to the characteristics of passive surface wave data and their use in H/V spectral ratio technique are discussed as additional information to be used independently or in conjunction with dispersion analysis. Finally, some recommendations for the use of surface wave methods are presented, while also outlining future trends in the research of this topic