Sensitivity analysis of optimization parameters on cooperative inversion of seismic reflection and gravity data

Reflections of seismic waves are strongly distorted by the presence of complex geological structures (e.g. salt bodies) and their vertical resolution is usually of the order of a few tens of meters, imposing limitations in the construction of subsurface models. One way to improve the reliability of such models is to integrate reflection seismic data with other types of geophysical data, such as gravimetric data, since the latter provide an additional link to map geological structures that exhibit density contrasts with respect to their surroundings. In a previous study, we developed a cooperative inversion method of 2D post-stack and migrated reflection seismic data, and gravimetric data. Using that inversion method, we minimize two problems: (1) the problem of the distortion of reflection seismic data due to the presence of complex geological bodies and (2) the problem of the greater ambiguity and the commonly lower resolution of the models obtained only from gravimetric anomalies. The method incorporates a technique to decrease the number of variables and is solved by optimization of the gravity inverse problem, thus reducing computing time. The objective function of cooperative inversion was minimized using three different methods of optimization: (1) simplex, (2) simulated annealing, and (3) genetic algorithm. However, these optimization methods have internal parameters which affect the convergence rate and objective function values. These parameters are usually chosen accordingly to previous references. Although the usage of these standard values is widely accepted, the best values to assure effectiveness and stability of convergence are case-dependent. In the present study, we propose a sensitivity analysis on the internal parameters of the optimization methods for the previously presented cooperative inversion. First, we developed the standard case, which is an inversion performed using all parameters at their standard values. Then, the sensitivity analysis is performed by running multiple inversions, each one with a set of parameters. Each set is obtained by modifying the value of a single parameter either for a lower or for a higher value, keeping all other values at their standard values. The results obtained by each setting are compared to the results of the standard case. The compared results are both the number of evaluations and the final value of the objective function. We then classify parameters accordingly to their relative influence on the optimization processes. The sensitivity analysis provides insight into the best practices to deal with object-based cooperative inversion schemes. The technique was tested using a synthetic model calculated from the Benchmark BP 2004, representing an offshore sedimentary basin containing salt bodies and small hydrocarbons reservoirs.