Joint inversion of Rayleigh wave fundamental and higher order mode phase velocity dispersion curves using multi‐objective grey wolf optimization

Abstract

Rayleigh wave dispersion curves can be inverted to retrieve subsurface seismic velocity profiles. The inverse problem is ill-posed, nonlinear and poorly conditioned, necessitating the application of global optimization methods. We present the application of the multi-objective grey wolf optimization algorithm to perform joint inversion of the phase velocity dispersion curves corresponding to the fundamental and higher order modes of Rayleigh waves to obtain shear (S-) and primary (P-) wave velocity profiles. Multi-objective grey wolf optimization is an extension of the grey wolf optimization algorithm for application to multi-objective optimization and can be adapted to solve joint inversion problems. We compare the joint inversion results obtained from the multi-objective grey wolf optimizer with those obtained from Markov chain Monte Carlo and fundamental mode inversion using the grey wolf optimizer on synthetic examples. The errors associated with phase velocity measurements are simulated by adding frequency-dependent noise, with a higher level of noise added to the phase velocities corresponding to lower frequencies as compared to the higher frequencies. In the multimode joint inversion problem, the multi-objective grey wolf optimizer gives a suite of solutions corresponding to each model parameter. The suite of S-wave and P-wave velocity profiles estimated from the multi-objective grey wolf optimizer matches closely with the true model for the synthetic case studies even in the presence of noise. However, the suite of solutions has a greater spread for the last few layers, qualitatively indicating a higher degree of uncertainty in the predicted model parameter. The uncertainty in the solution for the deeper layers is a consequence of the uncertainty in the phase velocity at lower frequencies. We demonstrate the efficacy of the algorithm on recorded data from a shallow seismic survey conducted at the Indian Institute of Technology Bombay. The results from the multi-objective grey wolf optimizer are in close agreement with those from Markov chain Monte Carlo, and the depth of investigation is found to be greater in comparison to results from refraction traveltime inversion.

Structural information derived from ambient noise tomography over a hydrocarbon‐producing region in the Cachar fold belt, lower Assam, northeast India

“…The dispersion curves can be inverted for shear velocities with either local (Cercato, 2009;Haney & Tsai, 2017;Xia et al, 1999) or global (Lehujeur et al, 2018;Pei et al, 2007) starting model to converge to the global solution. Global optimization techniques, for instance, simulated annealing (Pei et al, 2007), genetic algorithms (Moro et al, 2007), those based on particle swarm optimization (Vashisht et al, 2022) or Markov chain Monte Carlo methods (Lehujeur et al, 2018), explore a broad model space and are more likely to find the global solution. We invert the group velocity dispersion curves using the Metropolis algorithm implemented by Lehujeur et al (2018).…”

A Novel Combined Regularization for Identifying Random Load Sources on Coal‐Rock Structure

“…As a kind of inverse problem in system dynamics, load identifcation has become a hot research topic and has been widely used in coal mining machinery, aerospace, signal processing, and other felds [1,2]. Te main concern of many experts and scholars on the problem of load identifcation stems from the development of regularization theory for solving ubiquitous and ill-posed problems in practical engineering applications in recent years [3,4].…”