Multimutation Differential Evolution Algorithm and Its Application to Seismic Inversion

Seismic inversion problems often involve nonlinear relationships between data and model and usually have many local minima. Linearized inversion methods have been widely used to solve such problems. However, these kinds of methods often strongly depend on the initial model and are easily trapped in a local minimum. Global optimization methods, on the other hand, do not require a very good initial model and can approach a global minimum. However, global optimization methods are exhaustive search techniques that can be very time consuming. When the model dimension or the search space becomes large, these methods can be very slow to converge. In this paper, we propose a new global optimization algorithm by incorporating a new multimutation scheme into a differential evolution algorithm. Because mutation operation with the new multimutation scheme can generate better mutant vectors, the new global optimization algorithm has a very good ability of exploring the search space and can converge very fast. We apply the proposed algorithm to both synthetic and field data to test its performance. The results have clearly indicated that the new global optimization algorithm provides faster convergence and yields better results compared with the conventional global optimization methods in seismic inversion.