Способ определения положения координат глубин маркера при построении геологической модели месторождения

Изобретение относится к способу, устройству и машиночитаемому носителю данных, предназначенным для построения геологической модели нефтяного или иного месторождения, в частности, для определения коэффициентов корреляции для комплекса кривых ГИС и нахождения положений глубин маркера, для которых значение коэффициента корреляции является максимальным. Техническим результатом является повышение точности вычислений параметров, используемых при построении геологической модели расположения нефтяных или иных месторождений. Метод позволяет для маркера, уже имеющего отметки на некоторой, называемой опорной, группе скважин, вычислить их для скважин из другой группы. Для каждой скважины W, на которой ищется значение глубины маркера, выбираются скважины опорной группы, отстоящие от скважины W на заданном расстоянии, и среди них выбирается скважина с наибольшим значением коэффициента корреляции, при этом точка, в которой этот максимум достигается, назначается искомой отметки маркера. С помощью проверяющих тестов осуществляют поиск скважин, в которых функция корреляции меньше, чем максимальное значение коэффициента корреляции, а коэффициент качества корреляции больше, чем максимальное значение коэффициента корреляции. После чего добавляют найденную скважину к опорной группе скважин. 3 н. и 5 з.п. ф-лы, 9 ил.

 

Изобретение относится к способу, устройству и машиночитаемому носителю данных, предназначенным для построения геологической модели нефтяного или иного месторождения. В частности, изобретение относится к способу, устройству и машиночитаемому носителю данных, предназначенным для определения коэффициентов корреляции для комплекса кривых ГИС и нахождения положений глубин маркера, для которых значение коэффициента корреляции является максимальным.

Уровень техники

Из уровня техники известен источник ЕА 200600036 A1, E21B 7/04, 30.12.2008, описывающий пакет программ, хранимых в памяти на автоматизированном рабочем месте или в другой вычислительной компьютерной системе, предназначенных для построения модели одиночной скважины.

Известен источник US 2007/0276604 A1, G01V 1/00, 29.11.2007, который раскрывает способ визуализации и организации данных месторождений нефти и газа. Способ раскрывает обработку каротажных диаграмм скважин, которая заключается в использовании растровых изображений каротажных диаграмм, которые переводятся в цифровой формат для последующего нанесения на них маркера.

Источник US 2010/0004864 А1, G01V 9/00, 07.01.2010 раскрывает способ корреляции каротажных диаграмм, который заключается в том, что автоматически коррелируют данные множества каротажных диаграмм, описывающих информацию на различных участках местности.

Все из вышеуказанных источников в той или иной степени раскрывают общие принципы моделирования геологической модели месторождения полезных ископаемых, однако ни один из вышеприведенных источников не обеспечивает и не подразумевает генерирование высокоточной геологической модели на основании коэффициентов корреляции между большими группами скважин относительно положений глубинных маркеров, что в свою очередь позволяет повысить точность определения расположения месторождений полезных ископаемых.

Раскрытие изобретения

Задачей, на которое направлено заявленное изобретение, является построение геологической модели с точным определением положений нефтяных или иных месторождений.

Техническим результатом является повышение точности вычислений параметров, используемых при построении геологической модели расположения нефтяных или иных месторождений.

Задача корреляции комплексов кривых ГИС предполагает наличие группы скважин, для которых проведены исследования ГИС методами, достаточными для стратиграфического и литологического анализа. Без ограничения общности считается, что в результате этих исследований для каждой скважины группы получен набор кривых ГИС, по одной для каждого метода. Для построения геологической модели залежи важным этапом является прослеживание границ стратиграфических комплексов или литологических свойств. Эти границы могут быть выявлены вдоль скважин и продолжены на исследуемую площадь путем интерполяции. Такие границы, называемые маркерами, имеют на скважинах глубинные отметки, которые располагаются в местах совокупных особенностей кривых ГИС.

В одном варианте осуществления изобретения способ включает в себя определение коэффициентов корреляции для комплекса кривых ГИС на парах скважин, расположенных не дальше заданного расстояния друг от друга, и нахождение положений глубин маркера, для которых этот коэффициент является максимальным. Глубины маркера ищутся для скважин, не входящих в опорную группу, для которой эти глубины предварительно заданы. Метод состоит также в многократном повторении описанного поиска с включением вновь найденных скважин в опорную группу на каждой итерации. Поиск фильтруется набором тестов, повышающих точность вычислений. В способе могут быть задействованы трендовые маркеры, использование которых увеличивает его эффективность.

