We have developed a new iterative geostatistical seismic amplitude variation with angle (AVA) inversion algorithm that inverts prestack seismic data, sorted by angle gathers, directly for highresolution density, P-wave velocity, S-wave velocity, and facies models. This novel iterative geostatistical inverse procedure is based on two key main principles: the use of stochastic sequential simulation and cosimulation as the perturbation technique of the model parameter space and a global optimizer based on a crossover genetic algorithm to converge the simulated earth models toward an objective function, in this case, the mismatch between the recorded and synthetic prestack seismic data. As a geostatistical approach, all the elastic models simulated during the iterative pro cedure honors the well-log data at its own locations, the marginal prior distributions of P-wave velocity and S-wave velocity, and density estimated from the available well-log data, and the corresponding joint distributions between density versus P-wave velocity and P-wave versus S-wave velocity. We successfully tested and implemented this new algorithm on a synthetic prestack data set that mimicked the main properties of a real reservoir, and on a real seismic data set acquired over a deepwater turbidite oil reservoir. In both cases, the results showed a good convergence between the recorded and synthetic seismic. The synthetic example showed high-resolution inverted petroelastic models that reproduced the true petroelastic models. The inverted petroelastic models retrieved from the real case study found high resolution and do agree with previous seismic reservoir characterization studies.

No abstract

No abstract

Due to the nature of seismic inversion problems, there are multiple possible solutions that can equally fit the observed seismic data while diverging from the real subsurface model. Consequently, it is important to assess how inverse-impedance models are converging toward the real subsurface model. For this purpose, we evaluated a new methodology to combine the multidimensional scaling (MDS) technique with an iterative geostatistical elastic seismic inversion algorithm. The geostatistical inversion algorithm inverted partial angle stacks directly for acoustic and elastic impedance (AI and EI) models. It was based on a genetic algorithm in which the model perturbation at each iteration was performed recurring to stochastic sequential simulation. To assess the reliability and convergence of the inverted models at each step, the simulated models can be projected in a metric space computed by MDS. This projection allowed distinguishing similar from variable models and assessing the convergence of inverted models toward the real impedance ones. The geostatistical inversion results of a synthetic data set, in which the real AI and EI models are known, were plotted in this metric space along with the known impedance models. We applied the same principle to a real data set using a crossvalidation technique. These examples revealed that the MDS is a valuable tool to evaluate the convergence of the inverse methodology and the impedance model variability among each iteration of the inversion process. Particularly for the geostatistical inversion algorithm we evaluated, it retrieves reliable impedance models while still producing a set of simulated models with considerable variability.

Geostatistical seismic inversion methods are routinely used in reservoir characterisation studies because of their potential to infer the spatial distribution of the petro‐elastic properties of interest (e.g., density, elastic, and acoustic impedance) along with the associated spatial uncertainty. Within the geostatistical seismic inversion framework, the retrieved inverse elastic models are conditioned by a global probability distribution function and a global spatial continuity model as estimated from the available well‐log data for the entire inversion grid. However, the spatial distribution of the real subsurface elastic properties is complex, heterogeneous, and, in many cases, non‐stationary since they directly depend on the subsurface geology, i.e., the spatial distribution of the facies of interest. In these complex geological settings, the application of a single distribution function and a spatial continuity model is not enough to properly model the natural variability of the elastic properties of interest. In this study, we propose a three‐dimensional geostatistical inversion technique that is able to incorporate the reservoir's heterogeneities. This method uses a traditional geostatistical seismic inversion conditioned by local multi‐distribution functions and spatial continuity models under non‐stationary conditions. The procedure of the proposed methodology is based on a zonation criterion along the vertical direction of the reservoir grid. Each zone can be defined by conventional seismic interpretation, with the identification of the main seismic units and significant variations of seismic amplitudes. The proposed method was applied to a highly non‐stationary synthetic seismic dataset with different levels of noise. The results of this work clearly show the advantages of the proposed method against conventional geostatistical seismic inversion procedures. It is important to highlight the impact of this technique in terms of higher convergence between real and inverted reflection seismic data and the more realistic approximation towards the real subsurface geology comparing with traditional techniques.