Summary
Bulk-density (ρ) of soil is an important indicator of soil compaction and type. A knowledge of the spatial variability of in-situ soil density is important in geotechnical engineering, hydrology and agriculture. Surface geophysical methods have so far shown limited success in providing an accurate and high-resolution image of 3D soil-density distribution. In this pursuit, 3D seismic full-waveform inversion (FWI) is promising, provided the robustness and accuracy of density inversion via this approach can be established in the near-surface scale. However, simultaneous reconstruction of ρ and seismic wave velocities through multi-parameter FWI remains a challenging task. Near-surface seismic data are commonly dominated by dispersive surface waves whose velocities are controlled by the value and distribution of shear-wave velocity (VS). One major difficulty in estimating reliably ρ from near-surface seismic data is due to the relatively low sensitivity of the seismic wavefield to ρ compared to seismic velocities. Additionally, the accuracy of the estimated ρ decreases due to error in VS – an issue known as parameter coupling. Parameter coupling makes it difficult to estimate accurately ρ within the framework of conventional gradient-based FWI. More sophisticated optimization approaches (e.g., truncated Newton) can reduce the effect of parameter coupling, but these approaches are commonly not affordable in near-surface applications due to heavy computational burden. In this research, we have investigated how choosing correctly the force direction of the seismic source can contribute to a higher accuracy of ρ estimates through 3D FWI. Using scattered wavefields, the Hessian, and inversion tests, an in-depth and systematic investigation of datasets corresponding to different force directions has been carried out. A comparison of the scattered wavefields due to a point-localized ρ perturbation for different force directions shows the robustness of the horizontal-force dataset to noise compared to the vertical-force dataset. Furthermore, for a point-scatterer model, an analysis of the gradients of the misfit function using the Hessian shows that utilizing a horizontal-force source enables one to reconstruct the high-resolution gradient with relatively small parameter coupling. Finally, inversion tests for two different subsoil models demonstrate that 3D FWI on a horizontal-force-source seismic dataset is capable of providing a more accurate 3D ρ distribution in soil compared to a vertical-force-source dataset. Our results show that the use of a horizontal-force source might allow avoiding computationally demanding, costly optimization approaches in 3D FWI.