Метод позволяет для маркера, уже имеющего отметки на некоторой, называемой опорной, группе скважин, вычислить их для скважин из другой группы. Для каждой скважины W, на которой ищется значение глубины маркера, выбираются скважины опорной группы, отстоящие от скважины W на заданном расстоянии, и среди них выбирается скважина с наибольшим значением коэффициента корреляции. Тогда точка на скважине W, в которой этот максимум достигается, назначается в качестве искомой отметки маркера. Полагая далее результат работы описанного алгоритма итерацией обобщенного алгоритма, и, добавляя в опорную группу после каждой итерации, скважины его решения, процесс повторяют до тех пор, пока алгоритм выполнения способа не даст пустой результат, т.е. не сообщит о том, что все скважины, для которых были обнаружены отметки маркера, проанализированы и выявлены в заданной области. Такое “зацикливание” основного алгоритма позволяет получить решение для большого количества скважин. Для повышения надежности расчета в метод добавлен ряд тестов, не допускающих установки отметок маркера в “плохих” местах. Эти тесты следующие: порог корреляции, качество корреляции, уровень транзитивности, ограничения по значениям функций на скважине.

Во втором варианте осуществления изобретения способ включает в себя вычисление функционалов, представляющих собой сумму коэффициентов корреляции для комплекса кривых ГИС в точках маркеров, являющихся начальным приближением, на парах скважин, расположенных не дальше заданного расстояния друг от друга, и нахождение их максимумов с помощью адаптированного к задаче алгоритма градиентного спуска.

Такой способ позволяет для маркеров, выбранных в качестве начального решения, вычислить такие глубины маркера на каждой скважине, которые обеспечивают наилучшую суммарную корреляцию. Для каждого маркера, входящего в набор, определяется функционал, представляющий собой сумму коэффициентов корреляции комплекса методов ГИС для пар скважин, расположенных не далее заданного расстояния друг от друга. Для этого функционала вычисляются частные производные и полученный таким образом вектор сглаживается и используется для нахождения большего значения функционала на некотором отрезке вдоль этого вектора. Если большего значения не найдено, то последнее положение отметок маркера считается решением задачи, а если найдено, то производится сглаживание точки решения, и процесс повторяется снова. На каждой итерации алгоритма производится сортировка глубин маркеров.

В третьем варианте настоящее изобретение обеспечивает устройство для определения значения отметки глубины маркера на скважинах W при построении геологической модели месторождения. Устройство может представлять собой, но не ограничиваться: суперкомпьютер, персональный компьютер, портативный компьютер, планшетный компьютер, карманный компьютер, смартфон и тому подобное. Устройство обязательно содержит один или более процессоров, которые предназначены для выполнения компьютерных команд или кодов, хранящихся в памяти устройства с целью обеспечения выполнения способов по первому или второму вариантам осуществления настоящего изобретения, машиночитаемый носитель данных (память) и модули ввода/вывода (I/O). Модули I/O представляют собой, но не ограничиваются, типичные и известные из уровня техники средства управления устройством: манипулятор типа «мышь», клавиатура, джойстик, тачпад, трекбол, электронное перо, стилус, сенсорный дисплей и тому подобное. Также модули I/O представляют собой, но не ограничиваются, типичные и известные из уровня техники средства демонстрирования информации: монитор, проектор, принтер, графопостроитель и тому подобное. В качестве примера, а не ограничения, машиночитаемый носитель данных может включать в себя оперативную память (RAM); постоянное запоминающее устройство (ROM); электрически-стираемое программируемое постоянное запоминающее устройство (EEPROM); флэш-память или другие технологии памяти; CDROM, цифровой универсальный диск (DVD) или другие оптические или голографические носители данных; магнитные кассеты, магнитную пленку, запоминающее устройство на магнитных дисках или другие магнитные запоминающие устройства, несущие волны или другой носитель данных, который может быть использован для кодирования требуемой информации и к которому может быть осуществлен доступ посредством устройства, описанного выше.

