Optimized double kriging algorithm for the geological environment parameters modeling

Abstract

Параметры геологической среды имеют широкое применение при разведке и исследовании углеводородных месторождений. Они хранят в себе информацию о физических характеристиках горных пород, местоположении коллекторов и границах разделения пластов. Численные значения рассматриваемых параметров получают путем проведения геофизических исследований скважин (ГИС) либо путем пересчета уже имеющихся каротажных кривых на основании известных петрофизических закономерностей. Поскольку ГИС могут быть проведены лишь при условии наличия пробуренной скважины, параметры среды зачастую известны лишь в небольшом околоскважинном пространстве. Это делает актуальной задачу прогноза параметров среды в пределах всего месторождения. Существующие методы решения данной задачи условно можно разделить на две группы: использующие только данные ГИС и комплексно использующие данные 3Dсейсморазведки и ГИС. Первая группа эффективна при наличии густой сетки скважинных измерений. Вторая группа показывает лучшие результаты в условиях редкой нерегулярной сетки, однако требует наличия сейсмических данных в пределах всего исследуемого месторождения. В работе приводятся описание нового метода прогноза параметров геологической среды, базирующегося на совместном использовании данных, ГИС и 3D сейсморазведки, алгоритмическая реализация нового метода и рассматриваются способы оптимизации алгоритмической реализации. The geological environment parameters are widely used in hydrocarbon deposits exploration and study. They store information about the rocks' physical characteristics, the location of the reservoirs, the boundaries of the separation of the layers, etc. The numerical values of the parameters under consideration are obtained by conducting geophysical well logging (GWL) or recalculating existing logging curves based on known petrophysical laws. Because GWL can be carried out only if a drilled well exists, the environment parameters are often known only in a small near-wellbore space. That makes predicting the geological environment parameters within the entire field relevant. The existing methods for solving this task can be conditionally divided into two groups: using only GWL data; complexly using 3D seismic and GWL data. The first group is impactful when a dense grid of downhole measurements presents. The second group shows the best results in a sparse irregular grid. However, it requires seismic data availability within the entire studied field. This paper describes a new method for predicting geological environment parameters based on the GWL data and 3D seismic prospecting combined use. The study also provides an algorithmic implementation of the new procedure. Methods for optimizing algorithmic implementation are considered.