Document Type : Original Research Paper

Authors

School of Mining Engineering, College of Engineering, University of Tehran, Tehran, Iran.

Abstract

Seismic inversion is a critical technique for estimating the spatial distribution of petro-elastic properties in the subsurface, based on the seismic reflection data. This work introduces an iterative geostatistical seismic inversion method, designed to address challenges in complex geological settings by incorporating self-updating local variogram models. Unlike the conventional approaches that rely on a single global variogram or fixed local variograms, the proposed method dynamically updates the spatial continuity models at each iteration using automatic variogram modeling and clustering of variogram parameters. The optimal number of clusters is determined using three cluster validity indices: Silhouette Index (SI), Davies-Bouldin Index (DB), and Calinski-Harabasz Index (CH). The method’s effectiveness was evaluated using a three-dimensional non-stationary synthetic dataset, demonstrating robust convergence when employing the SI and CH indices, with both achieving a high global correlation coefficient of 0.9 between the predicted and true seismic data. Among these, the CH index provided the best balance between the computational efficiency and inversion accuracy. The results highlight the method’s ability to effectively capture local spatial variability, while maintaining a reasonable computational cost, making it a promising approach for seismic inversion in complex sub-surface environments.

Keywords

Main Subjects

[1] Bosch, M., Mukerji, T., & Gonzalez, E. F. (2010). Seismic inversion for reservoir properties combining statistical rock physics and geostatistics: A review. Geophysics75(5), 75A165-75A176.
[2] Mohammadi, A. K., Mohebian, R., & Moradzadeh, A. (2021). High-resolution seismic impedance inversion using improved ceemd with adaptive noise. Journal of Seismic Exploration30(5), 481-504.
[3] Tarantola, A. (2005). Inverse problem theory and methods for model parameter estimation. Society for industrial and applied mathematics.
[4] Grana, D., Azevedo, L., De Figueiredo, L., Connolly, P., & Mukerji, T. (2022). Probabilistic inversion of seismic data for reservoir petrophysical characterization: Review and examples. Geophysics87(5), M199-M216.
[5] Buland, A., & Omre, H. (2003). Bayesian linearized AVO inversion. Geophysics68(1), 185-198.
[6] Grana, D., & Della Rossa, E. (2010). Probabilistic petrophysical-properties estimation integrating statistical rock physics with seismic inversion. Geophysics75(3), O21-O37.
[7] Grana, D. (2016). Bayesian linearized rock-physics inversion. Geophysics81(6), D625-D641.
[8] Grana, D., Fjeldstad, T., & Omre, H. (2017). Bayesian Gaussian mixture linear inversion for geophysical inverse problems. Mathematical Geosciences49, 493-515.
[9] Azevedo, L., Narciso, J., Nunes, R., & Soares, A. (2021). Geostatistical seismic inversion with self-updating of local probability distributions. Mathematical Geosciences53(5), 1073-1093.
[10] Azevedo, L., & Soares, A. (2017). Geostatistical methods for reservoir geophysics (Vol. 143). Salmon Tower Building, NY: Springer International Publishing.
[11] Bortoli, L. J., Alabert, F., Haas, A., & Journel, A. (1993). Constraining stochastic images to seismic data: Stochastic simulation of synthetic seismograms. In Geostatistics Tróia’92: Volume 1 (pp. 325-337). Dordrecht: Springer Netherlands.
[12] Durrani, M. Z., Talib, M., Ali, A., Sarosh, B., & Rahman, S. A. (2021). Characterization of carbonate reservoir using post-stack global geostatistical acoustic inversion approach: A case study from a mature gas field, onshore Pakistan. Journal of Applied Geophysics188, 104313.
[13] Haas, A., & Dubrule, O. (1994). Geostatistical inversion-a sequential method of stochastic reservoir modelling constrained by seismic data. First break12(11).
[14]. Pereira, Â., Azevedo, L., & Soares, A. (2023). Updating local anisotropies with template matching during geostatistical seismic inversion. Mathematical Geosciences55(4), 497-519.
[15] Soares, A., Diet, J. D., & Guerreiro, L. (2007). Stochastic inversion with a global perturbation method. Petroleum geostatistics. Cascais, Portugal: EAGE.
[16] Nunes, R., Soares, A., Azevedo, L., & Pereira, P. (2017). Geostatistical seismic inversion with direct sequential simulation and co-simulation with multi-local distribution functions. Mathematical Geosciences49, 583-601.
[17] Sabeti, H., Moradzadeh, A., Ardejani, F. D., Azevedo, L., Soares, A., Pereira, P., & Nunes, R. (2017). Geostatistical seismic inversion for non‐stationary patterns using direct sequential simulation and co‐simulation. Geophysical Prospecting65, 25-48.
[18] Sevilla-Villanueva, B., Gibert, K., & Sànchez-Marrè, M. (2016, September). Using CVI for understanding class topology in unsupervised scenarios. In Conference of the Spanish Association for Artificial Intelligence (pp. 135-149). Cham: Springer International Publishing.
[19] Abdalameer, A. K., Alswaitti, M., Alsudani, A. A., & Isa, N. A. M. (2022). A new validity clustering index-based on finding new centroid positions using the mean of clustered data to determine the optimum number of clusters. Expert Systems with Applications191, 116329.
[20] Arbelaitz, O., Gurrutxaga, I., Muguerza, J., Pérez, J. M., & Perona, I. (2013). An extensive comparative study of cluster validity indices. Pattern recognition46(1), 243-256.
[21] José-García, A., & Gómez-Flores, W. (2021). A survey of cluster validity indices for automatic data clustering using differential evolution. In Proceedings of the genetic and evolutionary computation conference (pp. 314-322).
[22] Ikotun, A. M., Habyarimana, F., & Ezugwu, A. E. (2025). Cluster validity indices for automatic clustering: A comprehensive review. Heliyon 11: e41953.
[23] Soares, A. (2001). Direct sequential simulation and cosimulation. Mathematical Geology33, 911-926.