В четвертом варианте настоящее изобретение обеспечивает машиночитаемый носитель данных, содержащий код программы, который побуждает процессор и/или процессоры выполнять действия по способу, описанному в первом или втором вариантах осуществления изобретения и которые, соответственно, дополнительно не описываются. В качестве примера, а не ограничения, машиночитаемый носитель данных может включать в себя оперативную память (RAM); постоянное запоминающее устройство (ROM); электрически стираемое программируемое постоянное запоминающее устройство (EEPROM); флэш-память или другие технологии памяти; CDROM, цифровой универсальный диск (DVD) или другие оптические или голографические носители данных; магнитные кассеты, магнитную пленку, запоминающее устройство на магнитных дисках или другие магнитные запоминающие устройства, несущие волны или другой носитель данных, который может быть использован для кодирования требуемой информации, и к нему может быть осуществлен доступ посредством устройства, описанного в третьем варианте осуществления изобретения, и которое соответственно дополнительно не описывается.

Краткое описание чертежей

Фиг.1 иллюстрирует общую схему выполнения способа по первому варианту осуществления изобретения.

Фиг.2 иллюстрирует кривые метода ГИС на скважинах W и W₀ соответственно и функцию корреляции метода ГИС.

Фиг.3 иллюстрирует этап определения расстояния между группами скважин.

Фиг.4 иллюстрирует этап определения глубины z_max, соответствующей максимальному значению корреляции.

Фиг.5 иллюстрирует этап добавления к опорной группе скважин, глубина которых была определена на этапе определения глубины z_max.

Фиг.6 иллюстрирует общую схему выполнения способа по второму варианту осуществления изобретения.

Фиг.7 иллюстрирует этап составления вектора градиента для текущей точки решения.

Фиг.8 иллюстрирует этап определения в направлении вектора градиента значения функционала большего предыдущего.

Фиг.9 иллюстрирует этап сглаживания найденного решения.

Детальное описание чертежей

Метод решения задачи, описываемый далее в отношении первого варианта осуществления изобретения, состоит в нахождении глубин маркера путем вычисления коэффициента корреляции комплекса кривых ГИС на парах скважин, для одной из которых эта задача уже решена. А именно, пусть для некоторой подгруппы скважин, которую мы называем опорной, глубины маркера уже заданы. Тогда выбрав одну из опорных скважин W₀ с глубиной маркера z₀, для любой скважины W, не принадлежащей опорной группе, мы вычисляем функцию корреляции C(z), значениями которой являются коэффициенты корреляции комплекса кривых ГИС на скважине W в точке z, и на скважине W₀ в точке z₀. Обозначим z_max максимум этой функции. Если для скважины W найдены скважины опорной группы, расположенные в заданной ее окрестности, то глубина z_max, соответствующая опорной скважине с наибольшим значением максимума функции C(z), считается возможной глубиной маркера на скважине W. Функция корреляции принимает значения от -1 до 1 и для одного метода ГИС F определяется следующим образом:

$$C(z)=\frac{{\displaystyle \int F_W(T(z+x))F_{W_0}(z_0+x)dx}}{\sqrt{{\displaystyle \int F_W^2(T(z+x))dx}}\sqrt{{\displaystyle \int F_{W_0}^2(z_0+x)dx}}}$$ ,

где F_W и $$F_{W_0}$$ - кривые метода F на скважинах W и W₀ соответственно. Функция T имеет вид T(z)=az+b, причем а отлично от 1, a b от 0 только в том случае, когда заданы один или два трендовых маркера, соответственно. В этих случаях а и b вычисляются как решение системы очевидных линейных уравнений. Для комплекса методов $$\left\{F_i\right\}$$ ,i=0,…,n, функция корреляции определяется следующим образом:

$$C(z)=\frac{{\displaystyle \sum _{i=0}^n{w}_i{C}_i(z)}}{{\displaystyle \sum _{i=0}^n{w}_i}}$$ ,

где C_i(z) - функция корреляции для метода F_i, a w_i - весовые коэффициенты (Фиг.2).

В соответствии с фиг.1 далее будут подробно изложены этапы выполнения способа.

На первом этапе метода находятся скважины, расположенные в заданной окрестности опорной группы скважин, то есть расположенные не далее заданного расстояния R от какой-нибудь скважины опорной группы. Расстояние между скважинами может быть определено несколькими способами: как расстояние между точками скважины на фиксированной глубине, как расстояние между точками заданного маркера или как расстояние между точками пересечения скважин с заданной поверхностью (Фиг.3).

Второй этап метода состоит в вычислении функций корреляции и выборе для каждой, найденной на первом этапе, скважины W, соответствующей ей скважины из опорной группы W₀ с максимальным коэффициентом корреляции и соответствующей ей точке z_max на скважине W.

На третьем этапе каждая выбранная на втором этапе скважина W подвергается проверке на ряд тестов и, если эта проверка оказывается успешной, выбор скважины считается утвержденным, а глубина z_max, соответствующая максимальному значению корреляции, назначается глубиной маркера на скважине W (Фиг.4).

На четвертом этапе метода скважины с утвержденным на третьем этапе выбором добавляются к опорной группе и происходит переход снова к первому этапу (Фиг.5).

Этот процесс повторяется до тех пор, пока на первом, втором либо третьем этапе не будет найдено ни одной скважины. Тестами, входящими в проверку на третьем этапе, являются:

1) превышение функцией корреляции в точке z_max заданного порогового значения. Большие пороговые значения корреляции обеспечивают большую точность алгоритма, но уменьшают число найденных скважин. Например, если задано пороговое значение 0.9, а максимум функции корреляции в точке z_max равен 0.88, то скважина считается не прошедшей тест.

2) превышение заданного порогового значения коэффициентом качества корреляции в точке. Коэффициент качества корреляции определяется как коэффициент отклонения ближайшего по значению локального максимума функции корреляции от ее наибольшего локального максимума в точке z_max. Большие значения этого коэффициента означают, что найденный максимум коэффициента корреляции существенно больше других локальных максимумов функции корреляции. Например, если максимальное значение функции корреляции равно 0.9 и ближайший по величине максимум равен 0.89, то значение коэффициента качества корреляции будет равно (0.9-0.89)/0.9, что приблизительно равно 0.01. Тогда при пороговом значении коэффициента качества, равном, например, 0.5, скважина будет считаться не прошедшей тест.

3) превышение заданного порогового значения степенью транзитивности. Степень транзитивности определяется как количество предыдущих итераций алгоритма, для которых максимум функции корреляции найденной опорной скважины со скважиной текущей итерации проходит тест 1. Этот тест повышает степень надежности метода. Например, если задан порог для значений функции корреляции из теста 1, равный 0.9, и скважина А прошла этот тест при корреляции со скважиной В из опорной группы, а скважина В в свою очередь была утверждена при корреляции со скважиной С на предыдущей итерации, то значение порога степени транзитивности, равное 2, требует получения максимума функции корреляции не ниже 0.9 для скважины А и С как опорной для А. Пороговое значение степени транзитивности полагается равным числу предыдущих итераций алгоритма, если это число меньше, чем заданное пороговое значение.

4) принадлежность значений заданной функции на скважине в точке максимума корреляции заданному интервалу значений. Этот тест позволяет отбрасывать заведомо плохие места на основании значений какой-либо функции, содержащей эту информацию. В качестве такой функции метод может использовать, например, функцию когерентности или функцию отклонений экстремумов вейвлет-преобразования с увеличением периода.

Метод решения задачи, описываемый далее в отношении второго варианта осуществления изобретения, состоит в нахождении глубин маркеров путем вычисления максимумов функционалов, характеризующих степень похожести комплексов кривых ГИС в точках маркеров на скважинах. А именно, пусть $$\left\{z_i\right\}$$ ,i=0,…,n глубинные отметки на скважинах W_i. Определим функционал на n-мерном пространстве следующим образом:

$$C(z_0,\dots z_n)={\displaystyle \sum _{i=0}^n{\displaystyle \sum _{j=i+1}^{n}B(i,j)\cdot C(z_i,z_j)}}$$ ,

Здесь через C_k(z_i, z_j) обозначен коэффициент корреляции для k-того метода ГИС в точках z_i и z_j на скважинах i и j, соответственно. Этот коэффициент принимает значения от -1 до 1 и вычисляется по следующей формуле:

$$C_k(z_i,{\text{ z}}_{\text{j}})=\frac{{\displaystyle \int F_{k,i}(z_i+x)F_{k,j}(z_j+x)dx}}{\sqrt{{\displaystyle \int F_{k,i}^2(z_i+x)dx}}\sqrt{{\displaystyle \int F_{k,j}^2(z_j+x)dx}}}$$

Максимумы этих функционалов находятся с использованием метода градиентного спуска (в данном случае, подъема), который состоит в вычислении вектора градиента, координатами которого являются частные производные функционала, и последующем поиске наибольших значений вдоль этого вектора. Для нейтрализации типичных проблем применения метода градиентного спуска реализуется специальный прием "встряхивания" промежуточных решений алгоритма, сводящийся к сглаживанию текущей точки решения и текущего вектора градиента. Такое сглаживание выполняется с заданным, уменьшающимся на каждой итерации, коэффициентом и на последних итерациях алгоритма практически не происходит. Сглаживание производится, например, методом скользящего среднего в окне, причем коэффициентом сглаживания в этом случае является радиус окна. Поскольку при изменении текущей точки решения ее глубины на некоторых скважинах могут сильно отклониться от глубин на других скважинах и оказаться ближе по глубине к другим маркерам набора, то производится сортировка глубин маркеров набора на каждой скважине. В результате этой сортировки глубины, соответствующие одному маркеру набора, могут быть присвоены другому. Этот простой прием позволяет уменьшить разброс глубин маркеров на разных скважинах.

На первом этапе вычисляется значение функционала С в начальной точке решения $$\left\{z_i\right\}$$ .

На втором этапе метода вычисляются частные производные функционала и составляется вектор градиента для текущей точки решения $$\left\{z_i\right\}$$ (Фиг.7).

На третьем этапе производится сглаживание вектора градиента, то есть каждая компонента вектора на скважине W заменяется на среднее значение компонент на скважинах, попавших в окрестность скважины W заданного радиуса R.

На четвертом этапе для отрезка с начальной точкой $$\left\{z_i\right\}$$ и заданной длины в направлении вектора градиента с заданным шагом ищется значение функционала, большее предыдущего. Если такое значение не найдено, то алгоритм прекращает работу и текущее решение $$\left\{z_i\right\}$$ принимается за окончательное (Фиг.8).

На пятом этапе найденное решение уточняется за счет поиска большего значения в пределах шага, для которого обнаружен максимум, с более мелким шагом.

На шестом этапе текущим решением $$\left\{z_i\right\}$$ полагается точка, в которой достигнут максимум функционала и производится сглаживание этой точки тем же способом, что вектора градиента на третьем этапе, а радиус сглаживания R уменьшается на заданную величину (Фиг.9).

На седьмом этапе выполняется сортировка отметок маркеров набора по глубине и производится переход ко второму этапу алгоритма.

1. Способ определения положения координаты глубины маркера скважины W при построении геологической модели месторождения полезных ископаемых, содержащий этапы на которых:
1) определяют опорную группу скважин;
2) определяют скважины W, расположенные в заданной окрестности опорной группы скважин, при этом окрестность не превышает значение расстояния R от любой из скважин опорной группы;
3) определяют функцию корреляции C(z) и выбирают для каждой найденной на предыдущем этапе скважины W соответствующую ей скважину из опорной группы с максимальным значением коэффициента корреляции и соответствующую ей точку Z_max на скважине W;
4) проверяют скважину W с помощью проверяющих тестов и, если эта проверка оказывается успешной, то выбор скважины считается утвержденным, а глубина, соответствующая максимальному значению корреляции Z_max, назначается глубиной маркера на скважине W, причем по меньшей мере одним проверяющим тестом является проверяющий тест, на котором осуществляют проверку на предмет превышения максимумом функции корреляции C(z) в точке Z_max заданного порогового значения, причем значение максимума функции корреляции во время выполнения проверяющего теста округляют до десятых долей в сторону большего значения, причем если значение максимума корреляции при округлении оказывается равным или большим, чем заданное пороговое значение, то тестируемая скважина считается не прошедшей тест, причем по меньшей мере одним дополнительным проверяющим тестом является дополнительный проверяющий тест, на котором осуществляют проверку на предмет превышения коэффициентом качества корреляции заданного порогового значения, причем значение коэффициента качества корреляции определяется как коэффициент отклонения ближайшего по значению локального максимума функции корреляции от ее наибольшего локального максимума в точке Z_max, причем если значение коэффициента качества корреляции оказывается меньшим, чем заданное пороговое значение, то тестируемая скважина считается не прошедшей тест;
5) добавляют найденную скважину W к опорной группе скважин; и
6) итеративно повторяют все этапы способа до тех пор, пока на этапах 1) или 2) или 3) способа не будет определено ни одной скважины W.

2. Способ по п.1, отличающийся тем, что дополнительным проверяющим тестом осуществляют проверку на предмет превышения заданного порогового значения степенью транзитивности, которая определяется как количество предыдущих итераций алгоритма, для которых максимум функции корреляции найденной опорной скважины со скважиной текущей итерации не превышает максимум функции корреляции C(z) в точке Z_max.

3. Способ по п.2, отличающийся тем, что дополнительным проверяющим тестом устанавливают заданный интервал значений и осуществляют проверку на наличие скважин, которые не входят в указанный интервал, для их последующего отсеивания.

4. Устройство для определения положения координаты глубины маркера скважины W при построении геологической модели месторождения полезных ископаемых, содержащее:
один или множество процессоров;
модули ввода/вывода (I/O);
машиночитаемый носитель данных (память), который содержит код программы, который при выполнении побуждает процессор или процессоры выполнять этапы, на которых:
1) определяют опорную группу скважин;
2) определяют скважины W, расположенные в заданной окрестности опорной группы скважин, при этом окрестность не превышает значение расстояния R от любой из скважин опорной группы;
3) определяют функцию корреляции C(z) и выбирают для каждой найденной на предыдущем этапе скважины W соответствующую ей скважину из опорной группы с максимальным значением коэффициента корреляции и соответствующую ей точку Z_max на скважине W;
4) проверяют скважину W с помощью проверяющих тестов и, если эта проверка оказывается успешной, то выбор скважины считается утвержденным, а глубина, соответствующая максимальному значению корреляции Z_max, назначается глубиной маркера на скважине W, причем по меньшей мере одним проверяющим тестом является проверяющий тест, на котором осуществляют проверку на предмет превышения максимумом функции корреляции C(z) в точке Z_max заданного порогового значения, причем значение максимума функции корреляции во время выполнения проверяющего теста округляют до десятых долей в сторону большего значения, причем если значение максимума корреляции при округлении оказывается равным или большим, чем заданное пороговое значение, то тестируемая скважина считается не прошедшей тест, причем по меньшей мере одним дополнительным проверяющим тестом является дополнительный проверяющий тест, на котором осуществляют проверку на предмет превышения коэффициентом качества корреляции заданного порогового значения, причем значение коэффициента качества корреляции определяется как коэффициент отклонения ближайшего по значению локального максимума функции корреляции от ее наибольшего локального максимума в точке Z_max, причем если значение коэффициента качества корреляции оказывается меньшим, чем заданное пороговое значение, то тестируемая скважина считается не прошедшей тест;
5) добавляют найденную скважину W к опорной группе скважин; и
6) итеративно повторяют все этапы способа до тех пор, пока на этапах 1) или 2) или 3) способа не будет определено ни одной скважины W.

5. Устройство по п.4, отличающееся тем, что дополнительным проверяющим тестом осуществляют проверку на предмет превышения заданного порогового значения степенью транзитивности, которая определяется как количество предыдущих итераций алгоритма, для которых максимум функции корреляции найденной опорной скважины со скважиной текущей итерации не превышает максимум функции корреляции C(z) в точке Z_max.

6. Устройство по п.5, отличающееся тем, что дополнительным проверяющим тестом устанавливают заданный интервал значений и осуществляют проверку на наличие скважин, которые не входят в указанный интервал, для их последующего отсеивания.

7. Устройство по п.4, отличающееся тем, что дополнительным проверяющим тестом устанавливают заданный интервал значений и осуществляют проверку на наличие скважин, которые не входят в указанный интервал, для их последующего отсеивания.

8. Машиночитаемый носитель данных, пригодный для использования в устройстве по любому из пп.4-7, который содержит код программы, который при выполнении побуждает процессор или процессоры выполнять действия способов по любому из пп.1-3.